Search Documents

By M. Verri

The paper describes the design options available for the simultaneous production of liquid sulphur dioxide, via cryogenic condensation, and sulphuric acid in a sulphur-burning sulphuric acid plant. The impact of the operating conditions of the cryogenic condensation on plant capital and operating costs are discussed, with the main focus on the optimization of the most important economic and operative drivers of the non-ferrous mining industry, such as energy efficiency, reliability, and availability. A case history relevant to an integrated sulphurburning plant producing liquid sulphur dioxide and sulphuric acid, recently started up in the Democratic Republic of Congo, is described. Keywords sulphur burning, sulphuric acid, liquid SO2, cryogenic SO2 condensation, liquid sulphur dioxide, sulphur dioxide.

Aug 1, 2013

By J. R. de Wet

The Skorpion Zinc Project demonstrated the viability of the production of zinc from non-sulphide (commonly called ?zinc-oxide?) ores, using the leach/solvent extraction/electrowinning process route. The success of Skorpion raised the profile of a number of other zinc oxide deposits, considered to be similar. However, a study of three such deposits demonstrates that the application of whole-ore-leaching may be limited and that the concentration of zinc, or at least the separation from acid-consuming gangue minerals, prior to leaching is essential to the successful recovery of zinc from these deposits. The recovery of zinc is dependent on an understanding of the geology of the deposit, the mineralogical association of zinc minerals with gangue minerals (Boni, 2005) and limitations of mining practice. Flowsheet development and the associated metallurgical testwork to derive an economically viable treatment route from two of these deposits are discussed. Key parameters affecting zinc recovery by gravity separation and flotation are presented, demonstrating that even minor differences in mineralogy have a profound impact on zinc recovery and reagent consumption. The study also deals with the determination of appropriate design parameters through metallurgical mapping of the dense media separation process.

Jan 1, 2008

By C. -J. Rick

The ferroalloy market is characterized by a very long communication distance between the consumer and the producer. This has created a situation where it has been difficult to change ferroalloy properties. In this paper different FeNi and FeCr-compositions are discussed and their merits for stainless steel-makers in different situations are compared using the process models applied at Outokumpus Avesta works and Acerinox's Columbus operation Keywords Ferroalloy design, ferrochrome, ferronickel, stainless steel making, UTCAS software.

Jan 1, 2010

By J. A. L. Napier

The evaluation of stress distributions and deformations in the vicinity of tabular mine layouts has been carried out routinely for many years using the displacement discontinuity boundary element method (DDM). In this approach, mine layouts are approximated as irregular shaped planar cracks (or slits) where the ?width? of the crack, corresponding to the excavation height, is assumed to be negligible compared to the in-plane, lateral dimensions. Early implementations of this approach have employed square element tessellations with constant displacement discontinuity component values. This is generally restrictive when considering the analysis of irregular pillar layouts or the intersection of fault planes with excavations. This paper describes a general revision of this approach in which triangular or quadrilateral element shapes are introduced in conjunction with higher order variations of the displacement discontinuity shape functions. This facilitates an accurate evaluation of detailed stress and displacement components close to excavation surfaces and allows the investigation of local damage and failure behaviour associated with particular geological structures such as weak parting bands or weak layers. A method to embed these features in an efficient and accurate multi-level solution scheme is described. This formulation can be extended to allow the analysis of shallow-depth pillar layout problems, facilitating the assessment of the interactive response and stability of shallow depth pillar layouts. A simple deformation model is introduced to facilitate the assessment of pillar failure and bursting potential in shallow mine layouts. The computational techniques described here have been implemented in a computer code called TEXAN (Tabular EXcavation ANalyser) and are illustrated by the simulation of a particular multiple reef pillar mining problem in a shallow depth platinum mine.

Jan 1, 2007

By H. Gumede

