Search Documents

By J. H. Swisher, E. O. Fuchs

Rate equations were derived to describe the kinetics of internal oxidation of cylinders and spheres. The derived equations for cylinders were checked experimentally by means of sub scale thickness and electrical conductivity measurements on Cu-Cr alloy wires. The properties of the internally oxidized samples were examined with conductivity applications in mind. It was possible to produce uniform dispersions of Cr2O3 in copper with an initial chromium content as high as 3 wt pct. While electrical conductivities only a few pct less than that of OFHC copper were obtained, the Cr2Os particle size and spacing were too large for effective dispersion hardening. T.HE process of internal oxidation has been used widely in basic studies of the permeability of gases in metals. In a review article, Rapp1 has discussed the principles of internal oxidation in considerable detail. From a technological standpoint, internal oxidation is often considered undesirable, since it is a means by which inclusions can be introduced into an otherwise clean material. Another important aspect of internal oxidation is its use as a means of dispersion hardening a material. Broutman and Krock2 discuss this and other methods for making dispersion hardened alloys. The only internally oxidized material known to the authors which is commercially available is a Cu-BeO alloy.3'4 This alloy is made from Cu-Be alloy powder, using a so-called Rhines pack. It has a tensile strength of 80,000 psi and retains its strength at relatively high temperatures. The objectives of the present study were to derive rate equations for the internal oxidation of cylinders and spheres, to check the derived equations for cylinders experimentally, and to examine the structure and properties of internally oxidized Cu-Cr alloys. The Cu-Cr system was chosen for this study because uniform dispersions are obtainable at high alloy contents, which is a desirable characteristic in dispersion hardened materials. RATE EQUATIONS FOR VARIOUS GEOMETRIES A number of authors5--9 have derived equations to describe internal oxidation kinetics. These derivations differ somewhat in mathematical assumptions and approximations, and all except one of the derivations deal exclusively with the internal oxidation of plates. The exception is a brief treatment of cylindrical and spherical geometries given by Meijering and Druy-vesteyn9 as a part of a comprehensive paper on the general subject of internal oxidation. These authors did not obtain rate data to check their derivations, although they did show that the hardness profile across an internally oxidized sample is directly related to the rate of interface movement. For cylindrical and spherical geometries, a quasi-steady-state approximation is needed to circumvent mathematical complications in obtaining a solution to the basic differential equations. In using this approximation, we consider the concentration gradient of dissolved oxygen in the internally oxidized zone or sub-scale to be the same as the gradient which would be present if there were no movement of the subscale interface. The steady-state approximation introduces an error of about 1 pct in computing the rate of internal oxidation of an Fe-1.0 pct Mn alloy plate, if the present method is compared to the more exact method of Wagner.7'10 The details of the derivations of the rate equations for cylinders and spheres are given in the Appendix, and only the results of these derivations are given below. The final equations obtained by Meijering and Druyvesteyn9 can be shown to be equivalent to our Eqs. [1] and [2], although the two approaches are somewhat different. Cylindrical Geometry. [2] where r1 is the outer radius of the cylinder or sphere, cm, r2 is the radius of the unreacted core, cm, see Fig. l(a), D is the diffusion coefficient of oxygen in copper, cm2 per sec, %O is the concentration of dissolved oxygen at the surface of the specimen, wt pct, %Cr is the initial chromium concentration in the alloy, wt pct, and t is the reaction time, sec. Plate Geometry. The analogous rate equation for a plate has been derived previously for internal oxidation of Fe-Al alloys.8'11 For Cu-Cr alloys, we may write the same equation as follows: [3] where r1 is the half-thickness of the plate, cm, and r2 is the distance from the mid-plane to the subscale intherate is An analysis of Eqs. [1], [2], and [3] shows that for a plate the rate is completely parabolic. The initial

Jan 1, 1970

By Frank E. Woolley, John F. Elliott

Mathematical models of the transient thermal behavior of a high-temperature solution calorimeter1-3 have been developed. The thermal behavior of the calorimeter is appoxirrzated by linear lumped-parameter models, and hence is described by sets of linear ordinary differential equations with constant coefficients The response of the models to various inputs is shown to agree with the response of the real system. Application of the modeling to experimental design and analysis of data illustrates the usefulness of simple models of complex systems. The early eperiments1,2 with the high-temperature solution calorimeter indicated that the change in the temperature of the bath resulting from the addition of a solute sample to the bath involved not only the direct effect due to the solution process but also possibly a secondary effect arising from the change in coupling between the bath and the induction heating coil. Consequently, an extensive analysis of the calorimeter was carried out, and models of the transient thermal processes of the instrument were developed to aid in improving the design and interpreting the behavior of the system. This paper describes the dynamic modeling; the use of it in treating experimental results has been reported earlier.3 The high-temperature solution calorimeter was constructed to measure directly the partial molar heats of solution of solute elements in a variety of liquid metal solvents.1-3 The calorimeter consists of an induction-heated liquid metal bath into which small samples of a solute element can be dropped. The bath temperature is recorded continuously, and the change in the measured bath temperature with time, dTm = f(t), resulting from the solute addition are the raw data from which the enthalpy change caused by the addition is determined. To extract the rmodynamic results from the data, the temperature change must be compared with that resulting from calibration additions of known enthalpy change. Accordingly, it is necessary to understand the transient thermal processes arising as a result of the addition to the bath. Neither modeling nor experimentation alone could provide the required insight into the working of the calorimeter. The alternate use of both methods in conjunction greatly assisted the design of the equipment and experiments, and the interpretation of the data. THE PHYSICAL CHARACTER OF THE SYSTEM The essential parts of the calorimeter, Fig. 1, for model studies are the thermocouple, the liquid metal bath and the surrounding refractories. The system is the solvent metal bath and those refractories around it which undergo a temperature change as a result of an addition to the bath, and which determine the way the temperature of the bath responds to an input. The inputs are the combined transient thermal effects arising when an addition is made to the bath. They include the thermal effects of the addition itself and the results of changed coupling between the bath and the induction coil. The response is the variation in the measured bath temperature, dTm(t) = Tm(t) - Tm(O), from an initial steady state resulting from the inputs. It was assumed in this study that the physical properties of the various elements of the system are independent of the inputs and time, although these properties may vary as the result of changes in the composition and size of the bath during a series of additions. This separation of inputs and the system is equivalent to assuming that the system is linear, i.e., that its behavior can be described by linear differential equations with constant coefficients. Linear behavior can be expected whenever the departure of each portion of the system from its steady-state condition is small enough to cause negligible changes in the thermal properties of the materials and in the various heat-transfer coefficients. Radiative heat transfer is important in this system, so the assumption of linearity should be valid only for small temperature deviations. Several conclusions were drawn from operation of the calorimeter in earlier experimental studies: 1) Radiative heat transport from the top of the bath is a significant portion of the total heat lost from the bath. However, for small changes in the bath temperature the change in transport by this path could be assumed to be proportional to the change in the bath temperature. 2) A very small portion of the heat input is lost through the thermocouple to its water-cooled holder. The thermal resistance and thermal capacity of the thermocouple protection tube are small, so the temperature of the thermocouple should follow closely that of the bath. 3) The remainder of the total heat lost from the bath will pass by conduction through the crucible to, and through, the other refractories, eventually being absorbed by the water-cooled induction coil or by the water-cooled sides and bottom of the enclosure. 4) The thermal resistance between the bath and crucible is very small. Thus the thermal capacity of the crucible will affect the temperature of the bath very soon after an addition of heat to the bath. 5) The thermal resistance between the crucible and the silica sleeve is large, especially if a radiation shield is placed in the gap. The effect of the thermal capacity of the sleeve thus will be significant only at longer times. The thermal resistance through the packing below the crucible also is large, so the packing and the silica sleeve will have similar effects on the behavior of the system. 6) A large temperature drop exists across the gap containing the water-cooled induction coil. Thus for relatively small changes in the thermal input to the bath, the refractories beyond the sleeve

Jan 1, 1970

By J. Szekely, G. P. Martins

