Search Documents

By Peter Klaver

Deep tunnels have been used for decades for the temporary storage of large volumes of urban runoff, both storm water and combined sewer overflow (CSO). Tunnels are attractive because they reduce surface disruption and do double duty as storage and conveyance systems, collecting overflows from multiple locations and transporting them for treatment to various facilities. Like any hydraulic system, however, there are dynamic effects that must be fully understood during the design phase so that the final constructed product does not have unintended negative consequences. Preliminary sizing of storage tunnels is typically based on the total volume required. For example, if studies have determined that a volume of 56.8 ML (15 million gal) will provide the intended level of control, then that figure can be used, along with a preliminary idea of the desired tunnel length, to compute the required diameter. Of course, filling a long tube with water from various points along the way is not at all the same as pouring water into a bucket, and these dynamic effects require additional attention from design engineers. The rapid filling of storage tunnels can lead to the formation of bores, which are essentially moving walls of water. An open channel bore can build in height until it reaches the crown of the tunnel, at which point it becomes a pipe-filling bore. Figure 1 depicts both kinds of bores moving upwards in a rapidly filling tunnel. Bores can strike the ends of tun-nels with tremendous force, driving water and entrained air upwards into dropshafts and creating strong pressure surges in portions of the tunnel that are already filled. If the hydraulic surges are high enough to reach grade, they can cause overflows and/or damage structures at the surface. Trapped air pockets can occur through several phenomena, including the reflection of bores off dropshaft structures and the occurrence of ?premature pressurization.? The latter refers to the situation where the inflow to a tunnel is increasing more rapidly than the tunnel can convey it, such that the pipe becomes surcharged at an upstream point while points downstream are still flowing as an open channel. These air pockets are a great nuisance as they migrate toward points in the tunnel where they can rise to the surface, which are often shafts that are already filled with water (Fig. 2). The ?geyersing? that results can be hazardous and costly.

Jan 3, 2010

By Angelus S. C. Owili-Eger

Although the foundation of a unified treatment of fluid flow in variably saturated porous media was established more than seventy years ago, no serious attempt was made to develop it until the early 1970's. Ever since, a number of unsaturated-saturated numerical models have been reported in the literature. However, none of these models have treated multidimensional, multiphase flow in deforming media, parts of which may undergo extensive secondary fracturing. A situation such as this is not uncommon in coal mining operations in which the longwall method of extraction is employed. The present model has been structured to handle such physical systems. This paper presents the mathematical formulation, the numerical algorithm, and a further validation and evaluation of the modified model. The model treats a multiphase flow in a three- dimensional regime and handles mining excavations of complex and changing configurations.

Jan 1, 1980

By A. S. C. Owili-eger

The computer-aided simulation model re¬ported here is an extension of the one developed by the author while at The Pennsylvania State University. The present model, however, includes consideration of the deformability and secondary fracturing of the porous system. The modifications enable simulation of multiphase shallow unsaturated-saturated and/or artesian ground water systems in dynamic mining environments. The model has been validated, verified, and evaluated on a number of field case studies and the results have been excellent.

Jan 1, 1984

By Tsvetan Damyanov, Nikolay Hristov, Evgeni Kostadinov, Lubomir Kuzev

"The process of grinding in drum mills is directly dependent on the dynamic forces occuring in the free movement of the grinding bodies in the volume. The subject of discussion are the potential impacts between grinding bodies having different shape and equal weight in a drum mill that has a diameter of 1200 mm. The study utilised ball grinding bodies having a diameter of 120 mm and tetraballpebs having a conditional diameter of 120 mm. The study was carried out using the Solid Works Simulation software. The following cases were examined: • Impact between two identical grinding bodies; • Impact between a grinding body and the shell surface. The study provided numerical values for impact force, kinetic energy of impact, linear speed of impact, linear displacment of the center of gravity and duration of impact deformation. The impact force of the balls is approximately 179000 N. The impact force of the tetraballpebs with typical vertices, edges and spherical surfaces depending on the area of impact is 177000 to 185000 N. Upon impact between two spherical surfaces of tetraballpebs, the contact area is 263.55 mm2, which is more than twice the contact area upon impact between two balls, which is 114.39 mm2. This determines better conditions for the transmission of kinetic energy by reducing opportunities for overgrinding."

Jan 1, 2014