The challenges of noise emission in the mining industry are well known. These challenges are exacerbated by the increasing trend towards mechanized mining, and there is general consensus that an effective industry-wide Buy Quiet initiative (BQI) will significantly assist in reducing noise-induced hearing loss by minimizing noise at its source. In its endeavours to improve the health of workers, the MOSH Noise Team has facilitated development of a proposed strategy and structure for a BQI. Extensive consultation with stakeholders resulted in support for the initiative. The MOSH Noise Team was advised to solicit the contributions of consulting mechanical and electrical engineers (CM&EEs), whose support is viewed as critical for realization of the initiative, since elimination of noise at source is an engineering challenge, and most procurement standards and specifications are set and managed by engineers. A workshop with CM&EEs and their representatives, as well as group environmental engineers (GEEs), confirmed the urgent need for the initiative and provided valuable insights. The workshop used obstacle based planning techniques to identify the factors on which the success of this BQI will depend. Most potential obstacles to the implementation revolved around proper scoping, leadership commitment, and companies? adherence to the initiative. Lack of involvement by the relevant stakeholders and economic issues were also identified as potential obstacles. Delegates at the workshop felt that the successful implementation of the industry-wide BQI, will depend on the mining industry addressing the specific obstacles. The most critical and urgent issue was identified as the formation of a properly constituted industry-wide BQI task team that includes representatives from critical stakeholder groups and which is mandated and empowered to drive the initiative.

Jan 1, 2014

Hydrometallurgical processing of complex ores and concentrates is becoming increasingly important as the mining and metallurgical industry seeks to exploit mineral deposits that are difficult to treat using conventional mineral processing and pyrometallurgical technologies. Mineral processing is often challenged by the difficulty or inability to separate valuable minerals into marketable concentrates. Hydrometallurgical processing, using selective leaching technology, can often ?chemically beneficiate? such difficult deposits. Pyrometallurgical treatment of base metal concentrates is capital intensive and subject to ever more stringent environmental control. Hydrometallurgical processing is generally lower in capital cost for an equivalent metal production rate and avoids the gas and dust issues associated with pyrometallurgical processing. The possibility of by-product recovery may also increase with hydrometallurgical treatment. Two examples of new technology or flowsheet development for treatment of complex ores and concentrates are used as illustration. These include the El Boleo process of Baja Mining for recovery of copper, cobalt, zinc and manganese from a complex clayey ore and the PLATSOL? process for recovery of copper, nickel, cobalt, platinum, palladium and gold from a bulk sulphide concentrate.

Jan 1, 2009

By G. M. Dunn

Pressure leach autoclave circuits are employed in the leaching of ores, concentrates, mattes, alloys and intermediates for the recovery of metals into solution. Once the metals are extracted into solution, the value metals can be recovered by hydrometallurgical means such as by purification followed by electrowinning, hydrogen reduction, pyrohydrolysis, crystallization and other unit operations. In many of these integrated flowsheets the pressure leach step is pivotal to the recovery of the value metals from the host material. Furthermore, the autoclave circuit invariably is a high capital cost component of the plant and an area that is carefully scrutinized when debottlenecking and capacity increases are being considered. This paper identifies a unique proven way of increasing the capacity of existing or new exothermic pressure oxidative leach autoclave circuits by as much as two or three times.

Jan 1, 2009

By H. Moslemi

Metal sulphides are usually semiconductor and their electrical conductivity is low; therefore they are sources of electrons and promote electrochemical reactions. This phenomenon is of significant importance in sulphide flotation. The formation of a hydrophobic surface on the sulphide minerals is the result of various electrochemical reactions. Potential difference in the mineral-liquid interface plays an important role in the kinetics of the reactions, and consequently, on the flotation behaviour of sulphides. The value of the potential difference is influenced by different oxidation-reduction reactions at the interface area, and is measured by a two-electrode system in which an electrode is applied as reference electrode. In this study, we attempted to measure open-circuit potentials of pyrite, pyrrhotite, and steel electrodes with respect to the reference electrode. To verify the results, two reference electrodes, saturated calomel and gold, were used. Time-potential relations with respect to both reference electrodes show that the potential of both pyrite and pyrrhotite increases with time; however, this trend decreases for steel except when at pH 12. The results show a potential drop by an increase in pH value, except for steel, for which the opposite holds. This paper discusses galvanic interactions between minerals and grinding media in the grinding circuits, as well as among other minerals in the flotation processes.

Jan 1, 2012

By C. Chibwe, M. Tadie

