Search Documents

By A. Lebedev

Application of dynamic simulation for a base metal processing plant is discussed. The plant receives ore of variable quality from two pits and blending had to be applied to maximise recovery. The plant employs multiple stage crushing, screening, dense media separation, two-stage milling, flotation and further filtration. Stockpiles and bunkers are used in the plant to cater for breakdowns and operating delays, the latter being appropriately incorporated in the model. A detail dynamic simulation model was constructed from the rock face to the filtration plant to assist the plant design, notably with definition of the ore blending method, sizing storage capacity and in-plant conveyors. The plant model was broken down to capital equipment items and various feed types and mass balances were tested to ensure that the production target is met. The dynamic simulation tested the proposed blending method, identified bottlenecks, produced risk profiles for all strategic storage facilities and overall helped minimise the project risks and improve confidence in the plant design. Keywords: Dynamic simulation, ore grade, blending, capacity, risk profile, throughput, utilisation

Jan 1, 2020

By N. M. Skarbøvig

Safety nets have been used selectively in several South African mines over the years. Recently, however, an increased emphasis on mine safety has seen a resurgence of their utilization underground. This renewed interest has, however, also raised questions about the design and performance of these nets and their practical applications underground; these are the subjects of this paper. Generally, mine safety nets are used to protect personnel from falling rock in temporarily supported development ends, tunnels, ASGs and stopes. These diverse applications each require specific net designs, resulting in a wide product range. As a rule, however, mine safety nets need to carry specified static or dynamic loads with the least possible deformation. Because of the application-specific characteristics of mine safety nets, there is no South African national standard for them. Consequently, neither the SABS, nor the CSIR, performs standard tests on these components. In 2007, therefore, SIMRAC made available their Savuka facility for the testing of mine safety nets; the Savuka drop test became the accepted industry test standard for mine safety nets. From late 2007 to mid 2009, Norsenet, in collaboration with SRK Consulting, the designated operators of the site, conducted approximately 200 tests, which formed the core of Norsenet's database of safety net performance under various loading conditions. The test specimens comprised a wide range of net materials and config¬urations; the materials included: high density polyethylene, polypropylene, polyester, nylon, as well as a blend of polyester and polypropylene used primarily for reinforcement ropes. Net configurations encompassed the following parameters: net type, i.e., knotted nets or sewn nets; net strand type, i.e., braids, twines, or webbing; strand size; mesh size; mesh orientation, i.e., parallel or oblique to the net edge; net edge reinforcement, i.e., knotted, tied, with or without reinforcement ropes; internal reinforcement ropes, i.e., diagonal ropes and/or orthogonal ropes; as well as net attachment options, including loops and chainlinks. From the number of variables listed above, and the limited number of tests available, it can be appreciated that the findings of the test series required careful scrutiny. Nevertheless, principal performance factors clearly emerged: besides the obvious net mesh variables, such as net strand material, strand diameter and mesh size, it was the configuration of the reinforcement ropes and their attachment to the net panel which had the strongest influence on net performance. All the structural elements of a safety net need to work together to transfer the imposed load from the net panel, through the reinforcement ropes, to the net attachment points. In this load transference chain, the weakest link limits the strength of the structure. From this development work, Norsenet has gained valuable expertise to custom-design and manufacture nets to clients' specifications. The examples presented in the paper will also show the experience acquired in the installation of nets underground. Future work will include the establishment of a state-of-the-art net testing facility at Norsenet's factory in Vredenburg.. Keywords Mine safety net, safety net, protective device, mine safety, safety net design, safety net performance, safety net applications.

Jan 1, 2011

By D. M. Walters

To achieve safe and healthy working conditions in shaft sinking operations, ventilation, and where necessary cooling, must be provided. These requirements are considered for sinking vertical, sub-vertical and decline shafts. Regulatory requirements, determining air volumes and amount of cooling required, the various ventilation methods used such as force, exhaust, twin columns, air reversal using four gate systems and the equipment required to achieve this, are discussed. Consideration is also given to the sinking of single and twin declines using diesel equipment, calculating the air volumes required to dilute pollutants and the various methods for achieving this. Also discussed is the importance of having early ventilation design input to the overall design of the mine. In many cases ventilation can be the prime determinant for the size and configuration of the shafts or declines that have to be sunk.

Jan 1, 2009

