Search Documents

By George P. Demopoulos

In this paper the production and stability of crystalline scorodite (FeAsO4·2H2O) is discussed both on the basis of previously published work as well as recent work performed in our laboratory. The production of scorodite, especially from arsenic-rich and iron-deficient sulphate solutions offers a number of operational advantages such as high arsenic content, stoichiometriciron demand, and excellent dewatering characteristics. In addition to hydrothermal production of scorodite the work done over the last 10 years in our laboratory has determined that it is feasible to produce scorodite by step-wise lime neutralization at 90°C. The atmospheric scorodite possesses the same structural and solubility characteristics with the hydrothermally produced scorodite. Long term equilibration tests (up to 66 weeks) in the pH range 5 to 9 and temperatures from 22 °C to 75 °C revealed scorodite to undergo slow incongruent dissolutionyielding a highly metastable nano-sized 2-line ferrihydrite phase. The equilibrium arsenic concentration with this phase was 0.35, 0.96 and 5.87 mg/L at pH 5, 6, and 7 respectively at 22°C. In the presence of gypsum slightly lower concentrations were measured, this being more pronounced at pH>7. It was tentatively deduced that a Fe-Ca-AsO4 association (in the form of a phase resembling the mineral yukonite) contributes to the lowering of As concentration. Finally thermodynamic calculations determined that scorodite is stable in the presence of ferrihydriteunder oxic conditions up to pH 6.75 at 22 °C or higher pH at lower temperature and gypsum-saturated solutions. Keywords: Scorodite, Crystallization, Arsenic, Solubility, Stability, Precipitation, Disposal, Extraction, Processing, Metals, and Environment.

Jan 1, 2005

Arsenic occurs naturally in rocks and soils, water, air, plants and animals. It can be further released into the environment through natural activities such as volcanic action, erosion of rocks and through human actions. Higher levels of arsenic tend to be found more in ground water sources than in surface water sources of drinking water. Compared to the rest of USA, western states have more systems with arsenic levels greater than 10 parts per billion (ppb). The US Environmental Protection Agency (EPA) has historically regulated arsenic in drinking water at 50 ppb, but the agency has lowered the maximum contaminant level (MCL) for arsenic to 10 ppb-to be in compliance by January 2006.This paper reviews briefly a) aqueous arsenic chemistry and removal mechanisms and b) the available treatment processes to attain the drinking water standard of 10 ppb by 2006. The combining of an appropriate removal method with an existing system depends on several factors such as nature of arsenic, the water source, existing plant equipment and processes site conditions etc. The paper also includes one or two case studies from the State of Arizona. Key words: Arsenic, Drinking water, Treatment Technologies

Jan 1, 2005

By Suchun Zhang

The alkaline oxidation of refractory gold concentrates containing arsenopyrite and pyrite at ambient temperatures and pressures has been found to be accompanied by the simultaneous dissolution of gold without the addition of cyanide. A detailed study has been made of the kinetics of the oxidation of the pure minerals and of several concentrates. The reaction products of the oxidation of arsenic and sulphur have been monitored using chromatographic techniques. In the case of arsenopyrite, thiosulfate, monothioarsenate, arsenate and sulfite are the principal products while thiosulfate has been confirmed to be the dominant product of pyrite oxidation. Gold, either as added powder or as a constituent of the concentrates has been found to dissolve simultaneously during the alkaline oxidation of both arsenopyrite and pyrite, and this appears to be associated with the in-situ formation of thiosulfate as a lixiviant. The effects of temperature, stirring speed, oxygen partial pressure, alkalinity and pulp density on the oxidation of arsenopyrite and the dissolution of gold have also been investigated. Keywords: Gold, dissolution, arsenopyrite, pyrite, alkaline, oxidation, thiosulfate, thioarsenate

Jan 1, 2005

By C. G. Anderson

Throughout the world, there are many orebodies or materials which have significant value but also contain arsenic. As regulations on the transport, exposure, disposition and emission of arsenic have become more stringent; it has become increasingly more difficult to derive the values from these resources. This paper will outline the fundamentals of alkaline sulfide hydrometallurgy and its successful application to arsenic bearing ores, concentrates and materials.

Jan 1, 2005

By Yongfeng Jia

