Search Documents

By J. D. Steenkamp

In 2002 the commissioning of an ilmenite smelter on the North Coast of South Africa was extended by three months due to the failure and subsequent replacement of the magnesia-based refractory lining. The lining failed due to the hydration of magnesia caused by an unexpected source of water. The incident resulted in significant financial losses and a prolonged insurance claim which was settled in 2009. As magnesia-based refractories are used extensively in both ferrous and non-ferrous applications, the authors of the paper want to share the experience gained from this incident with others. The paper reviews the literature available on furnace start-up practices and explains the hydration of magnesia using available sources. The incident is studied in more detail, both technically and economically, and the costs incurred are quantified in terms of the cost of the original lining. The paper concludes with lessons learned and recommendations made for future work. The intention of the paper is to stimulate open debate regarding best practices in preheating of furnaces lined primarily with magnesia. Keywords Commissioning, hydration, start-up, magnesia.

Jan 1, 2011

By P. van der Westhuyzen

Gone are the days of stock-standard diamond beneficiation and final recovery circuits that could, with minor modifications, be adapted to treat any diamond-bearing ore source. A new era has dawned, offering the opportunity to streamline existing diamond processing operations or develop simpler, more efficient, and economical diamond processing plants, with the focus on more efficient comminution, beneficiation, and final recovery. This change in scene has been brought about by: ? The introduction, rapid development, and maturation of multiple comminution, sorting, and recovery technologies ? The need to adapt to a new standard of project approach post the commodity super-cycle phase, where optimizing existing operations and developing scalable, ?fit for purpose? new mines are fast becoming the norm in diamond processing plants in both primary and alluvial operations ? The quest for energy efficiency and lower labour costs ? More remote and inaccessible reserves (under lakes, tops of mountains etc.). This paper serves to identify some technology advances and demonstrates how these could be considered as replacements for or in combination with conventional technologies to arrive at an optimum techno-economic solution. To name a few applications/technologies: ? Comminution: conventional cone crusher, modified/specialized cone crusher, and the high-pressure grinding roll (HPGR). Although significant advances have been made in recent years, this paper only briefly covers comminution within the beneficiation circuit design ? Waste sorting: NIR (near-infrared) sorting, optical (colour) sorting, XRF (X-ray fluorescence) ? Primary concentration by combining dense media separation (DMS) with either XRT (X-ray transmissive), pulsed X-ray, jigs, or pan plants depending on the application, scale, and economics ? Final recovery of diamonds with either conventional X-ray technology, pulsed X-ray technology, or XRT. From recent studies, it can be concluded that there is no longer a standard solution, but rather the ?right? or appropriate solution. Through combining a sound knowledge of the ore source (ore dressing studies or onmine data gathering) and leveraging off the advances in technologies, one is able, through trade-off studies, to arrive at the ultimate techno-economic configuration, the ?right? solution. Emphasis is placed on maximizing diamond recoveries through appropriate technology selection and minimizing the associated costs in an effort to de-bottleneck or improve efficiencies of existing diamond processing plants, or to arrive at the ideal new diamond process plant design..

Jan 1, 2014

By S. W. Jacobsz, Y. Narainsamy

Several high-impact tailings dam failures around the world in recent years have placed a renewed focus on the stability of tailings dams and pointed to potential shortcomings in traditional drained design and safety evaluation procedures. A need to consider undrained shear strength in the design of tailings facilities has become apparent. However, there are specific requirements that need to be met before undrained shearing occurs. In South Africa, the last major failure was likely the Merriespruit disaster in 1994, which leads to the question of whether the conditions required for undrained shearing are readily applicable to South African tailings dams. This paper describes research at the University of Pretoria that has recently commenced to further investigate the conditions required for undrained failure to occur. The research includes laboratory and field testing to replicate these conditions in the laboratory and relate them to those found in an active tailings dam.

Jun 2, 2022

By G. Mahieu, L. Perez-Barnuevo, P. Barbeau-Verreault, O. Gravel

Raglan mine is a copper-nickel mining operation, part of Glencore, and located at the northern border of the province of Quebec, in Canada. The concentrator is fed by a mix of four different underground mines, thus it must accommodate periodic variations in run-of-mine mineralisation and composition. In order to concentrate nickel and copper sulphides, the flotation circuit also includes a regrind step at the cleaner stage. This leads to challenges in maintaining optimal recoveries and in diagnosing potential disturbances. This paper advocates for the combination of automated mineralogical analyses, reconciliated mass balances from diagnostic sampling campaigns and innovative multidimensional visualisation to assist in tackling these challenges. Through this approach, it has been possible to quantify the regrind step influence on the liberation of valuable mineral particles in the final flotation concentrate. The methodology also allows for virtual sieving of the products in order to track unit loss in the flotation process. The proposed analysis scheme presents an opportunity to easily diagnose flotation circuit performance based on mineralogical information. Keywords: Flotation diagnostics, multi-dimensional visualisation, mineral liberation analysis

