Search Documents

By Roberto Parreiras Tavares, Frederico da Costa Fernandes, Guilherme Antonio Defendi, Vangleik Ferreira da Cruz, Júlio César de Sousa Pena

"Near-net-shape casting is an important area of research in the iron and steel industry today. Among the near-net-shape casting processes, the single-belt caster seems to be the most promising, since it is the only one that can achieve levels of productivity similar to those obtained with conventional continuous casting machines. One of the most important aspects in the single-belt casting is the feeding system used to deliver liquid metal over the belt. This system determines the velocity distribution and the levels of turbulence of the liquid metal and affects the quality of the strips. In the present work, a metal delivery system having a simple geometry was proposed. Using the CFD code CFX 5.6, a mathematical model to simulate three-dimensional turbulent fluid flow and heat transfer in that system has been developed. This model was used to analyze different configurations of the metal delivery system and the free surface shape of the fluid was also determined. A physical model of the feeding system was built and used to validate the predictions of the mathematical model. The fluid flow calculations have been validated by dye injection experiments and laser doppler anemometry. Measurements of the free surface levels have also been performed. These experiments provided supporting evidence for the predictions of the mathematical model.IntroductionTo maintain its competitiveness in the present market, the metallurgical sector, particularly the steelmaking companies, are developing new techniques of production. Some of these new processes are focussed on the strip casting technologies. Among these technologies, the SingleBelt process is considered the most interesting option, reducing the number of rolling steps and, consequently, the final cost of the product, especially the part associated to the energy consumption (fuel and/or electric energy). This process can also reach levels of productivity that can not be matched by other strip casting methods(1). One of the most important aspects concerning these new processes is related to the liquid metal delivery system. This system can have a significant effect on the final quality of the strips, since it determines the velocity distribution and the levels of turbulence of the liquid metal."

Jan 1, 2004

By A. K. Kyllo

A study has been carried out on the dusts from Hudson Bay Mining and Smelting's Flin Flon copper smelter. The aim of the work was to determine the mechanisms of dust formation and the deportment of the various constituents. The work concentrated on the copper converters, but the electrostatic precipitators and baghouse also clean the reverberatory furnace off gas, so some mixing will have occurred. The dust composition was found to vary with distance from the furnace, and to be made up of three major components: solidified droplets, unreacted flux, and condensed particles. The solidified droplets were formed from material ejected from the converter during blowing, while the flux was blown out of the converters immediately after charging. The condensed particles contained the majority of the lead and zinc, as well as the volatile impurities. There was also a significant soluble content of the dusts, increasing with distance from the furnace. The condensed material may be separated from the remainder by a simple washing and sizing technique.

Jan 1, 1996

By R. D. Prengaman

Wrought lead-calcium-tin anodes developed by RSR Corporation have improved the performance compared with the other lead alloy anodes. These anodes have been in continuous service at various locations around the world for as long as fourteen years. Calcium provides mechanical strength to the anode to prevent the deformation while tin increases the mechanical strength, reduces the rate of corrosion, prevents the formation of non-conducting layers on the anode surface, and dramatically improves the conductivity of the anode. Tin is also instrumental in developing the unique rolled structure of the anode which leads to the reduced rates of corrosion, longer life, and improved oxygen evolution. Cobalt additions to the electrolyte further enhance the oxygen evolution and improve the anode life. Higher current densities and higher operating temperatures reduce the anode life. The unique patented method of attaching the rolled anode sheet to the bus bar substantially reduces the resistance of this important joint and offers the substantial savings in power over the life of anode.

Jan 1, 1999

By Andriy Burbelko, Daniel Gurgul, Wojciech Kapturkiewicz, Marcin Gómy

