Search Documents

By Baozhong Zhao

In-situ vitrification (ISV) is a new technology used for treating radioactive, organic and inorganic contaminated soils. In this process, electricity is applied through electrodes buried in the contaminated soil to melt it and form an environmentally stable glass-like solid. The organic contaminants and volatile metals are vaporized during heating up, the non-volatile inorganic elements are dissolved and incorporated into the melt. In present full-scale operation, the ISV process can treat 4 to 6 tons of soil per hour at consumption of about 1000 kwh per ton of soil. In this study, various configurations of electrodes were examined in a mathematical model. The results may be used to optimize the ISV process.

Jan 1, 1994

By M. E. Huckman

In this paper, we review the equations needed to construct an ion exchange column model. These models contain three major components: 1) a differential material balance around the colunm, 2) a differential material balance around an ion exchange particle, and 3) a constitutive equation for diffusive flux We then examine the ability of a fairly simple mathematical model to predict ion exchange column experimental data. This model contains the traditional absorption theory equations and is used to predict the behavior of a column filled with a hydrous crystalline silico-titanate and which is removing trace amounts of Cs+ from a high pH, complex electrolyte solution. Finally, we propose replacing Fick's law with the Nernst-Planck equation as the constitutive equation for diffusive flux We use a Nernst-Planck type model to predict batch experimental data, then evaluate it by comparing it's ability to predict experimental data to that of a more traditional Fick's law type model? The data from this comparison suggests that the Nemst-Planck type model yields better results.

Jan 1, 1996

By D. Gu, R. G. Anthony, C. V. Philip, M. E. Huckman

"In this paper, we review the equations needed to construct an ion exchange column model. These models contain three major components: 1) a differential material balance around the column, 2) a differential material balance around an ion exchange particle, and 3) a constitutive equation for diffusive flux. We then examine the ability of a fairly simple mathematical model to predict ion exchange column experimental data. This model contains the traditional absorption theory equations and is used to predict the behavior of a column filled with a hydrous crystalline silico-titanate and which is removing trace amounts of Cs+ from a high pH, complex electrolyte solution. Finally, we propose replacing Fick's law with the Nemst-Planck equation as the constitutive equation for diffusive flux. We use a NemstPlanck type model to predict batch experimental data, then evaluate it by comparing it's ability to predict experimental data to that of a more traditional Fick's law type model. The data from this comparison suggests that the Nemst-Planck type model yields better results.IntroductionMillions of gallons of strongly basic, radioactive waste are currently stored in underground storage tanks at Hanford and Savannah River and several locations within the United States. Due to the limited life expectancy of the tanks, there is a pressing need to remediate these wastes. This need has prompted the recent development of new ion-exchange materials for removing radioactive elements from aqueous solutions. A novel material called TAM-5 has been developed jointly by Texas A&M University and Sandia National Laboratories. TAM- 5, a crystalline silico-titanate, is highly selective for Cs + in the complex electrolytic solutions typical of these wastes. It is stable in very basic solutions and in radioactive environments [1, 2 ,3]. In laboratory testing to determine the selectivity of various materials, TAM-5 was superior to other candidate ion exchangers [4, 5]. The selectivities were determined by measuring a distribution coefficient, Kd, which is defined as the concentration, gig or mmols/g of Cs in the solid phase divided by the concentration, mg/L or mmols/mL, of Cs in the liquid phase. The units of Kd are ml/g or L/kg. Table 1 shows Kd values for TAM-5 compared to 9 other candidate materials in a typical waste solution [6]."

Jan 1, 1996

By Luís Marcelo Tavares, Henrique Dias Gatti Turrer, Luciana Pereira Alves

Desliming is a method of classifying mineral particles, which aims at removal of slimes so as to prepare the slurry to downstream concentration processes. It is a common item in mineral processing flowsheets, being mostly carried out in hydrocyclones, in which separation is result of complex fluid flow and its interaction with the mineral particles. When dealing with particles contained in size ranges relevant to desliming, the challenge in describing the process is even greater. The work deals with the fitting of the Narasimha-Mainza model of classification in hydrocyclones to data from desliming an itabirite iron ore, for prediction of feed slurry flowrate, corrected cut size and short-circuit to the underflow. As such, tests were conducted at both pilot and industrial plants in KREBS 2.7” and 4” cyclones operating under a variety of conditions. In the validation of the model, deviations between the predicted and measured values were lower than 34% for feed slurry flowrate, less than 22% for corrected cut size and less than 18% for short-circuit for underflow.