Current mine production planning, scheduling, and allocation of resources are based on mathematical programming models. In practice, the optimized solution cannot be attained without examining all possible combinations and permutations of the extraction sequence. Operations research methods have limited applications in large-scale surface mining operations because the number of variables becomes too large. The primary objective of this study is to develop and implement a hybrid simulation framework for the open pit scheduling problem. The paper investigates the dynamics of open pit geometry and the subsequent material movement as a continuous system described by time-dependent differential equations. The continuous open pit simulator (COPS) implemented in MATLAB, based on modified elliptical frustum is used to model the evolution of open pit geometry in time and space. Discrete open pit simulator (DOPS) mimics the periodic expansion of the open pit layouts. Function approximation of the discrete simulated push-backs provides the means to convert the set of partial differential equations (PDEs), capturing the dynamics of open pit layouts, to a system of ordinary differential equations (ODEs). Numerical integration with the Runge-Kutta scheme yields the trajectory of the pit geometry over time with the respective volume of materials and the net present value (NPV) of the mining operation. A case study of an iron ore mine with 114 000 blocks was carried out to verify and validate the model. The optimized pit limit was designed using Lerchs-Grossman?s algorithm. The best-case annual schedule, generated by the shells node in Whittle Four-X yielded an NPV of $449 million over a 21-year mine life at a discount rate of 10% per annum. DOPS best scenario out of 2 500 simulation iterations resulted in an NPV of $443 million and COPS yielded an NPV of $440 million over the same time span. The hybrid simulation model is the basis for future research using reinforcement learning based on goal-directed intelligent agents.

Jan 1, 2008

By D. F. Malan

This paper gives an overview of the difficulties associated with determining the strength of hard-rock pillars. Although a number of pillar design tools are available, pillar collapses still occur. Recent examples of large-scale pillar collapses in South Africa suggest that these were caused by weak partings that traversed the pillars. Currently two different methods are used to determine the strength of pillars, namely, empirical equations derived from back analyses of failed and stable cases, and numerical modelling tools using appropriate failure criteria. The paper illustrates that both techniques have their limitations and additional work is required to obtain a better understanding of pillar strength. Empirical methods based on observations of pillar behaviour in a given geotechnical setting are popular and easy to use, but care should be exercised that the results are not inappropriately extrapolated beyond the environment in which they are established. An example is the Hedley and Grant formula (derived for the Canadian uranium mines), which has been used for many years in the South African platinum and chrome mines (albeit with some adaptation of the K-value). Very few collapses have been reported in South Africa for layouts designed using this formula, suggesting that in some cases it might yield estimates of pillar strength that are too conservative. As an alternative, some engineers strongly advocate the use of numerical techniques to determine pillar strength. A close examination unfortunately reveals that these techniques also rely on many assumptions. An area where numerical modelling is invaluable, however, is in determining pillar stresses accurately and for studying specific pillar failure mechanisms, such as the influence of weak partings on pillar strength. In conclusion, it appears that neither empirical techniques nor numerical modelling can be used solely to provide a solid basis for conducting pillar design. It is therefore recommended that both these techniques should be utilized to obtain the best possible insight into a given design problem. Owing to the uncertainties regarding pillar strength and loading stiffness, monitoring in trial mining sections and in established mining areas is also an essential tool to test the stability of pillar layouts in particular geotechnical areas.

Jan 1, 2011

The effect of microwave radiation on the processing of Palabora open pit ore has been investigated. Significant reductions in Bond work index are demonstrated for microwave treated samples. Reductions in work index were found to be sensitive to applied power level and method of application. After single stage magnetic separation tests, improvements in recovery of FeO were noted. However, these were curtailed by oxidation of the magnetite species as a result of continued heating. The effects of microwave radiation upon the recovery of copper after froth flotation were also quantified. Small initial increases in recovery of copper were note dafter microwave treatment. However, with prolonged exposure, recovery was found to decrease. Overall recovery of copper was not affected except for samples treated for extended duration. The influence of microwave treatment upon the conventional Palabora grinding circuit has been quantified with the use of USIMPAC software. Conclusions are made regarding the economics of microwave treatment for this ore. Keywords : Microwave; Grinding; Magnetic Separation; Flotation; Economic Analysis

Jan 1, 2000

