Search Documents

By Leonard R. Weisberg

The general properties of diffusion in GaAs are reviewed. A total of .fourteen atoms have been studied to date, and activation energies for eleven reported are (in ev): Ga (5.6), As (lo), Zn (2.49), Cd (2.43), Sn (2.5). Mn (2.75). S (4.0). Se (4.2). Ag (0.33), Cu (0.52). and Li (1.0). The dijfusion of the first eight atoms, as measured by Goldstein using radiotracer techniques, established the principle of sub lattice diffusion for both host and impurity atoms, since atoms substituting for gallium have lower activation energies than those substituting for arsenic. The diffusion constant of interstitial copper has been determined for the first time by Hall and Racette using heavily doped p-type samples in which interstitial copper predominates. The anomalous behavior of zinc, which has an abrupt diffusion front and a strong dependence on zinc concentration, has recently been quantitatively explained in terms of an interstitial-substi-tutional equilibrium-diffusion mechanism. Lithium and silver diffuse mainly interstitially, while the diffusion of phosphorus is complicated by the continual formation ofGap. Junction-migration measurements indicate diiferent results than radiotracer techniques, and this is one problem that still remains, together with discrepancies in cadmium-diffusion results, and an unexplained dependence of impurity diffusion on arsenic pressure. BESIDES the usual insights into the behavior of atoms in solids, studies of diffusion in semiconductors have additional importance for device applications. In fact, the amount of diffusion studies performed on a semiconductor might serve as a simple barometer of its technological importance. Thus, the information accumulated for GaAs is not yet quite comparable to that for germanium or silicon, but considerable knowledge is available. The diffusion of at least fourteen atoms has been investigated for GaAs, and the diffusion constants have been determined for eleven of these. Interesting results were found for GaAs that elucidated problems not previously solved for germanium or silicon. For example, the interstitial-diffusion constant for the very fast diffusant, copper, was first measured in GaAs. Also, solutions for interstitial-substitu-tional equilibrium diffusion in extrinsic material were provided first for GaAs, although they are ap- plicable to germanium, silicon, and other semiconductors. Studies in GaAs firmly established the concept of sublattice diffusion, which should have application in many different compound systems. The purpose of this paper is to review briefly the general state of knowledge concerning diffusion in GaAs. First some general aspects of diffusion will be covered, and next attention will be given to the results mentioned above. Diffusion in GaAs has been measured both by radiotracer and junction-penetration techniques, and some discussion will be given of differences in the two sets of results. The comprehensive work of oldsteinl-' has served as the primary reference for the diffusion of well-behaved substitutional atoms in GaAs. Since a review has recently been presented of many of these reults, details of the diffusion data will not be repeated. Neither will any discussion be given of the experimental techniques of diffusion measurements, except to note here, in passing, the simplicity and effectiveness of the precision lapping device employed by Goldstein.' I) EXPERIMENTAL RESULTS AND DISCUSSION A) Summary of Experimental Results. Diffusion constants D = Do exp(E/kT) have been determined for eleven atoms in GaAs: a,' AS,' n,'" Cd,'" n,' n,8 s, e,' Ag,O CU," and i." These results are summarized in Table I. In addition, diffusion was investigated for phosphorus,4'12 germanium,I3 and ilicon,' but Do and E were not determined. Radiotracer techniques were used for most of the studies except for lithium, where a combination of chemical and electrical techniques was employed, and germanium and silicon, where junction-migration techniques were applied. The diffusion of zinc, manganese, sulfur, and tin was also measured using junction migration.13'14 Goldstein demon-

Jan 1, 1964

By C. S. Matthews, P. Hazebroek, F. Brons

A method has been developed for calculating the average pressure in a bounded reservoir. The reservoir is first divided into the individual drainage volumes of each well, by using the criterion that at steady state each individual drainage volume is proportional to a well's production rate. The average pressure in each drainage volume is then calculated by a method developed in the report. By volumetrically averaging these individual drainage volume pressures, the average pressure in the entire reservoir is obtained. To calculate the average pressure in each drainage volume, a correction is applied to the ordinary extrapolated pressure, i.e., the pressure obtained by extrapolating to infinite time the linear portion of the graph of closed-in pressure versus log[?t/(t + ?t)], where At is the closed-in time and t the production time. The correction, which is a function of the production time, is presented in graphical form for different shapes of the drainage area (horizontal cross section of the drainage volume). INTRODUCTION It is important to be able to find the volumetric average pressure in a reservoir so that the size of the reservoir may be determined from material balance calculations. It is also desirable to be able to find the approximate distribution of pressure within a reservoir for detection of fluid movement. 'The purpose of this paper is to present a method for calculating both the average reservoir pressure and the approximate distribution of pressure within a bounded reservoir — that is, a reservoir with no water drive. In reservoirs where the pressure builds up rapidly after wells are shut in, the determination of average pressure generally poses little problem, for one often need only average the final buildup pressures. It is when pressure buildup is slow that difficulties arise. For practical and economical reasons, the time allowable for closing in wells is limited. If at the maximum allowable closed-in time the pressure has not reached a con- stant value (and this is more often the case than is generally realized), calculation of average pressure presents difficulties. One method which has been used in the past for obtaining an "average" pressure involves the extrapolation to infinite closed-in time of curves of pressure as a function of closed-in time. Contour maps of these extrapolated pressures are used in arriving at an average pressure. In some cases wells are shut in for arbitrary Lengths of time, and contours of these arbitrary shut-in pressures are used for calculation. All these methods are somewhat arbitrary and suffer from lack of fundamental theoretical basis. For example, no precise meaning can be given either to the extrapolated or to the arbitrary shut-in pressure. Furthermore, there seems to be no reason why either set of pressures should be contourable. It was because of the vagueness of these old methods and their lack of a sound theoretical basis that, the present method for determining average reservoir pressure was developed. SYMBOLS A = drainage area of a well, cm' A, = drainage area of well j, cm2 A , = total drainage area of all wells in a reservoir, cm2 c = coefficient of compressibility of fluid at reservoir temperature, atm-1 f = hydrocarbin-filled porosity, fraction h = net thickness of formation, cm k = permeability of formation, darcy p = pressure at a well, atm p0 = initial reservoir pressure, atm p* = pressure obtained by straight-line extrapolation of linear portion of plot of p versus log [At/(t+ At)] p = volumetric average pressure inside drainage volume of a well q = volume flow rate at a well, cc/sec a: prevailing reservoir conditions q, = total volumetric flow rate from a reservoir, cc/sec at prevailing reservoir conditions rw = radius of well bore, cm T = kt/fµcA t = corrected time of production of a well, seconds, defined by

Jan 1, 1955

By J. N. Ong