Oct 3, 2019

By Barda K. Mishra, Surya P. Mehrotra

"INTRODUCTIONIn a jig, which is primarily a coarse mineral concentration device, the bed of mineral particles is pulsated with a stream of water resulting in the particles of different specific gravities. It operates in a cyclic manner – the jig cycle consists of four stages, namely, inlet, expansion, exhaust and compression. In the inlet stage the bed lifts up en masse. Near the end of the lift stroke the particles at the bottom of the bed start falling resulting in the loosing of the bed which, in turn, causes its expansion or dilation. During the third and fourth stages the particle resettle through the fluid, and the bed collapses back to its original volume. The pulsation and suction are repeated to bring about stratification with respect to specific gravity across the bed height. Various parameters affecting the stratification process include amplitude and frequency of pulsation, bed thickness, rate of hutch water flow, and feed characteristics. This paper first presents a review of various jigging theories and then describes a new approach to mathematical modeling of jigging employing the discrete element method.JIG MODEL BASED ON DISCRETE ELEMENT METHOD (DEM) APPROACHThe jig models based on DEM essentially involve tracking the position of each particle incrementally, applying Newton’s law of motion and a force displacement law at the contact. The model presented here is based on a three dimensional mathematical description of solid-solid and fluid-solid interactions. In this model, the force balance on a particle is extended to account for the effect of continuously varying fluid on the particle. The total force on each particle is made up of a known external body force, fluid particle interaction force and particle-particle and/or particle-wall interaction force. The fluid-particle interactive force is primarily made of drag and buoyancy. Di Felice’s correlation is used to evaluate the drag force. Solid-solid interaction is described by treating particles as discrete elements that are indestructible but interact with each other and any other elements present, such as walls, when they make physical contacts. A pair of spring and dashpot represents a contact The contact forces are computed from the amount and the rate of overlap of the particles at the contact. The motion of the particle is due to movement of the fluid, which transmits its energy through the particle-particle contacts. A computer program incorporating both the contact and the fluid model has been developed. It takes the operating parameters such as amplitude and frequency of the pulsating fluid, pertinent data concerning particle properties and the geometry of the jig bed as inputs."

Jan 1, 2003

By Edgar E. Vidal, Patrick R. Taylor

Standard textbook equations for heterogeneous kinetics typically utilize a constant driving force to solve for rates and times of reactions. In order to more closely model real systems, equations and solutions are presented for particle leaching kinetics in which the driving force is decreased as the reactions proceeds. In addition, models for continuous co-current and continuous counter current aqueous-solid reactions are presented.

Jan 1, 2005

By A. McLean, L. Gao, Z. Shi, Y. D. Yang, K. Chattopadhyay, G. F. Zhang

Levitation melting of metals is one of the emerging applications of magnetic and electric fields in liquid metal refining. The varying magnetic field generates induced current inside the metal droplet and leads to the Joule heating effect and Lorentz force against gravity. This contactless process can avoid contamination from the crucible and can be used to produce high-purity materials. Visualization of the levitation melting process is difficult, and thus generally investigated by mathematical modeling. In the present study, a three-dimensional mathematical model of metal levitation was developed using finite element method for proper understanding of the levitation process. The effects of harmonic magnetic field frequency, coils power input, and sample size on levitation nickel droplet were investigated. Temporal evolution of temperature fields were predicted directly though the model. Levitation and balance condition of the droplet were investigated as well.

Jan 1, 2015

By John J. Jonas, Sang-Hyun Cho, Ki-Bong Kang

