Search Documents

By Case JE

Complete Bouguer and isostatic anomaly maps compiled for Costa Rica provide new insights into the crustal structure and mineralization. Terrain corrections were approximated from the station to a radial distance of 167 km by computer using a 3- minute digital elevation model of Costa Rica and surrounding areas. Isostatic corrections, based on Airy-Heiskanen compensation, were also calculated from 0 to 167 km using the 3- minute elevation model and from 167 km to 180 degrees using combined topographic and isostatic reduction maps. Costa Rica can be divided into four main terranes based on gravity data: a Pacific coastal region, a Pacific inland area, a central region, and a Caribbean coastal region. Negative gravity anomalies along the Pacific coastal margin correlate with an offshore sedimentary basin, and large positive anomalies correlate with mafic and ultramafic igneous rocks on the Santa Elena, Nicoya, and Osa Peninsulas, and at Puerto Quepos. Pacific inland anomalies are separated from the Pacific coastal region by a steep gravity gradient, and over the Coast Ranges a major gravity low overlies a thick section of low- density sedimentary rocks. In the central region, a major gravity low occurs over a thick low-density Quaternary volcanic complex in the Cordillera Central, north of San Jose. A prominent feature of the Caribbean coastal region is an interrupted gravity low over the Limon Basin that reflects a thickening of low-density sedimentary rocks. A conspicuous gravity lineament, trending northeast separates northern and southern Costa Rica. The lineament is defined by the termination of several major gravity anomalies including those over the Coast Ranges, Cordillera Central, Puerto Quepos, and Limon Basin. In addition, the lineament correlates with other geological and geophysical data including: northeast-trending faults, the southern limit of large Quaternary volcanoes, the segmentation of the middle America volcanic belt, and a seismicity gap. This lineament may be the expression of a previously unrecognized major northeast- trending fracture zone in central Costa Rica. An important economic implication of the gravity lineament is a correlation with the termination of gold and porphyry copper deposits in Costa Rica.

Jan 1, 1987

By L O. Filippov, A Piçarra, I V. Filippova, N P. Lima, K Silva

The reverse cationic flotation using the starch as depressor of iron oxides is the most efficient upgrading process when is applied to the simple systems composed by hematite/magnetite and quartz. The fundamental studies performed with iron ore samples from different deposits over the world (Brazil, Russia, Mexico) confirmed the inefficiency of starch when the silica content in the concentrate is controlled by Fe/Mg-bearing minerals. The main approaches for solving these issues focused on using new collector formulations or new depressants. However, the alternative depressants tested inhibited flotation of quartz or were less efficient as iron oxide depressant. In this work, a comparison between new amidoamine collector with conventionally used etheramine and oleate collectors in the flotation of pure minerals and actual iron ore samples was performed. The electrophoretic mobility, contact angle and infrared spectroscopy measurements allowed to deduce the hindrance effect as main mechanism of selective adsorption of the amidoamine on the quartz surface compared to hematite and kaolinite. Molecular modelling studies revealed that this steric hindrance effect can be explained by the combination of the "shield" formed by the coordinating water molecules on the (001)-D2 hematite surface and the large volume of the amidoamine head group. The amidoamine tested as collector without use of a depressant showed excellent flotation selectivity allowing the total flotation of quartz from a mixture with hematite. High collector performance and selectivity were confirmed during reverse flotation of quartz using continuous flotation column (D = 500 mm; H = 4000 mm) and allowed to produce the iron concentrates with more than 60 per cent of iron from the rejected slimes samples of Vargem Grande 2 processing plant (Iron Quadrangle, MG, Brazil). The new collector allowed to develop flexible flotation flow sheets to process the iron ore with a complex gangue mineralogy and unlock the old altered tailings and slimes product as new resources

Aug 24, 2022

By J Euston

In the 1990s round flotation tanks were introduced as a result of the lower fabrication costs for larger cells. It was soon observed that round tanks held a number of further advantages over the traditional rectangular cells, such as better hydrodynamics, lower costs and reduced footprint per unit volume.In the last twenty years flotation cell sizes have increased almost exponentially and designs of 250 to 300m3 are a well established standard with plans for full scale pilot installations of 500m3. The option to install these large cells gives flexibility in circuit arrangements and allows efficient and economic matching of flotation lines with modern mill designs. The benefits of scale up are well known with reduced numbers of cells required, lower power costs per unit volume, etc. These benefits will be briefly reviewed with specific reference to the installation of 300m3 to 350m3 flotation cells at the Copperton Concentrator of Kennecott Copper in Utah, USA.The sophisticated techniques required to model these ‘supercells’ have led to a critical review of the fundamentals of flotation science and engineering. In addition to outlining the advantages of large flotation cells, this paper will review the science which supports, and occasionally challenges, the ongoing trend to larger flotation cells.