The oxidation of tantalum is assumed to occur by four simultaneous first-order chain reactions; solution of oxygen in the metal, nucleation and growth of a suboxide phase at the metal surface and two phase boundary processes which give rise to two different modifications of Ta2O5. A rate equation is then formulated which accounts for all principal details of oxidation behavior, unifies and correlatss the results of the principal investigations, and expresses the rate as a function of pressure, temperature, and time. For purposes of practical application, the expressions: may be used below and above 700T respectively to determine the regression rate of tantalum in cm per hr, where T is to be expressed in "K and P in atmospheres pressure of oxygm. The equations will reproduce within a factor of two all reported experimental rates in the temperature range 475" to 1400°C and pressure range 0.000026 to 40.8 atm 0. 1 HE oxidation kinetics of tantalum have been investigated thoroughly in the temperature range 50" to 1400 and pressure range 1.32 x l0-' to 40.8 atm pressure of 0,and as a result, the principal features of the oxidation of tantalum are known and have been described in considerable detail. The exposition of the characteristics of the behavior of tantalum in oxygen has been largely descriptive, with considerable attention being devoted to elucidating the structure, composition, and morphology of the resulting oxide products. The weight change vs time results have usually been reported graphically. With the exception of the linear rate below 700 C and the process of diffusion of interstitial oxygen into the metal,6'9 no effort has been made to characterize the reaction rates in functional form. It is the purpose of this paper to demonstrate that when the fundamental principles of chemical kinetics are rigorously followed, the results of the principal investigations will be correlated and unified and the oxidation rale may be expressed as an analytic function of pressure, temperature, and time. In addition to accounting for the temperature and pressure dependency of the oxidation rate, the following minor features will be accounted for as a result of the basic analysis: the initial period of nonline-arity in rate at low temperatures, which decreases with increasing temperature;3'5'7 the rate vs temperature cusps obtained at constant pressure"5 the concentration gradient of oxygen interstitially dissolved in the metal.5'7 The analysis will also be consistent with the overall chemistry of reaction and with all visual, microscopic, and X-ray observations. DEVELOPMENT OF RATE EQUATION As a fundamental basis for the treatment of metal-oxygen reaction kinetics, it will be regarded that the rate of reaction is determined by the concentration and not the activity or chemical potential. It is further assumed that all metal-oxygen reactions are first order complex chain reactions" whose rate is dependent on the concentration of an intermediate species which is a complex function of the oxygen concentration, or oxygen pressure, only. The basic rate equation is therefore: is the weight of oxygen consumed by the reaction, A is that portion of the area of the metal over which the reaction is occurring (cm2), t the time (sec), kf the rate constant for forward reaction, C the intermediate species concentration and f '([Oz]) and f (Po,) are complex functions of oxygen concentration and Oxygen pressure respectively. Reverse reactions may usually be neglected in metal-oxygen systems since a large negative free energy generally accompanies these reactions. Eq. [I] takes different forms when the rate of reaction is controlled by diffusion into or through a phase or at a phase boundary. For diffusion into the metal phase, the rate is given by1' Di is the diffusion coefficient of species diffusing into the metal phase (cm2/sec), and Q is the corre-

Jan 1, 1962

By P. G. Carpenter, R. G. Nisle, F. W. Bubb

It is the purpose of this paper to describe an electrical computer which has been constructed to solve the equations for vapor-liquid equilibrium in multi-component systems. The instrument consists of seven component-computing units each with proper indicating means and power supplies. Each unit is a resistance network with a voltage matching servomechanism, and each provides an output voltage proportional to the mol-fractions for vapor and liquid phases. These voltages are summed and matched with a reference voltage to provide the solutions. Any reasonable number of such units may be put together to make a computer. The theory and operation of the computer is discussed. A number of applications and examples of computer results are given. The computer yields the over-all vapor or liquid fraction to a probable error of 0.002. An interpolation method is described which reduces the probable error to 0.0002. INTRODUCTION The process known as Flash Vaporization may be described as follows: A mixture of known composition of relatively volatile components is allowed to come to thermodynamic equilibrium at some given temperature and pressure by any path whatsoever. In general, a vapor and a liquid phase will be present at equilbrium. The problem is to determine the fractions of the mixture in the vapor and liquid phases, and to determine the mol fractions of the various components in each phase. The relations that exist between these various quantities at equilibrium are well known and will be given later. These equations are difficult to solve, yielding only to trial-and-error methods. The virtues of the analog computer are its speed and the automatic character of its calculations. The device is operated by turning a crank which varies the value of v, the total vapor fraction. This actuates a number of servomechanisms which perform automatically all computations. This computer facilitates the solution of such problems as (to mention a few] : Analysis of separator operation, studies of changes in composition of reservoir fluids with pressure decline, and analysis of natural gasoline plant operation. Examples of some of these problems are given in detail later. STATEMENT OF THE PROBLEM The essential equations are xm = 1+(km-l)v . . (1) Kmzm ym =------------------------. . (2) 1 + (Km-l)v ; n n S S xm= --------------------=1 1 + (Km-l)v m = 1 m=1 ......(3) n n S = S Km Zm =1 S ym 1 + (Km-1)v m=1 m = l ......(4) where ym = Vm/V = Mol fraction of m-th component in vapor phase. xm = Lm/L = Mol fraction of m-th component in liquid phase. Zm = Fm/F = Mol fraction of m-th component in mixture. KM = ym/xm = Equilibrium constant for m-th component at the given temperature and pressure. v = V/F = Mol fraction of vapor in the mixture. F = Total mols of a mixture of n components. V - Total mols of vapor in the mixture. L = Total mols of liquid in the mixture. Fm = Total mols of m-th component. Lm = Mols of m-th component in liquid phase. Vm = Mols of m-th component in vapor phase. Assuming the total mols of the mixture F and its composition (all the Fm) to be known — either from a quantitative analysis or from the amounts put together to make the mixture — the zm = Fm/F may be regarded as known quantities. From the given temperature and pressure of the mixture the K., are known (principally from experimental data). The primary problem is then: having given all the zm and Km to calculate v and all the fractions xm and ym. Having calculated these ratios, it is a simple matter to calculate V. L. Vm, Lm from their defining equations. SOLUTION OF THE PROBLEM Basic Computer Unit The computer consists of n units, each of which provides voltages xEo and YE, proportional to x and y as in

Jan 1, 1950

By W. S. Hyler, W. L. Bruckart

The powder metallurgy and arc-cast types of wrought molybdenum stock were studied in rotating beam fatigue. Endurance ratios of unnotched specimens after 5x10 cycles were found to be 0.74 and 0.81, respectively. Fatigue notch factors (K,) of 1.36 and 2.05, respectively, were obtained on V-notched specimens having a theoretical stress-concentration factor Kt = 3.1. MOLYBDENUM and its alloys are rapidly gaining importance as engineering materials of construction for high temperature applications. AS a wrought material, molybdenum and its alloys are outstanding above 1800°F. It is necessary, however, with molybdenum to provide elevated temperature oxidation protection. This is currently being studied by a number of investigators. Thus far, there have been no published data on the fatigue behavior of molybdenum either at room temperature or at elevated temperature. The knowledge of fatigue behavior, particularly notched fatigue behavior, becomes increasingly important as the engineering importance of a material increases, especially from the design standpoint. The work described here was directed toward determination of the notched and unnotched rotating beam fatigue characteristics of molybdenum at room temperature. Also, it was desired to learn whether any difference exists between wrought molybdenum of the powder metallurgy and the are-cast types. Fabrication History of Bar Stock For this work, bar stock was procured from Climax Molybdenum Co. (are-cast type) and from Westing-house Electric Corp. (powder metallurgy type). The processing histories of these materials are listed in Table I. The first lot of are-cast molybdenum was deoxidized by aluminum additions. The second lot was deoxidized by carbon." Both lots of the powder Microstructure A microscopic examination of these materials, as received, revealed the powder metallurgy molybdenum to be free of inclusions and to be finely fibered. On the other hand, the are-cast C molybdenum possessed numerous carbide inclusions, and the arc-cast A appeared to have coarse grained fibers. These are shown in Fig. 1. Owing to the large apparent difference in fineness of fiber between the are-cast A and the powder product and are-cast C, a check was made of the recrystallized grain sizes to determine whether these differences persisted. A series of recrystallization anneals revealed mixed grain sizes in both the are-cast A and the powder product. Variations in grain diameter of as much as 50 times were found after full recrystallization at 1200°C for 30 min. There was no distinction between the two. Arc-cast C was not checked. Room Temperature Mechanical Properties Prior to subjecting the materials to fatigue tests, room temperature mechanical properties were determined. These data, together with the tensil: strength of the notched bars, data,togetherare listed in Table tensil: There were no significant differences in the tensile strengths of the two types of molybdenum. However, the powder metallurgy material was generally more ductile and yielded at a somewhat lower stress than the arc-cast material. Tensile strengths of the notched specimens indicated a slightly greater notch strengthening effect in the are-cast C than in the powder product, lot 2. Room Temperature Fatigue Properties Sample Preparation: Specimens were prepared for use in R. R. Moore rotating beam fatigue machines. These specimens had a continuous radius test section with a minimum diameter of 0.260 in. The specimens were prepared by grinding to shape and polishing to produce a smooth unworked surface. Polishing consisted of rotating a disk in the continuous radius section on which were attached, successively, emery belts of 240, 400, and 600 grit. The specimens were slowly rotated during the operation. The belt speed was high, and rotation was such that polishing scratches were longitudinal. Unnotched Fatigue Properties: Unnotched fatigue data were obtained from a number of specimens of each of the four lots of molybdenum. The results of the tests are presented in Fig. 2. From this figure,

