Imperial College London

Professor Jan Cilliers

Faculty of EngineeringDepartment of Earth Science & Engineering

Chair in Mineral Processing
 
 
 
//

Contact

 

+44 (0)20 7594 7360j.j.cilliers

 
 
//

Location

 

RSM 1.46BRoyal School of MinesSouth Kensington Campus

//

Summary

 

Publications

Publication Type
Year
to

155 results found

Vega D, Brito-Parada PR, Cilliers JJ, 2019, Small hydrocyclones for classiffication of particles in the micron range, Pages: 2398-2405

Small diameter (10 mm) hydrocyclones have been applied successfully for the separation of particle suspensions in the micron size range. These hydrocyclones are attractive because they show a bypass fraction larger that the water recovery, resulting in a high particle recovery to the underflow as well as low water recovery. However, this is a disadvantage when the purpose of the hydrocyclone is classification because of the large amount of fine particles that are misplaced in the underflow. The aim of this study is to explore, experimentally and computationally, the influence of design parameters on the classification process. In this work, a full factorial experimental design was defined to carry out comprehensive experimental tests using glass beads (0-20 µm) as the particulate system. We show that the dimensions of spigot and vortex finder diameter can be effectively manipulated to change the separation performance of the system and the energy consumption. A CFD model was developed that is able to predict particle size distribution. The numerical results for the partition curves showed very good agreement with the experimental data.

Conference paper

Morrison A, Brito-Parada P, Cilliers J, 2019, Developing a design modification for improved froth flotation performance through minimising turbulence at the pulp-froth interface, Pages: 1739-1747

The separation of valuable and gangue minerals according to their surface hydrophobicity in a froth flotation tank occurs in both its pulp and froth phases. In the pulp phase, hydrophobic particles must be brought into contact with rising bubbles introduced at or near the bottom of the tank. In a mechanical froth flotation tank, this is primarily achieved through the agitation of the pulp by an impeller or by a rotor-stator system. However, the turbulent and mixing effects of such an impeller or rotor-stator system, necessary to promote bubble-particle interactions in the pulp, are not confined to the pulp zone, but are transmitted into the froth zone at the pulp-froth interface. This impeller-induced turbulence at the pulp-froth interface is compounded by the impulse of buoyant bubbles entering the base of the froth phase, reducing the stability, and thus separation performance, of the resulting froth. Hence, turbulence improves separation performance if it occurs deep in the pulp phase, but diminishes it if it occurs close to the pulp-froth interface. For this work, a performance-enhancing pulp phase insert was designed to isolate the impeller-induced turbulence from the pulp-froth interface in a laboratory-scale, continuously-operated froth flotation tank. The result was a shifted grade-recovery curve that allows the insert to be designed to maximize the recovery from the system without sacrificing the grade of the concentrate, or vice versa, along a spectrum of improved performance. This paper will present the insert design and development process that resulted in this improved grade-recovery curve, as well as a set of heuristics that could be used to design a similar insert for pilot- and industrial-scale validation and deployment.

Conference paper

Reyes F, Cilliers JJ, Neethling SJ, 2019, Quantifying mineral liberation by grade and surface exposure using X-ray micro-tomography for flotation processes, Pages: 3985-3994

Liberation is a key driver in all mineral separation processes as it limits the maximum possible grade for a given recovery. In flotation, this is further complicated by the fact that it is surface exposure of the floatable minerals that determines the ultimate performance. Liberation, grade and surface exposure are commonly quantified using Scanning Electron Microscopy coupled to Energy Dispersive X-ray spectroscopy (SEM/EDX) analysis of polished sections. The intrinsically 2D nature of this technique can result in significant sampling errors and stereological effects that can affect the quantification of the ore's textural characteristics. X-ray micro-Tomography (XMT) is an imaging method that can non-invasively and non-destructively delineate ore fragments in 3D, thus providing an alternative method that eliminates the need for stereological corrections and readily provides surface exposure. A methodology and automated algorithm were designed for extracting this information from images of closely packed particles, thus allowing samples containing a large number of particles to be assessed. The main drawback of XMT is that, unlike SEM/EDX, it cannot directly measure the mineralogy of the sample, instead producing a 3D X-ray attenuation map. We therefore also present an algorithm for calibrating the thresholding of the XMT images based on SEM/EDX images of sections through the same sample, thus allowing the liberation analysis to be carried out on a 3D mineral map in which the uncertainty in the mineral assignment is small and statistically quantified. The methodology was tested on low grade porphyry copper ore as this is both an industrially relevant and traditionally difficult system to quantify using XMT due to the similarity in the X-ray attenuation of the sulphide species, especially the pyrite and chalcopyrite. As each 3D volume imaged contain 100-1000s of grains, large data sets can be readily produced. By dividing these particles into classes based on both