Jan 1, 2020

By M. C. Bekker, T. Govindasamy

The potential benefits of implementing the Last Planner System (LPS) on mining infrastructure projects were explored through two case studies in the South African mining industry. The first part of the study involved measurement of the impact of LPS through earned value techniques on secondary project data. The second part sought to establish the lean construction success factors leading to improvements in construction performance. This was done by means of a research questionnaire distributed to the project owner’s teams, contracted parties, and the LPS facilitators. The results revealed a positive correlation between LPS application and planned percentage complete. However, the performance achieved during the LPS pilot was not sustained. The success factors and benefits considered most evident in the case study pilot differed among the three stakeholder groups. The findings are expected to guide construction stakeholders to better define performance measures and focus on factors required to make LPS implementation more effective in the South African mining sector.

Feb 12, 2024

1. New books Extractionmetallurgy'81. London, Institution of Mining and Metallurgy, 1981. 441 pp. $40. Structures and properties of engineering alloys, by W. F. Smith. London, McGraw-Hill, 1981. 512 pp. R89,65. Rock and mineral analysis, by W. M. Johnson and J. A. Maxwell (editors). 2nd ed. New York, Wiley, 1981. 424 pp. $40. Mining methods and equipment, by Mining Information Services. London, McGraw-Hill, 1981. 218 pp. International directory of acronyms in library, information and computer sciences, by P. M. Vaillancourt. New York, Bowker Publishing Co., 1980.516 pp. 2. New journals Trends in analytical chemistry (TRAC), Amsterdam, Elsevier, 12 issues per year. Subs. for Volume I (16 issues) D.H 95.50. Engineering management international, Amsterdam, Elsevier, 4 issues per year. V.S. $32.50, D.H. 70.00. 3. Mineral policy sector publications Mineral Bulletin MR 188. Vanadium: an imported mineral commodity. $3.60. Mineral Bulletin MR 190. Canadian mines: 1980 perspective. $3.00. Canadian mineral survey 1980, by the staff of Mineral Policy and Energy Sectors. $1.00. 4. NIM reports Report 2090 Observations made during the dig-out of a 48 MV A ferrochromium furnace. Report 2110 The investigation of a fracture on one side of a skip bridle in a mine shaft. Report 2111 Pilot-plant tests for the recovery of chromite and the noble metals from the UG-2 Reef at Maandagshoek. Report 2112 The behaviour of UG-2 chromite concentrates during smeUing. Report 2113 M uUivariable control of a milling circuit at East Driefontein Gold Mine. Report 2115 Estimation of material variation in powdered materials Report 2116 Optimum operation of a direct-reading spectrometer with excitation by a 5 kW inductively coupled plasma torch. Report 2121 Five robust indicators of central value. Report 2124 The direct determination, by differential pulse anodicstripping voUammetry at the thin mercury-film electrode,of cadmium, lead, and copper. Report 2127 The preparation of a reference material of South African zirconium concentrate from Richards Bay. Report 2128 The determination, by differential pulse anodic-stripping voltammetry at the thin mercury-film electrode, of cadmium and thallium in six NIMROC reference materials Report 2130 The organizing of conferences. Research grant The Board of Directors of International Precious Metals Institute in conjunction with Gemini Industries, Inc., a California-based refiner of precious metals, recently announced a Gemini Research Grant of $3000 will be awarded annually to a graduate student who is pursuing studies in metallurgy, chemical processes as it pertains to precious metals, or materials sciences, and who is performing or planning research studies in the precious metals. Accidents and diseases The Third International Colloquium for the Prevention of Occupational Accidents and Diseases in the Iron and Metal Industry is to be held on 15th and 16th June, 1982, in Palma de Mallorca, Spain.

Jan 1, 1981

By I. C. R. Gilchrist

A principal requirement in the design of a backfill transportation pipeline is the accurate prediction of the values of key slurry flow parameters. Input variables, such as particle size and solids concentration, are usually determined by the load-bearing requirements of the placed fill. It is therefore necessary to identify the optimum conditions of transport for a particular combination of input variables. This paper critically reviews the methods available to predict the values of flow parameters in slurry pipeline systems in the light of experimental data generated at a test facility. Particular reference is made to rod milled waste and deslimed tailings transport systems and, in the case of the former, it is shown that both pumping and load-bearing requirements can be optimized.

