Search Documents

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 Shuji Owada, Mayumi Ito, Takehiko Nishimura, Toru Yamashita, Sadayuki Shinozaki, Takashi Ando, Tomoe Ota

Electrical disintegration (ED) was applied to Nantun coal (China) in order to improve liberation of coal substance for the pretreatment of producing “Hyper Coal” of which the combustion gas was planning to inject directly to turbine for power generation. It has been found that many cracks generated in the ED were existent at phase boundaries with mineral matters and liberation of coal substance became much higher than that in mechanical comminution methods such as roll crusher. Since the ED occurred more easily (with lower voltage applied) for the coals of higher mineral content, it is considered that this method could make a better contribution especially to comminute middling coals. We also found that the ED technique could apply not only single particles but also to coal particles assemblage, keeping almost the same or higher disintegration efficiency. In the mineral matters, Si and Al components were preferentially well liberated by the ED.

Jan 1, 2003

By C. Martinez, S. Helle, M. Pincheira

We studied the feasibility of automated mineralogy to discriminate between tennantite and enargite, the main sulfosalts found in the Rio Blanco deposit, Codelco-Chile, in order to predict the content of minor elements (As, Sb, Ag, Pb, Zn) that impact on the process and concentrate prices. The identification of these minerals is feasible using traditional microscopy, however not possible by automated mineralogy. This drawback prevents comparing the databases generated in the past by means of conventional microscopy with the information obtained at present by QEMSCAN. To differentiate enargite from tennantite the mineral composition of both sulfosalts was studied by electron microprobe. The presence of associated minor element contents of about 8% zinc in tennantites and not exceeding 0.07% in enargites were used to modify the SIP database of QEMSCAN to enable the quantification of these phases in samples of head, residue and concentrate coming from As-bearing domains of the deposits. Detection limits allowed to determine the mineral species in the head, concentrate and residues. The integration of methodologies such as automated mineralogy and microprobe can be used successfully to quantify the presence of some ore forming solid solution.

Jan 1, 2014

By K. Bru, D. B. Parvaz

"Electric pulse fragmentation (EPF) is an innovative technology that uses high-voltage pulsed power (HVPP) for the selective comminution of the treated material. Highly energetic electrical pulses are used to fracture rocks at mineral grain boundaries, resulting in extensive micro-cracks and/or complete fragmentation of the processed particles depending on the applied energy. In this study, the EPF technology was applied to a low-grade schist ore containing cassiterite, in order to investigate its influence on the selectivity of the fragmentation. The impact of EPF on subsequent comminution steps was studied by combining either the EPF technology or jaw crushing (used here as conventional crushing) with subsequent impact crushing. The results confirmed the potential of the EPF treatment for a selective fragmentation of the schist ore. This gave the potential advantages of reducing the total energy consumption in the comminution circuit and of improving the recovery of cassiterite during concentration processes. INTRODUCTION The rise in the world consumption of metals combined with the depletion of most of the high-grade ore reserves leads to an increased need to process ores of lower grade, with a more complex mineralogy. To recover the valuable minerals contained in such ores, fine grinding is required, but this process has high energy requirements and leads to the production of a large amounts of fines, which can impede the performance of later concentration processes such as gravimetric separation and flotation. Addressing this issue is a major challenge for the mining industry. Selective comminution methods have the potential to significantly contribute to the profitable exploitation of these low-grade ores by generating cracks specifically at mineral grain boundaries rather than randomly as in a conventional mechanical comminution system. This allows the liberation of valuable components from the feed material without overgrinding, i.e. this prevents unnecessary size reduction of already liberated grains, which leads to a reduction in production of fines and in energy consumption. In addition, rock not fragmented by the process are usually weakened by the dense fracture network generated, which increases grindability of the rock (Wang et al., 2011)."

Jan 1, 2018

By M. Evertsson, K. Bhadani, E. Hulthén, M. Bengtsson