Conference paper

Reyes F, Cilliers JJ, Neethling SJ, 2019, Quantifying mineral liberation by grade and surface exposure using X-ray micro-tomography for flotation processes, Pages: 1570-1579

Liberation is a key driver in all mineral separation processes as it limits the maximum possible grade for a given recovery. In flotation, this is further complicated by the fact that it is surface exposure of the floatable minerals that determines the ultimate performance. Liberation, grade and surface exposure are commonly quantified using Scanning Electron Microscopy coupled to Energy Dispersive X-ray spectroscopy (SEM/EDX) analysis of polished sections. The intrinsically 2D nature of this technique can result in significant sampling errors and stereological effects that can affect the quantification of the ore's textural characteristics. X-ray micro-Tomography (XMT) is an imaging method that can non-invasively and non-destructively delineate ore fragments in 3D, thus providing an alternative method that eliminates the need for stereological corrections and readily provides surface exposure. A methodology and automated algorithm were designed for extracting this information from images of closely packed particles, thus allowing samples containing a large number of particles to be assessed. The main drawback of XMT is that, unlike SEM/EDX, it cannot directly measure the mineralogy of the sample, instead producing a 3D X-ray attenuation map. We therefore also present an algorithm for calibrating the thresholding of the XMT images based on SEM/EDX images of sections through the same sample, thus allowing the liberation analysis to be carried out on a 3D mineral map in which the uncertainty in the mineral assignment is small and statistically quantified. The methodology was tested on low grade porphyry copper ore as this is both an industrially relevant and traditionally difficult system to quantify using XMT due to the similarity in the X-ray attenuation of the sulphide species, especially the pyrite and chalcopyrite. As each 3D volume imaged contain 100-1000s of grains, large data sets can be readily produced. By dividing these particles into classes based on both

Conference paper

Mackay I, Cilliers JJ, Videla AR, Brito-Parada PRet al., 2019, Optimising froth stability of copper flotation tailings, Pages: 1730-1738

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 optimise the froth stability of tailings flotation for enhanced metallurgical performance.

Conference paper

Brito Parada P, Dewes R, Vega Garcia D, Cilliers Jet al., 2018, The influence of design parameters on the separation of biomass in mini-hydrocyclones, Chemical Engineering and Technology, Vol: 41, Pages: 2323-2330, ISSN: 0930-7516

Small hydrocyclones are an attractive technology for biomass separation from fermentation processes. The interactive effect of design parameters on the performance of mini‐hydrocyclones is, however, not fully explored and studies are often limited by the challenges in manufacturing such small units. Here, 10‐mm mini‐hydrocyclones are produced by 3D printing and the impact of spigot diameter, vortex finder diameter and height on separation performance is studied. A central composite rotatable design was adopted to obtain information on the relation between the variables and their influence on concentration ratio and recovery of yeast from a highly diluted system. A Pareto front for separation performance was generated and shown to be suitable to select an optimal design for a set of process constraints.

Journal article

Cilliers J, Hutahaean J, Brito Parada PR, 2018, A multi‐criteria decision framework for the selection of biomass separation equipment, Chemical Engineering and Technology, Vol: 41, Pages: 2346-2357, ISSN: 0930-7516

For the first time, a two‐stage decision support framework for equipment selection, applied to biomass separation, is presented. In the first stage, the framework evaluates from a number of equipment based on the process requirements and outputs only those that offer a technically feasible separation. In the second stage, the analytic hierarchy process is applied for performing a multicriteria decision analysis to select amongst the feasible equipment based on separation performance and energy consumption criteria. This approach systematically considers the relative importance of those different alternatives and selection criteria by pairwise comparisons. The output of the framework is an overall ranking of equipment as well as a sensitivity analysis of the results for different weighting of the criteria. These results can be used to equip practitioners in the field of bioseparations with a tool for making more consistent and better‐informed equipment selection decisions.