A formulation is given and a general analytical solution is presented for the problem of spherical phase growth in multicomponent systems. The results are valid for describing the rate of bubble growth in the absence of translational motion and, therefore, should be relevant to heterogeneous nucleation processes. Typical worked examples are given, which include bubble growth due to the simultaneous diffusion of hydrogen and nitrogen in molten steel, and the growth of CO bubbles in steel, due to the reaction: C + 0 = CO. The calculations predict that the presence of additional components would enhance the growth rate and also increase the absolute rate at which the individual components are transferred into the bubble. These findings should be relevant to degassing processes. THE problem of spherical bubble growth in the absence of translational motion is appropriate to the description of the heterogeneous nucleation and subsequent growth of gas bubbles while attached to a solid surface. Such problems are of considerable practical relevance to various vacuum degassing and steelmaking processes where bubble nucleation and subsequent growth play an important part in the overall reaction mechanism. A definitive formulation of spherical phase growth has been given in the now classical paper by Scriven1 and the analogous problem for solid phase growth has been discussed by Horvay,2 furthermore a thorough review of the relevant literature has also become available.3 All this previous work was concerned with one component or binary systems. The purpose of this paper is to show that the technique developed by Scriven may be generalized to alow the consideration of multicomponent systems. The results of this generalization would be particularly relevant, as gas bubble nucleation in the metallurgical systems previously mentioned frequently involves the simultaneous diffusion of several components. In the following a general formulation is given and typical worked examples are presented. FORMULATION Consider the growth of a spherical gas bubble due to diffusion in a quiescent liquid of infinite extent, as sketched in Fig. 1. The formulation of such problems requires the statement of the following conservation equations relating to the liquid phase: Equation of continuity Equation of motion Conservation of the diffusing species—Diffusion Equation. The following assumptions are made: a) The fluid is incompressible, quiescent, and of infinite extent b) The bubble is stationary, i.e., it undergoes no translational motion c) Spherical symmetry is observed throughout d) The mass diffusivity is constant e) Chemical equilibrium exists at the surface of the bubble, and any reactions occurring there are fast compared to the diffusion process f) Inertia and surface tension effects are neglected. Equations derived on the basis of these assumptions by previous investigators were found to provide a reasonable representation of bubble growth while attached to a surface in organic and aqueous systems4,5 For gas bubble growth in molten metal systems assumptions a) to d) chould hold reasonably well. The exclusion of chemical kinetics as a possible rate limiting step may not be universally valid, but should hold for many systems operating at high temperatures. The effect of hydrodynamic inertia may be neglected for moderate growth rates; finally, the assumption neglecting the effect of interfacial tension on the bubble growth was made for the sake of mathematical convenience; its inclusion would have required the use of quite cumbersome numerical techniques. While surface tension effects could be appreciable on the initial values of the overall growth rate, especially at low absolute pressures and in case of small bubbles, they are not thought to modify the relative effects present in multicomponent systems, the study of which is the objective of this paper. Within the framework of the above assumptions

Jan 1, 1970

By D. de Fontaine

A convolution procedure is described for obtaining an analytical expression for the composition profiles of solute concentration:; in binary solid solutions. The free energy of solutions containing such compositional inhomogeneities can then be evaluated and the stability of spherically symmetric clusters (G.P. zones, for short) us a function of their radii and central compositions can be studied accordingly. Two characteristic zone radii are defined: the minimum radius (RM) corresponding to the smallest zone for which growth is energetically favored over dissolution, and the critical radius RK) corresponding to zones associated with the minimum work of formation. The values of RM and RK coincide at the two classical limits of nucleation and growth and spinodal decomposition. As the spinodal is approached, however, RM and RK diverge increasingly and it is proposed that the initial] decomposition morphology is governed more closely by the minimum, rather than by the critical radius since the former possesses no singularity within the m eta stable phase boundary. Reversion is briefly discussed in the light of the foregoing and it is shown that a set of reversion temperatures located between the coherent miscibility gap and the coherent spinodal must correspond to a given distribution of low temperature zones. SINCE the early works of Guinier1 and preston,2 the kinetics and morphology of GP zones have been studied extensively, see Ref. 3, for example, for a recent review. Nevertheless, a completely satisfying theoretical treatment of the decomposition of meta-stable or unstable binary solid solutions is still lacking due to the almost unsurmountable mathematical difficulties that one encounters in attempting to solve the pertinent diffusion equation which, in principle, describes the whole process. This equation, first derived by Cahn for the initial stages of spinodal decomposition,4 is a nonlinear partial differential equation of the fourth order in anisotropic space, a general analytical solution of which is completely out of the question. From this equation, some prediction can be made, however, concerning the later stages of spinodal decomposition5 and a recent numerical solution6 of the equation in the case of one- and two-dimensional composition modulations of a solid solution inside the spinodal has led to the following considerations: the sinusoidal morphology7 is to be expected only for solutions of average composition c, located near the center of the coherent miscibility gap.8 For solutions with co closer to the spinodal decomposition, depleted shell zones9 should constitute the first recognizable decomposition product, later evolving to iso- lated, enriched preprecipitates in a uniformly depleted matrix as this must, in all cases, constitute the final metastable morphology, barring, of course, loss of coherency and change in crystal structure. In agreement with the qualitative description of Bonfiglioli and Guinier,l0 it thus appears that the various models discussed by previous authors are, in fact, idealized successive stages of the general decomposition process. Unfortunately, corresponding calculations for a solid solution outside the spinodal are not yet available because of two additional difficulties: it is essential that the diffusion equation be solved in three-dimensional space, and furthermore, one must include in the equation an appropriate "source term" describing the statistics of the creation of composition fluctuations required to trigger the nucleation events. At present, a general treatment of this problem does not appear to be tractable. Nevertheless, by using a modelistic approach, i.e., by selecting a priori certain composition variation profiles and by evaluating their associated free energies, one can draw certain conclusions relative to the stability of these isolated composition modulations or "zones", both outside and inside the spinodal. This is the object of the present investigation. THEORY We shall now describe a standard four-parameter spherically symmetric GP zone, evaluate its free energy according to the Cahn and Hilliard11 formula (Eq. [5], below), and plot the corresponding energy surfaces as a function of zone amplitude and radius as

Jan 1, 1970

By P. H. Lindon, J. C. Billington

Fe-O melts containing 0.045 pct 0 were deoxidized with Mn-Si-A1 alloys. Product compositions were reluted to the melt and alloy compositions and were found to be most sensitive to the aluminum content of the alloy. Low residual oxygen contents could be obtained when aluminum oxide was present in the Products because of the reduction of silica and manganese oxide activities. Flotation of the Products from a quiescent melt was followed both by analysis of the oxygen content and metallographic measurement of inclusion concentration. MnO-SiO2-A12O3 products were found to float most rapidly when their composition was such that their viscosity may be expected to be low. Changes in the particle size distribution indicates that particle coalescence occurred and differences in the degree of coalescence are thought to be responsible for the different flotation rates observed between products 0f differing composition. Measured flotation rates were slower than those Predicted from a model based on Stoke's Law, although alumina flotation might be reasonably accounted for by this model. Interfacial effects between oxide particles and the melt are believed to be responsible for the discrepancy. It has been recognized that deoxidation products constitute a large proportion of the nonmetallic inclusions present in killed steel. The amount of oxide inclusions which originate as deoxidation products depends largely upon three factors. These may be summarized, according to P16ckinger1 as: 1) Amount of primary products remaining in the steel prior to cooling. 2) Residual dissolved oxygen content of the steel after deoxidation. 3) Amount of secondary products, formed during cooling and solidification, which remain entrapped in the solid steel. In a well-deoxidized steel containing residual aluminum and/or silicon, the equilibrium dissolved oxygen content is usually very low and so the maximum amount of oxide which may be produced as secondary deoxidation products is small in comparison with the amount of primary products. It may be seen, therefore, that the amount of indigenous nonmetallic inclusions may be minimized if a low dissolved oxygen content is achieved by deoxidation and if the primary deoxidation products are efficiently removed. Oxides which originate by reaction of the metal stream with the atmosphere during teeming are not considered in the present study. It is known that two or more deoxidizers may result in a lower equilibrium oxygen content when used in conjunction with one another than when any of the individual deoxidizers are used alone. Equilibrium studies by Hilty and crafts2 and by Bell3 have shown that manganese increases the effectiveness of silicon as a deoxidizer, and Walsh and Ramachandran4 relate this to a reduction in the activity of silica in the products as the manganese :silicon ratio in the steel increases. It was also shown by Herty's work on deoxidation of steel by silico manganese alloys,5 that there existed an optimum ratio of manganese to silicon which gave a minimum inclusion content. This ratio was in the range 4:l to 7:l and the (FeO-MnO-SiO2) products formed by such deoxidation practice were found to lie in a composition range having very low liquidus temperatures (1170 to 1250°C approx). The optimum manganese:silicon ratio was then explained by postulating that these fluid products were able to coalesce and that the larger particles formed floated out of the steel very quickly as predicted by Stoke's Law. The present work examines the effectiveness of various Mn-Si-A1 alloys as deoxidizers and their effects on the composition and removal of primary deoxidation products from a quiescent melt. EXPERIMENTAL TECHNIQUE Approximately 250 g of prepared Fe-O alloy, containing 0.045 to 0.055 pct O, were melted in an alumina crucible and deoxidized at 1550°C by plunging a thin steel cartridge containing the deoxidizer below the melt surface. A high frequency induction furnace supplying current at 8.5 kHz was used to heat a graphite susceptor, the interior of which had been machined to give a wall thickness of 0.85 in. to form a receptacle for the alumina crucible. The iron melt was essentially quiescent as the induced current was concentrated at the external surface of the graphite susceptor by the skin effect. A nonoxidizing atmosphere was maintained over the melt by passing a continuous stream of argon through the lid of the susceptor. The melt temperature was measured before deoxidation, and again at the end of an experiment by means