Lime neutralization and co-precipitation of As(V) with iron(III) is the most widely used and most economical method for the removal of arsenic from mineral processing solutions and effluents. However, the arsenate fixation mechanism in the coprecipitated solids is not fully understood. It was suggested by various researchers to be a ferrihydrite precipitation process during which arsenate is removed by adsorption on ferrihydrite. In this study, XRD, EXAFS and microprobe were employed for the speciation of arsenic and the effect of various factors on arsenic solubility was also investigated. The results indicate that coprecipitation is a more complicated process than simple adsorption resulting in lower arsenic solubilities. Poorly crystalline ferric arsenate was observed for the first time to be present in the coprecipitated solids. In case of lime as base, there exists some type of Ca-Fe-arsenate association. A general term ?coprecipitate? is suggested to describe the precipitated arsenic-bearing solids instead of those widely used specific terms ?arsenical ferrihydrite? or ?arsenic-bearing ferrihydrite?.

Jan 1, 2005

By W. F. Drinkard

Many in the environmental movement have focused on this element. But, nothing they or anyone does will change the amount of arsenic in the world. This report is on actual commercial plants that we developed and operated to recover or stabilize arsenic values from metallurgical dusts from gold and copper pyro operations and arsenical gold bearing ores. When they were available gold and copper values were recovered. Our experience in arsenic extends over thirty years. My first corporation, Mineral Research & Development Corporation was the first in North America to manufacture the then new CCA (copper-chrome-arsenic) wood preservatives. This class of wood preservatives proved to be so effective that it made wood into an entirely new material of construction. It is not uncommon for wood properly treated with arsenical CCA preservatives to last for thirty; forty, and sixty years. Over a hundred years have been projected in some field tests. Hundreds of thousands of trees have been saved. Millions more trees would have been saved had environmentalists not decided that their arsenic phobia was more important than trees. There are several variations of CCA. The most popular is CCA-C. Its formula is: CrO3 . . . . . . . . . . . . . . . . . . . 47.5% CuO . . . . . . . . . . . . . . . . . . . 18.5% As2O5 . . . . . . . . . . . . . . . . . 34.0% Until recent restrictions were placed on the consumption of arsenical wood preservatives, annual consumption in North America was about 150,000,000 lbs./yr. About 64,500,000 of arsenic as the pentoxide was safely consumed this way. World CCA consumption has been about double the above figures or 300,000,000 lbs. as 100% CCA and 129,000,000 lbs/yr as As2O5. Close to 100% of all the arsenic used in wood preservatives is derived from gold or copper processing. The valence of arsenic required for CCA wood preservatives and many other markets is the same valence that is required for stabilization. The economic challenge is how to get there. We have built and operated eleven plants to oxidize arsenic trioxide to its pentavalent form. Most of these have been classic in that they used a common oxidants such as hydrogen peroxide or nitric acid to achieve this.

Jan 1, 2005

A novel, low-pressure adaptation of Praxair?s coherent jet technology utilizing CoJet lances was developed for application at Inco&apos;s Copper Cliff Smelter where oxygen is available in large tonnage quantities but only at relatively low supply pressures (< 193 kPa header pressure). The project objectives were to evaluate the effect of coherent jet technology on controlling oxygen lance accretion (?horn?) formation and to increase oxygen efficiency in the Inco Oxygen Top-Blowing Nitrogen Bottom Stirring (TBBS) MK Reactor used for converting nickel- contaminated chalcocite to semi-blister copper. This paper describes the development and subsequent implementation and testing of two CoJet lances in this application.

Jan 1, 2005

By Lincoln McCrabb

Correct measurement and interpretation of viscosity data is critical for design of any thickening and pumping process that involves a non Newtonian slurry. The nickel laterite pressure acid leach (PAL) process relies on producing a slurry of maximum achievable solids concentration. Small variations of 1-2% solids concentration can alter economics significantly. For the Ravensthorpe Nickel Project, a rheology protocol has been developed which optimizes the solids concentration across a range of ore types. The vane yield stress technique has been utilized to accurately quantify the flow of thickened slurry in the lower shear range. A yield stress of 100Pa has been used as a benchmark to determine the optimum solids concentration achievable utilizing deep cone thickening. Statistical analysis of the data shows critical rheology distinctions between ore type, particle size, lithology and geochemistry. Linear and power law models have been applied to rheology flow curves to reconcile the material transport properties with the anticipated thickener underflow solids concentration. Keywords: Vane Yield Stress, Thickening, Rheograms