"Modelling and simulations for mineral processing plants have been successful in replicating and predicting predefined scenarios of an operating plant. However, there is a need to explore and increase the potential of such simulations to make them attractive for users. One of the tools to increase the attractiveness of the simulations is through applying optimization schemes. Optimization schemes, applied on mineral processing simulations, can identify non-intuitive solutions for a given problem. The problem definition itself is subjective in nature and is dependent on the purpose of the operating plant. The scope of this paper is to demonstrate two optimization schemes: Multi-Objective Optimization (MOO) using a Genetic Algorithm (GA) and Multi-Disciplinary Optimization (MDO) using an Individual Discipline Feasible (IDF) approach. A two stage coarse comminution plant is used as a case plant to demonstrate the applicability of the two optimization schemes. The two schemes are compared based on the problem formulations, types of result and computation time. Results show that the two optimization schemes are suitable in generating solutions to a defined problem and both schemes can be used together to produce complementary results. INTRODUCTION The comminution process for a mineral processing plant consists of multiple crushing and classification stages and aims to produce fine material to liberate ore from the minerals. The modelling and simulation of such processes are well established using static process models (King, 2001; Napier-Munn et al., 1996). A recent development in the dynamic modelling for such processes has shown further increase in fidelity compared to the static process models (Asbjörnsson, 2015). The developed simulations can be utilized to design, operate and control mineral processing plants based on the requirements of users (Bhadani et al., 2017). Typically, the requirements in the operation of the comminution process are conflicting in nature, as for instance, reduction of power draw by process units, maximization of production of fine materials, and increasing utilization of the process units. These conflicting requirements leads to a Multi-Objective Optimization (MOO) problem to operate a mineral processing plant. The MOO problems are normally referred to as non-linear problems and its fundamental understanding with respect to process behaviour is difficult to achieve."

Jan 1, 2018

By Kal Sastry, Anshuman Prasad

This paper presents a population balance model (PBM) for describing the agglomeration process and select a suitable coalescence kernel that can simultaneously describe nucleation, layering and coalescence mechanisms. The PBM is solved numerically by making use of a size discretization procedure in order to simulate iron ore agglomeration in rotating drums where moist feed is added to preformed seeds. Subsequently, a comparison of computer simulation results with the previously published experimental data is provided. Finally, the paper discusses how the current approach can be further extended to the simulation of industrial iron ore agglomeration circuits operated in closed circuit.

Jan 1, 2003

By Gaurav Soni, Roe-Hoan Yoon

A first principle flotation model has been derived from the basic mechanisms involved in flotation. The model consists of a series of analytical equations for bubble generation, bubble-particle collision, attachment, detachment, and froth phase recovery. The model can predict flotation from both hydrodynamic and surface chemistry parameters such as bubble size, particle size, energy dissipation rate, ?-potentials, contact angles (?), etc. The model predictions have been compared with the results of a series of laboratory flotation tests conducted using mono-sized glass beads of different particle sizes and contact angles under different operating conditions. Some of the experiments were conducted using a 1:9 mixtures if hydrophobic silica and hydrophilic magnetite, so that the product graded could be readily determined by simple magnetic separation. The model predictions are in good agreement with the experimental results both in recovery and grade. A computer simulator developed on the basis of the model is capable of designing a preliminary flotation circuit on the basis of feed characteristics, e.g., size-by-liberation data, and the hydrodynamic information on the flotation cells to be incorporated in the circuit.

Jan 1, 2014

By Katarina Theander, Robert J. Pugh