Jan 1, 1970

By S. N. Singh, M. C. Flemings

Results are presented of a detailed study on the combined influences of ingot dendrite am spacing and thermomechanical treatments on the structure and solution kinetics of high --purity cast and worked 7075 alloy. Solution kinetics were found to depend sensitively on ingot dendrite am spacing and on details of therrnomechanical processing, including amount of reduction and extent of&apos; solution treatment before reduction. An approximate analysis is given for rate of solution of nonequilibrium second phase in the cast and worked structres; results of the analysis are compared with experiment. MICROSEGREGATION in high strength aluminum alloys manifests itself as "coring" (composition differences within the primary aluminum-rich phase), and as interdendritic second phase. The mechanism of formation of the microsegregation is understood, and approximate prediction of the amount of second phase is possible for simple binary systems.1,2 When alloy elements or impurities are present in amounts less than their solid solubility at solution temperature, any phases forming from these elements are termed "nonequilibrium" and can be dissolved by appropriate solution treatment. The rate at which the nonequilibrium phases are removed depends sensitively on their spacing (dendrite arm spacing in the cast material, or band spacing in wrought material). When alloy elements or impurities are present in amounts in excess of their solubility at the solution temperature, second phase particles form an "equilibrium" second phase that does not dissolve in heat treatment and may, in fact, coarsen in such treatment. Usual commercial, high strength, wrought aluminum alloys contain nonequilibrium second phases that were not fully dissolved during ingot processing. They also contain equilibrium second phases resulting from impurities present in amounts greater than their solubility. As has been shown by Antes, Lipson, and Rosenthal,3 and will be demonstrated further in a subsequent paper by the authors,4 significant improvements in mechanical properties of high strength alloys can be achieved by reduction or elimination of these second phases. Methods of elimination are 1) to employ high purity materials to minimize amounts of equilibrium second phase, and 2) to employ suitable thermomechanical processing techniques to fully eliminate nonequilibrium second phases. Work reported herein comprises a study of selected thermomechani- cal processing treatments, and of their influence on solution kinetics of wrought high purity 7075 alloy. EXPERIMENTAL PROCEDURE Melting and Casting. The bulk of the work reported was performed on a single ingot of high purity 7075 alloy. The ingot was 4 in. by 4 in. by 8 in. high, uni-directionally solidified following a procedure previously described.5 The mold was heated to 1350°F before pouring the melt. The bottom chill was carbon coated stainless steel. Water was circulated through the chill after the melt was poured. The 7075 alloy was prepared from high purity virgin material (aluminum, zinc, magnesium) and from master alloys (Al-50 pct Cu, A1-15 pct Cr, A1-5 pct Ti). Final measured melt composition (wt pct) was: Zn Mg Cu Cr Ti Fe Si Al 5.70 2.28 1.35 0.18 0.15 <0.002 <0.012 bal Melting was done in a silicon carbide crucible; all tools were coated with zircon wash to minimize iron contamination; degassing was by bubbling chlorine through the melt. che-rmomechanical Treatments. Detailed studies were made on material taken from a location approximately 13 in. from the chill and 51/2 in. from the chill (i.e., from 1 in. thick slices taken between 1 and 2 in. from the chill and between 5 and 6 in. from the chill). Solution treatment was done at 860°F in an air-circulating furnace with a "bottom drop" arrangement to achieve minimum delay time between solution treatment and quench. Samples solution treated in this way were 2 in. by 2 in. by 1 in. Temperature of the quench water was approximately 10°C. Mechanical reduction was by cold rolling. Samples 11/2 and 51/2 in. from the chill were treated for 12 and 24 hr, respectively, before cold rolling. Reduction by cold rolling was then 4/1, 16/1, and 35/1. In each case, several intermediate anneals (1/2 hr at 860°F) were used to permit reaching the final thickness without cracking; two such anneals were used for the 4/1 reduction, five for 16/1, and six for 35/1. After working, materials were again solution treated for various lengths of time from 0 to 48 hr and quenched in water. Structural Measurements. Secondary dendrite arm spacings were measured using procedures previously described.&apos; For each measurement reported, five photomicrographs were first made at X75. Measurements were made of dendrite arm spacings in at least 20 different grains (grain structure was equiaxed). Grain size measurements were made by running a number of random traverses across photomicrographs of the samples and obtaining the mean lineal intercept. Measurement of the volume percent of second phase and porosity was done by quantitative metallography. A two-dimensional systematic point count was used

Jan 1, 1970

Jan 1, 1970

By A. Vignes, J. P. Sabatier

A study of ternary diffusion at 1315°C in the Fe-Co-Ni system covering the entire ternary phase diagram has been carried out. Experimental values of the four interdiffusion coefficients DFeNico,DFeCoCo, DFeCoNi were computed at 33 points of the ternary diagram. The limiting values of these coefficients on the three sides of the ternary diizgram have also been determined. This set of data is presented as a contour map on the ternary isotherm for each of the four coefficients. Assuming that the Fe-Co-Ni system can be considered as ideal, the phenomenological coefficients Lk ij (with k = Fe) have been computed and the Onsager reciprocal relutions verified. It is shown that, in the Fe-Co-Ni system, the diffusion process can also be described by only two interdiffusion coefficients DFeFe and DCoNiFe, the second cross diffusion coefficient DCoFeNi being equal to zero and the second direct difusion coefficient being given by The values of the cross interdiffusion coefficients and of the Onsager coefficients are explained in terms of atomic mobilities of the three components with the help of a kinetic calculation due to Lane and Kirkaldy.8 It is shown that all the results are explained if the tracer difusion coefficients of the three components are such that DFe > DNi ~ DCo in complete agreement with experimental data on the three binary systems. DIFFUSION in a ternary alloy system is basically and significantly different from diffusion in binary alloys. In binary diffusion, there is just one independent composition variable and only one coefficient Dij is needed to describe the interdiffusion process. In ternary diffusion, we must allow for cross interaction between two independent species, and four interdiffusion coefficients are necessary to describe the transport process. Kirkaldy et al.1,2 have determined interdiffusion coefficients in the systems, Cu-Mn-A1 and Zn-Cu-Al. Guy and Leroy4 have reported intrinsic diffusion coefficients for the r phase at 1300°C of Co-Ni-Cr alloys as a function of composition using the analysis of Philibert and Guy.3 Interdiffusion coefficients were measured at three different compositions at 850°C in the single phase copper-rich corner of the Cu-Zn-Mn system by Dayananda and race.5 Zie-bold and Ogilvie6 have determined interdiffusion coefficients in the single-phase field of the Cu-Ag-Au system at 725°C. Preliminary results on diffusion at 1136°C in the system Fe-Co-Ni have been reported.7 In this study, data are reported on a more complete determination of interdiffusion coefficients in the Fe-Co-Ni system in the entire ternary at 1315°C. There is complete solid solubility with the fcc structures at 1315ºC, and this system can be considered as nearly ideal. This approximation allows the computation of the Onsager phenomenological coefficients. These values can then be used to test a kinetic calculation of the Onsager L coefficients in ternary substitutional solid solutions due to Lane and Kirkaldy.8