By C. Pomeroy, J. Weyer, J. Brune

"Mining began in the Erzgebirge (“ore mountains”) region of Saxony, Germany more than 800 years ago, and extensive underground workings cut through the entire area. Some of these old workings are still used today: not for production, but for offering an underground venue for research, tourism and university teaching. One example is the Reiche Zeche mine at the Technische Universität Bergakadamie Freiberg which hosts one of the oldest mining engineering programs in the world. It is imperative to keep all active mine workings well ventilated to make the mine safe for students, visitors and researchers. Since the mine workings at Reiche Zeche are so vast, extensive and, to a large extent unexplored and unmapped, it has been difficult to provide a complete, functioning and adequate mine ventilation model for this mine that accounts for all known and unknown leakage air flows. Researchers at the Colorado School of Mines have cooperation with their colleagues at Freiberg to establish a mine ventilation network model based on a limited number of air readings in known and accessible drifts. In order to balance the measured air quantities, leakage flows were quantified and mapped by comparing input and output flows of mapped drifts, and modelled using Ventsim Visual, a powerful visualization software package. The resulting mine ventilation model is used to locate and quantify the leakage air flows and to examine how much air is lost in each area. With this model, ventilation changes and improvements underground can be documented reflected in the model. This enables the Bergakadamie Freiberg to project air changes required for scientific experiments and to set up mine ventilation in specific areas and for specific needs. The model also enables better ventilation control and monitoring than before. Finally, the model can be used to predict ventilation changes and escape way contamination in case of a fire, ground fall or inundation in the mine. The latter becomes important as the lower mine workings are often flooded during spring runoff, which has created uncertainty about mine ventilation in the past. HISTORICAL BACKGROUND The discovery of silver, lead, zinc, tin and other valuable metal ores in the Erzgebirge region of Saxony, Germany began in the12th century and started a mining boom than lasted for centuries. Surface mines quickly developed underground to access the higher grade veins. Explored by the early geologists and driven by the ingenuity of the miners, the underground works expanded vastly throughout Freiberg and the surrounding region. Development followed the highest grade veins, and engineers developed solutions to drain water and provide limited ventilation to the underground workings. Still, engineered and planned ventilation was usually not a priority in the development of a new drift, which today makes ventilation modeling a significant challenge."

Jan 1, 2015

By V. Bojarevics, K. Pericleous

"A number of different methods have been developed for noncontact electromagnetic treatment of liquid metal and to investigate the melting/solidification processes. Applying AC magnetic field in terrestrial conditions, along with the buoyancy and thermo-capillary forces, results in turbulent flow. The use of a homogenous DC magnetic field allows damping both of the turbulence and the large scale flow at different rates leading to surprising results at some stages. The dynamic interaction of the turbulent flow with the oscillating interface needs to be accounted if it is required to levitate liquid metal of up to few kilograms without the contact to container walls. At the high values of magnetic field some oscillation modes are damped quickly, while others are modified with a considerable shift of the oscillating droplet frequencies and the damping constants from the non-magnetic case. Numerical models are used to investigate behavior of liquid metal suspended in high DC magnetic fields providing microgravity-like conditions with laminar viscosity and heat transfer to the surroundings.IntroductionSince the early Electro Magnetic (EM) levitation experiments [1] in AC magnetic field it became apparent that the levitated liquid metal is prone to oscillation and instability. The EM and electrostatic levitation experiments with liquid metal droplets show difficulties related to confinement stability and a need for complex correction functions to establish a correlation between the measurements and the droplet material properties [2,3]. Intense internal fluid flow is visually observed, apparently being in the turbulent regime for earthbound conditions."

Jan 1, 2012