Copper electrowinning models are useful for prediction of process performance, optimisation and improved understanding of complex phenomena such as current distribution without the expense of physical experiments. Recently, finite element models such as those generated in COMSOL Multiphysics have been developed for the electrowinning process. The predictive capacity of these models is improved by incorporating physical and chemical phenomena. Since physicochemical properties of electrolytes can have a significant effect on the quality of copper produced, this work was conducted to investigate the effect of such properties on electrodeposition at varying levels of copper, sulphuric acid concentration, iron concentration and temperature, in the presence of a smoothing agent. Tests were conducted at conditions similar to those for industrial copper electrowinning. Mathematical correlations were developed for the following properties; density, conductivity and diffusivity which were compared to existing models in literature and thereafter applied in a COMSOL generated model. The model predicts, current distribution at the cathode surface, deposit thickness (morphology), concentration of species within the cell and current efficiency. Current distribution output shows the impact of several phenomena on the deposit generated including the relationship between increasing copper concentration and current density and lack of influence of acid concentration on this property. Keywords: Electrowinning, physicochemical properties, electrolyte composition, modelling

Jan 1, 2020

By P. T. L. Koh, M. P. Schwarz

A series of parametric studies were carried out for a flotation cell using computational fluid dynamic (CFD) modelling to determine the local response of hydrodynamic performance measures to change in operating conditions. As an example, an Outotec cell was used to carry out the response analysis. The hydrodynamics of the flotation cell has been previously investigated by laboratory measurements presented by Schwarz et al (2019), and a CFD model of the cell was developed and validated using these measurement data (Schwarz et al., 2019). In this work, parametric studies have been performed by modelling the flotation cell at different stirring speeds and gas flowrates, varying above and below standard conditions. Predictions of void fraction, bubble size, bubble surface area flux and particle–bubble collision rate are presented showing the effects of increasing and decreasing stirring speed and gas flowrate around the base case. Generally, the total gas holdup in the flotation cell was found to increase with increasing gas flowrate. The bubble diameter was found to decrease with increasing stirring speed and increase with increasing gas flowrate. The collision rate increased with both stirring speed and gas flowrate. The parametric results show that there is a correlation between the bubble surface area flux Sb which is affected by stirring speed and gas flowrate, and the particle–bubble collision rate. With the collision rate as an indication of flotation collection, it is also possible to show that cell performance can improve more by increasing stirring speed than by increasing gas flowrate in a situation where everything else remains constant. The CFD results can form the basis for a Response Surface Method (RSM) analysis which can then be used for automated optimization of cell performance. Keywords: Flotation kinetics, computational fluid dynamics, bubble surface area flux, particle-bubble collision rate

Jan 1, 2020

By A. S. Macfarlane

Western Deep Levels (North) is a deeplevel, longwall gold mine. Severe strata control problems and high rates of elastic closure, resulting in high rates of energy release, are associated with mining at depths of up to 3 500 metres. To assist in controlling these aspects of mine design and safety, a waste backfilling system has been designed and commissioned. The integration of such a system into an existing, producing longwall brings many problems of logistics, organization, planning and control. It is the purpose of this paper to examine these aspects of fill management, and to describe some of the practical problems involved with the integration of this project into the mining environment. Although many of the rock mechanics and soil mechanics aspects of backfilling systems are common there are specific problems involved which apply to longwall systems. Very strict guidelines are required for planning personnel to design their backfilling systems for a particular situation. These aspects are examined, and basic models developed which can be used without the necessity of comprehensive input from rock mechanics. A summary of the financial implications is included, and these are examined with regard to possible future applications and expansion of waste filling.

Jan 1, 1988

With the advent of the First World Conference on Sampling and Blending (WCSB1) in Denmark in 2003 and the subsequent biannual meetings, together with other sampling conferences held on a regular basis, the central tenets of sound sampling practice and theory are spreading throughout the minerals industry. End-users and project houses, or groups who execute the projects, are all now better-informed about the principles of ?correct? sampling and preparation. To support this, many large mining companies now have their own in-house champion who will advise and ensure that the supply of sampling equipment is correct in terms of design and application, which further contributes to improved metallurgical accounting. The presence of a sampling champion also ensures that consistent and acceptable standards in regard to sampling practice are maintained throughout the group. This has, in turn, placed a greater degree of responsibility on the shoulders of sampling equipment suppliers, in terms of their equipment designs and product offerings, as well as applications knowledge for effective installation of equipment on client plants. This paper attempts to highlight the significant challenges that sampling equipment manufacturers often face with regard to complying with the newly understood requirements of best practice and sampling equipment design. Manufacturers invariably face issues with end-users or project houses, including project cost control, restricted headroom for samplers in plant layouts, and a lack of knowledge in certain quarters. These all have a potentially negative effect on the final outcome of a proposed sampling solution or implemented project. In reality, what is intended by all parties and what is ultimately installed are at times very different, with undesirable consequences for all concerned. Against this challenging background, there are, however, a number of instances where intent and implementation are complementary, with very satisfactory results. Keywords sampling, mechanical samplers, design, installation, sampling theory.