Jan 1, 1970

By R. J. Fruehan, F. D. Richardson

Electrochemical measurements have been made of the activity of oxygen in copper and its alloys with silver and tin at 1100" and 1200°C. The galvanic cell used was Pt, Ni + NiO/solid ellectrolyte/[O] in metal, cermet, Pt The results do not support any of the equations so far designed for predicting the activities of dilute solutes in ternary solutions from activities in the corresponding binaries. If, however, a quasichemical equation is used with the coordination number set to unity, agreement between observed and calculated activities shows that this empirical relationship can be useful over a fair range of conditions. SEVERAL solution models have been proposed for predicting the activity coefficients of dilute solutes in ternary alloys from a knowledge of the three binary systems involved. Alcock and Richardson1 have shown that a regular model, and a quasichemical model,&apos; in which the dissolved oxygen is coordinated with eight or so metal atoms, can reasonably predict the behavior of both metal and nonmetal solutes when the heats of solution of the solute in the separate solvent metals are similar. But when this is not so, neither model gives useful predictions unless coordination numbers of one or two are assumed. Wada and Saito3 subsequently adopted a similar model to derive the interaction energies for two dilute solutes in a solvent metal. Belton and Tankins4 Rave proposed both regular and quasichemical type models in which the oxygen is bound into molecular species, such as NiO and CuO in mixtures of Cu + Ni + 0. However, their models have only been tested on systems in which the excess free energies of solution of the solute in the two separate metals differ by a few kilocalories. Ope of the reasons why more advanced models have not been proposed is because few complete sets of data exist for ternary systems in which the solute behaves very differently in the two separate metals. For this reason measurements have been made of the activities of oxygen dissolved in Cu + Ag and Cu + Sn. Measurements on both systems were made by means of the electrochemical cell, Pt, Ni + NiO/solid electrolyte/O(in alloy), cermet,Pt [1] The activity of oxygen was calculated from the electromotive force and the standard free energy of formation of NiO, which is accurately known.5 Before investigating the alloys, studies were made of oxygen in copper to test the reliability of the cell and to check the thermodynamics of the system. Of the previous studies those by Sano and Sakao,6 Tom-linson and Young,7 and Tankins et al.8,7 have been made with gas-metal equilibrium techniques; those by Diaz and Richardson,9 Osterwald,10 wilder," Plusch-kell and Engell,12 Rickert and wagner,13 and Fischer and Ackermann14 have been made by electrochemical methods. EXPERIMENTAL The apparatus employed was the same as described previously,9 apart from slight modification. The molten sample of approximately 40 g was held in a high grade alumina crucible 1.2 in. OD and 1.6 in. long. The solid electrolytes were ZrO2 + 7½ wt pct CaO and ZrO2 + 15 wt pct CaO; the tubes 4 in. OD and 6 in. long were supplied by the Zirconia Corp. of America. They were closed (flat) at one end. In one experiment with Cu + O, both electrolytes were used in the cell at the same time. The reference electrodes inside the electrolyte tubes consisted of a mixture of Ni + NiO. They were made by mixing the powdered materials and pressing them manually into the ends of the tubes, with a platinum lead embedded in them. The tubes were then sintered overnight in the electromotive force apparatus at 1100°C. By sintering the powders inside the tubes (instead of using a presintered pellet9) better contacts were obtained between the electrolyte, the powder, and the platinum lead. Troubles arising from polarization9 were thus much reduced. The electromotive force was measured by a Vibron Electrometer with an input impedence of 1017 ohm; the temperature was measured with a Pt:13 pct Rh + Pt thermocouple protected by an alumina sheath. The couple was calibrated against the melting point of copper. The cermet conducting lead of Cr + 28 pct Al2O3, previously found to be satisfactory9 for use with Cu + 0 was also found satisfactory with Cu + Ag + 0 and Cu + Sn + 0. Superficial oxidation was observed, but it did not interfere with the working of the cell. The reaction tube containing the cell was closed at each end with cooled brass heads and suspended in a platinum resistance furnace. The tube was electrically shielded by a Kanthal A-1 ribbon which was wound round it, and the ribbon was protected by a N2 atmosphere between the furnace and the reaction tube. The cell was protected by a stream of high purity argon which was dried and passed through copper gauze at 450°C and titanium chips at 900°C. All the metals used were of spectrographic standard. Procedure. In studies of the system Cu + 0, be-

Jan 1, 1970

By R. E. Cech, T. D. Tiemann

It has been shown that hydrogen may be made to serve as a rapid and eflicient reducing agent for Cu, Ni, Co, and Fe sulfides if a scavenging agent for hydrogen sulfide is intimately mixed with the sulfide particles being reduced. Accelerated reduction kinetics are demonstrated for nickel sulfide. Copper, nickel, and cobalt sulfides, when treated at certain temperatures in a combined reducing agent-scavenging agent system, are converted to voluminous masses of fibrous metal product. Studies have been carried out to determine the conditions which lead, on the one hand, to irregular poly crystalline fibers and, on the other, to long single crystal filaments a few microns in diameter. A mechanism is proposed to account for the formation of single crystal filuments. The sulfide minerals of Cu, Ni, Co, and Fe are an important source of these metals yet there has been comparatively little scientific effort devoted towards understanding reduction mechanisms of these minerals. This may be, in part, due to the fact that the most convenient reducing agents for carrying out such studies, viz., hydrogen and carbon, do not react appreciably with sulfides. We have found that the reaction of hydrogen with metal sulfides can be markedly accelerated by placing a scavenging agent for hydrogen sulfide in close proximity to the metal sulfide. A brief series of experiments demonstrating relative reduction rates is reported in this paper to illustrate the effect. With the reduction process thus accelerated we have observed an unusual type of reduction behavior on some of the sulfides investigated. Under certain conditions the metallic product of the reduction reaction takes the form of filaments growing outward from the sulfide particles. The present paper deals largely with efforts to classify the various types of growth forms observed. This study has shown that filamentary growths from sulfides take a much greater variety of forms than has heretofore been reported by Ercker,1 Hardy,2 and Nabarro and Jackson3 in their reviews of metallic growths from copper and silver sulfides. THERMODYNAMIC CONSIDERATIONS The thermodynamics for hydrogen reduction of metal sulfides is quite unfavorable. For the sulfides considered here equilibrium constants typically range from 10-3 to 10-5. These low equilibrium constants impose severe kinetic limitations on reduction since hydrogen sulfide must be transported out of the system at concentrations of only a few hundred ppm. Unless extremely high gas flow rates are employed the atmosphere surrounding any sulfide particle will always be essentially in equilibrium with the sulfide. If, however, one places an efficient scavenging agent for hydrogen sulfide in close proximity to the metal sulfide particles the concentration of H2S near the metal sulfide will be held to a very low value. This would permit the reduction reaction to proceed with little or no inhibition from a buildup of reaction product gas. It is well known that calcium oxide is capable of removing hydrogen sulfide from a hydrogen gas stream of low dew point.4 If a sufficient quantity of calcium oxide is mixed with the metal sulfide particles the reaction: CaO+H2S=CaS+ H2O [l] will substitute moisture in place of hydrogen sulfide in the gas stream and this will not affect, in a direct manner, the reaction: MeS +H2=Me + H2S [2] A convenient method of considering the thermodynamics of the combined reducing agent-scavenging agent system is to consider the atmosphere when the partial pressure of hydrogen sulfide is the same over both the metal sulfide and the scavenging agent, i.e., pH2S (1) =pH2S (2). As a consequence: pH2O (1) pH2(2) =K1K2 The chemical driving force for reduction will depend inversely upon the moisture content of the gas and will be 0 when, in the system, pH2O = pH2.K1K2. Table I lists values of the equilibrium constants for reduction and H2S scavenging reactions for a number of sulfides at several temperatures. Data are taken from Rosenqvist4,5 and Kelly.6 The equilibrium constant products calculated from this data show that the limiting level of gaseous reaction product has been increased by a factor of 10&apos; to l04 as a result of substituting a reducing agent-scavenging agent system for a simple reducing agent system. One possible side effect which must be considered is the possibility that the moisture evolved in the scavenging reaction might cause the atmosphere in the system to be sufficiently oxidizing to favor the formation of oxide rather than metal. This possibility was examined by comparing the equilibrium constant products listed in Table I with equilibrium constants for hydrogen reduction of the respective metal oxides. It was found that for copper, nickel, and cobalt the combined reduction-scavenging reactions could not develop a sufficiently high oxidizing potential in the