"In the conceptual stages of engineering, benchmarking and empirical ratios are commonly the only references for estimating the production rate of greenfield projects. For underground mining projects, the lack of similar cases often makes the comparison difficult due to each project’s unique geology, rock quality or other singular characteristic. This paper discusses an analytical alternative for determining the maximum achievable production rate for underground stoping mining by modeling production as a continuous process. The results show that long-term maximum capacity depends exclusively on the achievable development rate and the development intensity of the design. In operations constrained by backfill, the development rate could be limited by the backfilling rate. Extraction rate and stope size only have impacts on capacity estimates during the ramp-up period.These findings allow shaping of the maximum capacity profile of a project. They can be used as a tool for quick assessment of conceptual capacities, and sensitivities of mining rates. Benchmarking tests compared favorably with projects, mines in operation and scheduling exercises. IntroductionA key consideration in the value proposition of a mining project is its capacity or production rate (Tatman, 2001). At a conceptual level, it is common practice to refer to empirical references (Taylor, 1986; Long, 2009), rules of thumb (de la Vergne, 2014), benchmark ing studies or consultant experience to estimate or validate a target value. There is, however, a lack of analytical methodology to guide an early estimate of the maximum capacity.Araneda and Gaete (2004) presented an interesting mathematical approach to estimate the capacity of panel caving mines. Using a continuous model, they confirmed that the mine capacity is essentially determined by the rate of area incorporation, in terms of area per period, and ore column height. This simple relationship has been confirmed by multiple detailed scheduling works"

May 1, 2016

By F. M. Doyle

Metal removal kinetics and metal ion selectivity in ion flotation were studied from a thermodynamic perspective. Surface tension data could be used with the Gibbs adsorption equation to quantitatively predict the metal-ion removal kinetics in alkylsulfate-Cu systems. However, surface tensions predicted cation selectivity coefficients in the sodium dodecylsulfate (SDS)-Cu-Ca and SDS-Cu-Pb systems rather poorly, although the predicted order of the selectivity (Cu2+ < Ca2+ < Pb2+) was correct. Another thermodynamic selectivity model, which used the Grahame adsorption equation and a geometrical analysis, predicted cation selectivity coefficients that agreed very well with experimentally measured values.

Jan 1, 2002

By F. M. Doyle

Metal removal kinetics and metal ion selectivity in ion flotation were studied from a thermodynamic perspective. Surface tension data could be used with the Gibbs adsorption equation to quantitatively predict the metal ion removal kinetics in alkylsulfate-Cu systems. However, surface tensions predicted cation selectivity coefficients in the sodium dodecyl-sulfate (SDS)-Cu-Ca and SDS-Cu-Pb systems rather poorly, although the predicted order of the selectivity (Cu2+<Ca2+<Pb2+) was correct. Another thermodynamic selectivity model, which used the Grahame adsorption equation and a geometrical analysis, predicted cation selectivity coefficients that agreed very well with experimentally measured values.

Jan 1, 2000

By K. D. Luxbacher

Understanding methane emissions in underground coal mines is critical for a safe and productive mine. In addition to reasonable estimation of initial coalbed reservoir parameters, it is also crucial that changes in effective stress due to mining and pore pressure reduction are taken into account due to their effects on porosity and permeability. Primary parameters for estimation of emissions or modeling of the mining environment for this purpose are porosity and permeability which can change dramatically as a result of stress redistribution associated with mining and gas desorption from a large coal volume. These parameters affect the emission rates and ventilation requirements, as well as water inflow into the working environment. Stopping leakage, on the other hand, is a secondary stress dependent factor in estimation of emissions, as convergence of the roof and floor strata, compromising the integrity of the stopping, may result in leakage, making prediction of ventilation requirements difficult. This paper aims to examine the effects of porosity and permeability changes of the coal seam on methane emissions in an underground continuous miner section. The models were developed and executed in a dynamic fashion to simulate an advancing section. Through this process, the changes of effective stress in coal, particularly their change paths, on porosity and permeability were incorporated into the models and methane emissions, concentrations, air requirements, water inflow and possible leakage from the stoppings were investigated using a conventional coalbed methane reservoir model.

Jan 1, 2009

By Vaughan R. Voller

"Over the last ten years an increasing number of heat and mass transfer models of solidification processes have been presented in the literature. A key component in these models is the coupling between thermal and solute fields. This is achieved by developing relationships based on the behavior at the microscopic scale in the casting (typically a secondary dendrite arm space). An important class of models in this respect are models that describe the local scale redistribution of solute alloy components-collectively referred to as microsegregation models. This paper explores this class of models.IntroductionModeling takes a key role in processing materials for properties. This is particularly time in solidification processes where current computer techniques allow for a detailed treatment of the critical macroscopic phenomena. Typically these macro scale models consist ·of three main coupled components.1. A model of flow in the melt and solid + liquid mushy region.2. A model of heat transfer.3. A model of solute transport.The focus of this paper is the solute transport component, in particular the manner in which process scale transport is coupled to transport ·at the local scale of the solid-liquid interface. This step requires a modeling of the redistribution of solutes at the scale of the secondary arm spaces in the dendritic mushy region-a process referred to as microsegregation. Basic microsegregation models use the limit assumptions of no mass diffusion in the solid (Gulliver-Scheil [l]) or complete diffusion in the solid (equilibrium lever rule) in a fixed and space. Recently, however, more sophisticated microsegregation models that account for finite mass diffusion and coarsening of the arm space have been developed [2-19]. In this paper governing equations for a general microsegregation model for a multi-component alloy in a coarsening and spacing _are introduced. From these general equations various models for a binary system are developed and tested"

