Search Documents

By S. Pule, Z. Ndabezimbi, A. Hefer, E. Cupido

This article outlines various projects implemented at Rand Refinery to reduce effluent generation and freshwater consumption. The refinery generates about 2800 m³ of effluent annually and has initiated projects such as granulation water recycle and silver leach optimisation to address these issues. The granulation water recycle project has successfully reduced effluent volumes, while the silver leach optimisation has further decreased effluent by implementing stoichiometric addition of HCl and NaClO3. The vibratory shear enhanced processing (VSEP) trial showed promising results for gold recovery from wastewater, though further trials are recommended. The water savings of the effluent reduction projects are monitored by measuring the amount of effluent pumped to service yard effluent storage facilities with the use of flow meters connected on SCADA. Since the effluent reduction initiative of granulation effluent recycle in FY23, the effluent reduced to 2676 m³. In FY24, the stoichiometric addition of HCl and NaClO3 was implemented in the silver leach process, further reducing the effluent to 2083 m³.

Mar 5, 2025

By K. B. Le Bron, T. van Aard

The empirical pillar strength formulae have been widely, and in most cases successfully, applied to design underground pillars for mines in the Bushveld Igneous Complex (BIC). This paper presents the case for combining the empirical method with numerical modelling when it comes to crush pillar stability analyses. The benefits include an understanding of potential failure mechanisms, which may then be accounted for in the designs. In addition, the impact of weak layers, such as shear zones, the impact of the dip of the orebody, differences in layouts, such as holing widths and the presence of advanced strike gulleys (ASGs) immediately adjacent to the crush pillars, and stress regime (k-ratio) may be studied using numerical modelling, as it is not accounted for in the empirical method.

Mar 5, 2025

By L. Godfrey, Z. Nontso, F. Magweregwede, B. Maphalala, S. Khan

The aim of this paper was to establish if mineral resource scarcity is a driver for South Africa to transition to a circular economy. The project objectives were achieved by investigating economically viable, and critical or strategic minerals remaining in South Africa, as well as resource scarcity, as a driver for South Africa to transition to a circular economy. The study found that some minerals in South Africa, such as iron ore, lead, manganese, gold, and cobalt have less than 50 years of economically viable mining remaining. At least 18 minerals can be classified indicatively as critical or strategic: aggregates, bauxite, chromite, cobalt, copper, gold, graphite, iron, limestone, lithium, manganese, nickel, phosphorous, platinum group metals, rare earth elements, silver, titanium, and vanadium. The classification was based on the minerals’ economic importance, supply risk, and potential use in the development and manufacturing of emerging technologies. The findings suggest that the scarcity of mineral resources is a critical driver for the country to transition to a circular economy. Other identified drivers include socio-economic factors, climate change commitments, and business objectives. All minerals are non-renewable resources and at risk of depletion, therefore sustainable extraction and use are critical. Priority should be directed at investigating circular economy opportunities and substitutes for commodities at risk of being depleted within the next 50 years. There is an opportunity for South Africa to transition from a minerals export-based model to one that promotes local mineral beneficiation and manufacturing of end products.

Mar 5, 2025

By I. O. Ondicho, E. K. Mutinda, B. O. Alunda, E. Agyekum

Driven by the necessity to improve blast performance regarding fragment size distribution in limestone mines, this paper introduces the prediction and measurement of blast fragmentation distribution through the modified Kuz-Ram (MKR) and Kuznetsov-Cunningham-Ouchterlony (KCO) models. Precise prediction and measurement of rock fragmentation resulting from blasting play a critical role in hard rock ore extraction. This is due to the substantial influence of post-blast ore size distribution on the efficiency of all subsequent downstream rock handling and processing operations. Additionally, employing secondary breakage of boulders, when necessary, leads to an unjustified rise in production costs. In this paper, the MKR and KCO models were utilised to identify the most appropriate model for application in Bisil and Simba quarries situated in Kajiado County, Kenya. Split-Desktop was used for blast fragment measurement by processing digital images obtained after blasting and computing the fragment size distribution of the blasted rock. As per the investigation, the recorded percentages of boulders generated at the sites ranged from 7% to 12.3% for the Simba quarry and from 11.2% to 23.9% for the Bisil quarry, respectively. Regression, correlation, and root mean square error (RMSE) were employed as the primary performance indicators. The study concluded that KCO outperformed MKR in accurately predicting blast fragment sizes in the examined quarries.