Jan 1, 1970

By J. E. Gruzleski, O. Knacke, M. J. Meixner

The evolution of CO from small specimens of solidifying steel was investigated using a thermobalance technique. Gas evolution begins in the temperature range between 1400" to .1450°C. This behavior is independent of the initial carbon content and the cooling rate. At a constant temperature above 1220°C, the rate of gas evolution is temperature independent. Out-gassing of the specimen is thought to proceed through small channels formed by the coalescence of bubbles between dendrites. LIQUID metals have a higher solubility for gas than the solid. During freezing, these gases form bubbles which may either escape or remain trapped in the solid. Recent work1 on the evolution of oxygen from solidifying silver showed that the gas is not liberated until a large fraction of the melt has solidified, and that the gas is evolved in an almost explosive manner. Because of the commercial importance of CO evolution in solidifying steel, a similar investigation was undertaken to study the mechanism of gas liberation in steel. EXPERIMENTAL The thermogravimetric method of Knacke and eixner1 to follow gas evolution by weight changes was employed. Weight differences of 0.1 mg could be reproducibly measured on specimens weighing 2.5 g. These specimens were taken from alloys with varying carbon contents from <0.01 to 0.11 wt pct and <0.01 wt pct of each of silicon, manganese, and aluminum. Specimens contained in 8 mm ID magnesium alumi-nate crucibles were heated to 1550°C under 1 atm of very pure (99.995 pct) CO with a tungsten resistance furnace. The specimen temperature was measured with a Pt-Pt/Rh thermocouple located directly beside the crucible. It is felt that the temperature thus measured was within 5" to 10°C of the specimen temperature since on fusion the weight of the specimen increased due to uptake of CO, and this weight increase was always observed to begin within the temperature range of 1530" to 1540°C. Even after 3 hr there was a small weight increase with time indicating that the molten specimen was not completely saturated with CO. Therefore, complete saturation was not attempted, and all specimens were simply held for 3 hr at 1550°C before cooling. Two types of cooling experiments were performed. One consisted of continuous cooling at linear rates from 5°C to 20°C per min down to 1020°C. The second type of experiment involved continuous cooling at 10°C per min to a predetermined temperature, and then holding the specimen at this temperature for a certain time. During all of these experiments, the sample weight was recorded as a function of time and temperature. OBSERVATIONS Continuous Cooling Experiments. No gas evolution was detected above 1450°C. The initial weight decrease always occurred between 1400" and 1450°C with the average at 1430°C. With further cooling, gas evolution continued to about 1150°C. Between this temperature and 1020°C no further outgassing was observed. This behavior did not seem to depend on the initial carbon content nor on the cooling rate in the range investigated. Typical patterns of gas evolution are shown in Figs. 1 and 2 for different initial carbon levels and cooling rates. The total amount of gas liberated during solidification varied with the rate of cooling. The total weight change on continuous cooling at a rate of 10°C per min

Jan 1, 1970

By Robert A. Rapp, Ronald L. Pastorek

Solid-state electrochemical measurements by three alternative experimental procedures were made with the cell FeO, Fe3O4 |Zro.85Cao.15O1.85 |Cu| Zr0.85CaO.15O1.85 | FeO, Fe304 to establish the solubility and diffusivity of oxygen in solid copper in the temperature range 800" to 1030°C. The solubility of oxygen in solid copper and the diflusivity of oxygen in solid copper Dgu = 1.7 X 10-2 exp(-16,000/RT) Cm2/sec were determined and confirmed in alternative experiments. The enthalpy of solution of oxygen in solid copper equals —10,000 cal per mole; the partial excess entropy of the oxygen atoms in the Cu-O dilute solution is approximately the same as that found for interstitial atoms in other metals. The diffusivity of oxygen in solid copper is consistent with that expected for an interstitial atom. RELIABLE values for the saturation solubility N(s) and diffusivity DO of oxygen in solid copper have not been unambiguously established in the literature. Following three early determinations by others,1"3 Rhines and Mathewson4 reported that the solubility of oxygen in solid copper increased from 0.007 at. pct 0 at 600°C to about 0.015 pct at 1050°C. Phillips and skinner,, using essentiially the same analytical procedure, reported that the solid solubility increases from 0.0018 at. pct 0 at 550°C to about 0.0075 pct at 1050OC. The only previous value for the diffusivity of oxygen in solid copper was reported by Ransley.6 Ransley deoxidized Cu-Cu2O alloys in an atmosphere of carbon monoxide gas to yield a solubility-diffusivity product. He used the solubility data of Rhines and Mathewson to calculate the diffusivity values. Another method for obtaining the solubility-diffusivity product (N(s) DO) is by measuring the widths of internal-oxidation zones in copper alloys as reported by Verfurth and Rapp.7 However, the calculated N(S)Do products depend upon the alloy content of the specimen, so that the internal oxidation of copper alloys does not follow ideal internal oxidation kinetics. As a result, unequivocal values for the N(s) DO product were not obtained by this procedure. A solid-state coulometric titration technique similar to that employed in this work was introduced by C. Wagner8 to study the dependence on silver activity of the Ag/S ratio in silver sulfide in the temperature range of 160" to 300°C. Similar experiments have been carried out by C. Wagner and co-workers9-11 to study the stoichiometry range of silver and copper tellurides, cuprous sulfide, and cuprous selenide. Numerous authors have carried out electrochemical measurements with a solid oxygen-ion-conducting electrolyte to determine the solubility and/or diffusivity of dissolved oxygen in several liquid metals.12-l6 Rickert and Steiner17,18 have used solid-state electrochemical measurements to determine the diffusivity of oxygen in solid silver from 760" to 900°C. Two different cell geometries were used. In the cell of linear geometry Fe, FeO | ZrO2 + (CaO) | Ag + [0 (dissolved)] [1] oxygen diffused from the interior of the silver electrode to the silver/electrolyte interface where the oxygen activity had been lowered from a fixed initial value to practically zero by the application of voltage to the cell. The diffusivity of oxygen in solid silver was determined from the solution of the diffusion equation and the time dependence of the cell current. However, this determination of the diffusion coefficient depended upon a knowledge of the solubility of oxygen in solid silver. A cylindrical geometry was used for the cell Pt, O2(Po2 = 0.21 atm) | ZrO2 + (CaO) | Ag + [0 (dissolved)] [II] which also allowed the diffusivity of oxygen in solid silver to be determined. These values were in agreement with other available data.l9 Recently, Raleigh20,21 used a method involving the measurement of diffusion-limited currents in a cell involving the AgBr solid electrolyte to determine the diffusion coefficient of silver in Ag-Au alloys at 400°C. Diffusivity values on the order of l0-14 sq cm per sec were measured in the alloy composition range 10 to 60 at. pct Ag in a single experiment. From numerous electrical conductivity and galvanic cell measurements,9&apos;22"26 the solid solution Zr0.85 Ca0.15 O1.85 has been established as an electrolyte with predominant oxygen ion conduction over a wide range of intermediate and high oxygen activities. For interrelating the thermodynamics and the kinetics of the dissolution of oxygen in solid copper in this investigation, a galvanic cell was constructed with FeO-Fe3O4 as the reversible reference electrode, the Zr0.85Ca0.15 O1.85 electrolyte, and a pure copper specimen under-saturated in oxygen as the other electrode. THEORETICAL ANALYSIS Three variations of a high-temperature electrochemical technique were used in this study to provide two determinations each of the solubility and diffusivity

Jan 1, 1970

By G. L. F. Powell