"The mathematical model of the globular eutectic solidification in 2D was designed. The proposed model is based on the Cellular Automaton Finite Differences (CA-FD) calculation method. The model has been used for studies of the primary austenite and of globular eutectic grains growth during the ductile iron solidification in the thin wall casting. The model mentioned above takes into account, among other things, non-uniform temperature distribution in the casting wall cross-section, kinetics of the austenite and graphite grains nucleation, solute diffusion and non-equilibrium nature of the interphase boundary migration.Some aspects of stochastic nature of the solidification processes are taken into account. The results of the CA solidification modeling have been shown, that differences between cooling curves obtained in the same conditions of the casting solidification may be a result of the stochastic nature of the grain nucleation.The results of modeling were compared with the experimental microstructure. The model is well worked for thin walled ductile iron casting, for fast solidification process.IntroductionNodular graphite cast iron, also known as ductile iron (DI), has major applications in critical engineering parts due to its mechanical properties and castability. The mechanical and physical properties of this material depend on the shape and number of the graphite grains andmicrostructure of the metallic matrix.Solidification of DI was a subject of many computer modeling programs described in literature [ 1-5], in which the stationary conditions of carbon diffusion in austenite is pre-assumed. Recently, a tendency for production of thin-walled castings has been observed [ 6-8]. In this technology, the process of the fast solidification may be very far from equilibrium and steadystate conditions [9].The purpose of the present work is a two-dimension model development for simulation of the DI structure formation during the solidification in the condition of non steady-state temperature and diffusion fields in the thin-wall casting.Model of process The CA-FD is one of the known methods of the simulation ofmicrostructure formation during the solidification [ 10, 11]. In the CA microstructure modeling the outer grain shape is the result of the simulation and does not superimposed beforehand. The model development for a onephase microstructure evolution is a subject of the numerous researches [12 -23]. Model of the eutectic solidification of DI in the uniform temperature field and superimposed cooling rate is known [24]."

Jan 1, 2013

By D. C. Foguel, J. L. R. Bravo

For several years, Caraiba Metais S.A. (CMSA) had to stock the arsenic sludge generated from the electrolyte purification process of its tankhouse. The amount stocked · reached approximately 3,500 tonnes, as the effects of its recycling during the smelting process and consequently in the sulfuric acid plant were unknown. Besides, it had as limiting factors the high content of arsenic in the copper concentrates supplied by foreign countries as well as the reduced capacity of the tankhouse's purification plant. The choice for selling the mentioned residue, which generally contains more than 50% of Cu, was attractive, but the environmental policies and regulations had not allowed this option. Other alternatives, such as the production of bicupric arsenic by means of a solvent extraction process or arsenic trioxide, were considered but the economic feasibility was not attractive due to the absence of a domestic market and a poor international market. After a study of the arsenic partition in the metallurgical process, it was verified that the residue could be recycled without causing any operational problems or damage to the environment. To make this possible, the capacity of the purification plant was increased and the way of removing the arsenic sludge from the cells was modified from manual to mechanical, therefore avoiding the direct contact of the workers with this poisonous substance. Moreover, procedures for recycling in the smelter were established and tolerance limits were set for both the smelter and the sulfuric acid plant. The main amount of arsenic to be recycled is released with the smelting furnace off-gas, which is washed during the purification stage of the sulfuric acid plant. The acid water, rich in arsenic, is then treated in the waste plant and the arsenic, in its stabilized condition, is confined in an appropriate area.

Jan 1, 2000

By Ji-He Wei

Keywords: stainless steel, AOD refining, side and top combined blowing, mixing characteristics, water modeling The fluid mixing characteristics in the bath during the side and top combined blowing AOD (argon-oxygen decarburization) refining process of stainless steel have preliminarily been investigated on a water model unit of a 120 t AOD converter. The geometric similarity ratio between the model and its prototype (including the side tuyeres and the top lances) was 1:4. On the basis of the theoretical calculations for the parameters of the gas streams in the side tuyeres and the top lances, the gas blowing rates used for the model have more reasonably been determined. The influence of the tuyere number and position arrangement, and the gas flow rates for side and top blowing on the characteristics has been examined. The results demonstrated that the liquid in the bath underwent vigorous circulatory motion during gas blowing, without obvious dead zone in the bath, resulting in a well mixing effectiveness. The gas flow rate of the main tuyere had a governing role on the characteristics, a suitable increase in the gas flow rate of the subtuyere could improve mixing efficiency, and the gas jet from top lance could make the mixing time prolong. Corresponding to the oxygen top blowing rate specified by the technology, a roughly equivalent and good mixing effectiveness could be reached by using six side tuyeres with an angle of 27 degrees between each tuyere, and five side tuyeres with an angular separation of 22.5 or 27 degrees between each tuyere. The relationships of the mixing time with the gas blowing rates of main-tuyeres and sub-tuyeres and top lance, the angle between each tuyere, and the tuyere number have been obtained.

Jan 1, 2005

By Taiping Lou