Jan 1, 1956

By M. T. Hepworth, T. P. Papazoglou

Hydrogen diffusion in hcp and bcc titanium in the temperature range 610° to 900°C was measured by reacting cylindrical specimens of high-purity titanium at constant temperature with small constant pressures of hydrogen of the order of a few Torr. In the range 0 to 6.7 at. pct H in hcp titanium, the diffusion coefficient D in sq cm per sec is given by: where T is in °K and R in cal per g-atom OK. Only three measurements of diffusivity in bcc titanium were made: at 900' : the diffusivity in ß titanium was found to be 1.8 x 10-4 sq crn per sec. McQUILLAN1 has studied the constitution of the Ti-H system and reported the concentration of hydrogen as a function of pressure of hydrogen and temperature. He found hydrogen absorption was reversible and rapid at temperatures above 500°C in both crystallographic modifications, the a (hcp) phase and the ß (bcc) phase, with the latter phase stabilized by hydrogen additions. Wasilewski and Kehl 2 measured the diffusivity of hydrogen in hcp and bcc Ti-H alloys. For bcc alloys, they reacted cylindrical specimens of titanium with hydrogen gas at 1 atm pressure. They measured the rate of absorption of hydrogen isothermally on specimens ranging from zero hydrogen content when above the transformation temperature to hydrogen contents sufficient to stabilize the bcc phase when below the transformation temperature. By integrating Fick's second law for a cylindrical geometry, assuming a diffusion coefficient independent of composition, they calculated diffusivities as a function of temperature. In the range 650" to 1000°C their results in bcc titanium yielded the expression: where D is in sq cm per sec, T in OK, and R the gas constant in cal per g-atom-OK. Their studies on the hcp a phase were conducted by exposing long cylindrical titanium specimens to hydrogen at a sufficiently high pressure to establish a surface layer of fl phase 80 to 120 u thick for the duration of an isothermal experiment. Concentration profiles were then determined by hydrogen analysis from chips obtained upon machining quenched specimens. The diffusivity of hydrogen was then calculated yielding the expression: The present investigation was undertaken principally to redetermine the diffusivity data on a titanium by a direct absorption from the gas phase at sufficiently low hydrogen pressures to avoid the presence of a surface bcc phase and also to check the diffusivity data at low hydrogen concentrations in the bcc stmcture. EXPERIMENTAL PROCEDURE High-purity electrolytic sponge titanium (99.9 pct) of BHN 91 was arc-melted, forged, rolled, and machined into cylinders of final hardness, upon annealing, of BHN 92. Cylindrical specimens of two sets of dimensions were used: 1.524 by 2.931 cm and 1.905 by 3.663 cm diam and length, respectively. The experimental apparatus used for annealing specimens and for admitting known amounts of hydrogen at constant pressure is shown in Fig. 1. In a typical experiment, a specimen was lightly abraded with fine-grain silicon carbide, washed with water and detergent, and then thoroughly degreased with reagent-grade tri-chlorethylene and acetone. The specimen was suspended by a Chrome1 wire in the center of a quartz furnace tube, without contact with the tube, and heated to 925°C (as measured by a thermocouple adjacent to the specimen, outside the tube) for a period of at least 4 hr under a vacuum of from 10-6 to -7 Torr as measured by an ionization gage. The furnace temperature was then adjusted to the temperature chosen for the particular diffusion study and the specimen allowed to equilibrate thermally. Hydrogen was then ad-

Jan 1, 1969

By A. G. Guy

In a special type of electrodiffusion experiment the entire specimen is at constant temperature. The vacancies are then created at one end of the specimen and flow to the other end where they are destroyed. The energy required to operate the corresponding sources and sinks results in a super saturation of vacancies at one end of the specimen and an under -saturation at the other end. This gradient in vacancy concentration produces a chemical-diffusion flux, J(D), of vacancies in addition to the flux, J(E), caused by the electric current. Using the experimental data of Ch. Herzig,5 it is shown that a super saturation of 1 pct can cause J(D) to be about 25 pct of J(E), and thus account for the anomalous behavior of short specimens. THERE are advantages in a phenomenological treatment as contrasted to an atomistic one. For example, there are many reactions involving the element copper that are quantitatively and usefully understood from a phenomenological viewpoint, even though electron theory is not yet able to predict whether copper should be fcc or cph.' Similarly, the lack of a complete atomic theory of vacancies need not prevent the development of a treatment based on thermodynamics and chemical kinetics. In the following discussion, models will be used for the vacancies and for the mechanisms by which they move or are created or destroyed, but the phenomenological treatments are not dependent on the details of these models. The important product of a phenomenological analysis is an equation that describes the relation among the measurable quantities determining the behavior of the system. Based on such an equation, experiments can be performed to study these quantities in detail. Eventually atomic theory will supply an adequate theoretical base for the phenomena in question. A chemical species whose chemical potential was invariant during reactions in binary and ternary systems of varying composition would indeed be an anomaly. Yet this is the role that has frequently been assigned to the vacancies in pure metals, in dilute solutions, and even in ordinary alloys. The weight of available evidence supports the view that supersatu-rations of vacancies are on the order of 1 pct in usual diffusion processes.' Although it may seem strange that such a small deviation from equilibrium could significantly affect the thermodynamic properties of the vacancies, it is well known that the driving force for an ordinary chemical reaction is the departure of the components from their equilibrium state. Similarly, the reaction of vacancies that results in the formation of pores in various types of diffusion3 can reasonably be attributed to the supersaturation of vacancies. In the present paper, it will be shown that small degrees of supersaturation can significantly affect the kinetics of the diffusion process itself. ISOTHERMAL ELECTRODIFFUSION Two different experimental techniques are convenient for studying the behavior of vacancies in a pure metal in the presence of an electric field. Gilder and Lazarus4 studied self-electrodiffusion in a gold specimen in which only the central portion was at high temperatures. In this case the significant vacancy sources operate inside the specimen in the presence of a temperature gradient. In the laboratory of Th. Heumann a method has been devised for studying self-electrodiffusion at constant temperature.5 In principle, there is then no effect of thermodiffusion. It will now be shown why in this case the significant ATungsten disks \^^ ^^ ^Gold disks^^ ^^"T>.... "TTZA ? Up Iy'7 markers ///// ////// initial plane p I y/y/y of joining \I 7//// V7//A_________I_________V//// -I 0 +I Distance from center of specimen (a) -— Direction of flow of electrons —Direction of flw of vacancies %yy} Matano I z~irkendall /%W, vy/A inlerfacej "interface ///// -} HfiK- +/ 0 Distance from center of specimen (6) Fig. 1—-Schematic illustration of the essential features of isothermal self-electrodiffusion. (a) Initial condition of specimen. (b) Condition of specimen after electrodiffusion.

