Search Documents

By D. F. Malan, J. A. L. Napier

A number of shallow coal or hard rock mines employ pillar mining systems as a strategy for roof failure control. In certain platinum mine layouts, pillars are designed to ‘crush’ in a stable manner as they become loaded in the panel back area. The correct sizing of pillars demands some knowledge of the pillar strength and the overall layout stress distribution. It is particularly important to understand the impact of the layout geometry on the effective regional ‘stiffness’ of the rock mass around each pillar. An important design strategy is to model relatively detailed layout configurations which include a precise representation of the local pillar layout geometry and to analyse multiple mining scenarios and extraction sequences to select optimal pillar sizes and barrier pillar spacing. Although computational solution techniques are now impressive in terms of run time efficiency, a major difficulty is often encountered in assigning suitable material properties to the pillars and in devising an effective material description of the layered rock strata overlying the mine excavations. This paper outlines an efficient numerical strategy that can be used to assess large-scale pillar layout performance while retaining the ability to modify individual pillar constitutive behaviour. The proposed method is applied to selected layouts to compare estimated average pillar stress values against values determined by detailed modelling and against observed behaviour.

May 3, 2023

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 M. J. Sapko

The Pittsburgh Research Laboratory (PRL) of the National Institute for Occupational Safety and Health (NIOSH) recently conducted full-scale explosion experiments; these experiments evaluated the strength characteristics of various seal designs used for safely isolating worked-out areas in underground coal mines. Large-scale explosion tests were conducted within the multiple entry section of PRL's Lake Lynn Experimental Mine (LLEM) employing the only currently accepted test method endorsed by the Mine Safety and Health Administration (MSHA) for seal design in American mines. These explosion tests are labor-intensive, expensive to conduct, and can interfere with other critical underground safety and health research programs conducted by NIOSH. Therefore, the PRL has developed an alternative seal evaluation method, based on a hydrostatic pressure loading concept, that can facilitate the in-situ testing of seals in an operating mine. Two chambers within LLEM and one within the Safety Research Coal Mine (SRCM) were used for hydrostatic pressure loading various seal designs. The results from tests in these chambers compare favorably with those from the large-scale explosion tests in the multiple entries. Preliminary size-scaling relationships for predicting the strength of standard seal designs as a function of entry size is also presented. In addition to testing seal designs at the required 20 psi static pressure level, the new facilities also allow for the determination of the seal's ultimate failure pressure. This new approach shows promise as an alternative evaluation method for improving mine sealing technologies which, in turn, will enhance the safety of the underground personnel.

Jan 1, 2003

Primary roof support represents the first line of defence against rockrelated falls of ground in underground mines, and improper utilization or misunderstanding of the applicability and behaviour of primary support can be costly from a safety standpoint. This is a major concern for underground mines, as roof support is the single most costly expense from a mining operational perspective. This is further backed by the evidence that, in the USA, hundreds of injuries and fatalities still occur each year because of rib, roof, and massive roof falls. Additionally, the fully-grouted passive rebar, fully-grouted tension rebar, and resin-assisted mechanical anchor bolts, which constitute a large portion (89%) of the 68 million bolts installed each year in underground mines in the USA can vary in cost quite dramatically. To mitigate this concern a study was conducted in 2010 by the National Institute of Occupational Safety and Health, in conjunction with Southern Illinois University of Carbondale, to assess the performance of primary roofbolts in underground coal mines for improved safety and cost. This was accomplished using underground roofbolt monitoring solutions, field data, and numerical modelling to better understand the quasi-static behaviour of underground coal mine roofs and the response behaviour of the bolts. In particular, over 170 instrumented extensometers, closure meters, shear meters, fully-grouted passive rebar, fully-grouted tension rebar, and resin-assisted mechanical roofbolts were installed at three coal mines across the USA. Of these three mines, two used the room and pillar extraction method and the other used the longwall extraction method. There was no evidence to indicate a difference in performance of the active primary roofbolts compared with the passive primary roofbolts. Additionally, in the initial loading phase, the active bolts showed no difference in loading, indicating that tension bleed-off is of more of a concern than originally thought. Lastly, for the initial computer modelling studies, challenges still remain in obtaining a good match to the in situ bolt measurements and replicating the discontinuous roof rock and in situ bolt behaviour over time.

Jan 1, 2014

By N. Mangesana