Jan 1, 2014

By Q. G. Reynolds

The presence of dust and fume suspended in the freeboard region is a common feature of the operation of direct current (DC) plasma smelting furnaces. This occurs primarily as a result of the use of fine feed materials together with the open-arc, open-bath operation of such smelters, and is exacerbated by the high velocities and turbulent mixing of the gas in the vicinity of the arc jet. Dust and fume losses into the furnace off-gas system can be significant in some cases and may have economic, operational, and environmental impacts on the process. A computational modelling study is presented in which the concentration of dust material was considered as a continuous field subject to a governing partial differential equation. Settling behaviour was calculated as a function of particle size, local gas/plasma temperature, and other physical properties. Development of the coupling between the concentration field and a magneto-hydrodynamic description of the arc is shown, and the resulting models were used to compute various aspects of the behaviour of the concentration field in the arc region for a variety of furnace conditions. Time-averaged as well as transient models of the arc were used to generate the results presented. Qualitative case studies produced several practical suggestions for furnace operation, including increased dust capture by the bath when feed ports are located closer to the electrode, and the possible effects of feed segregation in the furnace freeboard based on dust particle size and density.

Jan 1, 2014

By D. Vogt

South African underground hard rock mines are typically managed using measurements made daily, weekly or even monthly of key parameters such as face advance, readiness to blast, blast success, temperature or dust levels. Safety and health can be greatly improved if a real-time measurement system can inform decision making. To enable this vision of widespread sensing, communication and decision support, CSIR has developed an open standard architecture for communication of sensor data, and a reference implementation using that standard. The standard is called AziSA, which means ?to inform? in isiZulu. AziSA provides an architecture that allows for connection of any type of sensor and that is particularly suited to wireless sensing. On the safety side, the CSIR is developing a suite of sensors specifically around the risk of rockfalls. Sensors such as an electronic replacement for the barring tool, or a thermal sensor that can detect loose rock, can already be used to identify potential hazards in the hangingwall. When these sensors are combined with location and time information in a single database, it becomes possible to build maps of risk and to extrapolate risk into unmined areas. It is also easy to confirm that routine safety procedures like barring are actually taking place. In the future, it is proposed that entry inspections will be made using remote techniques that will not place an individual miner at risk. While miners will still be required to bar or place support, they will do so with a clear idea of which portions of the hangingwall are safe, and which present risk. On the health side, if a network such as AziSA is deployed in the stope, it becomes much easier to monitor the condition of individual miners. The work strain experienced by miners can be monitored through heart-rate sensors, and their core body temperature can be monitored using various novel techniques discussed in the paper. Additional sensors can be added to quantify worker exposure to hazards such as noise or dust. Individual sensors allow for people who are at an unacceptably high risk of developing heat disorders to be treated timeously, and they also allow for team management based on objective measures. Over a longer period, position sensing and environment monitors, or personal dosimetry, provide a powerful management tool to prevent workers from being overexposed to hazards, and to confirm that management instructions are being followed. Cost-effective sensing in the stope is a major challenge, but one that can be overcome through technology like that used in AziSA. Sensing provides opportunities both to make immediate interventions when workers are exposed to unacceptable risk, and to manage the long-term exposure of workers to hazard.

Jan 1, 2010

The 114th Annual General Meeting of The Southern African Institute of Mining and Metallurgy was held at The Country Club, Johannesburg, Napier Road, Auckland Park on Thursday, 11 August 2011. Welcome The President extended a special welcome to the guests and representatives of our sister institutes and other associations, and also to recipients of awards, senior members of industry, Honorary Life Fellows, past-presidents, our members and other guests, among them the following: F. Cawood, President, Institute of Mine Surveyors of South Africa C. Roode, President, Mine Metallurgical Managers' Association of South Africa M. Cresswell, Chairman, South African Coal Processing Society A. Hay, Immediate Past President, South African Institute of Electrical Engineers L. Gardner, President, South African National Institute of Rock Engineers J. Krynauw, President, The Geological Society of South Africa (GSSA)