Jan 1, 1988

By R. W. O. Kersten

"The aim was to investigate stable pillar design procedures currently in use on chrome and platinum mines, subject them to a critical appraisal, and review some of the alternatives proposed in the literature. It is concluded that these pillar design methodologies suffer from drawbacks that can be detrimental to the mining industry and that an improvement is essential.It was decided that the increased availability of analytical models and failure criteria could be exploited to update the empirical to an analytical/empirical method. The proposed methodology adopted FLAC2D simulations using the Hoek-Brown failure criterion to calculate full stress-deformation curves for typical pillars. The mine stiffness concept was introduced to determine the pillar load, integrating the influence of the pillar and strata stiffness and mining geometry.The proposed methodology was calibrated by comparing experimental and predicted deformations and failure depth of pillars in one hard-rock mine. The conclusion was that the methodology is an improvement over the one currently in use. An improved factor of safety is given by the intersection point of curves for the stiffness of the system and the pillar strength. IntroductionThe paper describes part of the findings of a PhD thesis (Kersten, 2016) on an alternative stable pillar design methodology. This publication contains sufficient information to implement the method for further interpretation and testing. A detailed comparison of measured and predicted data is dealt with in the thesis, and the results show that there is promise in the methodology.The design of stable pillars in mining is of fundamental importance, not only in the Bushveld Complex mines but throughout the mining industry. Wherever mining occurs, pillars are formed at some stage and it is essential to predict their behaviour, i.e. whether they will burst, yield, or remain stable. Although of crucial importance, pillar design still suffers from major shortcomings and weaknesses."

May 1, 2019

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

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

Jan 1, 2020

By T. M. Louw, L. Auret, S. M. Bradshaw, F. F. Noelle

Modern industrial processes undergo frequent process state shifts due to changes in demand, feedstock or maintenance operations. Multimodal continuous process monitoring is therefore required for effective plant operation. Steady state analysis is a key component in effective multimodal monitoring. Various issues still remain in determining the starting and ending times of transition states. High dimensionality of process data is a characteristic that makes it difficult to select variables to consider when conducting steady state analysis. In order to address these problems, steady state detection on the principal component space is proposed. First, data undergoes principle component analysis to reduce dimensionality, then existing multivariate steady state detection algorithms are applied on the principle component scores. A novel algorithm, threshold PCA-SSD, is proposed and compared to an existing technique for steady state detection. The algorithms are investigated on high dimensional synthetic data Keywords: Multimode process, steady state detection, principle component analysis

Jan 1, 2020

By G. M. Evans, S. Mitra, G. J. Jameson, M. M. Hoque, E. Doroodchi, R. Moreno-Atanasio

In this study, the detachment behaviour of 314 μm diameter silica particles from a single 3 mm diameter bubble rotating in a confined cavity connected to a cross-flowing water channel was investigated. The liquid cross-flow generated a rotating flow inside the confined cavity, which was quantified using 2D particle image velocimetry (PIV) measurements. The fluctuating (x-y) velocity components inside the cavity for each cross-flow velocity, U, were then used to estimate the corresponding time/volume average turbulence intensity values, Ti. Turbulence intensity ranged from 9.3-16.2%. In the absence of a bubble, the rotating vortex within the cavity created a pressure field which decreased towards the centre of the vortex; wherein the lowest pressure zone was located at the centre of the cavity. The vorticity into the cavity was found to increase with increasing the turbulence intensity. A bubble, both without and with particles attached, was introduced into the cavity under different flow conditions and its motion was tracked using the shadowgraphy technique. Experimentally, it was found that the bubble-particle aggregate was trapped at the centre of the cavity (low-pressure region) following a spiral trajectory. An interface resolved 3D computational fluid dynamics (CFD) model was developed to predict the particle detachment conditions within the confined cavity. CFDpredicted vorticity in the absence of a bubble agreed well with experimental observation. The CFD analysis indicated that the vorticity inside the cavity was not significantly changed when a single bubble was introduced. On this basis, the computed vorticity was utilised in describing the force balance applied to the attached particle. Computationally, it was found that the critical turbulence intensity, beyond which the particle was detached from the bubble, was 14%. This value was consistent with the experimental observation of detachment at Ti at 15.7%, and above, but not at lower values. Keywords: Bubble, particle, turbulence intensity, vorticity, particle detachment, CFD, PIV

Jan 1, 2020

By E. K. Mutinda, D. K. Maina, R. M. Kasomo, B. O. Alunda