Stope support systems, typically consisting of props and packs, are used extensively in the gold and platinum mining industry to stabilize the rock mass in the excavation vicinity and to reduce the hazard associated with rockfalls and rockbursts. The design of stope support systems was historically based predominantly on past experience and practices, and cost considerations. The recent development of a new and improved support design methodology, has led to the potential for significant increases in worker safety and support cost savings. This paper reviews some of the fundamental rock mass and support criteria that form the basis of the improved support design methodology. The considered design criteria are the height of potential fall; quasi-static stope closure rates; dynamic stope closure rates; compressive hangingwall stresses; discontinuity spacing, orientation and interface properties; effect of support length (stoping width); effect of compression rate; consistency of support performance; a real coverage; support spacing; and zone of support influence. The proposed support design methodology combines both the zone of influence and key block stability theories, resulting in optimized support systems for rockfall and rockburst conditions. The site-specific methodology consists of two stages: (i) a tributary area analysis, and (ii) a zone of support influence and a stability analysis, considering hangingwall failure due to buckling, shear and block rotation, which gives maximum safe spacing of individual support units. The methodology is suitable for designing support systems for both blocky hangingwalls and hangingwalls fragmented by face-parallel extension and shear fractures.

Jan 1, 2001

By L. Andrews

Electron microbeam analysis of metallurgical samples covers abroad field and has applications in the areas of minerals processing, hydrometallurgy, pyrometallurgy, physical metallurgy and corrosion. Some idea of the possibilities of microbeam techniques can be illustrated by examining the following example of a run of mine gold ore sample. Milling and bulk chemical analysis of the ore will provide major and minor element analysis as well as the gold grade. More detailed chemistry can also reveal the presence of trace elements that may be advantageous or deleterious to downstream processing. Electron microbeam analysis can provide additional information in the form of size, liberation, association and composition of the gold, the percentage of refractory gold, the size, liberation, association and composition of enclosing sulphide or arsenide phases and the levels and association of the trace elements that may affect processing. In most cases electron microbeam analysis complements bulk chemistry, and does not replace it. Both techniques should be used in parallel for maximum effect. This paper touches on the origin and history of electron microbeam techniques, as well as the theory behind microbeam analysis and processing. The more common types of instruments and their applications will also be dealt with briefly.

Jan 1, 2007

By J. C. van Dyk

The objectives of this study were firstly to characterize and quantify the inorganic material in the Secunda gasification feedstock. Secondly, the effect of washing on coal properties and sulphur content was investigated. Coal from a single source and a coal blend from 6 different mines were tested. Composite samples were visually sorted by a geologist in the groups containing roof, floor, C-shale, pyrite, carbonates and siltstone material. Inorganic mineral types from the single source and coal blends were characterized separately. Characterization results indicate that the average mineral matter distribution (mass %) of the coal tested was as follows: carbonates (6%), carbon-shale (19%), pyrite (12%), roof (18%), floor (46%) and siltstone (22%). The roof, floor, pyrite, carbon-shale, carbonates and siltstone in the single source and coal blend have similar chemical properties and characteristics, except for the ash melting temperatures of pyrite and kaolinite in the siltstone. By implementing selective mining and increasing mining management practices, the percentage discard in the run of mine coal can be significantly reduced, by up to 45%, if the mining of roof and floor material can be avoided. No significant impact of mining practices on sulphur content is expected in this case due to the low concentrations of sulphur in the roof and floor material. In summary, selective mining practices will probably result in almost45% decrease in the discard content in the run of mine to gasification, an increase in the carbon content from 55.0% to 58% (air dry basis) and a decrease in sulphur content from 1.1% to 0.9% (air dry basis). It is clear from this study that the main source of sulphur is not in the high ash-content mineral types (e.g. floor and roof), but in situ in the coal structure. Keywords: inorganic material, effect of washing, roof, floor, pyrite, C-shale, carbonates and siltstone

Jan 1, 2005

By A. D. Engelbrecht