Jan 1, 2000

By Jan M. Mutmansky, Y. J. Wang

In the mathematical analysis of electrical networks, the electrical devices are often modeled using idealized network elements such as current sources, voltage sources and resistors. An analogous approach is outlined for the modeling of mine ventilation networks, which can be considered nonlinear resistive networks. The proposed modeling technique uses the following five basic (nonreducible) elements: H, Q, D, R and N. These elements are characterized by their pressure gain or loss, their air quantity and their pressure-vs.-quantity relationships, each of which may be known or unknown. In addition to the definition of the basic elements, the necessary conditions for which there exists a unique solution are discussed. This is performed using theoretical developments derived from electrical circuit theory and illustrated in example networks. The proposed procedures are useful in the solution and characterization of mine-ventilation networks.

Jan 1, 1998

By Y. J. Wang

In the mathematical analysis of electrical networks, the electrical devices are often modeled using idealized network elements such as current sources, voltage sources, and resistors. An analogous approach is outlined here for the modeling of mine ventilation networks, which can be considered nonlinear resistive networks. The modeling technique proposed uses five basic (non-reducible) elements: the H element, the Q element, the D element, the R element, and the N element. These elements are characterized by their (a) pressure gain or loss, (b) air quantity, and (c) pressure-versus-quantity relationship, each of which may be known or unknown. In addition to the definition of the basic elements, the necessary conditions for which there exists a unique solution are discussed. This is performed using theoretical developments derived from electrical circuit theory and illustrated in example networks. The proposed procedures are useful in the solution and characterization of mine ventilation networks.

Jan 1, 1996

By Y. J. Wang

In the, mathematical analysis of electrical networks, the electrical devices are often modeled using idealized network elements such as current sources, voltage sources and resistors. An analogous approach is outlined for the modeling of mine ventilation networks, which can be considered nonlinear resistive networks. The proposed modeling technique uses the following five basic (nonreducible) elements: H, Q, D, R and N. These elements are characterized by their pressure gain or loss, their air quantity and their pressure-vs. -quantity relationships, each of which may be known or unknown. In addition to the definition of the basic elements, the necessary conditions for which there exists a unique solution are discussed This is performed using theoretical developments derived from electrical circuit theory and illustrated in example networks. The proposed procedures are useful in the solution and characterization of mine-ventilation networks.

Jan 1, 1997

By A. L. Ling, A. K. Janzen

"Design of effective mine water treatment plants depends on the ability to reliably predict treatment outcomes, but systems with multiple chemical and physical treatment processes and recycle loops can be difficult to model. GoldPHREEQC is a modeling framework that integrates mass and flow routing capabilities of GoldSim software and water quality calculations provided by PHREEQC, and is well-suited to modeling parameters of interest in the mining industry, especially metals and sulfate. Physical processes such as membrane separation and greensand filtration are modeled by defining mass rejections for each parameter and flow routing based on vendor projections and bench- or pilot-scale equipment testing. GoldPHREEQC models chemical precipitation processes such as HDS metals removal and gypsum precipitation by (1) performing surface-aging and solubility calculations in PHREEQC and (2) using the GoldSim framework to direct supernatant and precipitate to subsequent process steps. Users can apply this modeling framework to predict effluent water quality, chemical requirements, sludge production, and required equipment capacity to support design, troubleshooting, and cost estimating. INTRODUCTION Design of effective mine water treatment plants depends on the ability to reliably predict treatment outcomes, but systems with multiple chemical and physical treatment processes and recycle loops can be difficult to model. GoldPHREEQC is a modeling framework that integrates mass and flow routing capabilities of GoldSim software (1) and water quality calculations provided by PHREEQC (2), and is well-suited to modeling treatment of parameters of interest in the mining industry, especially metals and sulfate. GoldPHREEQC can also be applied to other water treatment scenarios using inorganic reactions and processes. The GoldPHREEQC model framework moves water and mass between specified model nodes and calculates specified reaction and water quality outcomes. It can be organized to provide flow, load, and concentration information for each location, including effluent and byproduct streams, and outputs can be adapted to include power and chemical costs."