g samples of Cu-20 wt pct Ag alloy have been mdercooled to a maximum of 197°C by melting under a slag of commercial soda-lime glass in a vitreous silica crucible. No grain refinement of the primary copper was observed in samples undercooled to the maximum of 197°C. When the samples contained a small amount of oxygen, the copper dendrites were partially recrystallized at undercoolings greater than 97°C. In previous papers&apos;-3 reporting the grain structure of undercooled silver and copper, it was observed that grain refinement was dependent on both undercooling and oxygen content. Grain refinement occurred in undercooled silver when the degree of undercooling exceeded the range 153" to 175"C, while in Ag-0 alloys (0.12 wt pct) fine equiaxed grains were exhibited when undercooling was greater than 50°C. Similarly, copper samples undercooled as much as 208°C displayed fan-shaped growth from a single nucleation site, while the grain structure of Cu-O alloys (0.08 wt pct) was fine and equiaxed at undercoolings larger than 150°C. Thus the presence of oxygen greatly reduced the undercooling at which grain refinement occurred. It was also observed that the change in grain size resulted from recrystallization and was not due to an enhanced nucleation rate in the liquid-solid transformation. It is possible that the influence of oxygen on recrystallization is due primarily to its presence as a solute element. walker4,&apos; reported that, although a grain size change did not occur in pure nickel until the undercooling exceeded 150°C, small grains were observed in samples of Ni-Cu alloy solidified at small and large degrees of undercooling. Jackson et al.6 suggested that the fine grained structure of the Ni-Cu alloy resulted from the melting off of dendrite arms during recalescence. This remelting process may occur in alloys as a result of segregation during freezing which causes a variation in liquidus temperature from point to point within a dendrite. It was therefore decided to undercool copper with a metallic alloying element to ascertain whether the presence of a metallic solute would have a similar effect to oxygen in inducing grain refinement. A Cu-Ag alloy was chosen, since both metals had been shown to behave similarly on undercooling. The alloy Cu-20 wt pct Ag was selected since the eutectic constituent outlines the initial growth form of the primary copper, so that the as-frozen grain structure is not obscured if subsequent recrystallization occurs. This paper describes the results of undercooling experiments carried out with Cu-20 pct Ag samples undercooled to a maximum of 197°C and the effect of oxygen content on the grain structure of the undercooled samples. EXPERIMENTAL Melting was carried out in a small cylindrical resistance furnace using "fine" silver granulate and oxygen-free high conductivity copper. The procedure adopted was to melt the required quantity of silver in air in a clean vitreous silica crucible for approximately 15 min, freeze, and add granulated commercial soda-lime glass to form a complete surface slag cover, after which the sample was melted and frozen several times to reduce the oxygen content. The glass slag cover was approximately 3 in. thick. Pieces of copper (=50 g) were added to the crucible until the required quantity to make 350-g samples of alloy had been charged. Each piece was added quickly to the crucible which was held at a temperature slightly above the melting point of silver. The piece was quickly pushed beneath the glass to minimize oxidation and any oxide coating usually decomposed before the piece had settled down into the silver. After the full quantity of copper had been added, the melt was stirred with a silica rod to hasten homogenization and a Pt/Pt 13 pct Rh thermocouple enclosed in a vitreous silica sheath inserted for temperature measurement. Heating and cooling curves were recorded on a potentiometric chart recorder fitted with a zero suppression unit. The milli-voltage range of the recorder was adjusted so that temperatures could be read to 1°C. Heating and cooling curves were taken every hour until three consecutive readings gave the same solidus-liquidus range, consistent with the solidus-liquidus range for this alloy composition by reference to Hansen and Anderko.7 Metallographic examination of samples frozen at this stage, failed to show any variation in composition from bottom to top of the ingot. Consequently, it was considered that the melt was homogeneous at this stage and undercooling experiments were then car-

Jan 1, 1970

By P. Shahinian, R. L. Stegman, M. R. Achter

Creep strength has been determined as a function of oxygen pressure for nickel at 510" and 600°C. Creep-rupture life drops from its highest value at the ultimate vacuum of 10-7 torr to a high pressure plateau of nearly constant life which, at 600°C, commences at torr. This behavior is described in terms of a model that relates surface gas adsorption and the generation of fresh crack surfaces. In a calculation based a this model the predicted pressure beyond which the creep strength should not decrease further has been found to agree well with the observed value. The introduction of oxygen during creep in vacuum was found to cause a prompt rise of the creep rate. The magnitude of the environmental effect was much less than generally observed in fatigue and this is attributed to the diffarent crack morphologies in creep and fatigue. GASEOUS environments are known to affect creep resistance adversely1-5 but the effects have received less attention than the much larger ones in fatigue. At low temperatures it has been shown for nickel that the creep rates are faster and the rupture lives are shorter in air than in vacuum and that these environmental effects can be correlated with faster crack growth rates in the reactive environment.&apos; In support of a mechanism of gas adsorption at the crack tip to explain the above effects, a calculation based on kinetic gas theory showed that surface coverage by gas was not as rapid as the exposure of fresh metal at the crack tip in the ultimate vacuum but was in air.7 In fatigue a number of investigations have demonstrated that the reduction of life with gas pressure usually follows a step curve with most of the reduction occurring in a transition range at intermediate pressures and little or no change at high and low pressure.8 There have been no studies, however, in which creep strength has been determined as a function of reactive gas pressure to see whether similar effects would be found. This investigation was undertaken to study the kinetics of the environmental effect in the creep of nickel. EXPERIMENTAL PROCEDURES The environmental creep apparatus was constructed to minimize many of the experimental conditions which can cause scatter in data. Dead weight loading is employed and weights are contained within the environmental chamber to eliminate bellows and the necessity for compensation for gas pressure changes. Elongation over the specimen gage length is recorded from a Riehle Model DV Controlled Environment Extensometer which provides readout to 2.5 x 10- in. Temperature is monitored and controlled from a thermocouple at- tached to the specimen gage length to +0.5oC or better. The freon trapped oil diffusion pump vacuum system included an auxiliary liquid nitrogen trap and, with a wide range variable conductance leak valve, pressures from the ultimate vacuum of 10"7 to l03 torr of any desired gas are obtainable. Pressure is measured by means of a hot filament ionization gage below 10"3 torr and a NRC Model 530 Alphatron Vacuum Gauge above 10-3 torr. The oxygen, which was leaked into the environmental chamber, was an Aviator&apos;s Grade of 99.5 pct purity with a dew point of -80°F. Analysis of the nominal 99.8 pct Ni showed the following percentages of major impurities: 0.02 C, 4.01 Mn, 0.005 Si, 0.01 Cu, 0.02 Fe, and 0.034 Mg. Cylindrical specimens were machined with a 0.180

Jan 1, 1970

By M. J. Mullikin, J. B. Conway

Several different equation forms were studied to determine the extent to which each particular equation type yielded an accurate representation of a given set of first-stage creep data. Specially developed least-squares computer programs were employed to obtain unique solutions for the creep-equation constants. In studies of creep data for lead at 14.5"C, nickel at 70OCC, and arc-cast tungsten at 2400°C consistently good representations were afforded by a third-degvee polynomial in t113. Other effective correlations were obtained with the Andrade, Cottrell-Aytekin, and de Lacombe equations along with a third-degree polynomial in t112. A newly proposed creep equation involving a hyperbolic-sine term in t1&apos;3 was also found to yield very satisfactory results. This appears to be the first use of such a functional relationship in the analysis of transient creep behavior. In general, the deformation behavior of a material subjected to a tensile load is described in a strain-time plot. Experiments have shown that, following the instantaneous elongation which results on loading, the instantaneous slope (i.e., creep rate) of the strain-time curve decreases gradually with time and approaches what appears to be a constant value. Equations, therefore, which are proposed to provide an analytical representation of such first-stage or transient creep behavior must be consistent with these observations. In other words, at time zero, these equations must yield a finite deformation value corresponding to the elongation on loading. Furthermore, the first derivative of these equations must exhibit a fairly high although finite value at time zero, and then values which decrease with time. It is also preferable, in general, that the derivative eventually assume a constant value at high time values. Although many studies of first-stage creep behavior have been made and numerous equation forms have been proposed for expressing this type of deformation behavior, it is still not possible to state with any assurance which functional relationship will be adequate for a given set of measurements. Parabolic relationships have been employed most extensively although logarithmic expressions have been shown to be applicable in some instances. Combined forms1 have been studied in considerable detail and have on occasion been found to be very effective. One such expression is that of iXndradea and involves a modified parabolic expression: where 1 is the instantaneous length, lo is the length immediately following load application, t is the time, and p and k are constants. Another combined form has been proposed by weaver3 as follows: where E is the instantaneous strain, is is the linear creep rate, t is the time, and a and b are constants.

Jan 1, 1967

By Raymond L. Orr