Aug 8, 2011

By Eggleton R. A, Taylor G

As the search for major orebodies continues, the regolith dominated areas of the Broken Hill region remain one of the major exploration frontiers. The past exploration philosophies emphasising bedrock related characteristics have been limited in their ability to effectively explore in the regolith dominated areas, that are a feature of most of the Broken Hill region (Figure 1). If mineral exploration is to be more effective across all of the Broken Hill region, particularly in the regolith dominated areas, an appreciation and increased understanding of the regolith and landscape features are a necessity. Ironically, there have been very few studies of the regional regolith and landscape features of this region. Owing to the previous lack of an established regolith-landform framework, this study initially concentrated on the general characterisation and variation in regolith types corresponding to the area of the Broken Hill Block, within NSW (Willis, 1989; Stevens, 1986). The approach has been to produce a regional regolith-landform map at 1:100 000 scale and conduct a complimentary program of regolith sampling and site descriptions. This was a basis to identify areas for more detailed regolith-landform mapping (1:25 000 scale) and for producing a regional stratigraphic framework and models for landscape and regolith evolution. This information can then be used to provide a basis for regional mineral exploration models and strategies. This paper highlights some of the regionally significant regolith-landform features and their general significance to mineral exploration. More detailed information on specific regolith-landform features and study sites may be found in Hill et al (1994; 1996; in press) and Hill (in prep).

Jan 1, 1997

By Braaksma M, Dudley K

The Girilambone Copper Mine is located in central New South Wales and has been operational since May 1993. Girilambone is operated by the Girilambone Copper Company (GCC) which is a joint venture between Straits Resources (60 per cent) and Nord Resources (40 per cent). The mine complex consists of an open cut pit, ore crushing and stacking, heap leaching, solvent extraction (SX) and electrowinning (EW) allowing for the recovery of 14 600 tpa of LME grade A cathode copper. The leaching and SX-EW plant have been described previously by Rear, Taylor and Busby (1994) and Miller (1994).

Jan 1, 1996

By D Macfarlane, T McMorran

  This paper describes the design and construction of the upper section of the Pike River Mine Access Road. The road extends 11 km from Logburn Road alongside Big River on the West Coast of the South Island of New Zealand up Pike Stream to the new mine portal on White Knight Stream (Figure 1). This new road starts on alluvial terraces, then partially follows an old logging road before entering virgin bush and traversing the true left bank of a steep sided gorge on Pike Stream some 50 - 100 m above the stream before dropping down to the mine entrance. Its development has been an extended process from route selection, consenting through to design and construction.   The paper focuses on the 1 km section of road between the Amenities Area (Figure 1) and the portal. This section of road is a steep section through the gorge, principally on Conservation estate, and has presented the most challenges and required the greatest effort to obtain the regulatory approvals, prepare the design and to construct. This paper outlines the history of the challenges faced and how the road was successfully developed.

Jan 1, 2007

By Edwards A. B

Re-opening of the old Proprietary mine has provided new exposures of the system of uralite dolerite dykes observed in 19110 by the Geological Sub-Committee of the Scientific Society of Broken Hill. The newly exposed dykes have been invaded and replaced by the massive ore, with development of manganese garnet and biotite adjacent to the sulphides. A similarly altered dyke in massive ore has been encountered in the North Broken Hill mine.The textures of the unaltered portions of these dykes indicate that they have not been affected by the intense regional metamorphism that has recrystallized the other rocks of the area. Moreover, in mineral composition and textures, in chemical composition, and in cross-cutting relations, they are so similar to the uralite dolerites of the district which post-date the folding and metamorphism of the Broken Hill area, that they must be regarded as typical members of this suite of post-metamorphic dykes.The invasion and replacement of the dykes by the massive sulphide ore, and their metasomatic alteration adjacent to the invading sulphides, establishes that the lode was formed after the dykes were intruded, so that the high temperature features of the ore cannot be attributed to metamorphism and recrystallization, but are primary features of high temperature hydrothermal deposition.INTRODUCTIONTwo cross dykes are recorded on the plan of the 300-ft. level of the Broken Hill lode accompanying the report on the geology of the Broken Hill lode published by the Geological Sub-Committee of the former Scientific Society of Broken Hill (1910). One dyke occurs in Block 12, and the other in Block 13, and both were stated to consist of a dark basic amphibole rock, called diorite for convenience. The more southerly dyke was traced from some distance out in the footwall, and even through the lode, by means of occasionalpatches of it still existing in the ore, in a direction...