Jan 1, 2018

By C. V. Deutsch, D. A. Silva

The geologic modeling of multiple domains is required for the prediction of mineral resources and reserves. In general, explicit methods have wide acceptance among modelers due to the ability to manually add deemed important trends and geological features into the models. As an alternative, implicit geologic modeling techniques are increasingly adopted to reduce professional time requirements. These methods aim for objective, repeatable and flexible models by replacing the manual time-consuming interpretation of sections with the interpolation of signed distance values. These techniques, however, have been developed for the binary case of two domains at a time. The current paper presents an algorithm that extends the binary case to multiple domains. The method is able to transfer large-scale geological features from data to geologic models. A case study illustrates the benefits of the algorithm and practical implementation considerations.

Jan 1, 2016

By M. C. Williams

The solution of the dynamic population balance model (PBM) of a closed grinding network consisting of a mill, sump, and cyclone is problematic. The fact that the dynamic PBM model equation solution to the Inverse Problem for grinding systems is degenerate or underspecified is demonstrated. Two numerical solution approaches to the Inverse Problem are used. These are: 1) providing additional constraints on breakage selection functions or 2) performing the Arbiter-Bhrany (or other) normalization of the selection functions. Actual experimental anthracite batch grinding data is used to demonstrate the non-unique functionality of the batch dynamic mill selection and breakage functions for a real physical system. The Levenberg-Marquardt algorithm for systems of constrained non-linear equations is used to solve the batch dynamic PBM grinding equations to obtain the grinding selection and breakage rate functions. Different solutions to the same PBM transport equations are provided. The mill, sump and hydrocyclone were modeled as a CSTR operating at various retention times. Batch dynamic PBM data was used to provide the mill kinetic and breakage selection function data. Two different solutions were obtained depending on the numerical solution approach. The severity of the non-uniqueness problem for dynamic grinding is demonstrated. Each solution approach to a dynamic PBM with transport, while giving the same prediction for a single batch grinding time, gives different solutions or predictions for mill composition for other grinding times. This fact makes dynamic nodal analysis and control problematic. The fact that the constraint solution approach gives a solution may suggest that normalization for closed networks is not necessary. Differences in solutions to the PBM cannot be excused away by inaccuracies in the data used to model the grinding phenomenon.

Jan 1, 1995

By B. Farouk, Y. G. Kim, D. Apelian

"Continuous strip casting has received considerable attention in recent years due to its economic advantages with increased energy savings and the elimination of intermediate processing steps. A promising concept in strip casting is the twin-belt (Hazelett) process which incorporates two cooled continuously moving thin metal belts wrapped around pulleys which tension and drive the belts at a predetermined speed. Molten metal is introduced at the entry end of the caster through a tundish and feeding nozzle into the mold. The twin-belt process has been successfully used with non-ferrous metals such as aluminum, copper and brass, and is a prime candidate for strip casting of ferrous materials. A numerical model has been developed which predicts the interface locations, the temperature fields in the solid and liquid regions of the cast and the velocity fields of the liquid metal at the entrance region of the twin-belt caster. The process variables considered are the velocity and cooling conditions of the belt and the nozzle shape. A body fitted coordinate system is used to track the front locations and to solve .the Navier-Stokes equations in the liquid region. Results obtained can help design efficient casters for thin strip steel castings. The development of the numerical model is presented, and sample results are discussed. I. IntroductionContinuous strip casting of steel has attracted world-wide concerns in recent years due to its economic advantages with increased energy savings and the elimination of intermediate processing steps. New methods of thin strip casting of steel are developing. One of the promising concepts is the Hazelett process, which is being used in nonferrous applications. The process incorporates two continuously moving, intensely cooled, thin metal belts. Metal blocks linked together are sandwiched between the moving belts and form moving side mold walls.In this work, the freezing front position, the temperature field and the liquid region flow field will be analyzed for the strip casting of metal using the twin belt caster (the Hazelett process). Up to now almost all theoretical analyses undertaken of the continuous casting of metal [1-3] have used a heat conduction equation, assuming slug flow in the liquid region. The effect of convection in the liquid pool during continuous casting were first investigated by Szekely and Stanek [4]. They considered potential flow between the feed metal orifice and the liquidus front, and slug flow in the mushy and solid regions of a vertical static continuous casting process for steel. Viscous effects and natural convection affects were not included Kroeger and Ostrach [5] studied the effect of natural convection flow in the liquid pool during solidification for a pure metal."

