Search Documents

By Amjad Asad, Boyd Davis, Kinnor Chattopadhyay, Rüdiger Schwarze, Christoph Kratzsch, Donghui Li, Lei GAO

Sodium (Na) and Lithium (Li) are produced using molten salt electrolysis. The electrochemistry of the electrolyte is well-researched; however, there are benefits to understanding the melt flow and implications on it for cell design modifications. The basic configuration of alkali metal cells is the Downs cell. This consists of a central anode surrounded by a cathode, and this geometry was the basis for this mathematical modeling study. The behavior of gas bubbles in molten electrolyte was studied in both Na and Li cells through the use of computational fluid dynamics (CFD) techniques. The distance between the anode and the cathode was varied in the CFD model to ascertain whether strong circulatory flows would change significantly in the cell. The standard k-ε turbulence model was used to mimic turbulent flow, and a two-way coupled Discrete Phase Model (DPM) was adopted to simulate flotation behavior of chlorine bubbles and liquid metal droplets. The liquid metal distribution on the free surface was predicted using the Volume of Fluid (VOF) multi-phase model.

Mar 1, 2017

By Han-Tao Hu

Keywords: RH refining process, flow of molten steel, two-phase flow, gas holdup, circulation rate, mathematical modeling A three-dimensional mathematical model for the flow of the molten steel in the whole degasser during the RH refining process has been proposed and developed with considering the physical characteristics of the process, particularly the behaviors of gas-liquid two-phase flow in the up-snorkel and the momentum exchange between the two phases. The fluid flow fields and the gas holdups of liquid phases and others respectively in a 90 t RH degasser and its water model unit with a 1/5 linear scale have been computed using this model. The results showed that the flow pattern of molten steel in the whole RH degasser could be well modeled by the model. The liquid can be fully mixed during the refining process except the area close to the free surface of liquid and zone between the two snorkels in the ladle, but there is a boundary layer between the descending liquid stream from the down-snorkel and its surrounding liquid, which is a typical liquid-liquid two phase flow, and the molten steel in the ladle is not in a perfect mixing state. The lifting gas blown is rising mostly near the up-snorkel wall, which is more obvious under the conditions of a practical RH degasser, and the flow pattern of the bubbles and liquid in the up-snorkel is closer to an annular flow. The calculated circulation rates for the model unit at different lifting gas rates are in good agreement with the determined values.

Jan 1, 2005

By Rodríguez M. Esperanza

The efficient performance of Pachuca tanks is strongly linked to the suspension of mineral particles, which results from the motion of the liquid caused by injected gas rising, in general, through a central draft tube. This study reports a computational investigation on the effect of operating and design parameters on the suspension of particles. By extending the classical 2-phase drift-flux model to compute the gas hold-up, the three-phase (water-air-particles) system was simulated as consisting of a variable-density liquid plus the solid particles. These two phases were treated as interpenetrated continua using an Eulerian approach to set a turbulent recirculating flow model in 2-D and transient state. The model predicted adequately the measured critical superficial air velocity, uc, needed for maintaining particle suspensions, in pulps with 10-50 wt% solids. Additionally, the model was used to investigate the operation of industrial size reactors, obtaining results that agree well with the general trends reported in the literature, and highlighting important differences with respect to laboratory-scale reactors.

Jan 1, 2009

By Michael Zinigrad

The quality of metallic materials depends on their composition and structure and these are determined by various physico-chemical and technological factors. To effectively prepare materials with required composition, structure and properties it is necessary to carry out a research in two parallel directions: Comprehensive analysis of thennodynamics, kineties and mechanisms of the processes taking place at the solid-liquid-gaseous phase interface during welding processes and development of mathematical models of the specific welding technologies. We have developed unique method of mathematical modeling of phase interaction at high temperatures. This method allows us to build models taking into account the thennodynamic characteristics of the processes, the influence of the initial composition and temperature on the equilibrium state of the reactions, the kinetics of heterogeneous processes, the influence of the temperature, composition, hydrodynamic and thennal factors on the velocity of the chemical and diffusion processes. The model can be implemented in the optimization of various technological processes in welding, surfacing, casting as well as in manufacturing of steels and non-ferrous alloys, materials refining, alloying with special additives, removing of non-metallic inclusions.

Jan 1, 2003

By Michael Zinigrad

The quality of metallic materials depends on their composition and structure and these are determined by various physico-chemical and technological factors. To effectively prepare materials with required composition, structure and properties it is necessary to carry out a research in two parallel directions: Comprehensive analysis of thermodynamics, kinetics and mechanisms of the processes taking place at the solid-liquid-gaseous phase interface during welding processes and development of mathematical models of the specific welding technologies. We have developed unique method of mathematical modeling of phase interaction at high temperatures. This method allows us to build models taking into account the thermodynamic characteristics of the processes, the influence of the initial composition and temperature on the equilibrium state of the reactions, the kinetics of heterogeneous processes, the influence of the temperature, composition, hydrodynamic and thermal factors on the velocity of the chemical and diffusion processes. The model can be implemented in the optimization of various technological processes in welding, surfacing, casting as well as in manufacturing of steels and non-ferrous alloys, materials refining, alloying with special additives, removing of non-metallic inclusions.