Mar 5, 2025

By J. R. Bredell, J. Branehog, C. Coetzee, J. J. van Tonder

A nutating grinding mill, exemplified by the commercial HICOM mill, achieves efficient particle size reduction through the accelerated motion of the material charge within the nutating chamber. This design, characterised by high particle velocities and rapid charge circulation, offers a compact solution compared to tumbling mills. There have been concerns regarding the mechanical reliability of kinetic joints for certain specialised mining applications. Recognising the need to model internal forces in such systems, this study focuses on the uncoupled simulation using multibody dynamics and the discrete element method. The investigation employed a small-scale experimental test mill, referred to as the NuMILL, featuring a cylindrical chamber (140 mm internal diameter, 245 mm height, 6° nutation angle). Instrumented with an array of sensors, the NuMILL provided internal load data for validating the simulations. Experiments were conducted with an empty chamber or filled with charge material. Charge materials included coarse sand or steel pellets with a charge-to-chamber mass ratio of 0.25. The mill speed ranged between 100 to 700 RPM. Multibody dynamics modelled resultant inertial and frictional loads within 5%, while the discrete element method simulated the internal material contact loads with accuracies ranging from 1% to 19%. The combined simulation approach predicted the resultant forces in the critical kinetic joint within a 20% margin. Representing the charge using rigid lumped mass elements was also studied. The presented simulation techniques, validated against experimental data, offer a robust foundation for analysing larger systems, including the industrial-scale HICOM.

Mar 5, 2025

By T. P. Chidakwa, T. R. Stacey

This study investigated the corrosion rates of various steel types and coatings used for friction rock stabilisers in a deep-level gold mine under diverse underground environmental conditions. Testing was conducted in situ within ramps, intake airways, and exhaust airways, and through simulated laboratory experiments. Results revealed that environmental variations significantly influence the corrosion performance of different materials and coatings. Duplex coatings, made up of galvanised steel with an organic paint layer, demonstrated the highest corrosion resistance. Zinc-based coatings provided moderate protection, while 3CR12 stainless steel and uncoated S420 carbon steel exhibited comparatively poor resistance under both aqueous and atmospheric conditions.

Mar 4, 2025

By B. S. Xakalashe, M. I. Khama, Q. G. Reynolds

The carbonaceous reduction of manganese ore is an energy-intensive process that releases large amounts of pollutant gases. In an attempt to circumvent the greenhouse gas emissions, this research investigates hydrogen as an alternative reductant. The reduction of manganese ore with hydrogen results in the structural changes in the intermediate products, which might limit the diffusion of the reactant gases through the pores. Computational fluid dynamics (CFD) models capture the structural changes and predict the rate controlling mechanism and this is key in designing the prereduction reactors. The current work uses CFD models to predict the degree of reduction in a shaft furnace under hydrogen reducing conditions. Manganese ore is described as a collection of scalar and tensor fields such as thermal conductivity of the porous medium, density, and viscous resistance. The results indicate that advection contributes significantly towards the overall transport of hydrogen to the reaction sites. Porosity of the final product was found to be higher than porosity of the raw material.

Mar 4, 2025

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 J. K. Pothal, S. Basu, R. Barik, P. K. Senapati, A. Dubey

The aim of this paper is to study the characteristics of limestone particles (size varies from 4 mm to 20 mm) during its transportation from open cast mines through pipes. The influence of four different size fractions of limestone samples (4–6 mm, 6–10 mm, 10–15 mm, and 15–20 mm) on pipe flow characteristics during vertical transport in a 100 mm NB hose pipe was investigated. The flow velocities in the vertically lifting pipe was varied from 0.95 m/s to 2.78 m/s and a rotary feeder of 2 to 4 tonne/hr. capacity was employed to feed the limestone samples to the lifting pipe. The pressure drop data for the different size fractions of limestone samples at different mixture velocities were analysed. The effect of input volume fraction and superficial velocity of limestone samples on hold-up ratio were discussed. It was indicated that there was no substantial increase in pressure drop of slurry samples in the studied range of concentrations (CW = 4.22% - 12.54%), particle size, and low to medium range mixture velocities with respect to water run head loss data. The hold-up for the larger size fractions (15–20 mm) of limestone samples was found to be significant as compared to smaller size fractions, which may be due to terminal settling velocities of individual fractions of the studied samples.