Jan 1, 1986

By Raj K. Rajamani, Nikhil Dhawan

"Agglomeration has become common pre-treatment step to improve heap leach percolation thereby recovery, the fundamental understanding of the agglomeration process for crushed ores is still lacking. Agglomeration is required for ores which either contain excessive amounts of clay or an excessive quantity of fines generated during mining and crushing. Numerous experiments have been carried out in laboratory scale batch drums using different ores. It has been found out that experimental agglomerate size distributions exhibit self-preserving spectra. The self-preserving/ similarity distribution of the agglomerates is found to be independent of operating conditions and is uniquely determined by the size enlargement mechanism. An empirical model based on the absolute difference between moisture applied and optimum moisture was developed. This model uses self-preserving spectra in the form of a Rosin-Rammler distribution function to predict size distribution at different operating conditions. The empirical model was considerably easier to work with plant scale agglomeration drums. Also, a test known as dip test is developed which improves fundamental understanding of the agglomerate formation and growth mechanisms. This test also provides the progeny size distribution of each agglomerate size. The progeny size distribution was analyzed with a model that partitions the particles into a host and guest category. The ensuing partition coefficient identifies the behaviour of a particle as a host for other particles or a guest layering on other particles. This concept helps in identifying the type of agglomerate that would result for a specific feed size distribution. INTRODUCTIONHeap leaching is one of the hydrometallurgical processes that have served as a green and economically viable process for the treatment of low grade complex ores of copper, uranium, nickel, silver and gold (Dhawan et al., 2013). Heap leaching appears to be a favorable candidate for the least energy consuming process for low grade ores. In heap leach operations nowadays a pretreatment step known as crushed ore agglomeration has been implemented (Dhawan et al., 2013). A bed of agglomerates improves the uniform percolation of solution through the heaps of ore. Therefore, it was possible to treat wastes and mill tailings via agglomeration and heap leaching. Agglomeration pretreatment is required for ores which either contain excessive amounts of clay or an excessive quantity of fines generated during mining and crushing (Bouffard, 2005; Dhawan et al., 2013).While agglomeration has become common practice pretreatment step in heap leach operations, the fundamental understanding of the crushed ore agglomeration process is still lacking. The author reviewed publications up to 2012 on the crushed ore agglomeration process and did not come up with modeling studies on the topic (Dhawan et al., 2013). Some of the possible reasons can be the stochastic nature of the process, experimental difficulties, and lack of instrumentation for fundamental studies."

Jan 1, 2016

By Mingming Zhang, Ramana G. Reddy

A 3-D mathematical model was developed for the batch reactor of low temperature aluminum electrowinning using ionic liquid electrolytes. This model describes the deposition process by incorporating the mass transport of participating ionic species, homogeneous chemical reactions within the diffusion layer, and the associated electrochemical kinetics. Processing parameters, current and potential distribution, species concentration, fluid flow distribution, and electrode spacing were evaluated for the optimal reactor performance. The results indicated that the electrode spacing significantly affects the electrolyte fluid flow and current density distribution. The parallel electrode configuration (in line with electrolyte inlet) improved the convection and resulted in uniform current density distribution and electrolyte fluid flow. However, for this electrode configuration, electroactive species distribution was most favorable between the electrodes. Perpendicular configuration of electrodes resulted in a more non-uniform fluid flow within electrolyte domain, and it has a potential to cause non-uniform deposits. Aluminum electrowinning experiments were conducted using batch reactor at 80 °C, electrolyte flow rate of 5–20 ml/min, and applied cell voltage of 3–3.5 V. Good agreement was obtained between the model and the batch aluminum electrowinning experimental results.

Mar 1, 2017