Jan 1, 1956

By L Nicholls, A Parbhakar-Fox

Reprocessing Queensland’s historical copper tailings for critical minerals is an opportunity to circularise the state’s mining industry by converting waste liabilities into assets. Mine waste valorisation enables economic rehabilitation of legacy mining sites, improves environmental outcomes, and can establish a secondary resources sector that produces a sustainable supply of critical minerals to supplement growing global demands (Australian Government, 2023). The North-West Queensland Minerals Province (NWQMP) hosts numerous copper mines with known cobalt mineralisation (Degeling, 2020). Cobalt is a critical metal essential for green energy and battery technologies, however its global supply is at risk due to geopolitical and ethical concerns (Rachidi et al, 2021). Historically mining in the NWQMP has focused on copper, leaving cobalt to theoretically accumulate in the mine waste. First-pass geometallurgical characterisation of waste repositories by the Sustainable Minerals Institute identified Capricorn Copper as a key site for further investigation into cobalt deportment and recovery potential.

Aug 2, 2023

In the traditional analysis of airblast and ground vibration, it is often assumed that only that single blasthole possessing the maximum charge weight (or a small group of blastholes initiating in an 8 ms time window) will determine the peak level of airblast or vibration. This assumption is examined in light of waveform superposition, and is found to be unsound. More rational alternatives are given in which all blastholes (or a selected number of blastholes) are assumed to contribute, in varying degrees, to the peak vibration or airblast. However, it is also shown that there is generally no precise solution for determining the number of blastholes that contribute to the peak level. The technique of waveform superposition does not rely on evaluating which group of holes produced the peak, and so is the preferred method for prediction. Some numerical and analytical models of wave superposition for airblast and vibration are also reviewed and recent approaches discussed. The directionality of airblast, in particular, is analysed in detail using a new superposition model that is constructed of two submodels, onenumerical, the other analytical. This model shows that the airblast from even a single blasthole can be highly directional. The detail of this model is commercially sensitive and will not be discussed. Only results and a broad discussion will be given.

Jan 1, 2004

By M Astorga, J Jackson, J Gaunt

The Ban Houayxai mine (BHX) is a relatively low-grade, low-cost, open pit gold-silver deposit in Laos operated by Phu Bia Mining, a subsidiary of PanAust. The ore production rate at BHX is 4.5 Mt/a, with direct tipping to a semi-autogenous grinding/ball mill circuit with a carbon-in-leach process plant. Approximately 100 000 oz of gold is produced per annum.The operation is located in mountainous terrain with minimal run-of-mine stockpiling. This results in limited capacity for blending from stockpile, with the mill instead reliant upon direct feed ex-pit. Since commissioning in 2012, the plant has seen significant variation in milling rates due to variability in the feed properties of the oxide, transition and primary ores. As part of an ongoing continuous improvement program, an integrated approach was initiated, focusing on maintaining and enhancing production in the future as the proportion of harder primary ore increases. This will specifically focus on direct blending from the loading face using the ore properties and blast fragmentation to maintain mill throughput.This approach was based on the concept of physical assets management, commencing with improving information and knowledge of the condition of the orebody through modelling the characteristics, variability and performance of the feed for processes relevant to throughput. These models are used to support both mine and process plant production planning.The key ore feed characteristics and parameters modelled for the life-of-mine are blastability index, powder factor, crushability or impact resistance and grindability. These predictive spatial models were based on the data from diamond drill holes used in resource definition and geotechnical drilling programs by integrating geotechnical, geological, geochemical and metallurgical data.At the early stages of implementation, the models are being utilised for ore blending decisions to provide guidance and support for budgeting, long-term mine planning, blast design for mill feed and providing the mill with an expectation of performance.CITATION:Jackson, J, Gaunt, J and Astorga, M, 2014. Predicting mill ore feed variability using integrated geotechnical/geometallurgical models, in Proceedings Orebody Modelling and Strategic Mine Planning Symposium 2014, pp 165–176 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Nov 24, 2014

By Chenggang W, Junhua D