Assessment of blast fragment size distribution is critical in mining operations because it is the initial step towards mineral extraction. Different empirical models and techniques are available for predicting and investigating the consequences of blasting, one of which is the Kuznetsov Cunningham-Ouchterlony (KCO) model. In this paper we summarize the advances in the empirical models from inception until now, and explore the improvements that have been made so far with particular emphasis is on the most recent KCO model. Utilization of the model and the errors that arise between expected and the actual outcomes are analysed. The results indicate that the KCO model remains useful for predicting the blast fragmentation at limestone mine sites, despite the availability of other advanced prediction models. It is also a valuable instrument for pre-surveying the impact of varying certain parameters of a blast plan.

Mar 1, 2021

By Pentti Niskanen

"The ventilation system of the Pyhasalmi Mine is described. The mine is located in an area where the annual atmospheric temperature may vary from +30°C to -30°C. Natural ventilation causes serious difficulties because the underground mine is connected by drifts and finger raises to the bottom of the 120 m-deep open pit.The natural ventilation tends to make the ventilation net unstable. The paper describes how the instability characteristics of the ventilation net can be studied and how they have been taken into account when automating the control of the ventilation. In addition, the software, hardware and costs of the system are described. The system consists of 12 automated ventilation doors, two mine fans, 8 elements for sensing the air flow and temperatures, a control panel and a process computer."

Jan 1, 2014

By D. A. Hayman, N. Rambuda, N. A. Sithole, C. Hockaday, J. D. Steenkamp

"Mogale Alloys operates two alternating current (AC) submerged arc furnaces (SAFs), and two direct current (DC) open bath furnaces. These furnaces are capable of producing two products, silicomanganese (SiMn) and ferrochrome (FeCr), with a total production capacity of 110 000 t/a. This paper focuses on the silicomanganese process specific to Mogale Alloys, with a two-furnace operation facility that produces 55 000 t of SiMn per annum.IntroductionMogale Alloys is situated in Krugersdorp in the West Rand District Municipality in South Africa’s Gauteng Province. The Krugersdorp facilities were initiated by Palmiet Chrome Corporation in February 1963 to produce charge chrome (Basson, Curr, and Gericke, 2007). Mogale Alloys bought the Palmiet Chrome facilities in 2005 from Samancor Chrome. The Afarak Group currently owns 90% of Mogale Alloys.SiMn is produced in two 20 MVA submerged arc furnaces (SAFs), which were originally built for chromite processing, and ferrochromium (FeCr) in two DC open bath furnaces rated at 40 MVA and 10 MVA. When this paper was compiled, the 10 MVA furnace was not operational due to the high cost of production, and one of the 20 MVA SAFs was recently converted to produce FeCr due to market constraints. Mogale Alloys has 290 permanent employees with four contracting companies on site. The paper presented here focuses on the production of SiMn at Mogale Alloys.The high-level process flow at Mogale Alloys is summarized in Figure 1, which details the raw material input and the products produced by the two 20 MVA furnaces. SiMn is produced by carbothermic reduction of manganese-bearing raw materials from four different sources, namely Wessels, Mamatwan, and UMK ores, and BHP sinter feed. The carbon sources are bituminous coal from Forzando and Msobo, local suppliers in Mpumalanga."

Nov 1, 2018

By B. Oladipo, E. Govender-Opitz, T. V. Ojumu

Bioleaching, one of the several important processes of biohydrometallurgy is an accepted technique for the extraction of desired metals from sulphide ores. Although the mechanism is well understood, the slow metal extraction rate due to passivating layer and ferric ion precipitation, especially in a typical chalcopyrite leaching, impact negatively on the overall efficiency of the process. Iron precipitates are common in many biohydrometallurgical operations and these precipitates represent the exit route for undesired iron, alkali ions or sulphate ion from the processing circuit. While there have been reports on approaches to slow the formation of these precipitates by manipulation of various operating conditions such as pH, redox potential and temperature, the precipitates become highly significant over the long term in the continuous operation of bioleaching and may render heap impermeable. Although there is evidence supporting the nature of the passivation to be precipitates of iron, their scavenging mode of reducing extraction efficiency is still not well understood. This paper thus discusses iron precipitation in biohydrometallurgical processes, including its formation and significance on bioleach processes. The paper also details the associated passivating species identified with heap bioleaching of sulphide minerals and provides some insight into sorption mechanism that limits metal recovery in a typical bioleach operation. Keywords: Iron precipitation, bioleaching, chalcopyrite, passivation, jarosite, metal sorption