The effect of high solids concentration and coarse particle sizes (d50 ranges 90?300 µm) on the viscosity of a suspension of water and silica sand was investigated. The experiments were designed to include conditions that have been tested by other authors and those encountered in the operation of tumbling mills. The rheological data was measured using a tube rheometer with operating pseudo shear rates of up to 1500 s-1. The rheograms obtained depicted dilatant behaviour. The Otswald-de Waele, Bingham, Herschel-Buckley and Casson models were fitted to the experimental data. The Herschel-Buckely model, which approximates the yield point and gives an indication of the shear thinning or shear thickening behaviour of the suspension, was found to provide the best description of the flow curves for all slurries. The Herschel-Buckely model was then used to determine the apparent viscosity of all the tests. The apparent viscosity and yield stress increased with solids concentration and particle size at the different pseudo shear rates. The increase in slurry viscosity with solids concentration was attributed to increased frequency of particle-particle interactions. The increase in slurry viscosity with particle size was attributed to increased inertial effects. At low shear rates the sand slurry exhibited shear thinning behaviour for all solids concentrations. At higher shear rates above 1000 s-1the high solids concentration suspensions exhibited a transition from shear thinning to a shear thickening behaviour.

Jan 1, 2008

By S. W. Jacobsz, E. P. Kearsley, N. K. Schoeman

A physical modelling approach was adopted to simulate the process of cave mining propagation at various horizontal to vertical stress state ratios, and image data was captured and used to characterize the resulting failure mechanisms. Image processing was integrated into the study to determine minor principal strains in models to identify where critical extensional strains had been exceeded relative to crack propagation. The study exhibited a scale factor magnitude of 119 using a calculated critical extensional strain value of 0.014% for the modelling material used. It was perceived that cave propagation occurred through successive parallel extensional ‘fracture bands’ in all cases, as noticed in models of prior studies, deviating from that described by the Duplancic model. Cracks tended to bisect regions of minor principal strain equal to the critical extensional strain, which suggests ‘fracture banding’ is an acceptable failure mode for brittle, rock-like materials. Moreover, it was found that models with lower horizontal to vertical stress ratios (K ratios) showed a greater degree of vertical development of the caving mechanism (at similar vertical stresses), while the vertical extent of the cave was suppressed in models with higher K ratios. Models subjected to larger vertical stresses experienced slower cave formation in latter time-steps. These models exhibited larger two-dimensional caved perimeters and areas, respectively.

Nov 30, 2022

By J. B. A. Solar Menéndez, M. A. Martínez Fidalgo, L. A. Alvarez-Santullano

1. ANTECEDENTS This present communication refers to the results of the research project entitled "Training in the correct use of chemical oxygen self rescue apparatus and also an analysis of apparatus autonomy", carried out under the auspices of the Specific Collaboration Agreement signed between the Industry Council of the Principality of Asturias (Spain) and INSALUD (National Silicosis Institute, Spain). The antecedents of this project lie in work undertaken by the National Silicosis Institute aimed at optimising efficiency in the use of chemical oxygen self rescuers through improvements in worker training. To be precise, by means of Agreements drawn up between the above-mentioned parties, a simulator was designed and incorporated as the basis of a training station within the " Fundación Barredo" of the Industry Council of the Principality of Asturias. This has been operative since the year 2001 and, up to the present, 2300 Asturian coal mine workers have been trained there. What prompted this activity was the observation of various errors in the use of the self rescuers, deriving from their innate characteristics, especially with respect to the sensation of discomfort - the high temperature and humidity of the breathed oxygen, resistance and inconvenience - and limited autonomy. It was precisely this insight that, in critical situations, lack of familiarity with the conditions of discomfort likely to be experienced during the evacuation process might lead to rejection of the equipment, with fatal consequences, that has justified the setting up of this training and simulation centre. The centre is able to provide all mine workers with the necessary training, and experience of the characteristic sensations experienced when using such equipment (Fig.. 1).

Jan 1, 2003

By S. A. Ziebold