Fragmentation variation laws in loading and non-loading areas are obtained by cement block model blasting experiments under various stemming lengths. A reasonable ratio between loading length and stemming length is worked out according to minimum average fragment size. The present continuous loading structure, for its very long length to be stemmed, leads to insufficient fragmentation in the upper part of bench. So a new loading structure with bottom air cushion is put forward. The author analyzes blasting effects and processes of two loading structures' and proves that there exists a optimal length of air cushion.\

Jan 1, 1990

By K Faure

At Sams Creek a peralkaline granite porphyry dyke hosts gold-bearing sulfide-quartz-siderite veins. The granite dyke intrudes Ordovician-Silurian metapelite and quartzite, and has thin lamprophyre dykes along its contacts. The mineralised veins and the granite and lamprophyre dykes have been deformed during several phases of folding, the youngest of which is Early Cretaceous. The granite has been hydrothermally altered, with early magnetite-siderite¦biotite alteration overprinted by sulfide and quartz-siderite veins and quartz-perthite-albitesiderite-pyrite-rutile¦sericite alteration assemblages. The sulfide veins are composed of arsenopyrite + pyrite ¦ gold ¦ galena ¦ sphalerite ¦ chalcopyrite ¦ pyrrhotite ¦ graphite, with minor quartz and siderite. The gold contains significant silver (up to 30.8 wt per cent) and ranges in fineness from 692 to 844. Immiscible liquid-rich fluid inclusions, in quartz protected from recrystallisation by enclosing arsenopyrite, contain carbonic liquid that homogenised at temperatures of 320 - 355¦C. The Sams Creek gold deposit is distinguished from the orogenic gold deposit type by its stockwork veining, magnetite-siderite¦biotite alteration, high sulfide content, ore mineralogy and the absence of mineralisation in the metapelite country rocks. Features of the Sams Creek gold deposit do not fit neatly into existing classifications of granite-related gold-bearing deposits. The peralkaline composition of the Sams Creek granite and the absence of molybdenite and paucity of chalcopyrite, are major points of difference with gold-rich porphyry copper deposits. The presence of a lowsulf idation-state ore mineral assemblage, carbonate alteration and carbonic fluid inclusions at Sams Creek are typical of many of the reduced granitic intrusion-related type of gold deposit. However, notable differences are: the reduced type granites are metaluminous and calc-alkaline; and Sams Creek only has high As in contrast to elevated Bi, W, As, Mo, Te, and Sb in the reduced type. The Sams Creek deposit also shows affinities with gold deposits associated with alkaline magmatism, and it appears to be a new variant of this type as an As-Au deposit associated with a peralkaline granite and lamprophyre.

Jan 1, 2004

By A Tsang, D La Rosa, L Julian, W Hunt, R Pratama

Mineral extraction is facing the twin challenges of finite resources and increasing environmental and economic costs. Society’s growth is matched by an increase in demand for the earth’s resources. It is incumbent on us as contemporary mining professionals to promote and facilitate ethical and efficient practices. Technological developments have and will continue to allow a decrease in viable cut-off grades, creating new economic opportunities from deposits that were previously uneconomic. Mining low-grade deposits however, even with perfect recoveries, will have a significant environmental cost for the metal won. The greatest opportunity to increase metal with the least environmental cost is to increase the metal yield from existing mines. The determination of mining polygons is a critical juncture where rock is typically subjected to the binary grading of either ore or waste. This is undertaken using a number of factors, such as the mining equipment being used and the recovery in the processing plant. Traditionally, mining polygons were generated on the in situ grade control model. This ignores the complex differential blast dynamics that occur during blasting. Two opportunities exist to improve metal yields: optimise the blast design to achieve the post-blast muck pile that will support the greatest outcome; and generate optimised polygons on the post-blast muck pile. A data set from a Canadian Gold deposit is used to demonstrate the potential benefit of iterative blast movement modelling and grade control polygon optimisation on metal yield. The impact of blasting on the total recoverable metal is presented, highlighting blast design as a critical process in the mining value chain.

Sep 1, 2024

AT Broken Hill, N.S.W., where flotation is tbe final stage in concentration, the gangue minerals are discharged as an aqueous pulp, and water recovery from the pulp being essential, provision must be made for the clisposal of the solids and for the return of the clear water to the cemcentrating mill circuit. This water also contains an appreciable quantity of flotation reagents, and is more effective for flotation than supply water.In the development of a scheme for the disposal of large quantities of residues, several of the more important factors for consideration are:(a) Tonnage to be handled.(b) Grain size of residues and pulp density.(c) Dewatering of residues and return of clear water to mill circuit.(d) Possible uses of residues.-A varying proportion of sand residues is used for underground filling, and to a limited extent for yoncrete making and other engineering work.(e) Choice of a suitable dumping area. - Before deciding upon this area consideration will be given to possible future operations so that the dump will not be likely to interfere with any development which may be necessary. Again, residues dumped at the present time as worthless or...