The effects of slag composition, additive agents, and heat treatment conditions on the perovskite (CaO ? TiO2) crystal morphology, the precipitation behavior, the crystal growth and the volume fraction (VF) in the Ti-bearing blast furnace slags were investigated. The morphology of perovskite is dispersed dendriform in molten slags. The crystal growth mechanism of the remelting of fine dendrites and the coarsening of larger grains exists throughout the solidification of molten slags. With the temperature decrease, the VF was increased and Perovskite transformed from fine dendrites to coarsening ones or equiaxed grains. The lower cooling rate can increase the VF and coarsen the grain obviously, while higher cooling rate will decrease the VF and reduce the average equivalent diameter (AED) of Perovskite. When .increasing the content of CaO and Fe20}, the VF can increase obviosly. However, high basicity will give rise to the increasing of viscosity of slag, which result in the reducing of the AED. The experimental results showed that additive CaF2 and MnO can efficiently decreased the viscosity of the slags, and obviously improved the morphology of perovskite and promote growth.

Jan 1, 2000

By Jagdish C. Agarwal

Although copper ranks very low in the order of abundance of elements in the crust of the earth, it has been one of the first metals utilized by mar, at his emergence from the Stone Age, about 10,000 years ago. Use of copper for weapons, tools, utensils and ornaments dates back to the twilight of history, 6,000-5,000 B.C. Native copper outcrops and the malleability of the metal explain its very early use by the cave-dwellers. The working, melting and extraction of copper from its minerals are among the oldest metallurgical achievements of the human race. The oldest complete smelting installations for the extraction of copper from its minerals, dated about the fourth millenium B.C., were recently discovered in the Timna Valley, Wadi Arabah, Israel. (1) The history of the developments of the extractive metallurgy of copper has been adequately reviewed recently. (2,3) The center of copper production migrated from the Eastern Mediterranean to Northern Europe (Mansfeld, Germany, Swansea, Wales) towards the end of the Middle Ages (1500 A.D.). The production of the metal was boosted tremendously in North America since the end of the nineteenth century (Fig. 1).

Jan 1, 1976

By Ji-He Wei

Keywords: stainless steel, AOD refining, side and top combined blowing, mixing characteristics, water modeling The fluid mixing characteristics in the bath during the side and top combined blowing AOD (argon-oxygen decarburization) refining process of stainless steel have preliminarily been studied on a water model unit of a 120 t AOD converter. The geometric similarity ratio between the model and its prototype (including the side tuyeres and the top lances) was 1:4. On the basis of the theoretical calculations for the parameters of the gas streams in the side tuyeres and the top lances, the gas blowing rates used for the model have more reasonably been determined. The influence of the tuyere number and position arrangement, and the gas flow rates for side and top blowing on the characteristics has been examined. The results demonstrated that the liquid in the bath underwent vigorous circulatory motion during gas blowing, without obvious dead zone in the bath, resulting in a well mixing effectiveness. The gas flow rate of the main tuyere had a governing role on the characteristics, a suitable increase in the gas flow rate of the subtuyere could improve mixing efficiency, and the gas jet from top lance could make the mixing time prolong. Corresponding to the oxygen top blowing rate specified by the technology, a roughly equivalent and good mixing effectiveness could be reached by using six side tuyeres with an angle of 27° between each tuyere, and five side tuyeres with an angular separation of 22.5° or 27° between each tuyere. The relationships of the mixing time with the gas blowing rates of main-tuyeres and sub-tuyeres and top lance, the angle between each tuyere, and the tuyere number have been obtained.

Jan 1, 2005

By Donald O. Voit

A BASIC language program has been written that simulates the performance of an integrated solvent extraction circuit consisting of an extractor, a reflux fed scrubber, and a stripper. The program is designed to simulate the performance of a circuit having an aqueous feed containing each of the lanthanides as well as yttrium. The Kremser equation is used to determine the separation occuring in each section of the circuit. The required input variables are the feed composition, the separation factors, the light key extraction factors and extractor feed zone distribution coefficient, the number of stages, and the reflux ratios. The program calculates the composition of the streams at each node in the circuit, the total loading, and the remaining distribution coefficients. User interaction with the program is essential. The program has no capability to determine if the calculated values are consistent with various real restraints. Knowledge of the physical, chemical, and equilibrium behavior is essential to successfully utilize the program. The number of iterations required to achieve steady-state provides insight to the circuit response times. With suitable supervision the program can be used as a simulator for operator training.

Jan 1, 1989

By Y. Fautrelle, F. Felten, O. Metais, Y. Du Terrail