Jan 1, 2004

By H. Conrad, K. Jung

"The effect of an external electric field E=5(+) kV/cm applied during the solution heat treatment at 500°-550°C of three commercial Al-Mg-Si alloys (AA6061, AA6022 and AA61l1) on the asquenched resistivity Pw and hardness H.,NULL, w was investigated. It was found that the field gave a significant increase in Pw and in Hv,w, for AA6022 and AA6111, but had no clear effect on AA6061. Thermodynamic analysis gave that the effect of the field on AA6022 and AA6111 was to reduce both the enthalpy ?H, and entropy ?H, of solution. Considering the chemical composition of the three alloys it appears that the influence of the field may depend on the amount of Cu plus excess Si present in the alloy.IntroductionKlypin and Soloviev [1, 2] were the first to report that the application of an external electric field during solution heat treatment (SHT) of Al alloys gave an increase in the as-quenched resistivity and hardness. They attributed these effects of the field to an increase in solubility of the alloying constituents during the SHT. More recently Liu and coworkers [3, 4] also reported an increase in the as-quenched resistivity of Al-Li alloys with application of an electric field during SHT. Further, the present authors [5] found that the application of an electric field during the SHT of the commercial Al-Mg-Si alloys AA6022 and AA61i1 increased their as-quenched resistivity and hardness. Their thermodynamic analysis of the results indicated that the field decreased both the enthalpy ?H, and entropy ?S, of solution."

Jan 1, 2004

By T. -M. Guo, M. J. Braun

"This paper presents an experimental and numerical study of melt flow near the solid-liquid interface and its effect on the interface curvature in a horizontal unidirectional solidification system with succinonitrile. Melt flow is introduced by increasing the ampoule channel thickness. Visualization techniques are employed to observe the flow field and the solid-liquid interface curvature. A numerical model is also developed to simulate the heat transfer, melt flow and phase change in the system. Quantitative information on the flow and temperature fields is generated. A good agreement is achieved between the experimental observations and model predictions. The results could be used for the study of the effect of the melt flow on interface morphology during solidification.IntroductionClassic unidirectional solidification experiments of transparent alloys are conducted with thin film samples sandwiched between two glass slides to investigate the development of interface morphology during solidification [1]. Samples with small channel thickness are usually used to suppress any melt flow due to natural convection and solidification in the channel would then be in the diffusive regime so that the solidification theory can be significantly simplified [2, 3]. Effects of melt flow, in particular due to natural convection, on solidification have also been studied extensively in those systems with large geometrical dimensions [4-6]. These studies indicate a strong effect of the melt flow on interface curvature during solidification [5, 6]. Only limited efforts were made to address the effect of the melt flow on interface morphology in horizontal unidirectional solidification, in particular for the cases when the channel thickness is not thin enough to suppress completely the natural convection [7]."

Jan 1, 2004

By D. J. Donaldson

The PAL process for recovering nickel from laterite ore requires large volumes of thick abrasive slurry to be raised lo high temperature (245°C+) and pressure (40 atm.+). Flow diagrams of the process typically result in facilities that have serious operating and reliability problems. Most of the problems can be attributed to using multiple pump stations to raise the slurry to autoclave pressure, large amounts of flash and live injection steam to heat the slurry to autoclave temperature and several steps of flash cooling to lower the temperature and pressure of the abrasive autoclave slurry down to atmospheric conditions. It is proposed that operating and reliability problems of PAL facilities can be mitigated by using Tube-in-tube heat exchangers and "float type" positive displacement pumps. The subject paper describes an innovative application of such equipment in the facility design of the PAL process. Keywords: Nickel, Laterite, PAL, Heating, Slurry

Jan 1, 2004

By John Frederick Colloghon, O&apos

Bulong remains the only hydrometallurgical process plant in the world that recovers nickel and cobalt directly from high pressure acid leach (HPAL) laterite leach liquors without need for precipitation of an intermediate product (e.g. sulphide, hydroxide or carbonate). Direct solvent extraction (DSX) of both nickel and cobalt in combination with electrowinning of nickel has been operation at Bulong since early 1999. Focus in the first years has been on ramp-up and resolution of commissioning problems. More recently the focus has been on de-bottlenecking and process optimisation. A number of significant process problems have been resolved, the most difficult being the formation of gypsum in the nickel solvent extraction circuit. Sound control of process chemistry is essential to ensure the production of high quality nickel and cobalt products at the lowest possible cost. This paper describes the numerous changes and enhancements to the process chemistry and plant equipment at Bulong and discusses the importance of ore leach chemistry on downstream impurity metal removal. Potential additional enhancements to further reduce unit costs are also presented.