Jan 1, 1970

By B. Trafton, H. L. Hartman

An earlier paper1 introduced a rapid solution for mine ventilation network problems, employing the digital computer. In this paper, refinements and additions to the computer program which expand its scope of application and a solution to a ventilation network having a coal-mine configuration are presented. The revised program-written in FORTRAN for the IBM 7074 computer in Penn State's Computation Center-is applicable to any mine ventilation network employing natural splitting and having its fan located externally. It accomplishes the following: 1) Distributes an assumed total quantity of air flow among all airways of the network, performing iterations until an acceptable level of error is attained and the flow in each airway is known. 2) Computes the pressure drop (head loss) in all airways and the overall network and mine head losses. 3) Determines the operating point for the mine and a given fan, adjusting all air flows and head losses to the actual values. The basis of analysis used is the Hardy Cross method of successive approximation for solution of any fluid network, explained in detail in the earlier article. In brief, corrections to assumed air flows are computed for all the loops of a network; following distribution to the airways of each loop, the process is repeated until acceptable accuracy is attained. When a hand or desk-calculator solution is used, there is advantage to close approximation of the actual flows in the original assumption of air-flow distribution in the network. Since publication of the previous paper, in which assumption of initial quantities was also necessary in the computer solution, the program has been modified to make the first distribution itself, assigning the total quantity to one specified path in a split and then proceeding to correct the flows. The number of iterations required for a high-accuracy solution is little affected, at least in networks of few loops, 10 to 15 being sufficient to reduce the error to 1 pct or less. The principal limitation of the original program was that it contained no provision for handling a fan. Further, it did not determine the overall head loss in the network and was not designed for multiple parallel openings, such as occur in coal mines. The revised program was written to overcome these disadvantages and also to print out the adjusted flows for all airways, corrected for the system operating point. At the present time, the program can handle only an external fan. The procedure employed in solving a ventilation network is now as follows: 1) Input data consisting of airway characteristics and loop designation along with assumed direction of flow, for each branch of the network and any trunk airways, are fed into the computer. 2) Employing the resistance relationship, where R is resistance, K is friction factor, P is perimeter, L is length, A is area, and n is number of airways (where multiple intakes and returns are employed); values are computed for all airways. 3) An arbitrary, continuous route of flow through the circuit—from beginning to end—is selected and the path identified by airway numbers and read in. 4) Based on a total assumed quantity of flow Q (usually 100,000 cfm.) the computer employs the procedure described in the previous article in distributing quantities of flow among all the airways until a satisfactory level of accuracy has been attained (in the current program, the machine is turned off when the cumulative corrective error in cfm, is reduced to five times the number of airways). 5) Airway head losses H are calculated by the relation H = R Q2 and cumulated along the original designated flow route to find the overall head loss for the network. Trunk airways are then included to obtain the mine head.

Jan 1, 1964

By S. C. Miller, R. D. Carter, H. G. Riley

A merhod based on theoretical considerations is presented for estimating srabilized performance of gas wells from data obtained from short-term flow tests. Calculated results are presented for II gas wells evaluated with four-to five-hour flow tests. For seven wells on which comparison war available, satisfactory agreement between measured and calculated performance was obfained. Although the final results are favorable, diflerences in flow capacity k,h obtained from flow tests conducted at different rates were observed. This inconsistency did not seriously affect the calculated results in the test wells. INTRODUCTION In order to negotiate the sale of gas from newly drilled gas wells and to plan otherwise for the orderly exploitation of gas reserves, it is necessary to predict the long-term withdrawal rates physically possible for each well. Since low permeability gas wells do not stabilize except after a lengthy period of flow, the operator is unable, without wasting large amounts of gas, to measure directly the stabilized performance prior to contracting for sale of the gas and connecting the well to a pipeline. This problem is not new, but is one which has recently received considerable attention because of the rising importance of natural gas in the economy and new emphasis on the development of gas reserves. Thus, a method is needed for estimating stabilized performance from information obtained from short-flow tests in which stabilized flow is not attained. Such short tests would be suitable for wells not connected to a pipeline since the tests require only small amounts of time, labor and production of gas. Other requirements of such a method are short calculation time and use of a model of the gas well and its reservoir drainage area which is as consistent as possible with present theoretical, experimental and field knowledge. The method presented herein was developed in an effort to meet these requirements. It has been named the "two-Row method". In addition to the section on the evelopment and description of the method, a description of field application of the method and a discussion of well test results are in- cluded. Finally, there is a section on an additional consideration prompted by some of the test results. The Appendix contains certain details of the derivation of the two-flow method and a sample calculation for Test Well 1. MATHEMATICAL DEVELOPMENT AND DESCRIPTION OF THE METHOD The method is based upon the following assumptions: (1) The system can be described as a finite cylindrical gas reservoir produced by a single cylindrical well located at the center of the reservoir. Flow is strictly radial. (2) Flow obeys Darcy's law except very near the well where a quadratic flow law is assumed. (3) The system is also characterized by a laminar flow skin factor. The use of a quadratic or turbulent flow term is more fully justified in Refs. 1 and 2. When such a system is opened at the wellbore to constant rate flow from an initial condition of uniform pressure distribution, the radius or drainage quickly passes through the non-Darcy flow and "skin effect" region imrnediately surrounding the well. After the drainage radius or "transient" has passed beyond this region and before it reaches the external boundary, drawdown can be expressed by'." After the drainage radius reaches the external boundary, stabilized flow exists and drawdown is given by'," Eqs. 1 and 2 form the basis for the method. The validity of these equations will now be examined in the light of some well test data previously available. It has been found by actual measurements on gas wells that the relationship between drawdown and flow rate, that is the performance curve, is closely approximated by' Eq. 3 is a straight line on log-log graph paper. Eq. 2 and, for any lixed value of time, Eq. 1 plot as a curve which is slightly concave upward on log-log graph

By R. M. Rose, D. P. Ferriss, J. Wulff