Jan 1, 2011

By W. J. Kempson

As a result of the rising electrical energy costs in South Africa, a method was sought to reduce the overall electrical consumption of typical shaft systems. A typical shaft configuration was analysed and the primary energy consumers were identified. The ventilation fans for this system were found to consume 15% of the total energy of the shaft system. It was calculated that more than 50% of this energy is consumed by the shaft itself, more specifically, by the pressure losses that occur in the shaft as the ventilation air passes through it. In order to ensure that the theory being used for the evaluation of these shaft systems is accurate, a total of five shafts were instrumented and the actual pressure losses over the shafts plotted against time. These shafts were then analysed from a theoretical perspective. Finally, in order to ensure a thorough understanding of the behaviour of the ventilation air in shaft systems, the systems were simulated using computational fluid dynamics (CFD) techniques. There were significant discrepancies between the theoretical analysis and the CFD simulation during the initial comparisons. This discrepancy lessened as the complexity of the CFD models increased, until when the complete shaft was modelled using the full bunton sets, the pipes, and the flanges, the difference between the theoretical evaluation and the CFD simulation was small. This result demonstrates that the theory is insufficient and that the interrelated effect of the buntons and fittings has not been fully appreciated by current theory. The final phase of the work presented here was to evaluate the costeffectiveness of using different bunton shapes and shaft configurations. It is shown that the increase in the pressure losses and therefore the direct operating costs of the shaft can vary by as much as 80%, depending on the bunton configuration chosen. The placement of the piping in the shaft can increase the pressure losses, and therefore the direct operating costs of the shaft, by as much as 12%, depending on the placement of the piping in the shaft; this effect includes the use of flanges. The use of fairings on a large cage can reduce the resistance that the cage offers to the ventilation flow by as much as 30%. This, however, does not translate into a direct saving because as the cage moves through the shaft, the overall effect is transitory. These savings can be significant when the items highlighted in this work are applied correctly. Keywords shaft design, ventilation, energy consumption, modelling, CFD.

Dec 1, 2013

By L. Cortez, Á. Magalhães, F. Durão

Methods based on regionalized variables theory, commonly known as geostatistics, have been widely used since more than three decades to predict/estimate spatial distributed data, namely mineral resources, with quite good results. However, some restrictive theoretical assumptions, such as the need of some kind of stationarity of data and the linearity of the predictor/estimator, are often wrong options in many applications. An alternative approach based on a finite mixture modelling, or cluster weighted modelling (CWM) technique, is presented and compared, in its theoretical and practical aspects, with the geostatistical methodology. CWM is a general non-linear framework based on the estimation of joint density functions allowing the calculation of the conditional probability function of the target values given their locations. Validation tests are performed in a 2D case study, where testing data proceeds from a soil contaminated by a heavy metal. The results of ordinary kriging and CWM estimations show a very good agreement with actual data. Keywords: Non-linear regression, Geostatistics, Ordinary kriging, Finite mixture modelling, Cluster weighted modelling.

Jan 1, 2003

By C. A. Preece

Finsch Diamond Mine, part of De Beers Consolidated Mines (DBCM), is situated between Kimberley and Postmasburg in the Northern Cape, South Africa. The kimberlite pipe was discovered in 1960 and De Beers assumed full mining rights three years later. Production commenced in 1966 from an opencast mine, and the operation moved underground in 1990 by developing and implementing a new mining method of blast hole open stoping. A mechanized block cave was brought into production in 2005. While block caving was not new to DBCM, with productive caves having operated in Kimberley and Cullinan, from the outset, the Finsch caving operation was designed as the most advanced with a semi-automated ground handling system. The extraction level was laid out as a single sided herringbone system due to the very poor ground conditions on the south side of the pipe. These factors combined to create a unique environment, which necessitated a completely new approach to draw control. The required system was developed by a team consisting of DBCM personnel from head office and Finsch, GEMCOM and Sandvik-Tamrock, the equipment manufacturers. It was essential that the system could be run quickly, provide daily updates and issue calls directly to the LHDs. It was also important that the base data were obtained from the resource model, which was to be depleted on a monthly basis, thus forming a close link between mine planning and draw control. A further requirement of the system was that it would form a central database for the extraction operation where all data pertaining to delays and changes in conditions, be they production or ground condition related, could be recorded and analysed. The design brought together empirical rules and tables developed at Cullinan plus concepts from programming developed for the front cave at Koffiefontein so that the new system was based on sound practical experience. These facets were incorporated into PCBC (GEMCOM) and CMS and PCS (Sandvik-Tamrock programmes). The system has now been operational for just under a year and while modifications have had to be made, it has proved itself as a practical and reliable tool with which to manage the draw control at Finsch. The system is unforgiving and has forced the production team to focus on the specified drawpoints. Problems have been encountered but most of these have been related to re-educating a workforce, familiar with a flexible mining method, to the strict discipline required for a modern block cave.