"We deal with the prediction of electromagnetically-driven turbulent flows by means of a Large-Eddy- Simulation method (LES). The model is applied in the case of a liquid metal pool submitted to a polyphase linear electromagnetic stirrer. We investigate two cases: (i) the frequency of the magnetic field is equal to 50 Hz; in that case the effects of the pulsating part of the Lorentz forces may be neglected; (ii) the frequency of the magnetic field is low (f= 1 Hz), then the oscillating part of the electromagnetic forces is taken into account. The LES predictions agree well with the mean velocity measurements, as does the standard k-a model. However, as for the turbulent kinetic energy predictions, there is a large discrepancy between the two models. When the oscillating part of the Lorentz forces is taken into account, the computations show that the fluid flow is sensitive to the unsteady part of the forces. We find that the mean velocity is not affected by the fluctuating component of the force. As for the turbulence parameters, the presence of the pulsating part leads to a significant reduction of the turbulent kinetic energy, whilst the turbulence length scale decreases.IntroductionElectromagnetic stirring is commonly used in various metallurgical processes. In continuous casting for example, it may improve the structure of the solidified metal products. Fluid flows generated by electromagnetic stirrers have been extensively investigated (1 ). As for numerical predictions, many works have been published on the subject. In most of the previous works, onepoint closure models have been used. The most used model is the well-known k-i:: model. Those models yield fairly good results as for the mean velocity prediction. However, it has two main drawbacks:(i) the turbulence is not predicted accurately,(ii) it is not well fitted for unstationary flows as well as transitional turbulent motions, where theReynolds number is moderate."

Jan 1, 1999

By Miaoyong Zhu, Sen Luo

"A two-dimensional CAFD model was developed for simulating the solute diffusion-controlled dentritic growth in the low Peclet number regime. Mass conservation and local equilibrium at the solid/liquid interface were solved to determine the growth velocity of interface with consideration of constitutional and curvature undercooling. The modified decentered square growth proposed by Wang et al was adopted to capture the first eight neighboring cells. The comparison between the classic Lipton-Glicksman-Kurz (LGK) analytical model for steady state tip variables and the simulated values showed that the present model was capable of quantitatively determining the dendritic growth, and the simulation tests demonstrated that the effects of the mesh anisotropy in growth kinetics and crystallographic orientation were eliminated. The numerical simulations were performed to the dendritic growth ofFe-0.6%C binary alloy in undercooled melt, and the results showed that the present model could predict the dendrite morphology with artificial growth orientation in acceptable level.IntroductionDuring the solidification of metals and alloys, the dendrite grows from the substance or the undercooling melt and solute redistributes between the dendrite and the interdendritic liquid. The solidification microstructure is close related to the phase and solute evenness and has a close linkage with the properties of final production.In the past few decades, various kinds of deterministic and stochastic models have been developed to predict the evolution of microstructures in solidification, with the fast development of computer power. Among these methods, Phase Field (PF) model [1, 2] has been known as one of the most adequate deterministic models for directly simulating the dendritic growth. However, there is a vital drawback that the PF method requires considerable computation time and is limited to simulate very small domains with a few dendrites [3]. On the other hand, stochastic models, such as Monte Carlo (MC) method [4, 5] and Cellular Automaton (CA) method [6-8], have been usually applied to predict the solidification microstructure evolution. Comparing with the MC method, the dendritic growth kinetics is explicitly integrated in the CA model, while the MC model doesn't. Therefore, the CA method is more reasonable to predict the dendritic growth during the solidification of melt."

Jan 1, 2013

By Matthias Schnitzler, Herbert Pfeifer, Ralf Bolling

"To transform metals high temperatures are usually needed. A fast and efficient heating can reduce the C02 emitted. The fastest way to heat metals is called direct flame impingement (DFI). The flame impinges upon the item being heated. For this method it is important to know the shape and length of the flame. Laser-induced fluorescence (LIF) technology allows measurements of a two-dimensional field without influencing the reaction or the flow. Commonly used is the measurement of OH-radicals, an indicator of the reaction zone.In the presented study measurements of a turbulent diffusion flame are shown. The position of the reaction zone and the influence of the distance between the item being heated and the burner are investigated. Conventional measurements and LIP-measurements are compared and the potential of DFI to increase the efficiency of heating metals is discussed.IntroductionDuring industrial heating processes a fast and efficient heating technology is required. Two important aspects are the NOx emission and the energy efficiency. The direct flame impingement (DFI) technology has the highest heat transfer rate possible with a gas/air flame. The turbulent high velocity flame jet impinges upon the item being is heated [ 1]. Due to the high heat transfer rate furnace and exhaust gas temperatures are lower. Therefore the energy efficiency is increased and the amount of NOx emitted is reduced.Hydroxyl- (OH-) radicals are an intermediate product that is generated during the combustion of hydrocarbons [2]. They are commonly considered as an indicator for the reaction zone and position of heat release. Figure 1 shows the hydroxyl-mole fraction in dependence of the distance along the gas flow. The figure clearly shows how fast the OH-concentration increases in the reaction zone. Figure 1 also shows that OH-radicals are present in the post combustion region, however in smaller concentration [3].The laser-induced fluorescence technology serves to measure the spatial distribution of OH-radicals in a two dimensional cross-section area of the flame. Therefore the flame dimensions and location can be determined with high spatial resolution without influencing the flow or the reaction [ 4,5]."