Mar 4, 2025

By J. Moema, L. Mampuru, M. Phasha, D. Mkhonto

The global liquified natural gas (LNG) storage market size is projected to reach USD17.5 billion by 2025. Currently, owed to its high toughness, tensile strength, and excellent weldability, the 9% nickel steel is the most used in building infrastructure for cryogenic applications such as the construction of tanks and pipes for storage and transportation of LNG. However, this alloy is not readily available because of its complex production process and high cost. This provides an opportunity for suitable alternative materials to build infrastructure for cryogenic applications. Among strong contenders is the Fe-Mn based alloys, which have sparked global interest due to their desirable properties such as relatively lower density, low cost, high toughness and strength, due to a high manganese (Mn) content. These attractive properties render Fe-Mn based alloys preferred candidates to replace traditional steels in engineering applications in which strengthweight ratio is critical. In response, as a build-up to designing and developing austenitic Fe- Mn alloys, the present study employed density functional theory (DFT) based first-principles computational materials modelling technique to investigate the structural, thermodynamic, and magnetic properties of binary Fe-Mn alloy composition in competing FCC, HCP, and BCC crystal structures. Using this approach, it was possible to unravel the key underlying elastic properties that are directly correlated to experimental tensile strength and high toughness in binary Fe-Mn alloys. Besides successfully validating the existing experimental data, the current predicted properties and the deployed approach will serve as a benchmark from which further alloying required to improve other properties will be conducted.

Mar 4, 2025

By B. S. Xakalashe, I. J. Geldenhuys, B. Bowman, M. W. Erwee, G. A. Venter, J. Zietsman, C. Sandrock, Q. G. Reynolds

Understanding the behaviour of plasma arcs is an important part of the design and operation of direct current electric arc furnace smelting processes, which are used in the industrial production of many metallurgical commodities. In recent years numerical and computational modelling techniques have begun to illuminate the complexity of arc behaviour and how it is deeply connected to the thermodynamic and physical properties of the arc plasma. Plasma properties may be estimated from first principles using statistical mechanics methods, but this requires as input some knowledge of the composition of the gas phase that occurs over the process during smelting. In this paper the authors describe a workflow using thermochemistry software to predict the expected gas phase composition through a given process, followed by plasma property calculations using an open-source software tool that has been in development for several years. This tool, 'minplascalc' has been in development for several years and is also described in detail here. The workflow is then applied to the calculation of plasma properties of interest for a variety of historical and current pyrometallurgical processes, and the results are compared. The plasma property data are published in an open-access database for general use.

Mar 4, 2025

By R. J. Durrheim, Y. Yabe, M. S. D. Manzi, DSeis team, N. Z. Nkosi, S. B. Mngadi, Jr. H. Ogasawara

On 5 August 2014, the Orkney M5.5 earthquake, the largest in South African gold mining districts, occurred with an unusual strike-slip mechanism at great depth of 4.78 km below the surface. In a rare case, the ML 5.5 earthquake occurred within the area covered by the legacy 2D and 3D reflection seismic data acquired in the 1990s and 2000s for gold exploration and mine development. In addition, the earthquake ML 5.5 rupture plane was recorded and accurately delineated by the underground in-mine seismic network near the source region. The integration of the legacy 2D reflection seismic data and mine seismicity data allowed us to identify a near-vertical structure, striking northnorthwest-south southeast (NNW-SSE). The ICDP-DSeis team drilled three holes (Hole A, Hole B, and Hole C), and two holes (Hole B and Hole C), intersected the upper edges of the ML 5.5 rupture. These holes recovered metasediments, metabasalts, intrusive rocks including dolerite sills and lamprophyre dykes adjacent to the fault zone, and the fault gouge. Late Prof. Tullis Onstott and geomicrobiologists installed a packer in Hole A and successfully recovered saline water and detected gas (~10 MPa). Slip weakening and rupture propagation are significantly influenced by the existence of fault gouges and the production of wear material between two sliding rock surfaces. Using the fault gouge material recovered from Hole C, we conducted friction experiments under dry and room temperature conditions at high slip velocity of ~100 mm/s and normal stress of 2 MPa. The resulting steady-state frictional strength was ~0.66 over a slip weakening distance of ~9.1 m. The steady-state frictional strength was high, which may be caused by large gouge thickness, and the rate of wear generation. As a result, it is proposed that the Moab Khotsong ML 5.5 seismogenic zone is complex and could be controlled by three main processes: a) the complex structural architecture of the seismogenic zone (e.g., intersection of fault, lamprophyre dykes, and dolerite sills); b) the mechanical process induced by tectonic and/or mining related stresses; and c) the mechanical and chemical processes caused by the water and rock interaction.