"The recrystallization behavior of Nb micro alloyed steels containing Cr, Mo, Ni and Ti. was studied using hot torsion testing with the aim of modeling the recrystallization processes and the mean flow stress taking into account the microstructural behavior during hot rolling. The kinetics of static and metadynamic recrystallization were characterized and appropriate expressions were formulated for the recrystallization kinetics. These are shown that the recrystallization kinetics depend on steel composition and the processing conditions. The stress-strain curves were determined in order to derive new equations for the peak stress, peak strain, mean flow stress and softening kinetics. The peak strain is influenced by the presence of alloying elements; their addition, which has a solute drag effect, retards the rate of grain boundary motion, shifting the peak to the right. The addition of Nb to the steel results in a significant increment in the activation energy for hot deformation, and Mo has smaller effect than Nb, but Cr has the opposite effect. It was also found that the deformation activation energy in these steels was not altered by the addition of Ni. The rate of metadynamic recrystallization increases with strain rate and temperature and is observed to be independent of strain, in contrast to the observations for static recrystallization.The mean flow stress (MFS) calculations, based on the Sims analysis, are then compared with the predictions of a model based on an improved Misaka mean flow stress equation in which solute effects, strain accumulation, as well as the kinetics of static and metadynamic recrystallization are accounted for. Good agreement between the measured mean flow stresses, which is obtained from the mill logs of POSCO Kwangyang #3 HSM (hot strip mill), and 'predicted ones is obtained over the whole range of rolling temperatures. The results demonstrate that dynamic recrystallization followed by metadynamic recrystallization, supported by the unexpected load drops, sometimes occurs for the rolling schedules studied in the present analysis."

Jan 1, 2001

By J. R. Bhamidipati, Nagy EI-Kaddah

"Over the years, mathematical modeling has been established as a viable tool for design and analysis of melting and solidification processes. Modeling of electromagnetic casting and melting systems involves three coupled problems -- determination of the containment field, the shape of the meniscus, and the temperature field -- which have to be solved simultaneously. One of the difficulties in numerical simulation of these systems is in defining and gridding the solution domain as the meniscus shape of the liquid region is not known a priori. This paper describes a methodology for solving these coupled problems in a two-dimensional electromagnetic containment field. This method is based on the solution of the electromagnetic, free surface dynamic, and heat transfer equations in a fixed domain in computational space by mapping the physical space into an invariant computational space during the search for the equilibrium shape using a variable non-orthogonal coordinate transformation. This approach, also, eliminates the need for regridding the liquid and solid regions during phase change. This methodology is demonstrated for two electromagnetic confinement systems -- the Magnetic Suspension Melting process and the electromagnetic casting of aluminum. The model predictions are in good agreement with available data for these systems. It is suggested that the methodology presented here is likely to provide a useful tool in the design, analysis, and optimization of electromagnetic confinement systems. IntroductionDuring the past decade, mathematical modeling has emerged as a viable tool for the design, analysis, and optimization of casting and melting processes /1-2/. It is generally recognized that the usefulness of modeling to study the interplay between the various process parameters without resorting to extensive experimentation relies heavily on accurate representation of the physical phenomena in these systems. In modeling electromagnetic confinement casting and melting systems, it is important that the model address all aspects of electromagnetic field interaction with the metal, if accurate predictions are to be obtained /3/. The electromagnetic field affects the heat transfer phenomena primarily through joule heating. In addition, the electromagnetic forces play a key role in heat redistribution in the molten pool via melt stirring, and, in confinement systems, support the molten metal pool."

Jan 1, 1994

By Olusegun J. Ilegbusi

"A two-fluid model is used to calculate flow and heat transfer characteristics of a plasma jet. This model accounts for both gradient-diffusion mixing and unmixing that may result from pressure-gradient-body-force imbalance in the system. The unmixedness is considered to be essentially a Rayleigh-Taylor kind instability. Separate transport equations are solved for the plasma and ambient gas velocities, temperatures and volume fractions. The two fluids allowed to exchange mass, momentum and energy through empirical interaction parameters. The empirical coefficients are first established by comparison of predictions with available experimental data for shear flows. The model is then applied to an Argon plasma jet ejecting into stagnant air. The predicted results show the significant build-up of unmixed air within the plasma gas, even relatively far downstream of the plasma torch. By adjusting the inlet condition, the model adequately reproduces the available experimental data on the temperature profile.IntroductionPlasma jets have important applications in materials processing, some of which include spray deposition, melting and refining, heat treatment and materials synthesis. In a typical reactor, the hot plasma stream entrains a surrounding gas and the resulting heat and mass transfer and the condition of operation of the torch, determines to a large extent, the performance of the unit. The entrainment may prevent uniform mixing and produce regions of unmixed hot/cold gases and non-uniform deposits as a result of insufficient melting of deposit in cold regions. This phenomenon may also affect chemical reaction rate in plasma systems for NOx reduction in exhaust gases.The study of plasma phenomena is often conveniently divided into three parts namely, the plasma torch, the plasma jet and analysis of the particles carried in the jet. The present work concerns processes occuring in the near-field of the plasma jet i.e. just downstream of the torch."

Jan 1, 1994

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