Jan 1, 2003

By Y. M. Shneerson, V. M/ Shpaer, E. E. Zhmarin, A. Y. Lapin, E. M. Vigdorchik

This paper describes studies of zinc concentrate oxidizing pressure leaching using the method of mathematical modeling. The concept of a kinetic function was taken as a basis for mathematical modeling of the leaching process. The kinetic function describes dependence of the unreacted component share OJ on the relative time x (the relative time is the ratio of the test duration and the time of complete dissolution). This kinetic function is determined in accordance with the results of batch laboratory tests. The laboratory tests have been carried out on zinc concentrate produced from Kazakhstan zinc-copper ores. Zinc concentrate consisted mainly of three types of minerals: sphalerite (70%), pyrite (13%) and chalcopyrite (10%). Oxidation reactions of sphalerite, pyrite and chalcopyrite have been considered as principal ones. The kinetic characteristics of the process (kinetic function, reaction order, activation energy, duration of complete dissolution of the concentrate components) have been calculated in the course of laboratory tests. The mathematical model was used for estimation of the main technological leaching process parameters (such as process rate, degree of metals extraction, autoclave capacity, heat balance, etc.). Mathematical modeling based on the kinetic function concept allows safe designing of metallurgical equipment and scaling the technology to industrial capacity.

Jan 1, 2004

By Zden[e]k P. Bažant

Due to their heterogeneity, most rocks fracture with a zone of distributed cracking ahead of the fracture front. This makes linear elastic fracture mechanics inapplicable. The present study describes a mathematical model in which the progressive microcracking is taken into account by strain-softening stress-strain relation. In a simpler version of the model, the triaxial strain-softening stress- strain relation is obtained by direct adjustment of the elements of the compliance matrix. This formulation is rigorous only if the principal stress directions do not rotate during fracture formation within the fracture process zone. As a more sophisticated and generally applicable model, the so-called microplane model which is an analogue of the slip theory of plasticity is outlined. Comparison with fracture data for various rocks are given.

Jan 1, 1984

By Sylvie C. Bouffard

The need for a better understanding of the heap biooxidation pretreatment of pyritic refractory gold ores has prompted us to develop a rigorous mathematical model of the process. The hydrological structure of the model takes the form of lateral stagnant pores in contact with a uniformly distributed gas stream and connected to vertical plug flow channels, as suggested by the results of an independent study of column hydrology. A model of pyrite oxidation kinetics was derived based on electrochemical leaching theory and the results of potentiostatic stirred tank and column tests. The growth rates of attached and planktonic iron- and sulfur-oxidizing cells are modeled with separate Monod expressions over three separate temperature ranges corresponding to mesophiles, moderate thermophiles, and extreme thermophiles. Isothermal column tests validate the model simulations. Excellent fits of several leaching indicators (potential, extent of sulfide oxidation, iron concentration, cell numbers) reveal the rate-limiting step to shift from particle kinetics to oxygen gas/liquid mass transfer with increasing temperatures and pyrite head grades. Competition for oxygen between sulfur- and iron-oxidizing microorganisms lowers potentials and slows down pyrite oxidation. This modeling tool may assist engineers in the design and operation of heaps to give more complete gold liberation in a shorter time.

Jan 1, 2003

By B. Moudgil, M. Bogan

A mathematical model, based on first-order kinetics, is derived to describe batch settling data. The model is used to predict the settling behavior of solid suspensions under different conditions. The potential benefits of this model to solid-liquid separation processes are discussed.

Jan 1, 1994

Sedimentation is used to clarify process effluents or to thicken process slurries. The sedimentation cycle of a suspension is presented as a settling curve. This curve is simply a plot of the height of the interface between the settled suspension and clear liquor versus time. Figure 1 shows a series of graduate cylinders which contain a suspension at various stages in its settling cycle. Cylinder 1 shows the suspension completely mixed at time zero. After a brief period of time (see cylinder 2), the suspension is observed to settle from its initial height Ho to height H1, leaving behind a volume of clear (or turbid) liquor, called the supernatant. The interface between the supernatant and suspension is termed the mud line. It is also possible to distinguish at the bottom of the cylinder a layer of sludge representing solids which have settled from the suspension. The interface between it and the suspension is called the sludge line. The sludge line moves upward in the cylinder.

Jan 1, 1992

By Mike Trovant, Stavros A. Argyropoulos