Jan 1, 2012

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 Luiz R. P. De Andrade Lima

This paper presents a mathematical model for the oxide copper ore heap leaching process. The data used in the model include physical-chemical, .geometrical and operational data, as: leachable metals ore contents, .flow rate and acid concentration, parameters of passivity, ore size distribution, the average residence time of the solution in the heap, height, irrigated area and weight of ore in the heap. In this model the heap is divided into horizontal and plane layers of constant area. The flow of the liquid in the heap is considered unidirectional and the solid-fluid reaction model is considered under diffusive control. The average residence time of the solution into the heap and the apparent effective diffusivity of the acid through the ore particles have .a significant influence on the rate of the copper extraction and its concentration in the pregnant solution; therefore these parameters are used in the model calibration. When this algorithm was applied to a laboratory test, the results showed that the model is adequate to predict the process performance.

Jan 1, 1999

By C. P. Snyman

Three independent methods of coal classification can be distinguished, namely on the basis of grade, type and rank. The latter two can be determined accurately and unambiguously only by means of microscopical techniques. These three independent variables affect all the coal properties and therefore also the technological behaviour of a particular coal. The influence of the inorganic constituents of coal on the fusibility of the ash and on the length of the slag, and also on fouling and slagging, is discussed. The fractionation of the ash-forming constituents during combustion and gasification of pulverised coal shows that fouling and slagging are not related directly to the fusibility and composition of ash of the coal as a whole. The results of combustion experiments using a drop furnace technique indicate that the ignition temperatures of the light volatiles in the coal, and of the formed char, as well as the burn-out time of the char are more dependent on the rank than on the type or grade of the original coal. High-ash coals, however, apparently ignite at lower temperatures and are more reactive than low-ash coals of similar type and rank. The surface properties of coal, which are important in froth flotation and oil agglomeration, correlate well with the inherent moisture content (and by inference with the heat of wetting and the calorific value) and therefore with the type and rank of the coal. In coal exploration, variations of rank offer the possibility of predicting the presence of an intrusion near the seam, thus contributing to more effective exploration and mine planning.

Jan 1, 1991

By M. Nilmani, J. Sofra

The transverse distributions of the solids discharge velocity and mass flux at the discharge end of vertical and inclined lances were measured in order to detennine the transverse distributions of the axial momentum flux, which is a dominant factor influencing the penetration behaviour of particles into melts during submerged injection. The effects of gas velocity (10 - 30 m/s), solids loading (5 - 10), lance length (0.15 - 1.0 m), lance inclination (vertical downward - horizontal) and particle density (2500 - 8900 kg/m3 ) were investigated, while the lance diameter and particle size were fixed at 10 mm and 220 µm respectively. The particle velocities and mass fluxes were measured by image analysis of particles illuminated in a laser-generated light sheet aligned vertically through the axis of the lance at the discharge point. The particle discharge velocity profiles for vertical downward flows were synunetrical about the centre-line of the lance, with the maximum velocities occurring at the centre-line. The velocities increased with increasing gas velocity, increasing lance length and decreasing particle density in accordance with previous observations. The transverse mass flux profiles were also synunetrical for vertical downward flows, however at the lower gas velocities the peaks in the profiles were displaced from the centre-line of the lance by a fractional diameter of about 0.15. The particle discharge velocity profiles become increasingly asynunetric as the lance orientation was changed from vertical to horizontal, with higher particle velocities present near the top wall of the lance. The mass flux profile also became asyrrunetric as the horizontal orientation was approached, however in this case, the profile became biased towards the bottom of the lance. The effect was much more pronounced for a gas velocity of 10 m/s than for 30 m/s, showing the relatively stronger gravitational influence as the velocity was reduced.