Jan 1, 2020

By M. R. L. Blackwell

Over the period 1964 to 1969 Conzinc Riotinto of Australia undertook and completed an evaluation of a large low-grade copper orebody on Bougainville Island in the Territory of Papua and New Guinea. From the outset it was clear that the geographic, climatic and topographic characteristics of the area would necessitate a heavy capital outlay on production and support facilities. Both the physical problems and the magnitude of the investment amplified the need for a thorough and rigorous evaluation of the various strategies for exploiting the orebody. In view of the number of independent and inter-related variables and the flow of new data, a computer-based evaluation process extended the scope and depth of the analysis considerably. The system used permitted an iterative assessment and re-assessment of the various strategies with the results becoming increasingly refined with time. This process encompassed the modelling of the orebody, the simulation of alternative operating strategies and the determination of the financial profitability of the project,

Jan 1, 2014

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 K. Barani, F. Moradpouri, M. Godarzi

A sampling campaign carried out on the lead flotation circuit at Lakan lead-zinc processing plant demonstrated that the Pb recovery was only about 68-69%. Characterization of the Pb flotation tailing showed that 61.4% of the Pb is in the –37 μm size fraction, which is difficult to recover by flotation. Batch flotation experiments were carried out on a plant sample of the lead flotation feed and high-intensity conditioning (HIC) was used to improve flotation recovery. The results showed that increasing the collector conditioning time from 1 to 3 minutes significantly increases the Pb recovery from 42.5% to 57.3%. Further increases in conditioning time to 5 to 7 minutes decreased the Pb recovery. With a collector conditioning time of more than 3 minutes, due to the surface cleaning effect, increasing collector and frother dosage improved the flotation selectivity and recovery, even at higher collector dosages. The Pb recovery increased from 69.1% to 79.8% with increasing impeller speed from 600 to 1250 r/min. HIC increases the flotation kinetics and greatly improves the flotation recovery and concentrate grade.

Apr 2, 2022

By K. Esmeli, A. Ozkan

The enhancement of shear flocculation of a galena suspension by ultrasound in the presence of sodium isopropyl xanthate and potassium ethyl xanthate was investigated by studying the effects of surfactant concentration on flocculation efficiency, contact angle, and zeta potential. Strong flocculation was obtained under alkaline conditions (pH > 9) and a high ultrasound power level (150 W). Although the ultrasonic treatment increased the negativity of the zeta potential of the galena particles, this did not decrease the flocculation efficiency, thereby showing that the hydrophobic interactions between the particle were stronger than the electrical double-layer repulsion resulting from xanthate adsorption. The findings indicate that the ultrasonic treatment promoted the adsorption of surfactants onto the mineral surfaces.

Jul 5, 2023

By M. A. El-Sayed

"Double oxide films (bifilms) are significant defects in light alloy castings that occur when, during melt handling, the oxidized surface of the melt is folded over onto itself and then entrained into the liquid metal, encapsulating a portion of the local atmosphere. These defects create crevices in the final castings and have been shown to decrease tensile and fatigue properties, and also to increase their scatter, making casting properties unreproducible and unreliable. This paper reports on two different experiments performed to study the behaviour of such defects in Al alloys. In the first experiment and as a proxy for a bifilm, air bubbles were held in an Al-5Mg alloy melt for varying periods of time before the melt was allowed to solidify. The composition of the gas content of the bubbles was subsequently analysed using mass spectrometry to ascertain how it changed with time. In the other experiment the influence of the age of bifilms on the mechanical properties of Al-5Mg castings was assessed by holding castings in the liquid state for different times before solidification, followed by mechanical property testing and SEM and EDX studies of the fracture surfaces of tensile test bars. This research could provide a better understanding of the behaviour of double oxide film defects, and may lead to the development of new techniques by which these defects could be eliminated, or at least deactivated, in light alloy castings. Introduction One of the major structural defects affecting the reproducibility of mechanical properties of aluminium castings is the double oxide film defect (El-Sayed and Essa, 2017; El-Sayed and Ghazy, 2017). This type of defect is usually created as a result of poor filling system design. Owing to the very high reactivity of liquid Al-Mg alloy with oxygen in the air, the surface of the liquid metal can be considered to be permanently coated with an oxide film (Weigel and Fromm, 1990). As a consequence of surface turbulence, the oxidized metal surface can become folded over onto itself, but not fused, and then submerged into the bulk liquid, with a portion of the mould atmosphere being entrapped within it (Campbell, 2003; Al- Sayed, 2016; El-Sayed, Hassanin, and Essa, 2016)."

Nov 1, 2018