Single crystals of tantalum, tungsten, and the binary alloys thereof were grown by electron-beam zone melting and tested in tension between 77° and 373°K. The ductile-to-brittle transition temperature increased sharply with solute concentration, for both tungsten-rich and tantalum-rich alloys; apparently, the maximum transition temperature exists near the 50 pet concentration. Twinning was observed to be a consequence of cleavage-crack propagation. Cleavage occurred on(100)planes, and extensive tearing and "river markings" were present. The onset of brittleness could not be associated with any change in deformation mode, the yield point phenomenon, or the large rise in yield stress with decreasing temperature. Cleavage was observed to initiate at either the specimen surface, or subgrain boundaries. The observation was made that the grain boundaries of Polycrystalline materials are probably even better as nucleation sites for fracture, and are responsible for brittleness at high temperatures and lower solute concentrations than those of the corresponding single crystals. TUNGSTEN and tantalum, the two high-melting neighbors in the third transition period of the periodic table, present striking contrasts. Although both are bee, tungsten has far higher elastic moduli, and the highest ductile-to-brittle transition (DBT) of all the pure metals.' Tantalum, however, is ductile at liquid helium temperatures. Large strain aging and yield point phenomena are observed in tantalum,2,3 but only to a lesser extent in tungsten.4 Tungsten work-hardens at a high rate,5 whereas this effect is negligible in tantalum.' Since the mutual solid solubility of these two metals is complete,7 and no ordering has yet been reported in the binary system, the possibility of depressing tungsten's DBT by tantalum additions was considered. Also, tungsten additions to tantalum could produce hard, high-melting alloys which might be ductile at room temperature. Such alloys offered excellent bases for complex alloy development, particularly with regard to the effect of rhenium on Group VI metals.8"10 In addition, there was the opportunity to examine flow and fracture modes in the vicinity of the DBT, with particular emphasis on the correlation between brittleness and fundamental deformation mechanisms. Myers11 has reported that recrystallized alloys with up to 10 pet W in tantalum are ductile at room temperature, but alloys with 20 pct W could not be swaged at temperatures as high as 500°C. The whole spectrum of Ta-W alloys has been surveyed by Braun, Sedlatschek, and Kiefferl2 using sintered powder compacts. In addition to some physical properties they reported a maximum hardness at about 50 at. pet W, dropping to much lower hardness levels for the pure components. They were able to cold work tantalum-base alloys with as much as 30 pet W and to warm work tungsten-base alloys with up to 20 pet Ta by procedures similar to those used for pure tungsten. EXPERIMENTAL In order to study deformation modes and to keep the impurity concentrations at a minimum, electron-beam zone melting was used to prepare single crystals of high purity. The electron-beam zone melting apparatus is similar to that described by Calverley, Davis, and ever" and incorporates a grid control rectifier to maintain control of emission current. Pure metals were refined by making three zone passes at 4 mm per min up the rod. Alloys were prepared by melting together wire of one constituent on a previously melted rod of the major constituent. Subsequently, two zone passes were made in opposite directions to level the composition after the manner suggested by pfann.14 Purification of the starting materials was accomplished by evaporation of the several impurities, and the degree to which this was accomplished is indicated in the chemical analyses of Table I. The uniformity of composition along a zone levelled rod was within the i 1 pet accuracy of the fluorescent X-ray technique used to determine composition. Crystallographic orientation of the

Jan 1, 1962

By R. J. Teitel

The Pb-U system has been investigated by X-ray, thermal, and microscopic analyses. Two pyrophoric intermetallic compounds were found; UPb3 and UPb. The crystallographic structure of UPb3 is reported. A phase diagram was drawn for the system showing three eutectic and one syntectic reaction. AS a result of the interest in uranium metallurgy, the Pb-U system has been investigated. Practically no quantitative data appears in the literature concerning this system. Chips were machined from bars of high-purity uranium obtained from Ames Laboratory and the lead was 99.9 pct high-purity grade from Bel-mont Smelting & Refining Co. The chips were pickled in a dilute solution of nitric acid, washed in water and then alcohol, and dried prior to use for alloy preparation. The lead was simply sheared from the pig with no chemical treatment. Lead and uranium were placed in graphite crucibles and loaded into an induction furnace. In order to insure good alloy formation, the lead was loaded above the uranium to wash the uranium chips into the melt. After the charge was held at the desired preparation temperature for 15 to 30 min, it was allowed to cool in the crucible. The entire melting operation was accomplished in an atmosphere of helium (double charcoal refined grade with no further treatment). The procedure was to outgas the furnace by first evacuating it to 10-% m Hg, admitting helium to 8 psi gage and sealing the furnace off. An optical pyrometer focused on the charge through a hole in the crucible cover was employed for thermal measurements. Usually a 50° to 70°C correction had to be added to the instrument value. The most suitable alloy preparation temperature range was found between 1220° and 1280°C. The reason for this is explained by the phase diagram, Fig. 7. Above this temperature range two liquids separate preventing alloy formation upon cooling at a reasonable rate. On the other hand, below this range the uranium was coated with a layer of both compounds and the rate of reaction was diminished by diffusion processes. The fastest reaction rate occurred in this narrow range of temperature because only UPb coats the uranium and it evidently does not protect the uranium as well as UPb,. Apparatus and Procedures The development of apparatus and techniques was determined by the peculiar properties of the alloys. All the alloys were pyrophoric, and, consequently, the alloys had to be protected from oxidation by handling them in a dry box or inert atmosphere furnaces. Other operations performed outside of this equipment had to be done as quickly as possible. Thermal analyses were conducted in a furnace and arrangement similar to that described by the author.' Graphite crucibles having a thermocouple well in the bottom were used to contain the prepared alloy, while a helium atmosphere protected the melt from oxidation. The furnace was outgassed by pumping a 10-3 mm Hg vacuum until the furnace had been heated to 400°C. Then helium was admitted and the furnace sealed under 8 psi gage pressure. Chromel-alumel thermocouples were used for thermal analysis and the accuracy of thermal measurements was judged to be k 5°C. At first, an attempt was made to use platinum-platinum-rhodium thermocouples, but evidently lead vapors attack it at high temperatures forming low melting alloys. These lead vapors had direct access to the thermocouple because the graphite crucible was porous. A differential temperature control was used to regulate thermal analyses. The thermal gradient across an alundum tube placed between the windings and the crucible was maintained constant by a controller.

Jan 1, 1953

By W. C. Ellis

IN an ingot of germanium solidified progressively from the bottom by a method' devised by J. H. Scaff and H. C. Theuerer, examination discloses ex- IH. C. Torrey and C. A. Whitmer, Crystal Rectifiers. p. 364. 1948. New York. McGraw-Hill Pub. Co. tensive twin relationships. A slice taken about midway in a cylindrical ingot approximately 3/4 in. in diameter and 1 1/2 in. in length, after etching and polishing to delineate different orientations, is shown in fig. la. The orientations of the regions designated with letters were determined by Laue X-ray back-reflection methods and plotted on stereograms. The orientation relationships of the designated components were determined by comparing the values for the nine angles between the three respective cube poles of each orientation with the values calculated in a formal treatment of octahedral multiple twinning in cubic crystals. Germanium, diamond cubic in structure, has the octahedral twinning habit. I am indebted to K. H. Storks for the calculation of the angular relationships of cube poles in octahedral twinning to four orders. The set of values for the nine angles is distinctive of the order of twinning. Confirmation for a multiple twin relationship is obtained from the determined indices of the mirror planes. These exist for all orders of twinning through the third and for some of the fourth." If the orientations are not those for twinning, then the nine angles for the cube poles do not correspond in value and distribution to those formally calculated. Since all nine angles must correspond in values and distribution of values to the calculated ones, the Laue method with a precision of one or two degrees is adequate to establish or deny with confidence the existence of multiple twins. Regions of the same orientation can also frequently be identified by the optical reflections. For example, the narrow twin bands in A at 2: 00 o'clock have orientation B1. The interesting twin relationships shown schematically in fig. lb were found for the section in fig. la. Most of this slice was made up of orientation, A, appearing in widely separated regions. Of the four possible first order twins of A, three were present in this slice and are designated by B with subscripts. Four second order twins of A are designated by C with subscripts. It is worth emphasizing that two of these, C, and C,, are second orders of A through the first order missing in this section. Two regions, designated by D with subscripts are third order twins of A. D2 is a sixth order twin of D,. All orientations found were ascribable to multiple twinning, as if the orientations in this section originated earlier in solidification from a single nucleus. The presence of extensive multiple twinning to- gether with the absence of other orientations appears to be a characteristic result in many instances of solidification of germanium. Fertile conditions for twinning probably exist in the array of atoms in the diamond cubic structure. The high order twins in the section shown, it is reasonable to presume, are the result of progressive octahedral twinning through several generations in the material freezing earlier. Some of these generations do not appear in this section, it is suggested, because their orientations with respect to the growth direction were unfavorable for survival. The extensive twinning is probably promoted during solidification by constraint imposed by volume increase, and to an irregular growth front. Twinning may also be due, in part, to the attempt of the solidifying region to attain the orientation most favorable to the growth direction (maximum thermal gradient). The nucleus for a twinned orientation could have its origin in a "stacking fault" on the close-packed octahedral plane. Further work is being done in which the course of twinning throughout an ingot is being studied. The results appear to be important in the general theory of solidification.