Jan 1, 2007

By G. C. S. Burns

"A system has been developed to simulate the mining of the Mufulira orebody. The main output from the system is a picture-gram from which the start date and effective life of all proposed mining faces can be measured. This information provides schedule dates for the completion of major capital development projects and stope preparation activities.The orebody must be divided horizontally into main levels and vertically into blocks of equal interval. Ore reserves are defined by specifying the tonnage contained in each block and main level. Mining faces are defined by their co-ordinates, and may themselves be split into sublevels. The simulation maintains an echelon between contiguous sublevels. Other face parameters control the speed of mining, the availability and relative priority of the face.The required production from each crusher/shaft system is specified. The technique of simulating the mining process in order to achieve the required production is described."

Jan 1, 2014

By J-N. Oyourou, R. McCrindle, C. J. S. Fourie

"Many mines in South Africa were inappropriately closed and left unattended during the last century. These old mines are potential sources of environmental pollution and may pose a health risk to local populations, since the surrounds contain elevated levels of toxic elements. Soil from the old Edendale lead mine property in Mamelodi Extension 11, east of Pretoria, which was suspected to be contaminated with lead, was analysed. The mine, which primarily extracted galena, was decommissioned in 1938. During this study, the old mine property and surrounds were screened for zinc and lead using a field-portable X-ray fluorescence (pXRF) spectrometer. The metal concentrations, which ranged from 18 to 7 300 mg/kg for zinc and from 50 to 21 000 mg/kg, for lead, were plotted together with the GPS coordinates to produce a concentration distribution map for lead and zinc. This map revealed three anomalies of unacceptably high concentrations, possibly corresponding to areas where ore had been piled or dumped to waste from the old mine. Thirty-six of the soil samples that were analysed by pXRF on site were selected and analysed in the laboratory using inductively coupled plasmaoptical emission spectrometry (ICP-OES) for verification purposes. The zinc and lead concentrations obtained from the two analytical methods were compared using linear regression analysis. Excellent correlation was obtained between the results from the two methods (r2 = 0.99 for zinc and r2 = 0.95 for lead). IntroductionEarly economic development in South Africa was based mainly on mining. However, these activities left a legacy of approximately 6 000 improperly closed mines scattered throughout the country (SAI, 2009; van Schie, 2012). Due to their proximity to settlements, these abandoned mines may pose a potential health risk to humans and an environmental hazard. Mining activities generate large amounts of waste, mostly in the form of overburden and process tailings, which often contain high concentrations of toxic elements. Bioavailable toxic elements can contaminate crops growing in the polluted area (Kala, 2010; Marwa, Meharg, and Rice, 2012) or eventually migrate and detrimentally affect the groundwater (Rösner and van Schalkwyk, 2000). Moreover, sulfides accumulated in mine waste generate acid minewater (Rösner and van Schalkwyk 2000). Unfortunately, owners of mines, closed before the promulgation of the Minerals Act more than twenty years ago (Act 50 of 1991), avoided the compulsory rehabilitation of the environment enforced by this legislation. Toxic elements, such as lead, which were confined within ores, may become exposed through mining and are therefore more easily mobilized within the environment. Elements essential to human and plant metabolism, such as zinc, when exposed by mining activities, may become hazardous due to elevated concentrations (Aslibekan and Moles, 2003; Wuana and Okiemen, 2011). Humans exposed to potentially toxic metals may experience harmful effects as discussed elsewhere (Nriagu and Pacynar, 1988; Wayne and Ming-Ho, 2005)."

Jan 1, 2019