Journal article

Vega-Garcia D, Brito Parada P, Cilliers JJ, 2018, Optimising small hydrocyclone design using 3D printing and CFD simulations, Chemical Engineering Journal, Vol: 350, Pages: 653-659, ISSN: 1385-8947

The use of small hydrocyclones for the separation of particles in the micron range is of growing interest. However, these hydrocyclones are typically limited to conventional shapes or restricted to specific outlet sizes, which can lead to sub-optimal performance. The aim of this study is to present a method for the optimisation of small hydrocyclone design. This method consists of four steps that combine designing, Computational Fluid Dynamics (CFD) simulations, 3D printing and experimental testing. A 3D printed 10 mm hydrocyclone was shown first to match the performance of a ceramic equivalent, followed by factorial experiments with a set of printed hydrocyclones of different spigot and vortex finder diameters. A CFD model for small hydrocyclones was implemented and, following validation with the experimental data, used to simulate small hydrocyclone designs with parabolic walls. The model predicted improved separation performance compared to the conventional conic wall designs. In a novel development, a 10 mm hydrocyclone with parabolic walls was 3D printed and the prediction confirmed experimentally. The solids recovery and concentration ratio were increased by 10 percentage points and 0.2, respectively, for a 0.5 g/L yeast suspension and at an equivalent pressure drop. The use of 3D printing to manufacture small hydrocyclones of various designs has been proven in this study to be practical and to allow rapid prototyping design informed by CFD simulations. This is a significant improvement in the cost, time and versatility associated to hydrocyclone design and can lead to enhanced separation performance.

Journal article

Shean B, Hadler K, Neethling S, Cilliers JJet al., 2018, A dynamic model for level prediction in aerated tanks, Minerals Engineering, Vol: 125, Pages: 140-149, ISSN: 0892-6875

Stirred aerated tanks are a key unit operation in many industries, including froth flotation. Reliable and robust level control is of great importance in maintaining steady operation for successful implementation of higher level optimising control strategies, particularly when such tanks are arranged in series. When changes are made to the rate of aeration, there is a corresponding change in the pulp bubble size and gas holdup (the volume fraction of air in the tank), and consequently the pulp height. Stable operation of flotation tanks must, therefore, include the effect of air rate on pulp height in level control systems, especially if air rate is being actively controlled. In this paper, a model is developed from first principles to link the change in gas holdup with variation in air rate under dynamic conditions, accounting for the variability in gas holdup with height that results from differences in gas compressibility. This is validated experimentally. In order to test the model, experiments were carried out using a 70 L laboratory tank comprising water and reagent systems. For both simple and complex changes in air rate, the model showed good agreement with the experimental results when predicting the change in pulp height at steady state. Under dynamic conditions, the experimental system exhibited a slightly slower response than is predicted by the model; this is likely to be due to the well mixed assumption not being adequately met. This model provides a method to improve the operating stability of aerated tanks through better modelling of the dynamic pulp height changes that result from changes in air flowrate. In flotation tanks, this will enable greater control over froth height, which has been found to affect significantly mass pull, froth stability and flotation performance.

Journal article

Reyes F, Lin Q, Cilliers JJ, Neethling SJet al., 2018, Quantifying mineral liberation by particle grade and surface exposure using X-ray microCT, Minerals Engineering, Vol: 125, Pages: 75-82, ISSN: 0892-6875