Jan 1, 1951

By C. E. Birchenall, G. Derge

IN studying the kinetics of slag-metal reactions, it has become increasingly apparent that a complete knowledge of all aspects of interface phenomena will be required to clarify these processes adequately. A part of the data required for all oxidizing slags used in the refining of impure iron to steel is the rate of exchange of iron between the slag and metal phases. The simplest equilibrium of this type is that between liquid iron and iron silicate slag when both phases are equilibrated with the silica crucible containing them. In this study, the exchange of iron between the two liquid phases has been measured by using radioactive Fe as tracer. The observed high exchange rate indicates the ease with which such interfaces may be crossed in refining processes. The experiments required for this study involved relatively minor and normal modifications of well established techniques. Melting was accomplished by induction heating in heavy walled, fused silica, rotating crucibles." Rotation was used to minimize dilution of the slag by crucible erosion. The melting stock was ingot iron. Slags were prefused in iron crucibles using chemically pure grade oxide reagents. The special conditions applying to the active iron additions will be described for individual experiments. Suitable precautions were observed with regard to protection against radiation hazards.' Identification of Activity This particular study was initiated several years ago, but held in abeyance because the radioactive iron available at that time contained a foreign activity which concentrated in the slag so as to mask the iron activity and render the results unin-terpretable. With the benefit of experiences from other studies, radioiron of suitable purity can now be prepared. As received from the Isotope Div. of A.E.C., the unit appears to contain Fe + Fe50 and a high activity cobalt. The cobalt and iron are separated by repeated ether extractions, and the purified iron stored for sufficient time for the Fe'7 44 day half-life) to decay to negligible activity, leaving only Fe as the active isotope. This was confirmed by the absorption curves for slag and metal samples from heat No. 1. In Fig. 1 the experimentally measured points for the absorption curves of slag and metal are fitted by the calculated slope of 5.8 mg per sq cm in aluminum, which corresponds to the Mn K doublet emission resulting from the capture of an orbital electron of Fe" by the nucleus. This shows that the radioactive species are the same in both slag and metal and identifies it as Fe55 Description of Individual Heats For the first heat, the radioactive iron was contained in the prefused slag. The active Fe55 had been precipitated with carrier and was added as Fe2O3 during fusion in the iron crucible. The charge of 500 grams of ingot iron was melted, the protective iron sleeve added, rotation started, and a prefused, inactive slag was added to allow the system to approach equilibrium. After about 5 min at temperature, the wash slag was skimmed and the active slag added. Initial slag and metal samples were taken as soon as the slag was molten. Additional samples were taken at indicated intervals thereafter. Slag samples were dipped with a small iron spoon and metal samples were sucked into 1/4 in. ID silica tubes. The slag and metal samples were pulverized and mounted for radioactive counting. The slag samples were analyzed chemically for FeO. Temperature was read after the last slag sample by immersion of a Pt-Pt-Rh thermocouple protected by a fused silica sheath. Since no significant variation of rate with temperature was observed

Jan 1, 1954

By Kehsin Kuo

DIE steel with 1.47 pct C, 0.42 pct Mn, and 8.22 pct W contains Fe,C and WC in the annealed state and WC embedded in a martensitic matrix in the hardened state (quenched from 800°C). The presence of a large amount of undissolved angular WC carbide particles makes this steel extremely wear-resistant. However, the WC carbide begins to dissolve in the austenite when the heating temperature is raised above 1000°C, but this solution process is still not complete at 1250°C although partial melting has already taken place at the austenite grain boundaries. At 1325 C, this steel contains austenite and liquid metal but no WC carbide, Fig. 1. On quenching into brine, the liquid metal solidifies into a constituent which etches differently from the austenite-martensite matrix. This constituent appears lighter than the matrix at the 1250°C stage when etched with 2 pct nital but darker in Fig. 1 when etched with Murakami's reagent. Such a constituent has also been found in burnt tungsten high speed steels. This constituent has been isolated from the matrix by anodic dissolution in a 5 pct HC1 electrolyte and was found by the X-ray powder method to be iso-morphous with &manganese, Table I. The agreement in line intensity between these two phases indicates that there is no ordering in this new phase. The calculated cube edge of this phase is 6.395A, somewhat larger than that of ß-manganese (6.30A). Besides this phase, the anodic residue also contains free carbon from the matrix and amorphous WO,, although the X-ray powder pattern of the residue shows the presence of only this phase. This makes it impossible to determine directly the composition of this phase from the residue. Attempts to prepare an alloy containing only this phase have hitherto been fruitless. The following experiments serve to show that this phase can only be obtained by drastic quenching from the liquid state and that its tungsten and carbon contents must be lower than those in the Fe,W,C carbide. On air cooling from 1325oC, the liquid metal solidifies into a eutectic consisting of austenite and the Fe3W3C carbide. The nature of this dark-etching Fe3W3C has also been ascertained by X-rays. This eutectic becomes coarser on furnace cooling. The solidification of a specimen cooled for 5 min from 1325" to 1200°C and water quenched is still not complete and the fine eutectic is formed during the quenching. Since this new phase can be replaced by a eutectic of austenite and the Fe,W,C carbide, its composition must lie somewhere between the latter two phases. The special features of this phase are: 1—it can only form directly from the liquid state by drastic quenching and 2—it is a ternary phase of iron, tungsten, and carbon. The ß-manganese structure of this phase is also interesting. In the Periodic Table, iron and tungsten lie on opposite sides of manganese; a certain combination of these two elements, in the presence of carbon, may have average properties similar to those of manganese. It has already been established that the x phase, a ternary one of iron, chromium, and molybdenum, has the a-manganese structure."' A similar phase has also been found in the ternary system Fe-Cr-W and in the following binary systems: Re-Cb, Re-Ta, Re-Mo, and Re-W.6 References 1 G. D. Preston: Philosovhical Magazine 11928) 5. 7th Series, p. 1207. - K. W. Andres: Nature (1949) 164, p. 1015. 35. G. McMullin, S. F. Reiter, and D. G. Eberling: Trans. ASM 11954) 46, p. 799. ' J. S. Kasuer: Acta Metallurgica (19541 2, p. 456. • H. J. Goldschmidt: Trans. ASM (19541 46, p. 807. P. Greenfield and P. A. Beck: Private communication.

Jan 1, 1957

By Larry Kaufman