"A fixed grid numerical model for phase change problems is proposed. The model solves for heat transfer and fluid flow, via the coupled solution of the energy, continuity and momentum equations. In addition, the numerical scheme accounts for temperature varying thermophysical properties and determines the shrinkage profile resulting from phase and density change. Simulations are run to determine the effect of varying model parameters and wall cooling rates on the location and shape of the shrinkage void formation in cylindrical casts. Results are shown for both aluminum and magnesium.IntroductionThe mathematical modeling of phase change processes (be they solidification or melting problems) entails solving for many interrelated phenomena. A comprehensive numerical simulation of the casting of a metal, for example, would have to account for: (a) heat transfer in both the solid and liquid portions of the metal and in the mold, (b) the release of latent heat during phase change, (c) fluid flow behavior in the liquid portion of the metal, (d) temperature induced shrinkage caused by liquid contraction and solidification contraction, and (e) the development of thermal stresses with accompanying shrinkage in the solid portion of the cast and mold. A schematic of the relationships that exist between the coupled equations which govern solidification is shown in Figure 1."

Jan 1, 1994

By John C. Edwards

To have the capability to predict the development of localized spontaneous heating within a porous coalbed that is subjected to forced air ventilation or in an otherwise quiescent environment in which buoyancy develops, the Bureau of Mines developed three time-dependent mathematical models, which were used to calculate the temperature increase associated with chemisorption of oxygen by the coal. In each model, spontaneous heating is driven by an Arrhenius first-order reaction between the oxygen and coal. Two models are two-dimensional, and one is one-dimensional. In the first two-dimensional model, a constant-velocity forced convection airflow is specified; and in the other, buoyant flow is allowed to develop in the absence of forced convection. The third model evaluates the airflow from Darcy's law and a specification of the pressure at the surface of a one-dimensional porous coalbed. Numerical computations demonstrated how each model could be used to predict the onset of spontaneous heating when the porous coalbed was subjected to constraints of an imposed internal heat source or a high- temperature airflow. The effects of particle size and coalbed compaction upon spontaneous heating were examined with the third model.

Jan 1, 1990

By Xiaolin Wu

Most coal bump occurrences in the U.S. are directly related to mining operations under strong roofs. An understanding of the behavior of bridging or cantilevering strong roofs over working longwall faces is a key to a safe and productive longwall operation. In this paper, an attempt has been made to study the behavior of strong roof beds using the elastic beam theories. Four structural models are formulated to simulate the locations of strong roofs and their distinct weighting stages. For each model, analytical solutions for the deflection, bending moment, stress, and the strain energy stored in both the coal scam and the strong roof are obtained. Solutions for the critical spans, which are controlled by the caveability of the strong roof beds, are developed. The effects of rocks' differing elastic moduli, under tension and compression, on the actual flexural rigidity, and stress redistribution are analyzed.

Jan 1, 1993

By G. F. Friel, J. C. Edwards, C. C. Hwang

In the implementation of a mine fire alarm and location strategy, the Interpretation of the fire signature Is Important. Smoke and CO are the two dominant fire signatures, and many mine fire detection systems utilize sensors which respond to smoke or CO through a diffusion mode or a mechanical pump mode. The fire signature may be modified during its transport along airways and at crosscuts. The modification occurs through dilution of the contaminants. In this Investigation, computer codes were employed to compute the effects of contaminant dilution on the contaminant signature when a straight flow passage has one dead-end crosscut. The model predicted the airflow within always and the entrainment of contaminants into the crosscut from the airway. One application of the model is to determine the most suitable location for CO detectors for early mine fire detection.

Jan 1, 1995

By F. Chabchoub

A mathematical model has been developed to study the influence of process variables on the horizontal Ohno continuous casting process. The model simulates the heat transfer and fluid flow in a wide plate ingot of pure tin. A two dimensional numerical analysis based on a control-volume finite difference approach, the Semi-Implicit Method for Pressure-Linked Equations algorithm (SIMPLER), and an Enthalpy Transforming Model was used. The shape and position of the solid-liquid interface were calculated at different casting speeds, ingot thicknesses, and cooling device characteristics, as well as various mould temperatures. Results of the mathematical model were compared with

Jan 1, 1992

By Robert C. Bates

Background is given on the various equations used to describe diffusion and convective flow of radon in porous media. Bureau-developed equations for modeling diffusion and Darcy flow through multilayered porous media are described briefly, and examples of their use are given. Analysis results discussed include: The transient effect of pressure changes in different length models, time dependent flux changes in multi-layered models, and radon flux through pinholes in sealants

Jan 1, 1980

By B. Abedian, L. M. Racz, S. R. Berry, R. W. Hyers