An investigation was undertaken to ascertain the suitability of four selected low grade, South African coals for gasification in a bubbling fluidized bed for production of synthesis gas and for the development of integrated gasification combined cycles (IGCC). This study consisted of the characterization of the coals, laboratory evaluation of intrinsic reactivity, and experimentation with a pilot-plant fluidized bed gasifier. Results of the characterization experiments show that the selected coals are high in ash, rich in inertinites, very dense (low porosity), with low caking indices and high ash fusion temperatures. Reactivity measurements with a thermogravimetric analyser (TGA), under reaction rate controlling conditions with carbon dioxide, showed that the reactivity of the coal chars decreases with increase in rank of the coal, as expressed by the vitrinite random reflectance and carbon content of the parent coals. The fixed carbon conversion achieved in the fluidized bed gasifier also correlates well with the rank parameter of the coal, with higher conversions being obtained with the lower rank coals. Thermal shattering and attrition of the coal particles produce significant amounts of fines, which correlate with the grindability indices, and no agglomeration (non-caking) was observed. It is concluded that fluidized bed gasifiers are able to utilize typical low grade, high ash South African coals for synthesis gas production and for inclusion in integrated gasification combined cycles for power generation.

Jan 1, 2010

By Y. M. Zhao

The separating performance of a pilot dry beneficiation system using a gas-solid fluidized bed was investigated with coal from South Africa. The coal used for the study has a low inherent moisture content, a moderate ash content, a moderate volatile content, a low sulphur content and a high calorific value. Its washability is moderate. Experimental results show that the separating quality of the fluidized bed drops gradually as the feed-coal particle size decreases. The probable error, E, values for 50?25mm, 25?13 mm and 13?6 mm coals were 0.04 g/cm3, 0.06 g/cm3 and 0.09 g/cm3, respectively. The cause of the differences in separating characteristics was analysed by particle dynamics and numerical modelling. Furthermore, the pilot fluidized bed was employed to separate 50-6mm coal. In this experiment the coal ash content was reduced from 23.74% to 11.79%, with a probable error, E, value of 0.07 g/cm3and a recovery efficiency of 98.26%. This indicates that the fluidized bed is applicable to the separation of coal from South Africa and has good separating performance and wide applicability. Keywords: feed-coal particle size; gas-solid fluidized bed; dry coal beneficiation; numerical simulation.

Jan 1, 2010

Abrasivity and wear resistance are not intrinsic properties of materials, but rather parameters describing the interaction between the wear part and the material exerting the wearing action. Moreover, wear or, at least, that part of the wear which is due to the cumulative effect of minute scratches, is linked to the hardness disparity concept. Abrasivity, in this context, is a relative concept, pertaining to a metal-rock pair, rather than to a rock. Based on this assumption, a methodology has been developed and checked against actual cases, enabling description of the abrasiveness of a rock or the abrasion resistance of a metal through the ratio of comparable quantitative indicators. The paper describes the methods and provides examples of applications to practical cases (namely: rock cutting, excavation by means of tunnel boring machines (TBMs), rock crushing, and hydro-mill cutting).

Jan 1, 2012

By R. Rampersadh

Paper written on project work carried out in partial fulfilment of B Eng (Metallurgical Engineering) degree Steel was sampled after calcium treatment, to follow changes in inclusion composition during the processing of alumium-killed steel into continuously-cast billets, at Mittal Vereeniging. There is some shift in the distribution of inclusion compositions towards alumina, as the liquid steel moves from the ladle to the tundish. This shift may reflect some reoxidation. However, the effect on inclusion modification is small, since only a small fraction of the tundish oxide inclusions is fully solid?which agrees with the plant observation that these heats cast without clogging. Magnesium oxide in the inclusions contributes to liquid formation, substantially increasing the range of lime: alumina ratios over which the inclusions are partially or fully liquid.

Jan 1, 2006

Abrasivity and wear resistance are not intrinsic properties of materials, but rather parameters describing the interaction between the wear part and the material exerting the wearing action. Moreover, wear or, at least, that part of the wear which is due to the cumulative effect of minute scratches, is linked to the hardness disparity concept. Abrasivity, in this context, is a relative concept, pertaining to a metal-rock pair, rather than to a rock. Based on this assumption, a methodology has been developed and checked against actual cases, enabling description of the abrasiveness of arock or the abrasion resistance of a metal through the ratio of comparable quantitative indicators. The paper describes the methods and provides examples of applications to practical cases(namely: rock cutting, excavation by means of tunnel boring machines (TBMs), rock crushing, and hydro-mill cutting).

Jan 1, 2012

By O. Nyandoro, K. Khoza-Shangase, M. Madahana, N. Moroe