Jan 1, 2004

By Viktor. P. Popov, Oleksandr. V. Derev'yanko, Oleksandr. I. Raychenko

"This paper is devoted to scrutiny problem of obtaining electro-contacts by electro-discharge sintering technology. A new unconventional method of electro contacts production (on Aggraphite base) can yield essential optimization resource.Review of three electro-technological methods of processing contacts materials is presented:- using graphite powder cladded by Ag;- mixture of graphite powder(which is a non-cladded powder) with Ag powder and addition of Fe powder particles;- method of infiltration of melted Ag into graphite piece without deformation and rupture of object.This technological approach can give the possibility to increase of compression strength of electro-contacts materials up to 340 MPa and more under conditions of additional reinforcing components (for example, by Fe).IntroductionThe electric current can be an active ""participator"" of the powder metallurgy technology [l]. At the passage of electric current through powders or some powder mixtures the phenomena including phase transformations arise. Here such effects, which are observed in the Ag-graphite and Ag-Fe-graphite powder mixtures and in the pair Ag-porous graphite, are described.The Ag-graphite composition materials (CM) are used, for example, as electric contacts [2, 3]. Since a material being produced is destined for the contacts it is necessary to take into account the following aspects: commutating current, its character and voltage, active or reactive character of working load, geometrical dimensions of contacts, quantity and frequency of switching on and switching off, pressure on contacts, distance between contacts, method of arc control, surrounding medium, its content, pressure and temperature. Thus, these and, maybe, choice of some other circumstances influence design, materials and production technology of contacts. The range of problems concerned dictates such steps as alloying of contact materials, optimization of structure, selection of kind of graphite, its connection with metal [4-10]. It is expedient to recall the main procedures of the first stage of the conventional technology: dosage of the initial graphite powder and binder-phenolic resin plastic varnish (for example, Bakelite), mixing, densification, thermal treatment, crushing, grinding, sieving, drying, dosage of grained graphite, dosage of the initial silver powder, obtaining of the silver salt solution, settling of salt, filtration, washing, drying of the cladded graphite powder, thermal treatment, obtaining of the graphite powder of particles cladded by silver. The second stage of the convectional technology starts by obtaining of the cladded graphite powder and finishes by annealing of the pressed blanks (Fig. 1)."

Jan 1, 2004

By B. Q. Li, Y. Shu, K Lynn

"This paper discusses the effects of magnetic field and internal radiation on the melt flow and solidification morphology during solidification processing of oxide crystals for optoelectronic and remote sensing applications. The numerical solution of the integral differential equation characterizing the internal radiation and the magnetohydrodynamic equations describing the magnetic and transport phenomena is obtained by applying the combined discontinuous and continuous finite element method. The results show that both internal radiation and an externally applied magnetic field can have important effects on the melt flow and solidification behavior during the melt processing of oxide materials.IntroductionOptical single crystals are widely used in optical and electronic devices. These crystals are often grown from high temperature oxide melts under strictly controlled conditions. As in the well-known single crystal growth processes for semiconductor wafers, thermal gradient in the melt is a major contributing fuctor to the oxide melt pool. Unlike the semiconductor single crystal growth processes, however, the oxide melt is semitransparent and internal thermal radiation plays an important role in redistributing the thermal energy in the melt pool. Our recent experience with pilot-scale plant measurements also shows that the oxide melt is electrically conducting, although its electrical conductivity is considered much lower than the liquid metal or semiconductor melt.The fact that the oxide melt is absorbing and emitting in the thermal radiation frequency range makes it difficult to study the thermal energy transport in the melt pool. This is because thermal radiation is governed by a hyperbolic differential equation, which is different from those describing the transport processes in liquids. Also, our understanding of how thermal radiation interacts with other transport mechanisms such as convection and conduction is rather limited at point in time, in particular in high temperature oxide melts for optical single crystal processing."

Jan 1, 2004