Calorimetric measurements have been made of the heats of solution of gold axd indium in a number of liquid tin-rich alloys at a temperature of 705°K. Relative partial molar enthalpies of gold were determined as a function of dilute solute concentrations in the binary Au-Sn system and in the ternary Ag-Au-Sn and Au-In-Sn systems. Relative partial molar enthalpies of indium were determined for dilute solute concentrations in the In-Sn, Ag-In-Sn, and Au-In-Sn systems. Through application of the interaction coefficient concept introduced by Wagner and more recently extended by Lupis and Elliott, ualues have been obtained for the Au-Au, Ag-Au, Au-In, and Ag-In enthalpy interaction coefficients in liquid tin. The results are inte@reted in terms of solute atom distributions through comparisons with the predictions of dilute solution models. INTEREST in the thermodynamic behavior of dilute liquid alloys stems from two primary sources. From a practical viewpoint, it is often of importance to know or to be able to predict the effect that one solute will have on the thermodynamic properties of the other solutes in a multicomponent system. Studies of dilute solutions may also be rewarding from the theoretical point of view. Some of the difficulties arising from the more complex interactions possible in concentrated solutions are avoided, leading to easier interpretations in terms of solutibn models and bonding energies. For example, in a binary alloy the limiting values of the partial molar properties of the solute represent the case for which each solute atom is completely surrounded by atoms of the solvent, and no other interactions are possible. Solute - solute interactions in dilute solutions are conveniently treated by the interaction coefficient concept of ~agner.&apos; Using a Taylor series expansion for the logarithm of the activity coefficient, In yi, of a component, i, in a solution consisting of dilute solutes with atomic fractions xi, xj, xk, and so forth, in a solvent, s, and neglecting the second- and higher-order terms, Wagner obtained the expression: Where Wheyz is the limiting value of yi in the pure solvent, s. The interaction coefficients, e:, e:, and so forth, are defined by: Since yi is related to the excess partial molar Gibbs energy of component i by GfS = RT In yi, E: is re- ferred to as the "Gibbs energy interaction coefficient". Lupis and Elliott~ extended Wagner&apos;s treatment to the relative partial molar enthalpy, aHi, and the excess partial molar entropy, qS. They defined the enthalpy interaction coefficient as: and the excess entropy interaction coefficient as: leading to expressions similar to Eq. [I]: The three interaction coefficients are related by? The self-interaction coefficients, E:,qf,and at, represent the effects of interactions between atoms of component i in solvent s; €3, q{, and 03 represent the effects of interactions between i and j atoms in solvent s; and so forth. From the Maxwell-type relationships between partial molar quantities, it may be shown1&apos;2 that e{ = ej, r}\ =tj), and ai - crj. Experimental determination of the interaction coefficients requires extremely precise measurements of the appropriate properties as functions of solute concentrations in very dilute regions. Few such data exist for intermetallic alloys, even for binary systems, because of experimental difficulties. This is especially true with respect to enthalpy data for dilute alloys. Examination of the compilation of data for binary alloys by Hultgren, Orr, Anderson, and ~elle~~ reveals that determination of limiting values of &fii(Xi =0) often requires extrapolations to infinite dilution from xi = 0.05 to 0.10. Direct measurements of the enthalpy interaction coefficient, qi, for multicomponent systems are virtually nonexistent. This quantity is of course subject to direct calorimetric measurement, but its determination must be limited to cases where high experimental precision is possible. The heats of solution of gold and indium in liquid tin can be measured with relatively high precision, which makes determinations of 73 involving these metals as the measured solutes experimentally attractive. This paper presents the results of such determinations of the Au-Au, Ag-Au, Au-In, and Ag-In enthalpy interaction coefficients in liquid tin. MODELS FOR DILUTE-SOLUTION BEHAVIOR Random-Solution Behavior. In the quasi-chemical treatment of solutions, only nearest-neighbor bonds are considered, and a composition independent value

Jan 1, 1967

By G. R. Love, M. L. Picklesirner

THE use of color adds a new beauty, power, and versatility to metallography. This has been amply demonstrated in a number of public exhibits and on the walls of corporate, government, and university metallographic laboratories around the world. Indeed, in a number of cases, the color not only is of esthetic value but is scientifically essential to the competent examination of particular samples and systems. Optical strain analysis,&apos; determination of true grain size in U-A1 alloys: phase identification in Nb-base alloys,3 analysis of transformation rates in heat-treated Nb-Zr alloys,4 measurement of solidification segregation in cast Nb-Zr alloys,&apos; study of the intermediate phases in Nb/Sn diffusion couples,6 investigation of the oxidation kinetics at various surface orientations in single crystals,7&apos;8 and a number of other research efforts9&apos;10 would have been vastly more difficult without the extra clarity and distinction introduced by color—either the colors resulting from particular polarized light techniques or the interference colors obtained by oxidation. It is unfortunate that the beauty of colored micro-structures cannot be shared more widely, but the economics of publishing a technical magazine apparently prohibits color plates except under extraordinary circumstances. It is doubly unfortunate that no convenient, general method exists for communicating the information available in colored microstructures to the scientific community. True, in some instances the color changes observed are dramatic enough that direct black and white photography is adequate. In photoelastic stress analysis, several orders of interference fringes may be produced and it is frequently desirable to suppress the color effects (use monochromatic light) to simplify interpretation. On the other hand, relatively minute color variations may be readily detectable by eye, be of considerable assistance in interpreting a structure, and yet be enormously difficult to reduce to an unambiguous black and white print. In order to illustrate a possible compromise or solution to the above problems, we wish to show three black and white photographs of the same region of an Nb/Sn diffusion couple whose interpretation depends heavily upon use of color techniques. These micro-structures were obtained from a section through a rod consisting of commercial-purity niobium tubing filled with 99.995 pct Sn (melt-filled in 10"6 Torr vacuum). The composite rod, measuring 500 mm long, was bent into a hairpin shape and inserted into water-cooled copper grips in a vacuum high-temperature furnace. The rod was outgassed for 2 hr at a furnace temperature of 800°C and was further heated by self-resistance using a 2.2-kva ac power supply. The rod was raised to 1910°C over a 4+-hr period, held for 15 min, slow-cooled to 1500°C over a -hr period, and free-cooled to room temperature (all temperatures were measured with W5Re/W26Re thermocouple). Throughout the heat treatment the ambient pressure was kept below 8 x lo-&apos; Torr. Metallographic preparation of the sample was completely conventional except that it was necessary to extend the grinding steps to remove damage in the brittle NbsSn phase introduced by sectioning. No satisfactory chemical or electrochemical etch has been found for these diffusion couples; of those which have been tried, any that attack the niobium-rich phases at all attack the tin-rich phases disastrously. However, the final polishing steps did produce a slight relief effect and this, coupled with anodizing, yielded excellent clarity in the microstructures for direct observation. The anodizing voltage was chosen to give maximum change in color per unit change in composition. In this case, indeed normally, this maximum sensitivity occurred in the blue-to-red color range. The pure niobium tube wall anodized to a relatively bright blue

Jan 1, 1967

By E. Lifshin, R. E. Hanneman

QUANTITATIVE applications of the electron micro-beam probe frequently involve the evaluation of complex mathematical expressions and/or the analysis of large amounts of experimental data. The purpose of this communication is to describe briefly a versatile and useful computer program system that is applicable to analyze rapidly a wide variety of practical microprobe problems. This system consists of a group of ten FORTRAN programs that can be stored on tape, cards, or in the memory disc of the computer. These programs, or links, can be run individually or in any prespecified sequence without interrupting the operation of the computer or without destroying information which is being transfered from one link to another. For the program system described here a GE-235 computer with disc storage was used, so that the DCHAIN method of program linking was employed. Included in the library are programs to: 1) initiate analysis of a new set of data and transfer control between all other programs in any predetermined manner; 2) generate theoretical calibration curves of composition vs relative intensity; 3) generate empirical deviation parameters from least-square fits of experimental calibration data from standards of known composition; 4) convert raw X-ray data to corrected composition; 5) determine inter diffusion coefficients by Matano analysis of con centrat ion -distance data on a uniaxial diffusion couple; 6) determine activation energies and frequency factors of temperature-activated processes such as diffusion; and 7) generate calibration curves for determination of the thickness of thin films using microanalysis. A detailed description of these computer programs and their underlying principles is available on request from the authors."&apos; The first program to generate theoretical calibration curves of corrected relative intensities vs composition uses the Poole and Thomas atomic number correction&apos; and the Philibert absorption factor&apos; with a voltage-dependent mass absorption coefficient for electrons in the alloy. A modified Castaing fluores- cence correction is also used which includes the effects of both Ka and KO radiation and over voltage.&apos; Once the theoretical curves have been calculated in 1 wt pct intervals, these results are least-squares fit to obtain Ziebold deviation parameters&apos; which are stored in COMMON in the computer memory. The net discrepancies between the original theoretical calibration curve and the regenerated curve using the Ziebold parameter are computed. Although this link is explicitly written for K-K fluorescent interactions, it can be applied to K-L, L-K, and L-L interactions as previously disc~ssed."~ Similar programs have also been written to utilize the Wittry fluorescence correction, Birks combined corrections, and various other corrections.&apos; These modified programs have proved to be quite useful for quantitative comparisons of the results of the various theories. The program for conversion of raw X-ray data to corrected composition includes corrections for drift, backround, and instrumental dead time. The corrected intensities are converted to composition by use of the Ziebold equation"2 and parameter obtained from the program system. The results can be obtained for either atom or weight fractions. In addition to accurately computing interdiffusion coefficients the Matano analysis program calculates least-square smoothed values of concentration, concentration gradient, and curvature for each point on the raw input concentration profile. In order to obtain high accuracies a unique method of performing a least-square polynomial fit to incrementally advancing profile segments which overlap is used.&apos; This program has been successfully modified for use in ternary diffusion problems3 and can readily be modified to handle analysis of diffusion profiles which include phase boundary discontinuities. This link is generally applicable to analysis of interdiffusion data obtained by other techniques as well as by the microprobe. The primary function of the next program is to least-squares fit experimental diffusion data to the normal Arrhenius function: D = Doexp(-Q/RT), to obtain values of Do and Q. In addition the probable error and one, two, and three u statistical confidence limits of DO, Q, and log D are evaluated. This program is also directly useful for analysis of any other simple temperature-activated processes including conductivity, and certain deformation and chemical processes. The program to generate calibration curves for film-thickness determination using the microprobe is based on a numerical integration of the equation derived by Cockett and ~avis.&apos; Values of film thickness obtained by this program for copper on various substrates are in good agreement with measurements made by other techniques. Versions of the above program system have been prepared for use with or without a remote teletype connection to the computer for processing on either a real-time or time-share basis. The instrumentation coupling a microprobe to a teletype for automatic data collection and analysis by the presently described program system has been reported elsewhere by the authors.&apos; If teletype equipment is not used to communicate with the computer, standard methods of card reading and tape reading can be used. In either