"The presented work represents an effort to improve understanding and prediction of turbulent flow inside electromagnetically-levitated droplets. We show that the flow field in a test case, a nickel droplet levitated under microgravity conditions, is a low Reynolds number turbulent flow, i.e. in the transitional regime between laminar and turbulent. Past research efforts have used laminar, enhanced viscosity, and k - E turbulence models to describe these flows. In our study, we compare the RNG algorithm to laminar and k - E turbulence model results. We show that accurate description of the turbulent eddy viscosity µr is critical in order to obtain realistic velocity fields, and that µr cannot be uniform in levitated droplets. In the RNG method there are no characteristic length or time scales associated with the flow, thus allowing such anisotropic features to be captured. We perform calculations and analyses for two cases (1) a small, non-deforming, spherical nickel droplet and (2) a deforming, oscillating droplet. In addition to flow field calculations, we compare numerically-obtained surface tension and viscosity values to those obtained experimentally by the oscillating drop technique.IntroductionElectromagnetic levitation (EML) is a containerless processing technique for conductive materials. Its advantages are the avoidance of contamination or chemical reaction at container walls and the ability to perform long-term experiments in a metastable state. Because of its current limitation to small sample sizes, its primary use is in the laboratory. One application of interest is the measurement of thermophysical properties of undercooled liquid metals [1, 2, 3, 4]. In EML, a sample weighing approximately 1 gram is positioned inside an electromagnetic field generated by alternating currents flowing through an axisymmetric water-cooled copper coil. The coil must be designed appropriately so that thesample levitates in a stable manner.Traditionally, mathematical modeling has been performed in tandem with both earth-and spacebound electromagnetic levitation experiments, in order to determine the effect of the magnetic field on material and flow characteristics [5, 6, 7, 8]. These calculations are important because they allow design, planning, and correct interpretation of experiments. In addition, the calculations are of fundamental scientific interest, as they can be applied to other electromagnetic materials processing applications [9]."

Jan 1, 1999

By A. Navarra

Within conventional copper and nickel-copper sulfide smelters, ladles of molten matte are fed into converters. These converters blast oxygen-enriched air into the matte, thereby eliminating iron and sulfur. The converting reactions are exothermic and, indeed, the converter heat balance is often a limiting consideration on the smelter throughput. If a portion of the matte were fed in solidified (granulated) form, this would support higher oxygen enrichment, and lower volumes of converter offgas, allowing higher throughput. This approach is not applied in conventional smelters, partly because of the copious amounts of water that must be evaporated in typical granulators, as in the Kennecott-Outotec process. The current paper recalls a waterless matte granulator that had been pioneered in the 1990’s, and is applicable to conventional smelters. A mathematical formulation is presented that can be used in the context of simulation-based optimization, to estimate the size and impact of waterless matte granulation.

By Elliot Schwartz, Julian Szekely

"In this paper we describe an experiment that was recently performed on the IML-2 Space Shuttle mission in which the thermophysical properties of molten metals were measured using electromagnetic levitation. In the planning of these experiments calculations had to be performed to predict the levitation and stirring forces, heating rates, sample deformation, and internal flows within the sample. These predictive calculations were essential for both the planning of the experiments and for the re-planning of experimental runs during the actual Space Shuttle mission.Introduction""Materials Processing in the Computer Age"" involves the use of mathematical modeling to obtain a quantitative understanding of the phenomena that govern materials processing operations. Microgravity, or space shuttle experimentation represents a specific case in point here; the very high cost of the space shuttle flights, the very limited flight opportunities that are available, and the tight resources aboard a shuttle clearly mandate that these experiments be very precisely planned and that all possible opportunities for mathematical modeling be fully exploited.The purpose of this paper is to present the results of calculations performed for the planning of a space shuttle experiment involving electromagnetic levitation and wherever possible present a comparison of the computational predictions with both experimental measurements and asymptotic analyses. Finally we shall seek to discuss the practical implications of both the computational methodologies and the potential experimental findings."

Jan 1, 1994

By Edmond Donovan, O&apos

With the increase in complexity of conveyor systems and the inexorable demand for more capacity and availability, the ability of designers to accurately predict performance is becoming increasingly important. Dynamic analysis techniques for modelling wave mechanics in conveyors are finding increased applications. Unfortunately, the reliability of these modelling techniques is dependent on a thorough understanding of the performance characteristics of the components that make up the system. These characteristics include such things as response rate, hysteresis and efficiency. Some parameters such as idler rolling resistance can be easily quantified by measurement. More complicated drive parameters are difficult or impossible to measure, particularly in large horse-power equipment. One method of quantifying or at least bracketing the performance envelope of complicated pieces of equipment is to construct a mathematical model. This can be done quickly and economically, and the influence of various parameters can be determined. This paper sets out to demonstrate by example the usefulness of this technique in gaining an understanding of complicated systems. The model is developed and the influence of a number of variables are examined.

Jan 1, 1996