"Occupational noise-induced hearing loss (ONIHL) is classified as the leading work-related disability in the mining industry. ONIHL has a negative impact, on not only the health and occupational productivity of affected individuals, but arguably also on the country’s mining industry and economy. Hearing conservation programmes (HCPs) are an effective strategy in the management of ONIHL. However, current literature indicates that HCPs are not achieving the anticipated and desired outcomes in the South African mining sector despite the efforts focused on the management of ONIHL. Therefore, the purpose of this study is to propose the use of a feedback-based noise monitoring model as a tool for monitoring and managing ONIHL in the South Africa’s mining sector. This model is a basic static feedback model with practical applications such as estimating, monitoring. and providing quantitative information to aid miners, mining administrators, and policy-makers in decision-making around HCPs. Additionally, the model could form part of an early intervention and management strategy for ONIHL in the workplace. The strength of this model, although currently static, is that it encompasses all the pillars of HCPs and takes into account the policies concerned with the management of ONIHL in the mining sector. IntroductionHazardous noise levels (above 85 dBA) in the workplace place exposed employees at risk of developing a disabling occupational noise-induced hearing loss (ONIHL). ONIHL is the ‘partial or complete hearing loss in one or both ears as a result of one’s employment’ (Nandi and Dhatrak, 2008, p. 1). This type of hearing loss develops gradually as a result of being exposed to continuous or intermittent high levels of noise over a long period of time (McBride, 2004; Patel et al., 2010; Rappaport and Provencal, 2001). A hearing threshold below 40 dBs is classified as a disabling hearing loss (Yadav et al., 2015), resulting in a hidden condition that does not readily manifest itself (Tye-Murray, 2009, p. 3). Dugan (2003) describes disabling hearing loss as the ‘most prevalent, least recognized and least understood physical disability’, while Copley and Frederichs (2010) and Hermanus (2007) argue that permanent disabling hearing loss is a major contributor to the global burden of disease on individuals, families, communities, and countries."

Jul 1, 2019

By C. G. Anderson

"As demand for copper and precious metals continues to grow, the need to develop new supplies increases. In the imminent future, the new large global ore bodies are more complex, with enargite, tennantite and tetrahedrite as the abundant contained mineralogy. These resources also contain considerable gold and silver. As such, direct treatment by conventional smelting is prohibitive due to the potential environmental impacts and copper quality issues. This paper outlines the effective application and associated economics of industrially proven Alkaline Sulfide Leaching (ASL) hydrometallurgical technology to an enargite concentrate for selective removal and recovery of arsenic, antimony, gold and silver, leaving a clean high-grade copper concentrate. The fundamentals behind the technology are described, along with those of hydrometallurgical arsenic fixation and stabilization. Flowsheets and unit operations are elucidated, along with an in-depth economic analysis. INTRODUCTION Hydrometallurgical methods can be employed for treatment of gold-, arsenic-, antimony-, tin- and mercury-containing materials, concentrates and ores, as well as complex ones containing any number of metals (Conner & Anderson, 2013). The alkaline sulfide leaching (ASL) system is one of these. It is essentially a mixture of sodium sulfide and sodium hydroxide with other meta-stable alkaline species as required. This is a unique hydrometallurgical system beause it is a very selective lixiviant for the distinct leaching of tin, gold, antimony, arsenic, tellurium, silver and mercury (Anderson, 2000, 2001a, 2001b, 2002, 2003, 2012, 2014; Anderson & Twidwell, 2008a, 2008b, 2010; Anderson et al., 2004, 2005; Glazkov & Tseft, 1961; Gnatyshenko & Polyvyanni, 1961; Nadkarni et al., 1975; Nadkarni & Kusik, 1988; Polyvyanni, 1963). Worldwide, ASL technology has been employed industrially in the former CIS, China and the United States for the production of antimony (Anderson et al., 1992, 1994; Anderson & Krys, 1993; Holmes, 1944; Kaloc, 1967; Nordwick & Anderson, 1993; van Stein, 1971)."

Jan 1, 2015

By M. D. G. Salamon