A useful method for estimating the entropy of NaCl type compounds has been developed by combining the Debye theory of specific heat with the Lindemann formula. This method pemnits calculation of the entropy from a knowledge of the atomic weight, atomic volume, and melting point of the compound. Comparison of the calculated and observed entropies of 35 NaCl type compounds over a wide temperature range yields good agreement thus establishing the usefulness of the method. Application of the method to other crystal types is also discussed. KUBASCHEWSKI~ has suggested a useful method of estimating the enthalpy of formation of inter-metallic compounds from a knowledge of the crystal structure, lattice parameter, and heats of sublimation of the component elements. This method is quite valuable in predicting the tendency for compound formation in a given system and is of particular value in considering systems for which no experimental data are available. In order to supplement the enthalpy calculation afforded by Kubschewski's method, it would be of value to be able to estimate the entropy of compounds so that the entropy and free energy of formation could be evaluated. Although a considerable body of such data determined by calorimetric and equilibration techniques exists at present, the difficulty of making such measurements (particularly at high temperatures) and the ever increasing number of new compounds which are of interest make the desirability of such a computational method clear. This paper deals with a method for calculating the entropy of a compound from a knowledge of the atomic weight, atomic volume, and melting point. Since these data are generally easier to determine and therefore more readily available than thermo-dynamic data, they should provide an additional tool in estimating the stability and reactivity of compounds at high temperatures. I) THE ENTROPY OF PURE METALS As a first approximation, the specific heat of a pure metal can be represented by a sum of terms where C,($/tI is the vibrational specific heat calculated on the basis of the Debye theory, lo4 TC, (0/T) is an approximation to the Cp - C, correction. T is the electronic specific heat, and Cp (p) is the magnetic specific heat. Although the Debye theory does not provide a complete description of the vibrational specific heat, i.e., the Debye temperature 0 is not a constant, it does provide a good approximation when used properly. Blackman6'7 has written several reviews on this subject discussing the different values of 0 determined by various methods (specific heat, X-rays, elastic constants, electrical resistivity, and so forth), as well as the observed temperature dependence of 0 with T. In the present case the 0 in Eq. [I] naturally refers to the specific heat 0. The variation of Debye 0 (as determined from specific heat measurements) with temperature is observed to exhibit its most radical behavior at very low temperatures. At higher temperatures the 8 values are observed to become more nearly constant. Blackman6'7 has presented the variation of 0 with temperature for several metals. His tabulation is reproduced in Table I. Reference to Table I shows that above a certain temperature level, 0 remains fairly constant. The value of 0 which is most useful in computing the specific heat, and hence the entropy, is the value of 0 in the range 0/6 « T « 0. In 1910, Lindemann' suggested a relation for estimating Debye 0 values. According to this formula, where V is the atomic volume, T is the melting point in 'K, M is the atomic weight, and CL is a constant. Mott and Jones9 present this formula and state that there is at present no satisfactory theoretical explanation, but that it does have an empirical basis. By using the valuesof 0 in Table I along with known1' values of M, T, and the molar volume VM, where the Lindemann constant can be evaluated. This

Jan 1, 1962

By T. J. Koppenaal

The surface slip line structure has been investigated by optical microscopy during the easy glide deformation of Cu-A1 single crystals as a function of composition, testing temperature and prior thermal treatment. The slip line structure in alloys of 2.5 or less at. pct Al is similar to pure copper, while alloys of 5 or more at. pct Al show slip tracings typical of a brass. Cross slip is observed in the "a-brass type" crystals and the amount of this type of slip increases rapidly with aluminum content. The cross slip tendency in 14 at. pct Al crystals is signicantly higher when tested at 150°C as compared to 4.2" and 77°K. In addition, cross slip is more pronounced in this alloy in the furnace cooled condition than when quenched from 450°C. The cross slip behavior is analyzed in terms of existing models for this type of deformation. In two previous investigations, reports of yield points1 and solid solution strengthening2 in Cu-A1 single crystals have been given. This paper enlarges upon the brief description of the slip line structure discussed in these reports. Since the easy glide slip line pattern of copper changes markedly with the addition of -30 pct Zn (-brass),3 a study of slip lines as a function of solute additions in Cu solid solutions should help to explain some of these differences. In addition to composition variables, slip lines have also been studied as a function of testing temperature and prior thermal treatment. PROCEDURE The slip line observations were made on single crystals of 0.5 to 14 at. pct Al using a light microscope. Details concerning the preparation and basic testing procedure are given elsewhere.' Prior to testing, each crystal was homogenized at 900°C for 24 hrs, furnace cooled, and electrolytically polished in a 60 pct phosphoric acid solution. The strain rate was-5 X 10 sec-I in each case. Unless otherwise stated, all observations reported upon here were made during the early stages of easy glide. RESULTS AND DISCUSSION The maximum magnification generally used was 1000X, and slip lines could be resolved down to about 105 cm (1000A). This resolution is appreciably less than obtainable with electron microscopy replica techniques where the replication can resolve slip lines down to about 40 A. Thus, while many qualitative features could be studied, the type of quantitative data obtainable with an electron microscopy replica study could not be made. Crystals containing 0.5 and 2.5 at. pct Al deform much like pure copper exhibiting a high density of fine slip lines that cover the entire gage length of the specimens. Fig. 1 shows the surface slip line pattern in a 2.5 at. pct A1 extended 0.01 shear strain at room temperature. Crystals having 5 to 14 at. pct A1 showed appreciably coarser slip lines that were generally concentrated at one end of the gage length. During easy glide these coarse slip lines propagate along the gage length in a manner similar to that observed in a-brass single crystals by Piercy, Cahn, and cottrel14 and refered to by them as a Lüders band propagation. Fig. 2 shows a typical area of these coarse slips in a 14 at. pct A1 crystal after a room temperature shear strain of 0.01. Two slip systems were apparently active. Using a method described by Cullity, the angles, between the specimen axis and the normals of the primary and secondary slip planes were determined to be 50 and 84 deg, respectively. 'From the Laue photograph, was 52, 81, 66, and 30 deg for the primary, cross, conjugate, and fourth slip planes, respectively. The secondary active plane in Fig. 2 is thus established as the cross slip plane, and the operative slip system is evidently the cross slip system, e.g., slip

Jan 1, 1963

By J. L. Moriarty, Irving Johnson, H. M. Feder

The crystallization of ternary alloys has been studied as a coprecipitation phenomenon by measuring the substitutional carrying of trace metals with CeCd11 precipitating from liquid cadmium. The Doerner-Hoskins logarithmic distribution law has been found to apply in all cases; the coprecipitation coefficients obtained were: 1.48 (La), 1.06 (Th), 0.63 (Pr), 0.24 (Gd), 0.17 (Sm), 0.13 (U), 0.10 (Er), 0.09 (Sr). Li, Na, K, Ba, Y, Sc, and Zr were not measurably carried. It has been observed that coprecipitation was most extensive for those tracers which form with cadmium an intermetallic compmdnd isostructural with CeCd11 or closely related thereto. For these tracers an inverse relation between Solubility in cadmium and coprecipitation coefficient was noted. THE partition of impurities between the liquid and solid phases of a metal or a congruently melting intermetallic compound has been extensively studied by the methods of zone melting.' However, when the composition of a liquid alloy and of the solid phase crystallizing from it differ grossly, partition of an impurity cannot be investigated in this manner. Partition of the impurity under equilibrium conditions may be studied by investigation of the appropriate phase diagram and partition under nonequi-librium conditions may be studied by observation of coring in the solidified alloys.2 The nonequilibrium case is of special interest in metallurgical practice because many complex alloys are cooled too rapidly to permit the uniform diffusion of impurities throughout the crystallizing phase(s). With this in mind, the present study was undertaken to obtain quantitative data on the phenomenon of coprecipitation via substitutional solid solution in certain metallic systems. It will be noticed that rather than adopting metallurgical ter- minology we have chosen to refer to the partition of impurities during alloy crystallization as coprecipitation.3 In the terminology of coprecipitation the impurity is referred to as the tracer and the precipitating phase as the carrier. From aqueous solution, coprecipitation may occur as a result of adsorption on the surface of the carrier, by occlusion of the liquid within the crystals, or by solid solution of the tracer in the carrier. By slow cooling of liquid metallic solutions it is relatively easy to form large, regular crystals of the carrier; the effects of adsorption and occlusion may be minimized and the influence of substitutional solid-solution formation may be made to predominate. Previous studies of the crystallization of some intermetallic compounds4 and of copper, silver, and gold5 have shown that the coprecipitation or partition behavior of trace metals is entirely analogous to that observed in aqueous systems when ionic substitution occurs, e.g., coprecipitation of Ra++ with BaSo4. The precipitation of CeCd11 from liquid cadmium was chosen for study. This choice was influenced by the following consideration. If the relative atomic sizes of the tracer atom and the atom for which it substitutes determine the tendency to coprecipi-tate, then a lattice having a substitutional site of well-defined size would be desirable. In CeCd116 (structure type cubic BaHg117) the cerium atom is completely surrounded by a cage of twelve nearly equidistant cadmium atoms. Thus, atoms much larger or much smaller than cerium should be poorly accommodated in this cage. The choice of a lanthanide central atom allowed the testing of the size hypothesis with a number of lanthanide tracers which differed in size but not in valence or electro-negativity. EXPERIMENTAL The general technique of this investigation was to prepare a slightly undersaturated solution of carrier in liquid cadmium to which a tracer was added. The concentration of tracer was always so low that saturation over the entire temperature range investigated was never exceeded. The solution was cooled in a stepwise fashion. After each step a por-