Jan 1, 1927

By A Sainoki, D Fukuda, K Shioya, C Zhang, M Sharifzadeh, T Yamaguchi, S Saydam, J Kodama, C A. Amagu, Y Fuji

In an open pit mining, ensuring stability of rock slope has been one of the major challenges. As the instability of rock slopes can results in slope failures; consequently, endanger the lives of workers, affect mining operations, and often lead to economic losses. Therefore, it is practical important to study the stability of rock slopes to ensure safe, smooth, and effective mining. The stability of rock slope of about 130 m height at Higashi-shikagoe limestone quarry, which has experienced slope failures four times since 1996, was assessed. The rock mass at the quarry is consists of mainly limestone, schalstein and slate rocks intersected by clay seams of about 70 m thick at footwall of the rock slope. Previously, the effects of limestone excavation at the foot of the rock slope and shear failure due to rainfall infiltration have been investigated. Although clay is widely known to plays a vital role in the stability of slopes in engineering projects including mines, but its impacts on slope deformation at the study quarry are not adequately examined. Therefore, the present study addressed the impact of presence of clay on slope deformation at Higashi-shikagoe limestone quarry by employing both field displacement measurements and numerical simulation. Firstly, surface displacements measured by automated polar system (APS) for over eight years were analysed to understand the overall characteristics of the slope deformation. The measured results show that the displacement gradually decreases with time. The maximum displacement of about 100 mm was observed at middle of the rock slope. Secondly, the effect of deterioration and plastic behaviour of the thick clay layer were simulated using two-dimensional finite element method (2D FEM). The simulated results reveal that the clay exhibited deterioration by water contact and plastic behaviour arising from mining activities at the footwall of the clay zone. Based on the characteristics of the measured and simulated results, it was found that the displacement induced by the effect of deterioration of clay is more noticeable at the north side of the quarry whereas displacement induced by excavation is quite significant at the south side of the quarry if clay behaves as an elasto-plastic material.

Nov 29, 2022

By T Ireland

The Rosedale Outfall Project utilised an earth pressure balance (EPB) tunnel boring machine (TBM) for only the second time in Auckland. Although the tunnel was excavated fully through East Coast Bays Formation (ECBF), the downhill gradient along with some unexpected areas of extreme water posed some significant challenges during construction.The Rosedale Waste Water Treatment Plant Outfall Tunnel project was commissioned in response to the growing needs of North Shore City, New Zealand. The outfall tunnel comprises a 3 km long segmentally lined tunnel and a 2.1 km buried marine pipeline, and discharges 2.7 km offshore. The capacity of the new outfall is 6 m3/s compared to the existing outfall capacity of 1 m3/s.Following previous project experience, and with the tunnel located along its full length in the bedrock weak sandstone/siltstone bedrock of the East Coast Bays Formation, tunnelling conditions were expected to be benign. However during construction six high inflow events were encountered, including a full face (3.4 m diameter) of Albany Conglomerate over a 48 m length of tunnel. This material comprised cobbles in an uncemented sand matrix, and inflows were recorded at over 25 L/s.This paper describes the design challenges on the project including the design of a steel fibre reinforced segmental lining for internal pressures. The specification of the machine is discussed including the modifications required to enhance the risk mitigation measures already in place. The ground conditions encountered are discussed and have improved the knowledge of the range of conditions that can be expected in Auckland.Despite the challenges presented during construction, excellent progress rates were achieved and the tunnelling was successfully completed ahead of program.

Mar 8, 2011

The constructionhunnelling industry in Hong Kong has taken an upturn during the last few years due to the commencement of works to relocate the International Airport to Chek Lap Kok Island. The Airport Core Program consists of ten individual projects including roads, landfill, bridges and tunnels to build the airport and link it to the mainland of Hong Kong. Part of one of these projects is Cheung Ching Tunnel located on Tsing Yi Island. The proposed tunnel consists of a 1.6 km twin tube road tunnel through the island's central ridge, taking three lanes of traffic each way. Each tunnel is 16 m wide - (the widest road tunnel in Hong Kong) and are 10 m high. ICI Explosives (China) saw an opportunity to supply bulk explosives to the Cheung Ching tunnel and contacted ICI Australia to determine availability of a bulk system for tunnel blasting, and assist in its introduction. This paper discusses the methodology behind developing that bulk system to suit a large tunnel drill and blast operation. It also looks at the problems encountered along the way, from the development of the product to the application of the system in the tunnel and the eventual success of introducing a bulk system into a large tunnelling project.