Jan 12, 2025

By I. James, R. J. Durrheim, M. S. D. Manzi, O. Valeshin, J. Lepine, L. V. Socco, C. Colombero, A. Malehmir, M. K. Rapetsoa, M. Sihoyiya

The FUTURE project aims to enhance seismic surveying in South African mining by integrating advanced acquisition technologies. In this paper, we detail surface and underground seismic surveys at the South Deep gold mine as part of the project. Surface seismic acquisition used co-located broadband micro-electromechanical sensors, one and three component recorders connected to geophones, distributed acoustic sensing using a straight fibre optic cable, along with a 6 ton broadband seismic vibrator operating with 4–180 Hz linear sweeps. In addition, underground surveys employed a combination of straight, and helically wound fibre optic cables, cabled-based and nodal sensors, along with a 500 kg drophammer mounted on a skid steer bobcat. The surface and underground seismic arrays were time synchronised using a GPS-time transmitter to allow for accurate seismic data recording in a GPS-denied environment. Results show the efficacy of these techniques in acquiring high-resolution seismic data, revealing detailed subsurface geological structures and mineralisation down to 3500 m below ground surface. Surface data revealed clear seismic reflections corresponding to the Black Reef and Ventersdorp Contact Reef. An underground P-wave velocity model exhibits low velocity anomalies that are associated with fracturing, while a P-wave reflection stack section shows significant reflections and a degree of faulting below the tunnel floor that was not observed on the surface reflection seismic data. We highlight the potential of these technologies to enhance safety, efficiency, and sustainability in mining. This study sets a benchmark for future reflection seismic surveys, paving the way for safer, more efficient, deep mineral exploration and extraction practices.

Jan 12, 2025

By R. J. Durrheim, B. Manzunzu, V. Midzi

Selecting ground motion prediction models (GMPM) is crucial in assessing seismic hazard levels in a specific region. This study investigates the influence of ground motion model (GMM) selection on seismic hazard assessment in Johannesburg, South Africa. Most seismic events affecting Johannesburg are related to mining activities and water ingress into the underground mine voids. The most seismically active regions in South Africa include areas around Klerksdorp, Welkom, and Carletonville. These regions experience a significant number of earthquakes, with mining activities being a major factor contributing to seismic events. Ten ground motion prediction equations were selected from regions of similar tectonics setting and were assessed to see if they correspond to the available local data. Six of these were also compared to the ground motions produced by the earthquake of 14 August 2014, which occurred in the Orkney area, one of the mining areas in South Africa. Sensitivity analyses were also carried out to assess the influence of ground motion prediction equations (GMPE) on the seismic hazard of Johannesburg. The results indicate that the choice of GMPM significantly affects the estimation of seismic hazard levels in Johannesburg. Different GMM can lead to varying ground shaking intensity predictions, impacting risk assessments and engineering design considerations. The findings highlight the importance of carefully selecting appropriate GMM tailored to Johannesburg’s specific geological and tectonic characteristics. Understanding how different models perform in this region is essential for accurate seismic risk mitigation strategies.

Jan 12, 2025

By S. Plaatjie, L. Linzer, M. S. D. Manzi, M. Sihoyiya