It has been noted in recent years that the pillars in some bord-and-pillar workings scale or spall. In due course this process may result in the collapse of the affected panels. This paper presents a methodology for assessing the long-term behaviour of such workings, employing a simple model of pillar deterioration by scaling. It is assumed that scaling causes a gradual reduction in pillar width, hence the safety factor of pillars reduces in time. It is also postulated that the process of scaling will be arrested if the coal rubble around the pillar reaches a critical height. The probability that the pillars will survive for a certain period of time, or indefinitely, can be estimated. This can be accomplished by using the log-normal distribution proposed by Salamon and Munro4todescribe the distribution of the apparent safety factor values at failure. If an estimate of the rate of scaling is known, it is also possible to quantify, using a Monte Carlo technique, the expected life of the pillars. The method also facilitates pillar design by specifying a probability of survival for a given number of years or, alternatively, for an indefinite period. Using the limited pool of available data, analyses show that the proposed model does not appear to be inconsistent with observed behaviour of pillars in the Vaal Basin. The crucial parameter is the rate of scaling. Back calculations from eleven cases suggest that the rate of spalling in the Vaal Basin is about 0.2 m/year. It is noteworthy that two further cases in the same area give significantly higher rates (0.85 m/year). Some underground observations in the Witbank area indicate considerably lower rates in the order of 0.02?0.04 m/year. Further areas for study are recommended. These include a study of old workings to determine the extent and rate of scaling in the various coal fields, an investigation into the possible relationship between scaling and the chemical composition of the coal. In addition, model studies will be required to assess the practical impact of the research findings.

Jan 1, 1998

Zincor operates four Lurgi fluidized-bed roasters at its roast-leach-electrowin refinery in Springs, South Africa. In 2006, Zincor?s four roasters had reached ages ranging between 26 and 38 years, while the expected life of this type of roaster is generally 20 to 25 years. Roaster 3 was the first to be rebuilt in 2006 because its condition had deteriorated the most. Roaster 4 is being rebuilt in 2007 and the last two roasters will be rebuilt over the next few years as necessary. The main damage to Roaster 3 was caused by buckling of the roaster steel shell. This caused the freeboard and dome areas to shift, resulting in cracks in the dome brick lining and crumbling of the bricks lining the freeboard. The steel shell buckled owing to the combined effects of excessively high operating temperatures and corrosion. These were caused by the following two main factors: ?Thicker external insulation was installed during a previous maintenance shut, which resulted in an increased steel shell temperature and exacerbated the normal metal fatigue ?SO2/SO3 containing gas that had permeated the brick lining passed through its dew point during frequent planned and unplanned shuts over the life of the roaster. This caused acid formation between the refractory brick lining and steel shell, resulting in corrosion of the shell. The initial plan was to replace only the damaged areas of Roaster 3, but when it was dismantled more damage was found, and the roaster was rebuilding from the tuyeres upwards. During there building the opportunity was taken to modify the design of some of the brickwork to extend its life and make future repairs easier. This paper discusses the mode of failure of Roaster 3 and the lessons learnt during its dismantling and rebuilding.

Jan 1, 2007

By M. Sanchez, C. Morales, J. Palacios

The slag produced in the Peirce Smith converters at the Ilo smelter of Southern Peru Copper Corporation has shown a wide range of magnetite content, which introduces an additional problem in the smelting units due to materials? recirculation. However, the operation data analysis and the application of thermodynamic concepts have made it possible to produce slag with magnetite content between 18 and 20%. The results have suggested that matte grade enhances magnetite formation while silica controls its precipitation, and they are also influenced by the increasing temperature of the slag. Thus, in order to obtain the required magnetite content in the slag, matte grade is gradually enriched in equilibrium with slag containing from 25 to 32% silica.

Jan 1, 2004

By W. A. Naismith

On 16 July 2004, at 07:10 a collapse of about 1.8 million m3 (about 4 ½ million tons) took place on the north wall on Nchanga Open Pit. This collapse was the result of a developing instability that was first identified when the benches in this area were exposed during September 2002. The NOP geotechnical and survey departments closely monitored the development of the failure. Control measures were implemented to reduce risk to men and equipment. These measures were continually improved and adjusted in reaction to observations and changing conditions. The purpose of this paper is to describe in detail the monitoring results that were obtained both by traditional methods and the novel slope stability radar system (SSR). Difficulties encountered in establishing critical deformation levels and deformation rates that can be used to predict the onset of collapse are described. Details of the monitoring and control measures, which led to the successful completion of mining in the area, are also described in this paper.

Jan 1, 2005