By Scott Samson, Weidong Wang, Shekhar Bhansali, Vandana Upadhyay

"We present an improved design for the two-arm polysilicon thermal actuator. The classic thermal actuator design produces a localized hot spot at approximately the center of the thin actuator arm. This phenomenon limits the actuator efficiency, lifetime, and reliability. Since the actuator deflection is based on thermal expansion of the material, advantage may be gained by making the temperature profile more uniform along the hot-arm length. In our design, the geometry of the thin arm is modified by forming a nonlinear structure near the approximate location of the original hot spot. Electromechanical simulations of this design show an increase in deflection, due to a flattening of the temperature profile along the hot arm. The solid model generation, finite element analysis and simulations are carried out using PEA software. Experimental devices are in the process of being fabricated by using a polysilicon rnicromachining process.1. IntroductionMost active MEMS devices employ some kind of actuating mechanism to increase functionality. Various actuating mechanisms like magnetic, electrostatic, thermal, magnetostatic and piezoelectric have been investigated and incorporated for actuating microstructures. The primary features required from many on-chip actuating mechanisms are large deflection (>10 µm), low actuation voltage, low power.consumption, reliability, fast response, ease of design and fabrication, and small size. Thermal actuation is an attractive mechanism for rnicromechanical devices requiring high actuation forces, large deflections, good linearity and low actuation voltage. Thermal actuators, also known as ""electrothermal actuators"", ""heatuators"" or ""U-shaped thermal actuators"", have been largely used in MEMS based devices to provide both in-plane and out-of-plane actuation [l] .Thermal actuators are also easily compatible with standard CMOS circuitry due to their low driving voltage. However, high power consumption remains a major drawback of the thermal actuation mechanism."

Jan 1, 2004

This paper describes developments in ferro-nickel furnace technology, specifically improvements to cooling methods, furnace controls, and high voltage furnace operating regime. The evolution of these technologies is briefly presented, as context for describing the current state of the art that has enabled ferro-nickel furnace operation at over 75 MW and specific energy consumption of less than 400 kWh/t dry ore. Examples from existing operations are used to illustrate present best practices, and potential future trends are discussed: The paper begins with an overview of the fundamental aspects of pyrometallurgical treatment of lateritic ores that drive ferro-nickel furnace design. Furnace wall cooling methods are related to the requirements resulting from specific furnace process conditions, including slag and metal compositions, as well as arc and bath power. The benefits of high productivity and low specific energy consumption resulting from high voltage (shielded-arc) operating practice are discussed, along with the furnace controls and electrical power train that enable such operation even with captive power generation. This paper is an update of that which the authors prepared for the Infacon X ferro-alloys conference.

Jan 1, 2004

By Li. Jiarong, Zuojian Liang, Qingyan Xu, Baicheng Liu, Shizhong Liu, Hailong Yuan

Single crystal solidification of turbine blades is a key technology for the production of reliable and high efficiency gas turbines. The solidification process of single crystal investment castings causes a number of problems that have to be solved such as the increase of casting defects and complexity of production process. Computer modeling and simulation is helpful to predict the defects and improve the process design. In the paper, a mathematical model for three dimensional simulation of solidification and single crystal growth of investment castings was developed based on CA, mechanism of crystal growth and basic transfer equations such as heat, mass transfer equations. Many factors including constitutional undercooling, curvature undercooling and anisotropy, which had vital influences on microstructure evolution, were considered in the model. Temperature, crystal growth interface and structure defects of investment castings of single crystal superalloy were investigated at withdrawal rates of 4.5mm/min and 7mm/min. Simulated results showed that microstructure evolution and structure defects were predicted reasonably and calculated results were coincident with experimental ones. The study indicated that three dimensional solidification simulation technique is a powerful tool for understanding the fundamental of solidification and the formation of structure defects in single crystal investment castings and finally for optimizing process design.

Jan 1, 2004

By Chris Doyle