Jan 1, 1994

The metals titanium and zirconium and their alloys occupy individually unique positions in the catalogue of metallic materials available for engineering applications. Because of the uniqueness of their usage in the aerospace industry, in the case of titanium, and in the nuclear industry, in the case of zirconium, the except- ional resistance of both the metals to the majority of corrosive media is often overlooked. The cost factor in the application of these more exotic metals has, in the past, been a major considerat- ion in the rejection of titanium and zirconium as constructional materials in areas of usage where their superior corrosion resistance strongly favoured them over more conventional metals and alloys. The cost differential, vis-a-vis austenitic stainless steels, for example, is not a trifling matter, being of the order of 6 : 1, when compared on a weight basis. However, if one considers the cost, for, say, sheet product on an area basis, the ratio is redgced to approximately 3.4:1. If one also takes into account the enhanced performance of the more expensive material and the sub- sequent increase in service life and decrease in plant down-time, the economics of the usage of materials such as titanium becomes much more attractive and, indeed, highly competitive,

Jan 1, 1969

By S Oliver, D Willingham

"Although a resource model is central to the mineral resource value estimation process (Glacken and Snowden, 2001), creating it is a labour intensive task and in the end it is a limited representation of the actual orebody. Creating this resource model is a labour intensive task that requires inputs from geology, mining, metallurgical and commercial disciplines. It requires thousands of samples from hundreds of drill holes to be verified, grouped in geological domain, interpolated and then valued. Even after all this effort, a model is only an estimation; therefore further significant effort has gone into quantifying uncertainty in the model and subsequently risk around any value estimates derived from the model. Further risks to model quality are introduced through poorly conditioned data or incorrect assumptions.This paper proposes the application of machine learning to automate the resource model development. Machine learning is applied to the traditionally manual tasks of geological formation, domain identification and validation of the block model mineralogy. Through automation the resource estimation process can be accelerated, allowing more drill holes or a larger resource body to be processed in a given time frame, or allowing the process to be more agile to changes in input data and assumptions. A case study based on drill hole data from a Western Australian iron ore deposit (Government of Western Australia, Department of Mines and Petroleum, 2015) is used to demonstrate the application of machine learning in this process.CITATION:Oliver, S and Willingham, D, 2016. Maximise Orebody Value through the Automation of Resource Model Development Using Machine Learning, in Proceedings The Third AusIMM International Geometallurgy Conference (GeoMet) 2016, pp 295–302 (The Australasian Institute of Mining and Metallurgy: Melbourne)."

Jun 15, 2016

By A Kondratiev, A E. Moiseev

Viscosity of oxide melts is one of the fundamental physicochemical properties that plays a crucial role in important technological and natural processes like slag flow, slag/metal separation, volcano eruptions etc. However, many oxides melt at extremely high temperatures and are highly corrosive in the liquid state, which makes experimental measurement of melt viscosity by traditional experimental techniques a difficult and expensive, if not unachievable, task. In this case it might be useful to employ other methods of viscosity determination such as modelling or simulation. In the present paper the shear viscosity coefficients of selected pure refractory oxides (CaO, MgO, Al2O3 and others) have been simulated at temperatures above their melting points via the classical molecular dynamics and compared to the available viscosity data (eg experimental viscosities, other model predictions) collected by the authors. The simulation has been carried out using the noncommercial molecular dynamics simulation package LAMMPS and the Born-Mayer-Huggins potential, which parameters for the selected oxides were taken from the literature. Several viscosity calculation techniques (the Stokes-Einstein and Green-Kubo equations, the so called Einstein relations, and the non-equilibrium molecular dynamics methods) have been applied to ensure a more reliable viscosity calculation. It has been demonstrated that the simulated viscosity of the selected melts agrees with the available experimental data. It has also been shown that the simulated viscosity of a number of the oxide melts investigated is close to that calculated by phenomenological viscosity models, which supports viscosity extrapolation from the corresponding binaries and ternaries, if no experimental data is present for unary systems. In general, it has been demonstrated that the molecular dynamics simulation could provide a reasonable estimation of viscosity if no other viscosity data is available.

Jun 19, 2024