Fiber-Optic Sensing and UAV-Platform Techniques for Innovative Mineral Exploration (FUTURE) is a joint European Union and South African consortium project funded under ERA-NET Cofund ERA-MIN3. The project aims to develop cost-effective and novel seismic methods to image geological structures ahead of the mine face, ultimately leading to improved decision-making and safer mining practices. In this study we present results from numerical simulations conducted at South Deep Gold Mine for active in-mine seismic surveys. The study uses a finite differencing code known as WAVE3D, which models models the propagation of a seismic wave generated by stress pertubations that are applied to a defined grid. The numerical simulation model included two major ore bodies (Black Reef and Ventersdorp Contact Reef) and mine tunnels. The model was constrained using laboratory physical property measurements from underground borehole samples such as seismic velocities (e.g., compressional and shear waves) and bulk densities. Two case studies were created, one modelled the active shots and receivers on surface, and the other modelled the shots at depth (in a mine tunnel at ~4000 m depth) with receivers on surface. In both cases, the objective was to investigate the wavefield propagation through the rockmass between the mine tunnel and surface. The simulated wavefield produced by the source located along the mine tunnel was successfully recorded on surface. The surface-recorded wavefield exhibits clear one-way P-wave arrival times at approximately 450 ms to 500 ms, which can be used to calculate the surface-tunnel P-wave refraction tomogram. The results from the study were used to constrain the design, and for the acquisition and processing of the real surface and in-mine seismic surveys acquired under the FUTURE project.

Jan 12, 2025

By P. Lenegan, R. Masethe, R. I. L. Ferreira

Underground mining in high rock-stress environments can induce unanticipated dynamic rockmass deformation, posing a serious risk to the workforce. Data collected by seismic networks are imperative for the design-as-you-mine process. Often an inappropriate assessment of the seismic hazard is done when unnecessary data are included in analytical techniques, particularly short-term time frames. Depending on the methodology or application of data analysis and interpretation, some of the information may not be relevant and even detract from the objective. By using the spatial event clusters associated with active workplaces to define the shape and size of the polygons for data selection for subsequent analysis, a problematic subjective component can be eliminated. The mining of tabular reefs in South Africa, be they gold- or platinum group metals, is very extensive, and underground workplaces are often widespread. Production stoppages are very costly, especially for marginal mines, so there is just cause not to evacuate an entire operation’s workforce following the occurrence of a large mining-induced seismic event. A rationale based on the damaging peak particle velocities (PPV) of the past seismicity can exclude the unaffected underground personnel and allow a quick orderly withdrawal where required. Micro-deformation may precede rockbursting, which can occur on-shift close to the workface. A mine-wide array of geophones suitable for the accurate location and analysis of seismicity may not be able to detect these higher-frequency emissions. Inexpensive but fit-for-purpose accelerometers, providing real-time, continuous measurement of acoustic rockmass deformation near the workplace, integrated through underground and surface communications networks, are shown to record ‘rock-talk’, clearly distinguishable from the instrument background noise, which can be in the order of hours leading up to an imminent damaging seismic event.

Jan 12, 2025

By T. Matlou, S. Myendeki, V. Midzi, T. Zulu, H. Mabotja, T. Mulabisana, B. S. Zulu, L. Hlatshwayo, T. Makhateng

The Council for Geoscience (CGS) performed a macroseismic investigation following the occurrence of an earthquake in the Boksburg area on 11 June 2023. The event had a magnitude of ML 4.4 and a focal mechanism that showed the event as having strike-slip motion along a northwest-southeast oriented fault. The event resulted in aftershocks, some of which occurred away from the causative fault in the east in an area without a mapped fault. However, the cluster of past events in that area suggests that there is an active fault. Due to the size of the main earthquake, a macroseismic survey was conducted to determine its impact on the community in the region. Analysis of observed data obtained from questionnaires resulted in the creation of 153 intensity data points. Their distribution showed high intensity values close to the earthquake source area. However, the investigation highlighted a number of anomalies, such as the higher than expected intensity values observed far from the event epicentre. It could be that these anomalies are caused by site conditions, for example, amplification of ground motion due to thick soil layers.

Jan 12, 2025

By S. Durapraj, R. T. Masethe, U. Maqina, A. C. Adoko