The generation of froth with a critical degree of stability and height is an important aspect of the froth flotation. In this paper we have studied the synergism and anti-synergism in frothing behaviour in a mixed flotation collector/frother mixture (a model fatty acid; sodium oleate and a model frother; nonionic (polyethyleneoxide alkyl C14E6). Values of the collector / frother molecular interaction parameters (the ß-parameters) were determined from the surface tension data using regular solution theory. The frothing characteristics of the system were evaluated using a glass column equipped with a series of conductivity electrodes, which enabled the water content of the froth to be determined. From these measurements, the dynamic frothability and static froth stability were determined. It was concluded that the synergism/anti-synergism in frothability was depended on the molecular composition of the surfactant system, the diffusion of monomer to the thin film interface and on the breakdown of the mixed nonionic/anionic micelles.

Jan 1, 2003

By Marina Nicolova, Irena Spasova, Plamen Georgiev, Stoyan Groudev

Different strains of "silicate" bacteria were used to improve the quality of kaolin by different ways of treatment. All strains formed large slimy capsules consisting of exopolysaccharides and were related to different species of the genus Bacillus, mainly to B. circulans and B. mucilaginosus. It was found that even a short contact of some hours of well developed cultures of these bacteria resulted in improvement of the bending strength and the other ceramic properties of the kaolin. The improvement was connected with the action of some secreted microbial metabolites such as soluble heteropolysccharides, monosugars and organic acids. The best results were achieved by continuous cultivation of bacteria in the presence of kaolin with a high relative humidity and subjected to periodic stirring.

Jan 1, 2014

By R. Shiva Kumar, G. Sreedhar, S. Sharma

Some of the Indian origin iron ore deposits located in Bellary, Karnataka province are having peculiar characteristics of Iron ore deposits (Hematite) having Iron above 64% and 30-40% lump (-30+6mm) recovery which are not suitable to use in blast furnace due to not meeting the specification of Tumbler Index and Abrasion Index. In general for high grade Iron ores, process flow sheet involves crushing and screening (dry/wet). In the contest of high grade Iron ores located at Bellary sector of Karnataka state in India, lump generated from this process could not attain the required physical properties (Tumbler Index, Abrasion Index) and metallurgical properties , even though the grade of the lump ore is above 64% Fe. Washing techniques like Drum scrubber and Log washer play major role in improving the metallurgical properties of Iron ore lumps by removing the light and hard layers of clay, silica coated over the lumps. In this study the effects were compared of drum scrubbing and log washing on iron ore lumps. It is found that by incorporating simple washing technique i.e. drum scrubbing or log washing, Tumbler Index of Iron ore Lumps could be increased from 74% to 78% which is acceptable for use in Blast furnace. From the Drum scrubber and Log washer tests, it is clear that Log washers have an edge over Drum scrubber in removing hard coated clay on surface of Lumps. This simple washing technique facilitate us to utilize the lumps directly in blast furnace and eliminate the unnecessary crushing and grinding of the ore so as to meet the specifications of Sinter fines and pellet feed.

Jan 1, 2014

By R. K. Dwari

The electrostatic beneficiation of coal is based on tribocharging characteristics of ash forming minerals and coal particles. The tribocharging of quartz, pyrite and coal particles contacted with various metals and polymeric materials have been measured and the charge acquisition is assessed through surface energy calculations from liquid contact angle data. The contact angles, before and after tribocharging of solids, are measured using Krüss tensiometer and Washburn?s method where the sample holders in tensiometer are specially constructed with tribocharger materials. The charge polarity and amount of charge acquired by quartz, pyrite and carbon with metal tribochargers were found in good agreement with the reported work functions of the contacting surfaces. The charge results with polymer materials differed from the work function values, presumably caused by surface contamination. The surface energy of quartz calculated by all the theoretical approaches in the literature showed that the tribocharging increases the surface energy. Both polar and non-polar components computed from Fowkes and Owens-Wendt-Rabel-Kaelble approaches dis-play that both these components increased after triboelectrification. However, the polar component divided into acid and base components, as in van Oss approach, manifest decreasing acid part and increasing base part. Since quartz charged negatively during tribocharging with metal surfaces and therefore accepted electrons, the determined acid and base surface energy components are consistent with charge transfer process. The results also elucidate an explicit correlation between the charges generated by powders and the surface acceptor (acid) and donor (base) electronic states and thereby the work functions. A novel fluidized bed tribocharger was designed and developed with internal baffle system for efficient contact electrification of mineral and maceral to produce high differential charging for efficient separation. The results achieved on coal with 45% ash content are presented and discussed.

Jan 1, 2014

By Anuttam Patra, Romain Bordes, Andreas Fredriksson, Tommy Karlkvist, Krister Holmberg, Kota Hanumantha Rao

The present investigation aims to develop and distinguish mineral specific reagents with two functional groups for use in flotation of calcium containing minerals. For this purpose, a series of dicarboxylate-based surfactants with varying length between the carboxylate groups (one, two or three methylene groups) were synthesized. As reference, a surfactant with the same alkyl chain length but with only one carboxylate group in the polar part was synthesized. The adsorption behavior of these new reagents on pure apatite and calcite mineral surfaces was studied using Hallimond tube flotation, FTIR and ? potential measurements. The relation between the adsorption behavior of a given surfactant on a specific mineral surface and its molecular structure over a range of concentration and pH values, as well as the region of maximum recovery were established. It was found that one of the reagents, with a specific distance between the carboxylate groups, was much more selective for a particular mineral surface than the other homologues synthesized. This selective adsorption of a given surfactant to a particular mineral surface relative to other mineral surfaces as evidenced in flotation studies is substantiated by ? potential and infra-red spectroscopy data.

Jan 1, 2014

By Seungwoo Park, Gaurav Soni, Roe-Hoan Yoon, Hyunsun Do, Kaiwu Huang

In a flotation froth (and foam), air bubbles become larger due to coalescence, causing bubble size to increase, bubble surface area to decrease, and hence casing less hydrophobic particles to drop off to the pulp phase below. Thus, bubble coarsening provides an important mechanism by which product grades are increased. On the other hand, excessive bubble coarsening results in low recoveries. In the present work, a model describing the process of bubbles becoming coarser in a foam as they rise to the top by deriving a mathematical relation between the Plateau border area, which controls film drainage rate, and the lamella film thickness, which controls bubble-coalescence rate. The model has been derived by assuming that liquid drainage rate increases due to slip at the air/water interface, and that the critical rupture thickness is determined by the capillary waves created by thermal motion. The model developed in the present work can predict the bubble size ratio between the top and bottom of a foam as functions of aeration rate, foam height, and surface tension (frother dosage). The model predictions are in good agreement with the changes in bubble sizes measured using a high-speed camera.

Jan 1, 2014

By Y. Foucaud, C. Eswaraiah, L. Filippov

Enhanced gravity separators such as Falcon concentrator which uses additional centrifugal force, are the most common gravity concentration techniques used for fine particles processing. Hence, understanding the kinetics and separation mechanisms within these separators is of paramount interest. Recent research yielded a predictive physical model for the Falcon Ultra Fine (UF) which however does not explain some results obtained with industrial ores. This study aims to identify the mechanisms involved in Falcon UF separation with the objective of confirming or rejecting the hypotheses on which the current model is built and suggesting potential points of improvement. The Falcon UF kinetics is investigated through the processing of complex low-grade and fine-grained ores from the Altenberg tin deposit (Germany), the Tabuaço tungsten deposit (Portugal) as well as a synthetic iron ore. Saturation tests are performed on the Falcon UF by controlling the separation efficiency in terms of recovery and enrichment ratio with time or feed mass. Results show an evolution of Falcon UF performance with time in contradiction with the stationary separation hypothesis on which the physical model is based. Observation of the concentrate bed suggests the presence of erosion figures in furrows which may play locally an active role in the separation. In terms of Falcon UF separation timing, three phases can be distinguished. First, upon initial feeding of the bowl, a bed quickly grows unselectively and ultra-fine particles are ejected from the bowl until the bed reaches its final profile. Next, actual separation occurs with selective concentration of dense particles. Finally, capture sites saturate and recovery drops. The evolution of partition curves over time confirms the existence of two separation mechanisms. First, differential particles settling within the flowing film responsible for the ejection of small size particles during the whole separation cycle which is already accounted for in the physical model. The second mechanism involved may be explained by the action of lift force, neglected in the physical model, which acts preferentially on coarse particles deposited at the surface of the bed resulting in coarser particles being rejected in the tailings. Hence, the addition of a lift force component to the existing model is suggested to account for this second separation mechanism.

Jan 1, 2018

By Morin D. H., Battaglia-Brunet F., D'Hugues P.

"In the mineral-processing field, bioleaching of metal sulphide concentrates is currently operated in very large mechanically agitated tanks. Such bioreactors have considerable requirements in terms of power consumption for mixing, aeration and heat transfer. The use of slurry bubble column, where the particles are suspended by gas motion as well as by liquid upward flow, could be an alternative to achieve bioleaching in more economical conditions. In order to compare these two techniques, a unit of stirred reactors in cascade and a bubble column were comparatively operated at laboratory scale in the bioleaching of a pyrite concentrate. For both pieces of equipment, experiments were performed not only in batch mode but also in continuous conditions.Bubble column performances were significantly affected by unpredictable technical difficulties met to maintain continuous feeding in the bubble column and by a significant settling of the heavy solid particles. Nonetheless, theoretical kinetic data, extrapolated from results in batch conditions, could be compared to those obtained in continuous tests, in bubble column and in mechanically agitated reactors. These kinetic data were used to evaluate the influence of reactor design and configuration choice on bioleaching performances."

Jan 1, 2003

By D. B. Mazzinghy, L. C. de R. Machado, H. Delboni, M. G. Bergerman

"The use of vertical stirred mills in the mining industry has increased remarkably over the past few decades, as a result of growing requirement for finer ore grinding. This equipment delivers higher energy efficiency in fine grinding operations when compared to conventional tubular mills. Methods for designing vertical stirred mills involve operational experience, pilot plants and bench tests. An important issue is that these tests require, at a laboratory scale, at least 10–20 kg of material, depending on its density, which is not available in many cases, particularly in the early stages of greenfield projects. For the regrind of flotation concentrates, several bench scale flotation cycles are required to generate such a sample. This paper describes the development of a laboratory jar mill test using a 6”x8” jar – smaller than that commonly used – based on the overall specific energy method, which requires significantly less mass compared with existing tests. The proposed test has proven to be more practical than tests normally used, requiring about a quarter of the mass, and also shows good adherence to industrial results for the studied iron ore processing plant. INTRODUCTION A significant demand for finer grinding has been noticed over the past few decades in several mining projects to achieve the mineral liberation required for the concentration of valuable minerals. However, the finer the product size, the greater the grinding energy consumption, resulting in an overall increase in operating costs, which may reach or even exceed 50% of the mineral processing cost (VIEIRA, PERES, 2013; DANIEL, LEWIS-GREY, 2011). Comminution technology can have a major impact on subsequent ore concentration stages (DANIEL, LEWIS-GREY, 2011; MARSDEN, 2009; NAPIER-MUNN, 2012; PEASE et al., 2010), which further strengthens the importance of proper technology and equipment selection. Adequate grinding circuit selection and design is thus a key aspect associated with mining projects."

Jan 1, 2018

By Frank van der Meer, Lea Quincke

"High Pressure Grinding is applied to achieve size reduction of the particles by inter-particle breakage and to generate micro-cracks in the product particles as a contributing factor in e.g. leaching or mineral liberation. As alternative to classification, a closed circuit operation for HPGR can be established by product recycling, where a proportion of the coarser edge material and parts of the center material are split-off to be recirculated to achieve a finer final (center) discharge product for subsequent processing steps as ball mill grinding or beneficiation. As a result of the product recycling, the content of fines in the feed material increases and may affect the grinding efficiency in terms of production of fines and the generation of micro-cracks in the particles. In the scope of this paper, the effects of the proportion of fines in HPGR feed on grinding efficiency and micro-crack generation were evaluated. Grinding tests were executed with different feed size distributions, from a fines-void situation (truncated feed) to a fines-abundant situation (high product recycle). The evaluation of grinding efficiency and generation of micro-cracks was done by particle size analysis, Bond Work Index tests, impact tests on HPGR feed and products, and X-ray computed tomography of single product particles. This publication summarizes the main aspects of the outcome of these investigations."

Jan 1, 2014

By Lucy van de Ruit, Megan Becker, Aubrey Mainza, Jenny Wiese

"UG2 chromitite ore is a complex, highly disseminated platinum ore from the Bushveld Complex in South Africa, with typical platinum group mineral (PGM) grain sizes below 15µm. The introduction of stirred mills to the platinum industry has enabled grinding to finer product sizes, which has improved liberation and recovery of the finely grained PGMs. However, reduction in particle size leads to an increase in chromite entrainment in the flotation concentrate, which has a negative impact on the subsequent smelting process. Factors such as particle size, froth structure, solids concentration and water recovery have all been identified as impacting entrainment, but whether particle shape also has a significant effect is unclear. The purpose of this study was to determine whether the recovery of chromite particles by entrainment during the flotation of UG2 ore is affected by particle morphology. A laboratory vertical stirred mill with inert ceramic media was used for grinding UG2 ore which was then floated in a batch flotation cell. Flotation feed and concentrate samples were analyzed using QEMSCAN to quantify the angularity and aspect ratio distributions of the particles. The stirred mill product was found to have relatively broad aspect ratio and angularity distributions, with the major mineral components (chromite and silicates) having very similar shape characteristics. Analysis of flotation feed and concentrate samples suggested that particle shape does not have a significant effect on chromite entrainment."

Jan 1, 2014

By Zhang Yuejun, Feng Tianran, Shen Zhengchang, Han Dengfeng

"With continuous mining industry development and dressing plant scale expansion, flotation cell is developed towards to large scale and intelligence. Based on flotation cell scale-up theory and industrial practice, KYF-320 air-forced mechanical flotation cell with effective volume of 320m3 was successfully developed by Beijing General Research Institute of Mining and Metallurgy (BGRIMM) and adopted by China Gold Inner Mongolia Mining Co., Ltd, Luming project of China Railway Resource Group Co., Ltd, and Toromocho project of Chinalco in Peru. In this paper, the structure and technical characteristics of KYF-320 flotation cell are introduced, and the application of that equipment in China Gold Inner Mongolia Mining Co., LTD is demonstrated. The tank of KYF-320 flotation cell is cylinder-shaped, 8200mm in diameter and 7500mm high, with double froth lauders and large inclined angle cone bottom. 16 sets of KYF-320 flotation cells of which the processing capacity is 26,000,000 tons/a, are used in the second phase project of China Gold Inner Mongolia Mining Co., LTD. Their stable running shows the long life of impeller and stator and the power intensity (power consumption divided by effective volume of flotation cell) of KYF-320 flotation cell is only 0.73kW/m3 which is significantly less than that of small scale flotation cells. The industrial production obtains good results that the Cu grade of final concentrate of entire flotation circuit is 22.90% at the recovery of 86.73%, when Cu grade of feed is 0.35%, which indicates that KYF-320 flotation cell has good performance and can be widely used in large scale mine."

Jan 1, 2014

By A. D. Ramirez, L. E. Gutierrez, M. A. Garces

"Scarcity of freshwater in mining operations in arid zones in different parts of the world has led to the mining companies to think about using seawater as an alternative source in mineral processing operations. Examples of this situation are Minera Centinela (AMSA) and Minera Sierra Gorda (KGHM), both plants located in northern Chile. Previous studies showed that molybdenite flotation is depressed at high pH values in seawater as a result of the action of hydroxy-complexes of magnesium and calcium that attach to the faces of the molybdenite particles, reducing its hydrophobicity and as a result flotation recovery. In this work, a study of the interactions between molybdenite and kaolinite is presented. The effect of sodium hexametaphosphate (SHMP) dispersant on the recovery of molybdenite in the presence of kaolinite is studied in different aqueous media and it is correlated with induction time data. The results allow to conclude that kaolinite has a stronger depressing effect on molybdenite in seawater compared to its effect in freshwater. Magnesium and calcium are the cations that have the strongest effect in enhancing the interactions between molybdenite and kaolinite. SHMP reduces the negative effect of kaolinite on the floatability of molybdenite both in seawater and freshwater. The induction time data correlate well with the micro-flotation results.INTRODUCTION The need for water in mining processes in arid and semi-arid zones has led to the search for new alternatives for this resource, among which seawater has gained an important space. In Chile, most of the mining operations are located in the northern part of the country, an area that has limited availability of freshwater for the processing of copper and molybdenum sulfides (Castro and Laskowski, 2011, Castro et al., 2012 a). Seawater is an aqueous solution with high concentrations of ions such as Na+, K+, Mg2+, Ca2+, So42-, Cl-, and other less concentrated species which generate important effects in the flotation process as they modify the properties of the solid/liquid, and liquid/gas interfaces, and also increase the consumption of some important reagents used to depress pyrite such as lime. It has also been proposed that flotation of copper sulfides is negatively affected in seawater because of the reduction of oxygen solubility compared to freshwater (Millero et al., 2002; Duan and Sun, 2003). In contrast, it has been shown that the presence of ions in solution in the process of flotation of natural hydrophobic resources improve recovery as result of the compression of the electric double layer (Laskowski, 1974; Klassen and Mokrousov, 1963, Laskowski and Castro, 2008."

Jan 1, 2018