From 1997-2000. PT INCO undertook the Fourth Line Expansion Project (FLEX) to increase annual production rates from 90-100 to 150 Mlbs Ni/year (68,000 MT Ni/yr). In 2003, PT INCO nickel production underwent a major sustainable step jump that finally exceeded the targeted rate. This paper presents the changes that led to this achievement. The focus of mining changed. Improving Ni grade was important but only if electric furnace feed met target chemistry so that on-specification slag could be produced. The intensity of diamond core drilling for exploration was dramatically upgraded, mine planning strategy was totally re- vamped and an ore chemistry quality assurance team was formed. In the processing plant, measurement systems were upgraded for better control. Teamwork between sections vastly improved. The final step jump in production occurred when the fifth hydro-electric power generator was added to the furnace grid to allow the power-starved furnaces to cruise in operation at 70+ MW. Keywords: Free cash flow, Ore Chemistry, Measure, Control, Improve

Jan 1, 2004

By Dieter Neuschiitz, Pengfu Tan

"In order to understand the sintering phenomena during iron ore sintering process, a CFD model of iron ore sintering has been developed, combining with the kinetics of 15 relevant chemical reactions with the aim of describing the temperature fields and the gas flow patterns in the sintering bed. This CFD model can be used to simulate the concentrations of the gaseous and solid species, temperatures of the gas and solid, melting and molten fractions of the mixture, gas velocities, temperature zones and residence times of the gas in each temperature zone. IntroductionThe off-gas analyses clearly revealed that the off-gas of a sinter strand does not have a composition corresponding to thermodynamic equilibrium, for example, the CO, CO2 and O2 concentrations are not the mutually consistent equilibrium values. Apparently the chemical reactions in the sinter bed are controlled kinetically. Therefore kinetic modeling of sintering reactions combined with CFD modeling is necessary to understand the sintering phenomena during iron ore sintering process.Several models of iron ore sintering have been developed since 1970 [l-l 4l, each being more refined than its predecessor. The models were formulated as a series of differential equations to describe chemical and physical changes occurring at a point in the bed. These equations were solved numerically in a vertical plane along the strand. The present CFD model, based on the Nath's model [14], has been published in the reference [15]. In this work, some new chemical kinetic models of gas-solid reactions are incorporated in the present CFD model. The concentrations of gaseous and solid species, temperatures of gas and solid, melting and molten fractions of solid, velocity of gas, and residence time of gas in each combustion zone can be calculated by this model."

Jan 1, 2004

By E. Saez, M. Bellet, T. Coupez, O. Jaouen

"The paper presents a 3D finite element solver for non-steady state fluid flow. The thermomechanical simulation package REM3D® has been initially developed to model Stokes' flow for polymer injection. The mechanical module is based on a Eulerian velocitypressure formulation. The spatial discretization uses P 1+IP1 tetrahedral elements. The front tracking module consists of the resolution of a transport equation, thanks to a space-time discontinuous Galerkin formulation. It includes mesh adaptation, which permits a dynamic refinement at the front, thus controling the numerical diffusion. The same approach is used to solve heat transfer. In this paper, we focus on the specific adaptations done to treat Navier Stokes flow. The treatment of inertia is detailed, as well as the implementation of sliding conditions at mould surface, using conservative normals. Industrial applications are presented.IntroductionThe study of the mould filling stage in casting processes requires the ability to simulate complex non steady-state three-dimensional flows. Most commercial software use the volume of fluid (VOF) method, associated with a Eulerian formulation, see for instance the series of conferences ""Modelling of Casting, Welding and Advanced Solidification Processes"" [1-3]. The VOF method [4] consists in solving the conservation equation dF / dt = 0 for a variable F - the volumic fraction of fluid - whose value is 1 in filled regions and 0 elsewhere. Two fluids are actually considered: the molten alloy and a gas in the empty regions. This method suffers generally from numerical diffusion in the resolution of the free surface tracking equation, which is of pure advective type. The precision of the computation is then strongly dependent on the mesh density. Since the mesh is fixed (Eulerian formulation) this implies that in case of complex flow the user must have a priori a fairly good idea of the flow characteristics in order to capture properly the critical features of the flow. In addition, the method can hardly handle the discontinuity of the material viscosity at the interface. To prevent this difficulty, the viscosity value is smoothed around the interfuce, adding once again some inaccuracy. Since the value F = ½ is supposed to represent the free surface, it results in a smeared fluid surface. An alternative to the fixed mesh methods are the methods in which the mesh covers the fluid domain only, as proposed by Navti and Lewis [5-6], or Gaston and Bellet [7]. Based on updated Lagrangian schemes, they must be associated with an ALE formulation (arbitrary LagrangianEulerian) in order to reduce the number of full remeshing operations."

Jan 1, 2004