Jan 1, 1964

By R. M. Fowler, H. L. Hamner

DURING the last 25 years, there appeared in the literature a number of papers describing equipment and operating techniques for the determination of total oxygen in iron and steel. In the early papers, comparatively little attention was paid to side reactions that occur when iron is saturated with carbon at a high temperature in a high vacuum. As more experience was gained with the new method, the limitations appeared. Metals, with a boiling point lower than iron, distilled out of the crucible, condensed in the cooler part of the furnace, and subsequently reacted with carbon monoxide to form oxides and carbon. Aluminum, manganese, and the alkaline earth metals proved troublesome in this respect. Modifications have been devised to overcome the difficulties due to aluminum and manganese, while the alkaline earth metals are not of too much concern as yet to the steel chemist. The reactions of sulphur have received less attention, although it is generally recognized that some carbon disulphide is probably formed. Liquid air traps have been used in some cases to prevent any carbon disulphide from entering the analytical train. Sulphur was not a particular problem at the Union Carbide and Carbon Research Laboratories until the analysis of the irons of higher than usual sulphur contents became necessary. In order to describe the procedure developed for such irons, it will be necessary to first describe the type of vacuum-fusion equipment employed because the operation of vacuum-fusion equipment is an art to a considerable extent. At the Laboratories, vacuum-fusion equipment has been in daily operation since about 1930. Over the years, the evolution of this apparatus has been governed by the versatility required to analyze a wide variety of materials. Since all types of materials from high purity metals to slags are analyzed and, since seldom more than a few samples of one type are analyzed in succession, the apparatus was designed for accuracy and flexibility rather than for speed. A general view of the equipment is shown in Fig. 1, while Fig. 2 shows the relationship of the various parts. In Fig. 2, A is a vertical double-glazed 24X27/8 in. OD silica tube heated with a 20 kva mercury-arc high frequency generator. The tube is closed with a water-cooled brass head attached with de Khotinsky cement. A turret cap fits on this head and is sealed with a 1/8x3 in. "0" ring. The head contains a sight glass and two 8 in. horizontal glass tubes, for sample storage, cemented in with de Khotinsky cement. B is a four-stage Gaede mercury diffusion pump connected to the water-cooled head with a 11/2 in. copper tube through a flange with an "0" ring. The output of the vacuum pump feeds into an automatic Toepler F pump which compresses the gas sample over mercury in the measuring burettes, G and H. There are bypasses in this line for a McLeod gage C and a Piranni gage D, a manometer E, and for initial pumping out of the system with a Cenco Megavac backing-up pump. The reduction crucible shown in Fig. 3 is machined of Tungar Bulb quality (A.G.K.T.E.-2) graphite. The sample is guided into the crucible with the graphite funnel depicted on this figure. The arrangement of the crucible, radiation screens, etc. in the furnace is as follows: A 1 15/16x61/2 in. Alundum thimble, porosity RA-84, with two 1/8 in. holes near the top and opposite each other which take a pair of special tongs used for lowering it in place, is prepared. A 3/4 in. layer of powdered graphite is placed in this Alundum thimble, and the graphite crucible containing 20 g of ingot iron and 0.5 g of carbon is carefully centered on this bed of graphite powder and graphite powder is packed around it to the top of the graphite cru-

Jan 1, 1953

By C. W. Sherman, John Chipman

IN the mathematical statement of the law of mass action, the activity of each substance consumed or produced in a reaction is used to obtain a numerical constant which is characteristic of the equilibrium of that reaction at the given temperature. The activity in many cases is equal to, or proportional to, the concentration. If this is the case, any units denoting concentrations, such as weight percentage or mol fraction, may be substituted in the equilibrium expression. It is even permissible to use mixed units provided only that each unit is that used in defining the activity of the corresponding substance. The indiscriminate use of concentration units in writing equilibrium constants sometimes leads to very erroneous results. On the other hand, what appear to be contradictory results from similar investigations often might be resolved if the activity is the quantity that is used to determine the equilibrium constant for the reaction in question. De-sulphurization serves as an important illustration. Sulphur in lron and Steel The control of the sulphur content of the finished product is one of the oldest and broadest problems in the metallurgy of iron and steel. In general, this control must be exercised through the action of a slag, and for this reason experimental studies of the distribution of sulphur between slag and metal have assumed considerable importance. Under slags of the type found in the basic open-hearth furnace the ratio (% S)/[% S] does not in general exceed about 8. Under blast-furnace conditions, however, de-sulphurization ratios of 50 to 100 are readily obtained, and under laboratory conditions using similar slags, the ratio may be as high as 400. The fact that blast-furnace metal is much more readily desulphurized than is molten steel rests upon two distinct phenomena. The first of these is connected with the slag reaction. In the open-hearth furnace, the equilibrium in the following reaction places a definite limit on the sulphur content of the slag: S + 0-- = 0 + S- - where the symbols underlined indicate elements dissolved in the liquid metal, while slag constituents are represented as ions. The high oxygen content of the open-hearth bath prevents this reaction from proceeding far enough to effect good desul-phurization. In the blast-furnace process, however, the oxygen activity is limited by the presence of carbon in the metal and in the coke bed and the controlling reaction becomes: S_+O- + C_=S- + CO(gas) [2] The second phenomenon which contributes to the easier desulphurization of blast-furnace metal is the effect of high concentrations of other elements. These elements, notably carbon, silicon, and phosphorus, increase the thermodynamic activity of sulphur several fold. As a result of this increase, the escaping tendency of sulphur from the metal phase is increased and its removal facilitated. Previous Studies It becomes, therefore, a matter of importance to determine the effects of various elements upon the activity of sulphur in the melt. The experimental approach to this problem has been made through the equilibrium between sulphur in the melt and gaseous atmospheres containing hydrogen and hydrogen sulphide. The reaction and its equilibrium constant are written as follows: S + H2(gas) =H2S(gas);K1 = -- [3] Pn2 • as where the symbol a, represents the activity of the sulphur in solution in the metal. The understanding of the chemical behavior of sulphur in steel and pig iron rests on a firm knowledge of the binary system Fe-S. A number of studies of the equilibrium have given closely agreeing results from which the activity and free energy of the dissolved sulphur may be found. The recent paper by Sherman, Elvander, and Chipman1 includes a discussion and quantitative comparison of the results of several of the significant researchesr-" on dilute solutions of sulphur in iron. The results of this study indicate that the relationship between the ratio Ph,s/p+ and the percentage of sulphur in solution is not a linear one. This

Jan 1, 1953

Jan 1, 1920