Liberation is a key driver in all mineral separation processes as it limits the maximum possible grade for a given recovery. In flotation, this is further complicated by the fact that it is surface exposure of the floatable minerals that determines the ultimate performance. Liberation, grade and surface exposure are commonly quantified using Scanning Electron Microscopy coupled to Energy Dispersive X-ray spectroscopy (SEM/EDX) analysis of polished sections. The intrinsically 2D nature of this technique can result in significant sampling errors and stereological effects that can affect the quantification of the ore's textural characteristics. X-ray microCT (XMT) is an imaging method that can non-invasively and non-destructively delineate ore fragments in 3D, thus providing an alternative method that eliminates the need for stereological corrections and readily provides surface exposure. A methodology and automated algorithm were developed for extracting this information from images of closely packed particles. By dividing these particles into classes based on both their surface exposure and grade, the extent to which there is preferential breakage of the particles can be assessed—an important consideration if sufficient surface liberation for good flotation performance is to be achieved at coarser particle sizes. Using low energy scanning simple 3D mineral maps can be obtained via XMT, allowing for the assessment of liberation and surface exposure for each mineral species. The methodology was tested on low grade porphyry copper ore as this is representative of the most commonly treated ore types for copper production.

Journal article

Mackay I, Mendez E, Molina I, Videla AR, Cilliers JJ, Brito Parada PRet al., 2018, Dynamic froth stability of copper flotation tailings, Minerals Engineering, Vol: 124, Pages: 103-107, ISSN: 0892-6875

In this work, dynamic froth stability is used for the first time to investigate the flotation behaviour of copper tailings. Reprocessing of material from tailings dams is not only environmentally desirable, but also increasingly economically feasible as head grades can be high compared to new deposits. Flotation tailings, however, usually contain a large proportion of fine (10–50 m) and ultra fine (<) material and the effect of these particle sizes on froth stability is not yet fully understood.For this study, samples were obtained from the overflow and underflow streams of the primary hydrocyclone at a concentrator that reprocesses copper flotation tailings. These samples were combined in different ratios to assess the dynamic froth stabilities at a wide range of particle size distributions and superficial gas velocities. The findings have shown that the effect of particle size on dynamic froth stability can be more complex than previously thought, with a local maximum in dynamic froth stability found at each air rate. Moreover, batch tests suggest that a local maximum in stability can be linked to improvements in flotation performance. Thus this work demonstrates that the dynamic froth stability can be used to find an optimum particle size distribution required to enhance flotation. This also has important implications for the reprocessing of copper tailings as it could inform the selection of the cut size for the hydrocyclones.

Journal article

Dobson KJ, Harrison STL, Lin Q, Bhreasail AN, Fagan-Endres MA, Neethling SJ, Lee PD, Cilliers JJet al., 2017, Insights into ferric leaching of low grade metal sulfide-containing ores in an unsaturated ore bed using x-ray computed tomography, Minerals, Vol: 7, ISSN: 2075-163X

The distribution of the metal-bearing mineral grains within a particulate ore prepared for leaching, and the impact of this spatial heterogeneity on overall extraction efficiency is of key importance to a mining industry that must continuously target ever-reducing grades and more complex ore bodies. If accessibility and recovery of the target minerals is to be improved, a more detailed understanding of the behaviour of the system must be developed. We present an in situ analysis using X-ray computed tomography to quantify the rates of volume reduction of sulfide mineral grains in low grade agglomerated copper bearing ores during a miniature laboratory scale column leaching experiment. The data shows the scale of the heterogeneity in the leaching behaviour, with an overall reduction of sulphide mineral grains of 50%, but that this value masks significant mm3 to cm3 scale variability in reduction. On the scale of individual ore fragments, leaching efficiency ranged from 22% to 99%. We use novel quantitative methods to determine the volume fraction of the sulfide that is accessible to the leachate solution.

Journal article

Norori-McCormac A, Brito Parada P, Hadler K, Cole K, Cilliers JJet al., 2017, The effect of particle size distribution on froth stability in flotation, Separation and Purification Technology, Vol: 184, Pages: 240-247, ISSN: 1873-3794