Jan 1, 1967

By E. R. Gilbert

Klopp et al. have reported data on tensile and creep properties of are-melted and electron-beam-melted tungsten. We would like to point out some similarities between their creep results and ours on are-melted and powder-metallurgy tungstenand offer some alternate interpretations of the data which lead to different rate-controlling mechanisms. Klopp et al. observed extensive transient creep in the 2250" to 3500°F temperature range where the acthetivation energy is 106,400 ± 3300 cal mole-1. They suggest cross slip as the possible rate-controlling mechanism since "cross slip may occur only by the movement of screw dislocations, and not edge dislocations, thus allowing only partial recovery of the strain hardening and resulting in a continuously decreasing creep rate, characteristic of transient creep.&apos;&apos; Our tests on are-melted tungsten and powder-metallurgy tungsten displayed an extensive transient creep stage in the same temperature range but they also displayed a steady-state creep stage. The stress levels and strain rates of the two investigations, 1.4 to 7 kg mm-2 and 10-7 to 10-4 sec-1 for that of Klopp et al. and 4.5 to 10 kg mm-2 and 10-8 to 10-4 sec-1 for those of Gilbert et al., are overlapping. The activation energy for creep of our powder-metallurgy tungsten is 105,000 cal mole-1. This close agreement with the value found by Klopp et al. suggests that the same mechanism may be common to both types of materials. It may be possible to gain some additional insight into the rate-controlling mechanism by applying the analysis made by conrad22 in terms of the activation volume. The value for the activation volume depends upon the dislocation mechanism and may be used to help differentiate between various dislocation mechanisms. The activation volume for cross slip is 10 to 100 b3. Values determined from our data were 130 to 220 b3 which are in better agreement with the range of 102 to 104 b3 for nonconservative motion of jogs. Values determined from data of Klopp et al. are 240 to 1040 b3. We suggest that the rate-controlling mechanism for creep in this temperature range may be nonconservative motion of jogs assisted by enhanced diffusion along grain boundaries or dislocations. This model requires an activation energy lower than that for self-diffusion as is observed. If this is the controlling mechanism, a higher activation energy which is closer to that for self-diffusion might be expected for tests conducted on specimens with reduced grain boundary area, i.e., larger grain size. In support of this premise we have conducted tests on speci- mens having one to three grains per cross section for which a higher activation energy of 130,000 cal mole-1 was observed. Also, in agreement with results of Klopp et al., a faster creep rate resulted from the increased grain size and probable decreased dislocation density. Further argument against the cross slip mechanism is based on the constant activation energy we obsserved over the range of 5.0 to 10 kg mm-2,21 whereas the activation energy for cross slip displays a significant stress dependence.23 In addition to this, the powe on the stress term for the strain-rate dependence appears too high for cross slip. The values for tungsten fall in the range of 4.5 to 7 (Klopp et al. and Ref. 21) whereas values of approximately 2 are observed in creep by cross slip for hep metals.24 At temperatures from 3300° to 4000°F, Klopp et al. observed an activation energy of 141,000 ± 4000 cal mole-&apos; which is less than that for self-diffusion, 153,000 cal mole-1, 25 and less than that reported for high-temperature creep of powder-metallurgy tungsten by Green:&apos; 160,000 cal mole-1. Klopp et al. suggested that recovery of strain hardening by dislocation climb is the rate-controlling mechanism in this temperature region. Actually an activation energy larger than that for self-diffusion is predicted by this model if the temperature dependence of the elastic modulus27 is taken into consideration. Klopp et al. attributed this lack of correlation to the compositional differences between powder-metallurgy tungsten used by Green and their own arc- and EB-melted tungsten. It is our position that the value of 141,000 ± 4000 cal mole-&apos; represents a transitional value between the higher and lower activation energies and does not represent a discrete mechanism of deformation. If the tests had been extended to higher temperatures the higher value of activation energy would have been observed as has been demonstrated by both our analy sis&apos;&apos; of isothermal test data by GE-NMPO28 and differential temperature tests28 on are-melted tungsten. The activation energies determined (170,000 and 160,000 cal mole-1) were greater than that for self-diffusion and are in better agreement with the value found by Green. The model of dislocation climb is associated with ; activation volume of 1 b3 according to Conrad.22 The activation volume calculated for the creep data of fou independent investigators in the temperature range where the activation energy for creep is near that for self-diffusion (Klopp et al. and Refs. 26, 28, and 30) is 110 to 2800 b3 and is in the range suggested for nonconservative motion of jogs. In the case of tungsten this mechanism may be preferred over that of dislocation climb. Sincere appreciation is expressed to Messrs. J. E. Flinn and F. L. Yaggee for the helpful discussions and suggestions in the preparation of this communication.

Jan 1, 1967

By J. E. Flinn, E. R. Gilbert

Recently, Green reported some high-temperature, 0.60 to 0.89 T3, creep results on unalloyed tantalum. From his study he determined an activation energy of 114 kcal per mole for steady-state creep using the T, differential-temperature method25 for testing. He accounted for the difference between this value and 110 kcal per mole for self-diffusion by the temperature dependence of the elastic modulus, 4 kcal per mole. This reported activation energy is in apparent disagreement with that which can be determined from Green&apos;s isothermal data, which we have plotted as log creep rate vs 1/T in Fig. 15. From this figure, it appears that there are two activation energies of -185 kcal per mole at the higher temperatures and -75 kcal per mole at the lower temperatures. It is not obvious that the data can be fitted to a straight line corresponding to an activation energy of 114 kcal per mole. One is led to two possible conclusions regarding Green&apos;s study: there are two discrete processes controlling the creep rate from 0.60 to 0.89 Tm or the apparent activation energy is not discrete but continuously increases with increasing temperature. The former speculation would be in closer agreement with the high-temperature creep behavior of tungsten.28 The apparent disagreement between the activation energies determined by the AT and isothermal methods certainly suggests that further study on the high-temperature creep of tantalum is needed.

Jan 1, 1967