There is often a need to increase capacity through existing mist eliminator installations without significant tower modifications. Often there is not enough extra space outside existing vessels unless more expensive changes can be justified. Also, plant downtime must be considered. The challenge in retrofitting is to overcome velocity limitations inherent in the design while at the same time extend performance capability. The solution most often is to create new ways to add extra filter area while remaining constrained to the existing space. This strategy, however, can dangerously backfire if basic design principles are not followed. Mist eliminator design is a complex integration of many conflicting performance parameters. Often the objective is to increase scrubbing capacity while maximizing collection efficiency and service life and minimizing operating pressure drop (energy cost), particle regeneration (re-entrainment) and overall maintenance costs. Results of plant retrofitting using mesh pads, impaction fibre beds and diffusion fibre beds are presented based on actual plant projects. Often the retrofit project can be more effective by increasing the scope and combining other new plant equipment with the upgrade such as new tower components. In some cases, the type of mist eliminators can be upgraded from impaction type (velocity dependent) to Brownian diffusion type (velocity independent) if additional downstream protection or reduced emissions are desirable, especially at turndown rates. Trade-offs between different fibre bed types and styles are discussed as they relate to performance and maintenance.

Jan 1, 2007

By L. Cong, C. Liu, N. Zhang, W. Cao, C. Zhou, X. Fan

"The modes of occurrence of arsenic and the effects of low-intensity leaching-flotation on arsenic removal from coal slime from the Tanggongta Plant, Inner Mongolia, China were investigated. The coal slime was examined using hydride generation inductively coupled plasma optical emission spectrometry and X-ray diffraction to obtain the content of elements and the major minerals. The modes of occurrence of arsenic in minerals were determined using selective leaching, float-and-sink analysis, and polarized light microscopy. The results indicate that pyrite is the dominant carrier of arsenic in the coal slime, which predominantly exists in association with clay. Significant proportion of the arsenic is removed by a low-intensity leaching-flotation process, consistent with the data from selective leaching and flotation of gangue. The results show that low intensity leaching-alkali washing-flotation is more efficient than direct flotation and low-intensity leaching-flotation.IntroductionUnderstanding the occurrence of arsenic in coal and its removal is significant in optimizing coal utilization, because arsenic in the environment is viewed as a potentially toxic trace element (Finkelman, Belkin, and Zhang, 1999; Hall, 2002; Smith et al., 1992; Smedley et al., 2003; Zheng et al., 1999). Many studies on the modes of occurrence and concentration trends of arsenic in different coal and washing products have been performed (Diehl, Goldhaber, and Hatch, 2004; Kolker et al., 2000; Quick and Irons, 2002; Wang et al., 2006; Fan et al., 2016). Trace elements are usually preferentially associated with certain minerals and the form that these minerals take will influence the efficiency of removal during cleaning (Quick and Irons, 2002). Arsenic is mainly associated with minerals such as pyrite, carbonate, and silicate (Demir et al., 1998; Diehl, Goldhaber, and Hatch, 2004; Zhou et al., 2014). The existence of organically bound arsenic was also confirmed using X-ray absorption (Kolker et al., 2000). However, inorganic arsenic is the principal source of pollution because of its leachability and high concentration compared with organic arsenic (Fujino et al., 2004; Gao, Lu, and Wang, 2010; Jiang et al., 2008; Kolker et al., 2000)."

Jan 1, 2017

By L. Thomas, D. Muteba, L. I. Nsenga, B. Mbao, G. Longo Banza

Owing to the high ore grades in the African Copperbelt and variable acid consumption, agitated leaching is common practice and is used more widely in Africa than in other regions. Two common challenges for agitated leach–solvent extraction circuits are higher presence of colloidal silica in the pregnant leach solution and elevated levels of total suspended solids. Both can negatively affect the physical and metallurgical performance of the solvent-extraction circuit, and lead to higher processing costs or reduced copper production. Operating in organic continuity is preferred to address such challenges, but is not always possible: additional treatment options may be advantageous. For this reason, Solvay recently introduced ACORGA® CR60LT, a new chemical additive to address these concerns. Continuous addition to the pregnant leach solution improves physical performance without negatively affecting Cu transfer or cathode quality. Results and benefits from commercial utilization of ACORGA® CR60LT at different operations in the Copperbelt are reviewed

Jan 1, 2023

By I. V. Filippova, Y. Foucaud, M. Badawi, O. Barrès, L. O. Filippov, S. Lebègue