Separation of particles of different surface properties using froth flotation is a widely-used industrial process, particularly in the minerals industry where it is used to concentrate minerals from ore. One of the key challenges in developing models to predict flotation performance is the interdependent nature of the process variables and operating parameters, which limits the application of optimising process control strategies at industrial scale. Froth stability, which can be quantified using air recovery (the fraction of air entering a flotation cell that overflows in the concentrate as unburst bubbles), has been shown to be linked to flotation separation performance, with stable froths yielding improved mineral recoveries. While it is widely acknowledged that there is an optimum particle size range for collection of particles in the pulp phase, the role of particle size on the measured air recovery and the resulting link to changes in flotation performance is less well understood. This is related to the difficulty in separating particle size and liberation effects.In this work, the effects of particle size distribution on air recovery are studied in a single species (silica) system using a continuous steady-state laboratory flotation cell. This allows an investigation into the effects of particle size distribution only on froth stability, using solids content and solids recovery as indicators of flotation performance. It is shown that, as the cell air rate is increased, the air recovery of the silica system passes through a peak, exhibiting the same froth behaviour as measured industrially. The air recovery profiles of systems with three different particle size distributions (d80 of 89.6, 103.5 and 157.1 μm) are compared. The results show that, at lower air rates, the intermediate particle size distribution (103.5 μm) yields the most stable froth, while at higher air rates, the finest particles (89.6 μm) result in higher air recoveries. This is subseq

Journal article

Tong M, Cole K, Brito-Parada PR, Neethling S, Cilliers JJet al., 2017, Geometry and Topology of Two-Dimensional Dry Foams: Computer Simulation and Experimental Characterization, LANGMUIR, Vol: 33, Pages: 3839-3846, ISSN: 0743-7463

Pseudo-two-dimensional (2D) foams are commonly used in foam studies as it is experimentally easier to measure the bubble size distribution and other geometric and topological properties of these foams than it is for a 3D foam. Despite the widespread use of 2D foams in both simulation and experimental studies, many important geometric and topological relationships are still not well understood. Film size, for example, is a key parameter in the stability of bubbles and the overall structure of foams. The relationship between the size distribution of the films in a foam and that of the bubbles themselves is thus a key relationship in the modeling and simulation of unstable foams. This work uses structural simulation from Surface Evolver to statistically analyze this relationship and to ultimately formulate a relationship for the film size in 2D foams that is shown to be valid across a wide range of different bubble polydispersities. These results and other topological features are then validated using digital image analysis of experimental pseudo-2D foams produced in a vertical Hele–Shaw cell, which contains a monolayer of bubbles between two plates. From both the experimental and computational results, it is shown that there is a distribution of sizes that a film can adopt and that this distribution is very strongly dependent on the sizes of the two bubbles to which the film is attached, especially the smaller one, but that it is virtually independent of the underlying polydispersity of the foam.

Journal article

Fagan-Endres MA, Cilliers JJ, Sederman AJ, Harrison STLet al., 2017, Spatial variations in leaching of a low-grade, low-porosity chalcopyrite ore identified using X-ray mu CT, Minerals Engineering, Vol: 105, Pages: 63-68, ISSN: 0892-6875

This study presents an investigation, using 3D X-ray micro computed tomography (μCT), into the effect of sulfide mineral position within an ore particle on leaching efficiency. Three sections of an unsaturated mini-leaching column that had been packed with agglomerated low-grade, low-porosity chalcopyrite ore and leached with an acidified ferric iron solution were imaged at different stages of a 102 day experiment. Image analysis was used to quantify changes in the mineral content and the influence on this of the mineral distance from the ore particle surface, local voidage and radial position within the column. The main factor affecting the mineral recovery was identified to be proximity of the mineral to the ore particle surface, with recovery decreasing with increasing distance from the ore surface. A maximum leaching penetration was observed to exist at 2 mm from the surface, beyond which no recovery was achieved. Higher recoveries at the column wall indicated that preferential flow in this higher voidage had an additional, albeit smaller, impact on leaching efficiency.

Journal article

Shean B, Hadler K, Cilliers JJ, 2017, A flotation control system to optimise performance using peak air recovery, Chemical Engineering Research and Design, Vol: 117, Pages: 57-65, ISSN: 0263-8762