This study investigates seismicity and rockburst control within the Merensky Reef and UG2 orebody of an intermediate-depth platinum mine located in South Africa's Bushveld Complex. Utilising a mine-wide seismic monitoring system, over 1,900 seismic events were recorded between April 2023 and May 2024, revealing complex seismic patterns driven by mining-induced stress changes. The largest event, with a magnitude of ML 2.0, underscores the seismic risks associated with mining at increasing depths. Key findings include the spatial and temporal variability of seismic events, with high seismic activity observed near pillars and geological structures. A significant portion of events correlated with pillar and foundation failures. Moment tensor analysis identified dominant failure modes, such as stress-driven fracture closures and dynamic pillar failures. To mitigate risks, the study emphasises advanced support system redesign and rockburst risk assessment using the rockburst damage potential (RDP) method. The research highlights the importance of monitoring, geotechnical compliance, and adaptive designs in enhancing safety and reducing seismic hazards in underground mining operations. The study applied the RDP framework, classifying damage across R1 to R5 scales, where most instances fell into low-risk categories (R1 and R3), indicating minor to moderate damage. However, data limitations and the tabular stoping method's complexity led to inconclusive results, highlighting the need for more robust data and customised risk assessment tools. Despite these challenges, the RDP approach provided valuable insights into rockburst mechanisms and the effectiveness of existing support systems.

Jan 12, 2025

By J. Petersen, S. W. Robertson

Heap leaching is used extensively for the processing of low-grade copper, gold, and uranium ores. The hydrodynamics is often the most important factor that limits the rate of metal recovery. Over the past decades the author has performed physical and hydraulic testing of a large number of ore samples considered for heap leaching. The data were reviewed to investigate the effect of physical properties such as particle size distribution (PSD) and bulk density on the hydrology. The unsaturated flow through the bed is traditionally modelled with soil science models such as Brooks-Corey (BC) and Van Genuchten-Mualem (VGM), which relate the hydraulic conductivity to the degree of saturation, also known as the hydraulic conductivity function (HCF). However, most authors found a disconnect between the saturated hydraulic conductivity (KS) and the rest of the HCF, which is explained as resulting from a transition from capillary-to gravity-controlled flow. However, the author found no such transition, and found the flow to be governed by capillary flow throughout. Traditional BC or VGM fits were found to be appropriate, however it was necessary to modify the model to include a point of inflection, above which the hydraulic conductivity remains constant while the moisture content keeps on increasing. It was proposed that this is the result of unconnected or poorly connected pores, which fill with solution but do not create additional flow channels. The point of discontinuity was also found to correspond with the air entry point, measured independently from the air conductivity curves. On average, the air entry point corresponds to the experiential rule that the void saturation should remain below 65% for aerated heaps.

Dec 11, 2024

By J. Petersen, A. E. Lewis, J. Chivavava

Rare earth elements (REEs) have special properties that prompted their extensive use in high-tech applications. The demand for REE has, therefore, increased over the past years, resulting in the supply risk of the materials. Extraction from non-conventional low-grade ores, like easy-to-mine ion adsorption clays (IACs), by hydrometallurgical methods is being explored to supplement the supply of REEs. The extraction of REEs from low tenor IAC leach solutions requires enrichment using cost-effective methods such as precipitation. Therefore, the aim of this study was to understand the recovery of REEs from IAC leach solutions using various reagents for precipitation. In the first part of the study, recoveries of REEs from IAC leach solutions through direct precipitation using H2C2O4 and NH4HCO3, as well as antisolvent crystallization using C2H5OH, were predicted via thermodynamic modelling using OLI Stream Analyzer. Verification experiments were then conducted through the use of an agitated reactor. Simulation results showed that high yields of REEs were possible using each of the reagents, but large quantities of C2H5OH were required. Experimental results confirmed the high yields predicted from simulations. Precipitation with H2C2O4 was found to selectively recover REEs while rejecting Al. However, utilization of NH4HCO3 and C2H5OH resulted in the co-precipitation of Al. Furthermore, the yield of Y from antisolvent crystallization was significantly lower than the theoretical value. The product recovered from antisolvent crystallization consisted of well faceted, large crystalline particles, which filtered faster than the fine, amorphous REE carbonates recovered using NH4HCO3. REE oxalates were crystalline, larger than the REE carbonates and, overall, showed the highest filtration rate.

Dec 11, 2024