Sodium silicate (Na2SiO3) is one of the most common reagents used in the froth flotation process, in which it is an efficient dispersant as well as a powerful gangue minerals depressant. Due to its high affinity with a lot of minerals, mostly silicates, it adsorbs on their surfaces and prevents the adsorption of the collectors, maintaining them hydrophilic. Hence, Na2SiO3 is widely used in the flotation separation of rare earth minerals, zinc minerals, iron minerals, scheelite, fluorite, apatite, and many other minerals. For decades, authors have demonstrated higher efficiency of Na2SiO3 when it is used in combination with sodium carbonate (Na2CO3), a common pH regulator. Among many well-known depressants, this system presented the best efficiency for many ores, including the Tabuaço scheelitebearing skarn studied herein. However, the mechanisms resulting in these strong positive synergistic effects are not well understood at the moment. Hence, the response surface methodology (RSM) was applied to this ore, considered as a case study, to optimize the depressing conditions and to study deeply the flotation behavior of this system. Strong synergistic effects between the two reagents were exhibited by the RSM, mostly affecting the scheelite, fluorite, and calcium-bearing silicates recoveries. This was attributed to surface carbonation of the minerals by Na2CO3, subsequently leading to a better depression by Na2SiO3. Overall, a 1:1 ratio between the two reagents displayed the best results in terms of WO3 grade and recovery. To gain understanding in the molecular mechanisms responsible for the aforementioned synergistic effects, FTIR studies were combined with state-of-the-art molecular modelling (Ab Initio Molecular Dynamics at 300 K) on fluorite, an archetype of calcium minerals. Keywords: Synergistic effects, scheelite, depressants, FTIR, DFT calculations

Jan 1, 2020

By M. S. Shishvan, J. R. van Duijvenbode, L. M. Cloete, M. W. N. Buxton

Geochemical and mineralogical datasets from Tropicana Gold Mine, Australia, have been used to define ore fingerprints. VNIR/SWIR spectral data were represented by four normalised wavelength regions and were clustered to form spectral classes. Sequentially, these spectral class proportions within a block and collocated XRF data were clustered to form material types (fingerprints). The material types were related to an Equotip-BWi correlation. These correlations can be used to extrapolate a hardness signature and generate a BWi proxy for different blocks. The combined fingerprints and BWi proxy can assist as a tool for enhancing the prediction of comminution behaviour. They can explain specific domain-related hardness variations. For example, one material type could be separated into a softer (~15-18 kWh/t), and harder (>20 kWh/t) material blend. This was accomplished using the commonly overlooked VNIR region at 605 nm. This outcome has significance for blending strategies.

Sep 1, 2021

By N. C. Joughin

A simplified economic model is derived for the financial return from a new gold mining project. The model is used to analyze the effects of changes in the mining rate, capital costs, working costs, value of the reef and the stoping width. The most effective avenues for improving the return from a new project are indicated. It is shown that those technological innovations which would improve the return significantly, and therefore make it possible to open more mines, are the use of smaller mines, means for improving the mining rate to capital costs ratio for smaller mines, planning and control of mines by computer simulation, an improved mine valuation procedure and means for mining at a decreased stope tramming width.

Jan 1, 2014

By L. M. Tavares, T. M. Campos, P. M. Faria, R. M. Carvalho

Recent advances in mechanistic ball mill modelling are expected to allow for more accurate predictions that overcome the well-known limitations of phenomenological approaches, including their inability to predict the effect of liner geometry and their wear, in grinding. Although able to describe continuous operation, the UFRJ model of ball mills has only been validated in the past in mills that operate in batch mode, either in laboratory or pilot scale. The present work demonstrates the application of the mechanistic UFRJ ball mill model to predict fine grinding of iron ore in a pelletising plant in Brazil, the first application of such a mechanistic modelling in continuous mode. It shows that the model provides accurate descriptions of flows after iterations in a closed-circuit with hydrocyclones. Keywords: Ball mill, iron ore, grinding, modelling, particle breakage

Jan 1, 2020

By A. L. Edwards

A retrospective study was conducted on the Rand Mutual Assurance (RMA) noise-induced hearing loss (NIHL) compensation claims from 1998 to 2008 to determine if the 2008 milestone agreed at the 2003 Mine Health and Safety Summit was achieved. The number and costs of NIHL compensation claims in different commodity sectors and workplaces were collated. A secondary analysis of the ages of employees compensated after 2008 was conducted. A complementary retrospective analysis of audiogram data investigated the percentage loss of hearing (PLH) shift in different homogeneously exposed groups and occupations at two gold mines. The compensation claims analysis indicated a significant decrease in NIHL claims from 1998 to 2008, but the milestone was not achieved. The reason may be either that claims have not been submitted timeously as required by Instruction 171 and that the current submissions are a result of pre-2003 noise exposure, or that employees who were baselined are still developing NIHL owing to ineffective hearing conservation programmes in place in the mining industry. On the basis of best practice for hearing conservation, recommendations are made for leading indicators in hearing conservation programmes and for reducing the risks of NIHL in order to achieve the 2013 milestone.