Automatic control of industrial flotation cells and circuits presents a set of significant challenges due to the number of variables, the sensitivity of flotation cells to variation in these variables and the complexity of predicting flotation performance and/or developing a strategy for optimisation. Air recovery, a measure of froth stability, has been shown to pass through a peak as flotation cell aeration increases. Furthermore, the air rate at which the peak air recovery (PAR) is obtained results in optimal flotation performance, whether improved concentrate grade, recovery or both grade and recovery. Peak air recovery, therefore, presents a clear optimising control strategy for the operation of flotation cells which is generic to all flotation cells regardless of position in the flotation circuit. In this study, a novel control system based on PAR is developed and demonstrated using a large continuous laboratory flotation cell.In this study, a direct search optimisation algorithm based on the GSS (generating set search) methodology was developed using a 70 l continuous flotation cell operating with a two-phase system (surfactant solution and air only). Characterisation of the laboratory system showed that it was stable for up to 6 h and exhibited a reproducible peak in air recovery. A dynamic model of the response of the system with regards to changes in air recovery was developed that allowed simulations of the proposed optimising control system to be carried out. The optimisation algorithm was then applied to the experimental system. The trialled GSS algorithm was shown to find the PAR air rate when starting above, below and at the PAR air rate, and additionally with a disturbance introduced into the system. While the direct search approach can be slow, it is simple and robust. This demonstrates an innovative approach to optimising control for froth flotation and is the first application of froth stability maximisation for flotation control.

Journal article

Kiziroglou M, boyle D, Yeatman E, Cilliers Jet al., 2016, Opportunities for sensing systems in mining, IEEE Transactions on Industrial Informatics, Vol: 13, Pages: 278-286, ISSN: 1551-3203

Pervasive sensing - the capability to deploy large numbers of sensors, to link them to communication networks, and to analyze their collective data - is transforming many industries. In mining, networked sensors are already used for remote operation, automation including driverless vehicles, health and safety, and exploration. In this paper, the state-of-the-art sensing and monitoring technologies are assessed as solutions against the main challenges and opportunities in the mining industry. Localization, mapping, remote operation, maintenance and health and safety are identified as the main beneficiaries, from rapidly developing technologies such as 3D visualization, augmented reality, energy autonomous sensor nodes, distributed sensing, smart network protocols and big data analytics. It is shown that the identification and management of ore grade in particular, which transcends each stage of the mining process, may critically benefit from certain arising sensing technologies, where major efficiency improvements are possible in exploration, extraction, haulage and processing activities.

Journal article

Boucher D, Jordens A, Sovechles J, Langlois R, Leadbeater TW, Rowson NA, Cilliers JJ, Waters KEet al., 2016, Direct mineral tracer activation in positron emission particle tracking of a flotation cell, Minerals Engineering, Vol: 100, Pages: 155-165, ISSN: 0892-6875

Understanding the complex interplay of physics and chemistry inside a flotation cell is the ultimate goal of most flotation research. Key to the development of a model of flotation is the ability to validate it from measurements of a real flotation system. This work uses positron emission particle tracking (PEPT) to track directly activated mineral particles, hydrophobic and hydrophilic, in a lab-scale flotation cell. In contrast to other particle activation methods the direct activation technique allows mineral particles with their original surface characteristics to be used in PEPT experiments. In this work the flotation separation investigated was the separation of hematite from quartz from a synthetic ore using a combination of an oleic acid collector and sodium silicate depressant. This work represents the first time in which particles of typical flotation size (−106 + 90 μm diameter) with real bulk mineral properties and surface chemistry have been tracked in a flotation cell. The results illustrate small particles flow behaviour in the cell for a hydrophilic particle. The trajectory and velocities of the tracer particle are shown as it is transported inside the flotation cell.

Journal article

Bhutani G, Brito Parada PR, Cilliers JJ, 2016, Polydispersed flow modelling using population balances in an adaptive mesh finite element framework, Computers and Chemical Engineering, Vol: 87, Pages: 208-225, ISSN: 1873-4375

An open-source finite element framework to model multiphase polydispersed flows is presented in this work. The Eulerian–Eulerian method was coupled to a population balance equation and solved using a highly-parallelised finite element code—Fluidity. The population balance equation was solved using DQMOM. A hybrid finite element–control volume method for solving the coupled system of equations was established. To enhance the efficiency of this solver, fully-unstructured non-homogeneous anisotropic mesh adaptivity was applied to systematically adapt the mesh based on the underlying physics of the problem. This is the first time mesh adaptivity has been applied to the external coordinates of the population balance equation for modelling polydispersed flows. Rigorous model verification and benchmarking were also performed to demonstrate the accuracy of this implementation. This finite element framework provides an efficient alternative to model polydispersed flow problems over the other available finite volumeCFD packages.