Jan 1, 1996

By Juha Järvi

Outokumpu Zinc Oy has two separate roasting lines with fluidized bed furnaces and waste heat boilers for calcine collection and heat recovery. During the last 30 years, zinc production in Kokkola plant has increased continuously and at the same time the raw materials have changed. Therefore, the waste heat boiler has been under continuous research and development. The waste heat boiler plays a significant role in the process and consequently the reliability and efficiency of the device is important for the effectiveness of the entire process. Many changes have been made to the boiler in the last decades. New radiation tube screens, tube bundles, and baffle-plate constructions have been added to the boiler and also rapping systems for the dust have been developed. The latest improvement to the boilers was a baffle wall, located in the middle of the radiation section roof in the first boiler, and baffle plates in the lower part of the convection section in the second boiler. The new boiler design and optimization were carried out using physical water and air-dust modeling as well as computational fluid dynamics. The modeling work agrees well with experiences in practice. The latest developments have increased the operating hours per year and the effectiveness of the waste heat boiler.

Jan 1, 2002

The entrainment of copper matte in the slag during the slag skimming stage of a Peirce-Smith converter was studied with a dam break model. Experiments were carried out with a range of fluid pairs where the lighter fluid was allowed to overflow a weir and the amount of entrained heavier fluid measured. The resulting entrainment is measured as the ratio of the heavy to light fluid withdrawn. The fluids include oil-water, water- glycerol and salt solutions-water pairs to cover a wide range of densities and viscosities. It was found that the Richardson number, Ri, correlates the entrainment well. For the case of the liquid-liquid interface initially at the same height as the weir, the entrainment was found to vary with Ri-1/3. The transient flow was simulated with a volume of fluid (VOF) code written in-house and used a stairstep approach to represent flow obstacles. The code was validated against the experimental results with good agreement for the variation of the heights of the two fluid with time. The simulation showed that the prevailing mechanism for entrainment is due to a pressure differential created across the liquid in the converter during pouring resulting in the heavier liquid being rotated upwards and over the weir.

Jan 1, 1999

By S. -H. Jeon, Y. -S. Choi

"The three dimensional current distribution on the rectangular cathode is predicted using a boundary element method with a linear element of trigonal type. The effect of cell geometry and anode shape on current distribution is simulated. Calculated results compare well with experimental results. The proper combination of anode size and distance between electrodes increase the current uniformity.The importance of electrodeposition as a fabrication technology in the electronics industry is large and growing. The main advantages of this wet processing method, in comparison to other dry processing methods, lie in its high yield, rapid processing and low cost. In nearly every application of electrodeposition the pursuit of deposit uniformity, especially thickness uniformity, has been technically imperative. The thickness uniformity of the deposit depends largely on the current distribution over the cathode. The current distribution is primarily determined by the geometrical characteristics of the electrodes and the cell. Also the polarization at the electrode surface and the mass transfer affect the current distribution.The current distribution is an important consideration for those involved in designing electrochemical systems and electroplating systems in particular. Although it is important, the common practice in industry is to use trial and error method to design the cell that optimize the uniformity of current distributions in electroplating. The purpose of this paper is to illustrate the useof boundary element method to predict three-dimensional current distribution in rectangular cell.The calculation of the current distribution is possible by various numerical techniques, which include the finite difference method (FDM)1, the finite element method (FEM)2, and the boundary element method (BEM)3,4. These methods of the numerical calculation have various problems, respectively, but the BEM is considered to be the most efficient. This is due to the fact that the BEM has the advantage of solving field equations only on the surface of a calculation domain. A good comparison between the BEM and the FEM for the secondary current distribution in a Hull cell was made by Matlosz et al.5The fundamentals on the theory of current distribution in electroplating are well-established in the literature by a number of authors.6, 7, 8 However, a small fraction of these works have been specially addressed to three-dimensional analysis, and, mostly thel·10 have focused on the problem qualitatively rather than quantitatively.Extending the previous work to three-dimension, we developed the program for the quantitative calculation of 3-dimensional current distribution using a boundary element method with a linear element of trigonal type. Integrals with singularity from the fundamental solution have been analytically computed and it is shown in detail in the Iiterature11. This paper deals with the three dimensional current distribution for the electrode position of copper in a cubic type cell and the effect of anode shape, size and the distance between electrodes on current distribution is studied."

Jan 1, 1997