Jan 1, 2012

By T. T. Matuludi

Kleinkopje Colliery has five Wemco drums that were installed and commissioned in 1978. Over the years, there have been remarkable changes in mining and geological conditions at the colliery. First, an increase in the percentage of fines was recorded, which resulted in an uneven split or an imbalance in silo levels between material feeding the cyclone (-12 mm) and Wemco drum (+12 mm) sections respectively. This problem was solved by doing modifications around the banana screens to allow finer material to the drums. The imbalance challenge was then followed by spontaneous combustion that takes place in situ and results in low product yield. The results obtained from efficiency tests indicate that although there is a limit in the percentage of -6 mm material in the feed, Wemco drums are still very efficient and applicable to the current mining conditions. These results also suggest that it is essential to lower the throughput when the drums are treating low yielding coal.

Jan 1, 2010

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 C. Musingwini, T. C. Sebutsoe

"The primary objective of any mining business unit is to make profit by extracting, processing, and selling minerals from a particular mineral deposit. It is important to optimize extraction of the mineral resource given time, capital, and space constraints. Mineral extraction is often associated with uncertainty due to variable technical and human factors. Technical factors such as grade distribution, ground conditions, and equipment reliability influence the performance of a mining production system (MPS). The performance of the MPS is also impacted by human factors such as employee skills, health, and attendance. Uncertainty associated with technical and human factors often leads to actual output differing from what was planned. Therefore an in-depth analysis of the significant causes of deviations from the planned outcomes is relevant.In this paper we investigate the empirical relationships between inputs and outputs in a MPS in order to assist management in directing efforts at key production drivers. Once a MPS has been empirically characterized, more effort and resources can be focused on the key decision-making variables (DMVs) in order to meet the planned outcomes.For a typical underground, tabular, narrow-reef conventional drill-andblast platinum mine the face advance, face length mined, number of teams, team efficiencies, and team size have a statistically significant relationship with the centares (m²) produced, which is a key performance indicator (KPI). The production function derived in this study can be used to align the physical, technical, and human factors together to predict the optimal output level. The most significant production lever of the MPS is the face advance. IntroductionA mining production system (MPS) is a result of an iterative process of design, planning, and optimization of mining input variables and decision-making variables (DMVs). The MPS exists within the mineral extraction link of the complete mining value chain. It represents the stage where mining companies have the opportunity, through production, to start generating returns on the investments undertaken. Returns on shareholders’ investments can be realized at this stage."

Jan 1, 2017

By D. C. Krige

The estimation of the tonnage and grade of payable ore in a mine and the correct policy of selective mining based on such estimates, is of vital importance to the mining engineer it is surprising, therefore, that more attention has not been devoted on the Witwatersrand to the scientific improvement of mine valuation methods. At present these methods consist almost entirely of the application of simple arithmetic and empirical formulae guided by practical experience, and ignore the many advantages to be gained from a careful statistical analysis of the behaviour of the gold values. The science of statistics has expanded rapidly during the last two decades and its value as a powerful and indispensable tool is now recognized not only by research workers and scientists but also, ever-increasingly, by the commercial and industrial world. This being the case it is noteworthy that in a mining field such as the Rand with its highly developed and advanced mining methods, singularly little attention has been paid to the analysis of mine valuation problems on a modern statistical basis. This omission is even more striking when cognizance is taken of the wealth of sampling data concerning the gold ore which is available and of the far-reaching decisions and deductions constantly being based on such data. Some contributions have however been made from time to time towards the application of statistics to mine valuation on the Rand (References Nos. 7, 8, 9, 11, and 12), but none of these have been generally applied in practice and a systematic practical approach based on clearly deemed fundamental concepts still appears to be lacking. The object of this paper is, therefore, to define these statistical concepts and to indicate briefly the lines along which statistical methods can be applied profitability in solving some of the existing problems and in improving the general standard of mine valuation on the Rand For this purpose a knowledge of mathematical statistics is essential, but as this paper is primarily intended to provide the mine valuator with a practical introduction to statistical application in his field, the mathematics is confined to a short annexure, together with references from which details can be obtained.

Jan 1, 1994

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