Journal article

Brito-Parada PR, Bhutani G, Cilliers JJ, 2016, A novel framework for modelling flotation column hydrodynamics

The physical processes occurring in flotation columns are very complex as they take place at different scales over a range of sizes and multiple phases. Column design and optimum operating parameters are usually obtained from empirical studies and although progress has been made on its numerical modelling, most column models rely on simplifications that do not fully capture the hydrodynamics. In particular, Computational Fluid Dynamics (CFD) models of the pulp phase in flotation usually assume the bubble population to be monodispersed. While this assumption may be reasonable in some situations, the polydispersed nature of the gas phase can have a significant impact on the predicted flotation rates. Techniques to model size distribution do exist in the literature and have been rigorously analysed in the past two decades for bubble columns, which are very similar to flotation columns from a hydrodynamics point of view. In fact, available CFD models of flotation columns are an extension of bubble column models, sometimes including additional equations for solids treated as passive scalars. This work focuses on the development of an efficient CFD model for flotation columns. The bubble size distribution in the column has been modelled using the Population Balance Equation (PBE), which is solved using the Direct Quadrature Method of Moments (DQMOM). Coupled with the PBE, the multiphase flow equations in the column have been modelled using the standard Eulerian–Eulerian approach. The solution of this coupled system has been made tractable in this work with the implementation of a highly parallelised adaptive-mesh finite element method. This powerful framework has been implemented in the open source CFD code Fluidity and it allows the carrying out of detailed simulations of flotation columns in 2D and 3D, with important implications for industrial column evaluation and design.

Conference paper

Morris G, Hadler K, Cilliers J, 2015, Particles in thin liquid films and at interfaces, Current Opinion in Colloid & Interface Science, Vol: 20, Pages: 98-104, ISSN: 1359-0294

The behaviour of small particles at interfaces and in thin liquid films has been studied for many years and recent advances in experimental and numerical techniques have allowed a wealth of new research to be conducted. This manuscript reviews the last five years of work investigating the effect of particle shape and packing density on their behaviour when attached to a thin liquid film or at an interface between two immiscible fluids. We discuss advances at the individual particle scale, covering shape and surface heterogeneity, as well as breakthroughs in experimental and numerical modelling of larger scale systems.

Journal article

Bhutani G, Brito-Parada PR, Cilliers J, 2015, Modelling polydispersed flotation columns using population balances in a finite element framework, Computational Modelling '15

Conference paper

Cole K, Brito-Parada PR, Morrison A, Govender I, Buffler A, Hadler K, Cilliers JJet al., 2014, Using positron emission tomography (PET) to determine liquid content in overflowing foam, Chemical Engineering Research & Design, Vol: 94, Pages: 721-725, ISSN: 1744-3563

Journal article

Morris GDM, Cilliers JJ, 2014, Behaviour of a galena particle in a thin film, revisiting dippenaar, INTERNATIONAL JOURNAL OF MINERAL PROCESSING, Vol: 131, Pages: 1-6, ISSN: 0301-7516

Journal article

Cole K, Buffler A, Cilliers JJ, Govender I, Heng JYY, Liu C, Parker DJ, Shah UV, van Heerden M, Fan Xet al., 2014, A surface coating method to modify tracers for positron emission particle tracking (PEPT) measurements of froth flotation, Powder Technology, Vol: 263, Pages: 26-30, ISSN: 0032-5910

Positron emission particle tracking (PEPT) is a technique by which particle behaviour can be measured in a system of flow. The quality of the measurement is related to the spatial and temporal precision of the PET scanner and the characteristics of the tracer, which must replicate physical and chemical properties of the system bulk. Tracer particles can be made from ion exchange resins which have a high capacity for the commonly used positron emitting radionuclides 18F or 68Ga. However, these resins have a polymer composition and are naturally hydrophilic, which limits their application in systems involving mineral particles. This work presents a method to modify ion exchange resins with a coating to change the physical properties of the tracer. Two types of tracer were fabricated in this way, with hydrophobic and hydrophilic surfaces, to investigate the behaviour of valuable and gangue minerals in froth flotation with PEPT. The PEPT data were used to determine the spatial occupancies of each tracer, showing that the hydrophobic tracer has the highest occupancy in the froth region and the hydrophilic tracer is rarely entrained.

Journal article

Morris GDM, Neethling SJ, Cilliers JJ, 2014, The stability of a thin liquid film supported by a double layer of spherical particles, COLLOIDS AND SURFACES A-PHYSICOCHEMICAL AND ENGINEERING ASPECTS, Vol: 443, Pages: 44-51, ISSN: 0927-7757

Journal article

Morris G, Neethling SJ, Cilliers JJ, 2014, Predicting the failure of a thin liquid film loaded with spherical particles, Langmuir: the ACS journal of surfaces and colloids, Vol: 30, Pages: 995-1003, ISSN: 0743-7463

A model is presented for predicting the failure of a thin liquid film stabilized by attached inert particles. A statistical analysis of roughly 3500 Surface Evolver1 simulations was used to identify the relationship between the packing density of the particles on the film, their contact angle distribution, and the capillary pressure required to rupture the film. The model presented allows a fast and simple method of calculating the range of pressures a thin film in a three-phase froth will fail at based upon three variables: the film loading, mean particle contact angle, and the standard distribution of contact angles round the mean. The predicted range of failure pressures can be used in simulations of bulk froth properties where bubble coalescence is an important factor governing the froth properties.

Journal article

Norori-McCormac A, Hadler K, Cilliers JJ, 2014, Peak Air Recovery: An investigation into the effect of particle size

Air recovery, the fraction of air entering a flotation cell that overflows the cell lip, is linked to flotation performance, and it has been shown through industrial testwork that operating flotation cells at air rates that yield the Peak Air Recovery (PAR) results in high mineral recovery. For operating parameters over which flotation plants have little control, such as particle size, investigating the effect on the position of PAR, and the subsequent effect on flotation performance, proves more challenging. To this end, we have developed a bench scale flotation system that runs continuously and exhibits froth behavior similar to that found industrially, allowing the study of such parameters as particle size in addition to air rate. The bench scale system comprises a 4 l baffled cylindrical flotation cell, stirred by a Rushton impeller. The concentrate is recycled back into the pulp, allowing for continuous operation. Glass beads are used as solid particles, allowing a wide range of particle sizes to be tested, with TTAB as collector and MIBC as frother. The system allows operation at air rates typical of industrial flotation systems (1-2 cms<sup>-1</sup> superficial gas velocity), in addition to yielding a froth that coalesces and bursts and exhibits a peak in air recovery as cell air rate in increased. Solids and water recovery is also measured. The results show that air recovery is sensitive to changes in particle size. The intermediate size distribution showed a peak in air recovery (PAR) that corresponded with the maximum solids flowrate, as found industrially. It is shown that finer particles do result in more stable froths and higher recoveries, particularly at low air rates. The effect of particle size on the position of PAR is not yet unequivocally determined from these results.

Conference paper

Randall EW, Wilkinson AJ, Cilliers JJ, Xie W, Neethling SJet al., 2014, Current pulse technique for electrical resistance tomography measurements, Pages: 493-501

Experiments were performed on a model system of a non-conducting object(s) in a conducting liquid. The performance of two reconstruction algorithms, linear back projection (LBP) and Newton-Raphson (NR), were compared. In the system where the conductivities are widely different, LBP can detect the position, but not the size of the object(s), while NR clearly delineates them.

Conference paper

Bhutani G, Brito-Parada PR, Hadler K, Cilliers JJet al., 2014, Modelling laminar multiphase dispersed flows using population balances in an adaptive mesh finite element framework, 6th European Conference on Computational Fluid Dynamics

Conference paper

This data is extracted from the Web of Science and reproduced under a licence from Thomson Reuters. You may not copy or re-distribute this data in whole or in part without the written consent of the Science business of Thomson Reuters.

Request URL: http://wlsprd.imperial.ac.uk:80/respub/WEB-INF/jsp/search-html.jsp Request URI: /respub/WEB-INF/jsp/search-html.jsp Query String: limit=30&id=00448171&person=true&page=2&respub-action=search.html