Search Documents

By Paul A. Beck, John Towers, W. D. Manly

Recent work on grain growth in high purity aluminum and in a solid solution type alloy of aluminum and magnesium1 showed that the isothermal increase of the average grain diameter D with time follows the equation D = K(tg + A)2 Eq 1 Here tg is that portion of the total annealing time 1, which is available for grain growth, that is the annealing time less the time necessary for complete recrystallization. The parameters A and n depend only on the temperature. A has the dimension of time, and is usually of the same order of magnitude as the time for recrystallization: R. In some instances R may be substituted for A in Eq I. This is the case for example at long periods of annealing or at high temperatures, where both A and R become small in comparison with t. In such instances Eq I is simplified to D = K, tn Eq 2 For high purity aluminum the exponent n increases in the ratio of about I to 5 in. the temperature range of 350' to 6w°C. 70-30 brass has been long known to exhibit grain growth of the continuous type, such as occurs in pure metals. For this reason, it was of interest to determine whether grain growth in brass also follows Eq I and 2 as it does in aluminum. Isotherma1 grain growth data for 70-30 brass were published in recent years by R. S. French2 and by H. L. Walker.' Fig I shows a log D vs. log t plot of French's data, together with the 500, 600 and 700°C data by Walker for 42.8 pct reduction by rolling. In such a log D vs. log t plot Eq 2 corresponds to a straight line with a slope of n. As seen in Fig I, the data of French and the 600°C data of Walker can be represented by straight lines. However, the 500 and 7c0°C data of Walker give curved lines and the average slope of his 500°C line is very different from the slope of the line representing French's data for the same temperature. Furthermore the curvature of the 700°C line suggests a tendency for grain growth to stop after 8 hr. If confirmed, this would constitute a definite deviation from Eq 2, not explainable on the basis of the specimen thickness effect found in high purity aluminum.' Walker's specimens were 0.125 in., or 3.18 mm thick, while the deviation occurs at a grain size of approximately 0.4 mm. A further point of divergence arose recently, when J. E. Burke suggested4 on the basis of Walker's data, that, instead of Eq 2, grain growth in 70-30 brass follows the equation D - D, = K . tgn Eq 3 Here D, is the grain size as recrystallized at the temperature in question. Although it was possible to show5 that, for pure aluminum, for which isothermal grain growth has been carefully investigated, Eq 3 is incompatible with the experimental facts, the situation remained somewhat obscure with respect to 70-30 brass. In view of the discrepancies between the published data on brass from different

Jan 1, 1949

By John Towers, Paul A. Beck, W. D. Manly

Recent work on grain growth in high purity aluminum and in a solid solution type alloy of aluminum and magnesium1 showed that the isothermal increase of the average grain diameter D with time follows the equation D = K(tg + A)2 Eq 1 Here tg is that portion of the total annealing time 1, which is available for grain growth, that is the annealing time less the time necessary for complete recrystallization. The parameters A and n depend only on the temperature. A has the dimension of time, and is usually of the same order of magnitude as the time for recrystallization: R. In some instances R may be substituted for A in Eq I. This is the case for example at long periods of annealing or at high temperatures, where both A and R become small in comparison with t. In such instances Eq I is simplified to D = K, tn Eq 2 For high purity aluminum the exponent n increases in the ratio of about I to 5 in. the temperature range of 350' to 6w°C. 70-30 brass has been long known to exhibit grain growth of the continuous type, such as occurs in pure metals. For this reason, it was of interest to determine whether grain growth in brass also follows Eq I and 2 as it does in aluminum. Isotherma1 grain growth data for 70-30 brass were published in recent years by R. S. French2 and by H. L. Walker.' Fig I shows a log D vs. log t plot of French's data, together with the 500, 600 and 700°C data by Walker for 42.8 pct reduction by rolling. In such a log D vs. log t plot Eq 2 corresponds to a straight line with a slope of n. As seen in Fig I, the data of French and the 600°C data of Walker can be represented by straight lines. However, the 500 and 7c0°C data of Walker give curved lines and the average slope of his 500°C line is very different from the slope of the line representing French's data for the same temperature. Furthermore the curvature of the 700°C line suggests a tendency for grain growth to stop after 8 hr. If confirmed, this would constitute a definite deviation from Eq 2, not explainable on the basis of the specimen thickness effect found in high purity aluminum.' Walker's specimens were 0.125 in., or 3.18 mm thick, while the deviation occurs at a grain size of approximately 0.4 mm. A further point of divergence arose recently, when J. E. Burke suggested4 on the basis of Walker's data, that, instead of Eq 2, grain growth in 70-30 brass follows the equation D - D, = K . tgn Eq 3 Here D, is the grain size as recrystallized at the temperature in question. Although it was possible to show5 that, for pure aluminum, for which isothermal grain growth has been carefully investigated, Eq 3 is incompatible with the experimental facts, the situation remained somewhat obscure with respect to 70-30 brass. In view of the discrepancies between the published data on brass from different

Jan 1, 1949

By L. K. Nemeth, H. T. Kennedy

The investigation presented in this paper was performd to develop a relationship between the dewpoint pressure of a hydrocarbon reservoir fluid and its composition, temperature and characteristics of the heptanes-plus fraction. This is part of an over-all program to develop correlations for the essential properties of reservoir fluids with the composition, temperature and pressure. In all, 579 dewpoint pressures, determined experimentally on 480 different condensate systems along with their compositions and temperatures, were obtained from the literature and from previously unpublished sources. Compositions were expressed in terms of mol fractions of methane through heptanes-plus, nitrogen, carbon dioxide and hydrogen sulfide, and the molecular weight and specific gravity of the heptanes-plus fraction. Distribution of the samples assured a good representation of most of the major oil-producing areas in the world. A multiple regression analysis program written for the IBhl 7094 digital computer was used to compare numerous equations and to calculate their coefficients. The final equation containing 11 coefficients resulted in an average absolute deviation of 7.4 percent for all the samples. Segregation of data on the basis of one or more variables did not materially improve the results. The correlation developed is useful in predicting the dewpoint pressure of condensate reservoir fluids from variables normally measured in the laboratory. INTRODUCTION In any discussion of the dewpoint of hydrocarbon mixtures, it should be noted that there are two kinds of dewpoints which must be distinguished. The first type, or normal dewpoint, usually occurs at low pressures as dry gas is compressed to the point where liquid first forms. The second type, sometimes called the retrograde dewpoint, occurs when a gas mixture containing heavy hydrocarbons in solution is depressured until a liquid is formed. The first type, as applied to condensate mixtures, would normally be below atmospheric pressure and would be of no interest in reservoir performance. The second type is of great significance in interpreting reservoir behavior and is the property referred to throughout this paper. Laboratory measurement of dewpoints is laborious, expensive and subject to many errors, and it is desirable to find a reasonably accurate method to predict dewpoint pressures from the composition and other readily available properties normally measured on the fluids. A possible method of calculating dewpoint pressures would be to employ a set of equilibrium ratios (K-values). 1 However, the K-values themselves depend on composition as well as temperature and pressure; presently their use adds nothing to the accuracy of predicting dewpoints. If future correlations provide a significantly better method of estimating K-values for heavy hydrocarbons, they may be of greater use. REVIEW OF LITERATURE In 1944, Olds, Sage and Lacey 2 studied the behavior of fluids from the Paloma field. Properties of seven mixtures were studied at temperatures ranging from 100 to 250 F and at pressures up to 5,000 psia. From the information obtained, the influence of composition and temperature upon the dewpoint was determined. It was found that removal of the intermediate molecular weight components from the mixture resulted in considerable increases in the dewpoint pressures. Results indicated that the effect of changes in temperature was relatively minor when compared with the effect of modifying the mixture by removing the intermediate components. Thus, removal of propane and butanes resulted in an increase in dewpoint pressure more than 10 times greater than a change in temperature from 100 to 250 F. The authors stated that althouah

By D. G. Westlake

NUMEROUS hydride habit planes in zirconium, Zircaloy-2, and Zircaloy-4 have been observed by various researchers and have been compiled by Louthan and Angerman.1 While only partial agreement exists on the identification of these habit planes, there is a growing consensus that a tensile stress applied on a polycrystalline rod during slow cooling will cause the hydride platelets to lie perpendicular to the tensile axis. walters2 agrees with Louthan and Angerman who reason that "The particular plane on which precipitation will occur is the plane that will result in the maximum relief of stress by precipitation of the low-density hydride." The experiment described below was initiated to test this hypothesis on a single crystal of zirconium. A single-crystalline tensile specimen of zone-refined zirconium, designated number 13 in a previous study,3 had a square cross section with one pair of parallel faces only 4 deg from the basal plane of the hcp structure. The tensile axis was about 86 deg from the (0001) pole, 14 deg from the (1010) pole, and 16 deg from the (1150) pole. The specimen was charged to 50 ppm H by allowing reaction with the correct amount of hydrogen at 1073°K, homogenizing at this temperature for 2 hr, and cooling in the furnace. Two- surface analysis revealed that hydride platelets had precipitated parallel to the three {1010} prism planes. The specimen was heated to 673°K, in air, for 16 hr to bring the hydrogen back into solid solution, and furnace-cooled to 573°K.* It was then stressed in ten- sion at 1.2 kg per sq mm. (This produced a resolved shear stress of 0.6 kg per sq mm for glide on the prism plane with the highest Schmid factor. A higher stress would have produced significant strain at this temperature because the critical resolved shear stress for prism-plane glide in this crystal was only 1.0 kg per sq mm at ambient temperature.) The poles of the three prism planes were at angles of 14, 46, and 74 deg to the tensile axis so that an applied stress would exert normal stresses on the three planes in the approximate ratios 10:5:1, respectively. As soon as the stress was applied, the sample was furnace-cooled to room temperature. It was ground on wet Sic paper and chemically polished in 50 parts HNO3, 50 parts water, and 7 parts HF before examination. The result of the experiment was negative; i.e., the three (1010) planes were of equal importance as hydride habit planes.* It must be admitted that the ex- periment is not entirely comparable to those performed by walters,2 Louthan,5 and Louthan and Marshall6 on polycrystalline materials. They have found that some critical stress (below the macroscopic yield stress) is required for stress-orienting the hydrides. The applied stress in the present experiment was, by necessity, smaller by an order of magnitude than the critical stresses they report for the various materials. If the critical stress was related to the microstructure, one might postulate that it should be decreased by increasing the degree of perfection of the metal lattice and by increasing the grain size. But the experiment above suggests that, for the extreme case of a zone-refined single crystal of zirconium, the critical stress lies above that which one can apply without introducing gross plastic strain as another variable.

Jan 1, 1967

By R. E. Collins

The existence of fluid migration across fixed boundaries in oil and gas reservoirs has been known for many years. Several techniques have been developed in the past for estimating The rate of migration across fixed boundaries as an aid in planninng field development and in the valuation of oilfied properties. The principal de/errc7nt to the use of these techniques lies in the rather extensive reservoir and field data required for ca1cula1ion.s of fluid migration. For this reason, a new, simplified procedure has been developed which makes possible the calculation of fluid migration with a minimum of field and resrvoir data. This new technique is based on certain. solutions of the differential equations describing flour in the reservoir which assume, for the portion of the reservoir of interest, that the formation can be approximated by a homogeneous rock of uniforrn thickness, that only a Single mobile fluid phase exists, and {hut fluid production at the well is solely a result of expun.rioti of rile resersoir fluds. The results of the, present work are compiled in a set of curves. These curves can he used to calculate both rule and cummulative fluid migration when the abovc assumptions ore ju.r/ifiril. The only data required for such calcu1ations are the production hisrories of all wells in the field, the permeability and porosity of the reservoir., the, compressibility and viscosity of the fluid. and the shape of the reservoir, INTRODUCTION The migration of fluids across fixed boundaries in oil and gas reservoirs has long been recognized as an engineering and economic problem, but no methods for estimating the extent of migration during the development stage of a reservoir have been published. The investigation reported here has been directed toward providing a simple technique for estimating the extent of such migration and has therefore been restricted to the simple case of a homogeneous reservoir of uniform thickncss containing compressible liquids. It is assumed that only a single mobile fluid phase exists and that fluid production at the wells is solely by expansion of the reservoir fluids. Obviously, few, if any. reservoirs conform to these assumptions during their entire productive life. On the other hand, most reservoirs approximate fluid expansion reservoirs during their initial stages of primary production. Consequently, the present work should be applicable during the field development of most reservoirs, regardless of whether or not they ultimately are gas cap, dissolved-gas, or water-drive reservoirs. The mathematical analysis of this problem has been outlined in the Appendix, while the results of the investigation are presented in the following discussion. DISCUSSION OF CALCULATION PROCEDURE SLNGLE Well, Ini-iimite Reservoir The procedure for calculating fluid migrations is based on solutions of the partial differentia1 equations describing the flow of a homogeneous compressible liquid in a homogeneous reservoir of uniform thickness. This differential equation has the form of the well-known diffusion equation. The equation, as applicable to our problem, is where p is fluid density, is porosity, EL is viscosity, c is compressibility,* and k is permeability; .r, y and t are the space coordinates and time. Derivation of this equation is given in the Appendix. The analytical solution of this equation for a single well of negligible radius producing at a constant rate. q, from a reservoir of infinite areal extent for a time, 1. can he readily obtained. From this solution the rate of flow and also the cumulative flow across a straight line segment of length, 1, at a perpendicular distance. 6. from the well can be computed. One end of the line must be on the perpendicular joining the well to the

By A. A. Bauer, M. S. Farkas, F. A. Rough

An investigation of the d-phase relationships between the uranium-zirconium and uranium-molybdenum systems was conducted. A ternary cut joining U-31.5 at. pct Mo to U-74 at. pct Zr is presented on the basis of differential thermal analysis, metallographic and X-ray diffraction data. HyDothetical isothermal ternary sections, inferred from the determined d-phase relationships, are presented for the uraniztm-zirconium-molybdenum system. A ternary eutectoid, at 454ºC, was found in zirconium-rich alloys. Its composition, based on metallography and electron-beam micro-analysis, is estimated as U-61.7 at. pct Zr-7.1 at. pct Mo. URANIUM binary systems have been studied extensively in the past to supply basic information needed in reactor-fuel development. Alloys of uranium with zirconium, molybdenum, and titanium are of particular interest because, in these binaries, extensive or complete solid solubility between the body-centered-cubic forms of these elements and y uranium exists. In addition, an intermediate phase of limited solid solubility forms in each of these systems by transformation of the high-temperature ?-uranium phase. Previous investigation of the phase relationships existing between the intermediate 6 phases in the uranium-zirconium-titanium and the uranium-molyb- denum-titanium ternaries have already been reported.1,2 The results of an investigation of the phase relationships existing between the 6 phases in the uranium-zirconium-molybdenum system is reported here. The phase equilibria that occur in the solid-state portions of the pertinent binary systems are shown in Fig. 1. Phase nomenclature appearing throughout the report is as shown in these diagrams. The crystal structure of the phases in this system are given in Table I. EXPERIMENTAL PROCEDURES Alloys were cast and fabricated to 1/4-in.-diam rod. Specimens were heat treated and examined by metallographic and X-ray techniques. Thermal-analysis data were obtained.

Jan 1, 1960

By J. Kucera, J. Vreštál

DIFFERENT authors1 ' have been engaged in meas-uring the vapor pressure of pure cobalt. The results of their measurements satisfy the expected temperature dependence of vapor pressure and are in mutual agreement, with the exception of the results by De Dyk Man and ~esmejanov.~ The latter authors, by using the effusion method, have achieved considerably higher values than other workers. By using the isotope exchange method they obtained much lower values than other authors. The determination of thermodynamic functions of the binary Ni-Co system was carried out, as far as we know, only by De Dyk Man and Nesmejanov.4 This determination was based on extremely low values of vapor pressure over pure cobalt obtained by the isotope exchange method. Therefore we considered it useful to repeat these measurements. Furthermore, the reliable values of thermodynamic functions should complement the migration characteristics, the inner magnetic field data obtained by Miissbauer effect measurements, specific heat data, and the K effect and positron annihilation characteristics that have been achieved for the Ni-Co system in the Institute of Metallurgy of the Czechoslovak Academy of Sciences in Brno.5-9 The measurements in this paper were performed by the isotope exchange method. In this method, as has been shown by analysis,10-'3 the vaporization rate, w, can be calculated from the relation w = k/1oSw [1] i0Sw) [1] where L, is the activity per unit mass of the source, S is the area of evaporation, k is the tangent at the point t = 0 to the curve giving the radioactivity de-pendence of the target vs annealing time, t, and wis the probability for the tracer atom condensation on the target surface. The value for cc. can be obtained from the equation w = a + 2P(l-a) [2] where P is Clausing's coefficient and a is the condensation coefficient. In the present work the cylinder shaped containers for annealing the samples were used, characterized by the ratio 1:d = 1.06. Here 1 is the length and d the diameter of the container. For this ratio P = 0.514 and w =0.5, and w is independent of the condensation coefficient. The vaporization rate is related to the vapor pressure p by the relation where R is the gas constant, T is the absolute temperature, and M is the molecular weight of the component investigated. Measurements of the vapor pressure of the i component over the pure metal, p!, and the alloy, pi, at the given temperature and composition enable one to calculate the thermodynamic activity of this component from the equation ai=pi/pi [4] is and the activity coefficient ri=ai/xi [5] xi where xi is the molar fraction of the investigated alloy component. EXPERIMENTAL Evaporization rate measurements were carried out with pure cobalt and eight different Ni-Co alloys, which were kindly delivered by the Centre National de Recherches Me'tallurgiques Liege, Belgium. The composition of the alloys is given in Table I. The cobalt used in these alloys contained the following impurities: 0.016 wt pct Fe, 0.014 wt pct Cu, 0.01 wt pct Mn. Other elements were present in much lower concentrations. From the melts, samples in the form of small cylinders of 12 mm in diam were ma-chined. The thickness of the targets and sources was 2.5 and 1.2 mm, respectively. The lower thickness of the sources was chosen to make their over-all ac-tivity as low as possible. Sources were activated in

Jan 1, 1970

By Charles Stein

The effects of several different prior treatments an the creep behavior of 99.9995 pct aluminum at 260°C and 1000 psi canstant stress are compared with annealed specimens. Quenching from 538oC, irradiation with 2 mev electrons, and tension-compression room temperature fatigue damage were used to change the substructure of the specimens prior to their creep testing. The quenched and the irradiated specimens showed a larger primary and transient stage contribution to the creep curve than those specimens which were only annealed prim to creep testing. The specimen receiving prior fatigue damage at room temperature showed no first stage or transient creep when tested under identical conditions as the above specimens and had an average creep rate seven orders of magnitude lower than that of either the annealed, the quenched or the irradiated creep specimens. In a previous electron transmission microscopy investigationl of the substructure developed during the creep of pure aluminum at 0.57 Tm, it was noted that a large concentration of vacancy loops and dislocation loops were present in and near the subboundaries while the interior of the subgrains contained few of these defects and had a low dislocation density, see Fig. 1. Yim and Grant2 and Hazlett3 have shown that the presence in nickel of a substructure developed by prior cold work effectively reduces the primary and transient contributions to the creep curve. Hultgren,4 McLean,5 Gervais, Norton, and Grant,6 Chang and Grant,7 and others have shown the same effect with the higher stacking fault energy material, aluminum. However, the specific defect responsible for the change in the creep behavior developed by the prior cold work was not established. Specifically, what effectively interfered with dislocation motion in these materials at elevated temperatures? Was it jogs on dislocations produced by their interaction with other dislocations or with vacancy or dipole loops, see Fig. 1, or did the subgrain walls determine the mean free path for glissile dislocations? This would be a realistic possibility only if subgrain boundaries are effective barriers to dislocation motion. The ability of a subgrain boundary to act effectively against glissile dislocations depends on the number and the arrangement of dislocations in the subboundary,8 which in turn is a function of the creep strain.' This paper compares the creep rate of specimens possessing a large vacancy concentration to that of annealed specimens and with specimens having a stable subboundary wall containing Lomer locks produced by prior fatigue damage. EXPERIMENTAL PROCEDURE Vacancy Loops. Aluminum creep specimens, 99.9995 pct pure, were machined from zone refined rod to a 2-in. gage length and a cross-sectional area of 0.025 in. These specimens were annealed at 538oC, 1000°F, for 10 min and quenched into ice water. They were then reheated to the creep test temperature of 260°C, 500°F, and held at this temperature for 18 min to al-low for vacancy condensation and loop formation. The specimens were subsequently creep tested at 1000 psi constant stress at 260°C. The temperature along the gage length was monitored by three Pt-Pt, 13 pct Rh thermocouples, two embedded in the upper and lower shoulders of the specimen and one attached to the mid-point of the gage length- The temperature gradient was held to ± l°C or less throughout the Creep test. Elongations were measured by two LVDT's 'connected in an averaging 'On-figuration, with quartz extension arms placed against the specimen at the extremes of the 2-in. gage length. i K 3J?--' * *H - >^> - ' ^L r 77 . *y J OJL, J Fig. 1—Transmission electron micrograph of pure aluminum strained 8.2 pct at 260°C and 1000 psi constant stress. Note the presence of numerous loops near the light dislocation tangle

Jan 1, 1970

By J. Buchwald, R. W. Heckel

THE DeHoff method' of analysis of ellipsoids of constant shape provides a quantitative metallo-graphic technique whereby the actual size distribution of the ellipsoids may be determined from polished and etched random planar sections. This method provides a useful metallographic tool since it appears to be the only method available for determining the size distribution of nonequiaxed second-phase particles. This method, however, assumes that the type of particle being studied is known to be either oblate or prolate of a fixed minor to major axis ratio, q. This, in general, means that one assumes either oblate or prolate ellipsoids of a given q and checks the total volume and surface area of the second phase as determined from the particle size distribution to values of these quantities determined by other techniques. Suitable check methods may be point counting, lineal fraction, and areal fraction methods for obtaining the total volume and the Smith and Guttman2 method of intercepts for obtaining the total surface area. If the initial guess at particle shape and q were incorrect, i.c., the volume and surface information obtained from the particle size distribution data does not agree with those values obtained by other techniques, another assumption must be made and the particle size dist distribution must be recalculated. A current study in this laboratory has involved a large number of these calculations. To facilitate the handling of the calculations, the DeHoff method was programmed in the FORTRAN language for an IBM 1620 computer. The purpose of this note is to briefly describe this program. The complete program along with a list of variable names, a sample input, and a sample output may be obtained from the senior author. Basically, the program provides for the calculation of particle size distributions for a series of q values for both oblate and prolate ellipsoids for each set of input data. The total volume of second-phase particles per unit volume are calculated for each calculated particle size distribution. The output of the program is arranged to allow a rapid comparison of the calculated volume and surface values to data obtained by independent check measurements. A condensed flow chart of the program in the notation of DeHoff1 is given in Fig. 1 The calculation of the number of particles per unit volume in the jth size class, Ni, according to the DeHoff method' is: where k(q) is the shape-correction factor given by DeHoff as a function of q for both oblate and prolate ellipsoids, A is the size-class increment, k is the number of size classes, ni is the number of particles in the ith class on the planar measurement surface, and B(j,i) are the saltykov3 coefficients which correct for measurement surfaces which do not pass through the major and minor axes of the particles. The values of /3(j,i) are included as a two-dimensional array in the program. Particular values of k(q) are included in the program in increments of q of 0.05. The first calculation is carried out for q = 1 (spheres), the second for q = 0.95 oblate, the third for q = 0.95 prolate, and so forth. The number of q values chosen for the calculation is determined by a control card in the input for each set of data.

Jan 1, 1965

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 H. I. Aaronson, H. A. Domian, G. M. Pound

The thernodyna,nics of I the Proeutectoid ferrite re-action ha1.e been investigated on the bases of three diifevent descviptions of the statistical thernzodynamics of interstitzal solid solutions. Especially at low ternperatures, substantial dijeretzces are found in the valltes of the extrapolated y/y + a and cr/a + y equilibriurn curls, and in the To- co7nposilion metastable eqliilibriunz curl-es thus co,npited. Although neither actiztity nor phase-diagram date are able to distingrish udeqiutely among the various descviptions, the statustics of Lacher, Fouller, and Guggenheinz provide the ?nost conplete epresentation az3ailable of the single carbon-ctrbon interaction energy (w) nzodel of an intevstitiril solid solution, and the relationships deduced frojrl then2 are accordingly pre.ferred. The findings tlzut w in austenite 'aries significu?ztly uith tenpevntrue and that phase-diagawr data arc not accurately con?palible usith eqcrctions deduced on the basis of actility infoation indicate, hou'e13er, that a wzore cornplex nlodel of the uustenite , und pvobably of the .feyrite, solirl solutions will have to forn the basis of f4tuve statistical t1.eat1nents. Some suggestions for such a ruodel arc advanced. Althoug11 w represents a vepltlsio energy in austozite, it was found to corvcsporzd in ferrite to a binding energy, in qualitatit.e accord zc,itlz the internal-friclion results of Keefer and wert. ThE transformation of austenite to proeutectoid fer-rite is becoming recognized as an exceptionally useful vehicle for quantitative studies of the kinetics of phase transformations.'" Reasons for this popularity include well-characterized microstructures which develop on a convenient scale with acceptable kinetics within a considerable temperature-composition region, the absence of great sensitivity to common impurities, and, especially. the availability of a unique wealth of the ancillary information indispensable for the quantitative interpretation of kinetic measurements, including phase diagram,3 crystallographies diffusivity,5 inter-facial energy,=17 elastiit," specific vlume, ''' and ativity'-' data. Development of considerable amounts of additional thermodynamic information from activity data is of particular importance with respect to kinetic studies, since much of this information is needed in temperature-composition regions not directly accessible to presently available techniques for the measurement of equilibrium properties. A generally useful approach to the thermodynamics of solid solutions, which is of especial value in this situation, is based upon the formulation of the free energy of such solutions as a function of temperature and composition. The free energy of solid solutions has been resolved into many components. Particularly in the case of interstitial solutions, with which we shall be concerned in the present study, however, only the positional entropy component has been intensively studied. The other components are conventionally taken into account by fitting equations for the activity of the interstitial species or the temperature-composition path of phase boundaries, derived principally from positional entropy considerations, to experimental information on these quantities through adjustment of one or more constants in the equations. A number of studies of the thermodynamics of interstitial solid solutions have been made primarily on this basis,'"25 but only Kaufman, Radcliffe, and cohenZ2 (hereafter KRC) and ener''% ave applied the results of such studies to develop relationships for the thermodynamic quantities of interest in the austenite —- proeutectoid ferrite transformation. KRC made use of a model developed by cheil'' (on the basis of a previous study by ohansson'), and subsequently rede-rived correctly by Speiser and retnak,' in which a carbon atom is considered to inhibit occupancy of a certain number of its nearest-neighboring interstitial sites. schei120 determined the number of sites in austenite thus excluded by matching calculated and experimentally measured curves for the variation of carbon activity with the mole fraction of carbon. zener17 used a different approach. He assumed that the positional entropy of carbon is ideal both in austenite and in ferrite (i.e., no excluded interstitial sites), and included all other contributions of carbon in a temperature-independent "free-energy change" accompanying the transfer of 1 mole of carbon between the two phases. This quantity was determined by "pinning" an equation derived for the boundaries of the a + y region to experimental data on the compositions of these boundaries at a single temperature. In the present investigation, three different approaches to the thermodynamics of the proeutectoid ferrite reaction are examined comparatively: I) the method of Zener; 2) the model employed by KRC and predecessors; and 3) the more detailed statistical thermodynamic treatment of the same model developed by acher and Fowler and uenheim, ' in which the problem of overlapping regions of influence of individual interstitial atoms is resolved. In order to make these comparisons, it is necessary to undertake an extension of one aspect of Zener's method; since the Lacher- Fowler -Guggenheim (LFG) treatment has not been previously used in the present context, the necessary thermodynamic "superstructure" for this treatment will be developed in parallel

Jan 1, 1967

By D. J. DeLazaro, W. Rostoker, R. E. Riley, M. Hansen

At very low concentrations, carbon dissolves interstitially in molybdenum resulting in a linear expansion of lattice parameter with increase of carbon in solid solution. Geometrical consideration of the relative size of carbon atom to size of interstice approximately predicts the observed volume expansion. THE element carbon, occurring in molybdenum, either as an unavoidable impurity, or purposely added, markedly affects the mechanical properties of molybdenum. Several investigators1-" have noted the occurrence of molybdenum carbide at the grain boundaries of cast molybdenum, and Fischer and Rengstorff' have shown that this intergranular constituent can induce intergranular brittleness. The appreciable change, with temperature, of the solid solubility of carbon in molybdenum, found by Few and Manning," suggests other important effects upon mechanical properties. Finally the relative atom diameters of carbon and molybdenum (the ratio of atom diameters is 0.59) are indicative, by analogy with iron (ratio of carbon atom diameter to iron atom diameter is 0.65), of a possible interstitial solid solution. If carbon dissolves intersti-tially in molybdenum, anelastic effects and age-hardening effects, similar to those observed in other body-centered cubic metals, could be expected. However, from the work of Sykes, Van Horn, and Tucker," it might be inferred that carbon forms only a substitutional solid solution with molybdenum. These considerations made it desirable to study the effect of carbon on the lattice parameter of molybdenum in greater detail than heretofore. All heat treating was done in a high temperature furnace7 built specifically for heat treatment of molybdenum. The temperatures and times of heat treatment were determined from the a solubility limit for the Mo-C system." These heat treatments were carried out in argon which was purified with respect to N, and H.!0." Purification of the argon was accomplished by passing the argon over hot (750°C) titanium and then through a magnesium perchlorate drying tower. As many as six samples were run concurrently and quenched at different times without interrupting the furnace cycle. Both the purity of the furnace atmosphere and ability to quench samples from any temperature up to 4000 °F were of prime importance in the successful heat treatment of samples for the lattice parameter investigation. Carbon Analysis The Leco Combustion apparatus was used for carbon analysis. This apparatus is capable of an accuracy of & 0.003 pct for a 1 g sample in the range below 0.10 pct. However, it is felt the degree of accuracy is somewhat better than this limit in view of the consistency evident by cross checking a number of determinations. The Battelle Analytical Laboratory consistently reported check results to within k 0.001 pct C. X-ray Measurements The molybdenum wire specimens (20 mils) were carefully filed down to about 12 mils and then etched with nitric acid to 8 mils. Carbon analyses were made on specimens before reduction in size; however, since there was no evidence of surface carburization, it is probable that the carbon was uniformly distributed throughout the sample. The resultant specimens were strain free and small enough to minimize the error in the lattice parameter determinations due to X-ray absorption. Diffraction patterns were obtained with a Debye-Scherrer powder camera 114.59 mm in diam utilizing the Straumanis type asymmetrical film mounting. Nickel-filtered copper. Ka radiation and zirconium-filtered molybdenum Ka radiation were used. Line positions on the film were measured with a traveling microscope equipped with a vernier scale permitting measurement to 0.001 cm. The average of five readings was used to determine the position of each line. The lattice parameter for each specimen was determined from the data by the extrapolation method of Bradley and Jay, Y.e., the intercept at 90" of the straight line obtained by plotting the apparent lat-

Jan 1, 1953

By Donald C. Barton

The theory of gravitation is based on Newton's law that any two bodies exert a mutual attraction which is proportional to the product of their masses and inversely proportional to the square of the distance between them; i. e., A = KMRm2 (1) where K is a constant, M and m are the respective masses of the two bodies, and R the distance between them. In the earth's gravitational system, the earth is one of the bodies; the other may be any body within or without the earth. In theoretical discussions and calculations,in regard to gravity, that other body usually is taken as a unit of mass: in the C:. C:. S. system, a body of one gram in mass. Gravity and Level Surfaces Gravity is defined in geophysics as the force of attraction exerted by the earth, towards itself, on a body of unit mass. The intensity of gravity is usually expressed in dynes of force or in centimeters of acceleration In geophysical work in America, most commonly it is expressed in dynes (per gram) with the "per gram" understood and not directly expressed. In geophysical work in Europe, most commonly it is expressed in centimeters per second per second (cm. /sec.2) and less commonly lrierely as centimeter-gram-second units (c. g. s). At sea level, the intensity of gravity amounts very nearly to 980 dynes. The vertical is the direction along which the attraction of gravity is exerted at any point. In most cases, it is very nearly, but not quite, perpendicular to the earth's surface. The path that would be taken by a body falling freely through space toward the earth is vertical at each point and may be spoken of also as ('the vertical," or as the line of the vertical, or as a line of force of gravity. At the earth's surface it comrrlonly is very nearly, but not quite, a straight line.

Jan 1, 1929

By A. E. Wheeler, H. Y. Eagle

The copper industry in the Province of Katanga, in the Belgian Congo, which is now controlled and operated by the Union Minére du Haut Katanga, had its inception many years ago in the vision of a Scotsman, Robert Williams, now Sir Robert Williams, Baronet. He it was who personally negotiated the mineral concession with King Leopold of Belgium, sent in the prospecting parties who located most of the present mines of the company, later arranged for the Belgian participation from which the present company was formed, and has been its active mentor ever since. The idea back of his efforts was the establishment of an industry on the Congo-Zambezi divide, of sufficient importance to draw up to it the head of the Cape to Cairo Railroad, the dream of Sir Cecil Rhodes, and the consummation of which he had left as a legacy to Mr. Williams, his successor. The reports of the prospecting expeditions showed that copper in profitable quantity was there, but many difficulties lay in the way of the establishment of the infant industry. The isolated location of the mines, various tribal wars, the Boer War, and many other factors contributed to making progress very slow. However, in spite of these difficulties the railroad was finally built and a blast-furnace plant erected at Lubumbashi, near Elisabethville, the capital of the Province. Operations of this plant were started in 1911, treating high-grade oxidized ores and producing pig copper by a reducing smelting. Production has been continuous since that time and the company has grown so that at present it ranks as the largest individual copper producer in the industry. It is a fitting monument to the foresight, optimism and driving force of its founder, and the progressive policies of his Belgian colleagues. In January, 1914, the Union Miniere du Haut Katanga engaged as its consulting engineer, Archer E. Wheeler, metallurgical engineer, at that time general superintendent of the Great Falls plant of the Anaconda Copper Mining Co. Mr. Wheeler was commissioned by the Board of Directors to proceed into Katanga, and make all necessary investigations to enable him to report on a comprehensive development program, and to recommend a scheme of treatment particularly for their low grade

Jan 1, 1934

By AIME AIME

HENRY DAYTON WILDE-"Date," informally- manager- of the technical and research department of the Humble Oil & Refining Co., was born at Aguascalientes, Mexico, Sept. 25, 1900. He is a graduate in chemical engineering: M.S. from the University of Texas and I Sc.D. fro111 the Massachusetts Institute of Technology. t For a year and a half in 1926 and 1927 he was a research : associate in the research laboratory of applied chemistry at M.I.T., under W. K. Lewis. In the latter year Dr. Wilde went with the Humble as head of the laboratory experimental division of its refining development department. In 1929 he was made director of production research. Dr. Wilde is the author of important papers in the TRANSACTIONS of the Institute, wherein are also mans dis-

Jan 1, 1938

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 H. I. Aaronson, H. A. Domian

ELectvon-probe analgsis has been used to in7:estigatr the partition ot alloying- elements between aus-tenite and proeutectoid ferrite (or the Jerritic component of bainite) in high-purity Fe-C-X alloys, lsuully at tzco levels of carbon content fca 0.1 and 0.4 pct in ntost steels) at temperatures between the M, and the Ae3 in Si, mn, Ni, mo, Co, Al, Cr, Cu, and PI steels. No partition ous found during the early stages of transformation at any temperature studied in the Si, ,mo. Co, Al, Cr, and Cu steels. Partition did take place, Jrowel'er, in Ihe Mn, ATi, and Pt steels, occurring in lewtperature. regions below the Ae3 of ca 130°, 65", UNDERSTANDING the mechanisms through which alloying elements influence the kinetics of austenite decomposition reactions which proceed by nucleation and 35c respectively. A break was discovered in the TTT curve for the initiation oj the proeutectoid ferrite reaction in the latter steels at the temperature at which partition began. The alloy content in ferrite at the a:y bondavy remained at constart lecels, shown experimentally to be significantly above those corresponding to eqiilib+iu??l, e'en at the longest reaction times investigated (ca 1o6sec). The partition and the TTT curve observations foere rationalized on the basis of the thermodynamic treatment developed in an accompany ing paper. and diffusional growth has been a major research objective of physical metallurgy since the latter part of the 19th century. The measurements made on the dif-fusivities of carbon and of substitutional alloying elements in austenite during the 1930's and 1940's seemed to provide a useful clue as to the nature of these mechanisms. These measurements showed that the diffu-sivity of carbon is several orders of magnitude higher than that of the alloying elements at a given tempera-

Jan 1, 1967

By Milton E. Wadsworth, R. Ted Wimber

Cobalt sulfate solutions containing ammonium acetate and chloroplatinic acid were reduced by hydrogen in a pyrex-glass lined autoclave in the temperature range of 170o to 232°C and hydrogen partial pressure range of 115 to 830 psia. The reduction rate was directly proportional to the hydrogen partial pressure and surface area of the pyrex glass and was independent of the quantity of chloroplatinic acid added initially. Experiments involving the variation of the relative concentration of ammonium acetate indicated that the reducible cobalt complex was probably the diacetate complex of cobalt, Co(AC)4H20, or a new mononcetate complex Co Ac, which was in solubility equilibrium with a pink precipitate of CO(AC)-4HzO. THE reaction in which a metal is dissolved by an acid to produce gaseous hydrogen and a salt solution was discovered early in the history of chemistry. In 1859 Beketoff found experimentally that this reaction could be reversed; i.e., a salt solution could be reduced by gaseous hydrogen to produce a metal and an acid. A review of work done on this phenomenon since that time may be found elsewhere., The hydrogen reduction of a cobalt salt solution is facilitated by complexing the cobalt ion. An ammonia complex of cobalt has been reduced commercially in the recovery of cobalt metal. A new reducible complex of cobalt was discovered5 when it was found that a co-baltous sulfate solution containing ammonium acetate could be reduced by hydrogen at temperatures in the region of 200°C. When a cobalt sulfate-ammonium acetate solution was heated to a temperature below its normal boiling point, a violet color became apparent, indicating complex formation. The nature of this complex was investigated5 by the addition of NH4Ac to a CoSO4 solution maintained at 85o. During the first additions of NH, Ac, the pH of the solution remained fairly constant at about 5.85. However, as the ratio of NH,Ac to CoSO, approached two, the pH rose and then leveled off at about 6.05. The absorption spectra of a Co(Ac), solution and a CoSO,, NH Ac solution were obtained at 85°C and were compared and found to be the same. These findings suggested that the diacetate complex of cobalt, Co(Ac),.4H20, was formed at 85°C. When a cobalt sulfate-ammonium acetate solution was heated to a temperature above about 165o, a finely divided pink precipitate appeared. The X-ray diffraction pattern of this precipitate indicated that it was Co(Ac), - 4H,O. In addition, it was discovered that when chloroplatinic acid, H,PtCl;, was added initially to the cobalt sulfate-ammonium acetate solution, a faster reduction was obtained. The roles of the solution complex, pink precipitate and chloroplatinic acid in the reduction process were then investigated. APPARATUS The experimental work was carried out in a two-liter stainless-steel autoclave. Adetailed description of the autoclave and the auxiliary equipment used in maintaining constant temperature and pressure may be found elsewhere.= Because the stainless steel was corroded, and also because it acted as a hydrogena-tion catalyst, an all-glass liner was fabricated such that the solution came only into contact with flame-polished pyrex glass. EXPERIMENTAL PROCEDURE The solutions used in the experimental work were prepared by dissolving reagent grade chemicals in distilled water. Although variation of the brand of ammonium acetate appeared to have no effect on the experimental results, CoSO, 7H O from the J. T. Baker Chemical Co. of Phillipsburg, N. J.,was found to give faster reductions than that prepared by Allied Chemical and Dye Corp., N. Y. The former was used throughout the course of the experimental work and was weighed up at the outset of each experiment. The ammonium acetate was dissolved to form a 6M stock solution, which was stored under refrigeration. A 10 pct solution of chloroplatinic acid (J. T. Baker Chemical Co.) was diluted to a 1.15 x 1Q2 M stock solution, which was standardized by precipitation of K,PtCl, as outlined by Scott. The appropriate volume of the chloroplatinic acid, H,PtCl,, solution was pipetted into the clean, dry glass liner. The cobalt sulfate-ammonium acetate solution, which had previously been saturated with

Jan 1, 1962

By G. W. Brock

Experiments in the form of aging of overstrained tension specimens and elevated temperature tension testing, have been carried out on recrystallized arc-cast molybdenum. The aging behavior of molybdenum was found to be strongly dependent on the degree of plastic strain present in the lattice: the reason for this being the low solid solubility for the common interstitial elements in molybdenum. GEORGE W. BROCK, formerly Research Associate, Department of Theoretical and Applied Mechanics, University of Illinois, Urbana, III., is Research Associate, Department of Engineering Sciences, IBM Research Center, Yorktown Heights, N. Y. Manuscript submitted January 19, 1961. IMD THE term strain aging is familiar in connection with the recovery of the yield point in mild steel after tensile overstraining, and there are other instances where it is known to exert an appreciable effect on mechanical properties. For example, steel exhibits a greater tensile strength at mildly elevated temperatures than at room temperature and serrated stress-strain flow curves occur under certain conditions. It has been known for a number of years that the phenomenon of strain aging is directly connected with the presence of carbon and nitrogen in the lattice; and extremely small amounts of either element are sufficient to produce the effect. In this study, the aging behavior of overstrained tensile specimens of recrystallized molybdenum has

Jan 1, 1962

By M. H. Pilkuhn, H. Rupprecht

The junction luminescence of GaAsxP1-x diodes containing up to 47 pct Gap was studied. Diodes were prepared by diffusing zinc into n-type material which was eithw boat- or vapor-gvown. Observations concerning the energy and the line shape of the spontaneous emission are presented and discussed; in particular the influence of inhomogeneous phosphorus distribution is considered. The role of reabsorption of the junction radiation is pointed out. With increasing current, the peak emission shifts to higher enwgy. This is alm true in cases where slow line narrowing occurs simultaneously. The length dependence of the threshold current density for lasing at 77 °K was examined for lasers made from the same material. It was found that longer lasers have lower threshold curvent densities, Estimates of the characteristic laser losses are given. Data concwning the range of Gap concentration in which stimulated emission was observed are presented. FORWARD-BIASED p-n junctions in mixed alloys of GaAsxPl-, emit light whose wave length depenfs on the composition of the alloy ranging from 8400A (x = 1) to 5350A (x = 0) at 77°K. Laser action in vapor-grown material was first reported by Holonyak and evacqua,' and in boat-grown material by Ainslie, Pilkuhn, and Rupprecht.' In both cases this was observed for low Gap concentration where GaAs,Pl-, is still a "direct" semiconductor with the lowest conductance-band minimum at k = 0 and a subsidiary minimum in the [loo] direction. As the Gap concentration is increased, the subsidiary [loo] minimum in the conduction band is lowered and finally becomes the lowest minimum for Gap concentrations above 50 pt,3 making the alloy an ''indirect" semiconductor. It is of interest to determine up to which Gap concentration stimulated emission can be observed. EXPERIMENTAL Boat-grown material was prepared by N. Ainslie of IBM by horizontal Bridgman technique. The material was n type and silicon-doped up to a concentration of 10'' cm-'. Vapor-grown material was made by E. Hull of IBM or received from Merck & Co. Diodes were fabricated in a similar manner as in the case of GaAs by zinc diffusion. Details about the diffusion can be found in a later paper.4 All diodes were of Fabry-Perot structure with two parallel sides cleaved and the other two sides sawed. Their width ranged from 50 to 150 p and their length from 50 to 700 p. RESULTS AND DISCUSSION Spectral Distributions. A typical display of the junction luminescence of a GaAsxPl-, diode is shown in Fig. 1. The liquid-nitrogen and room-temperature curves which were taken at a very low current (25 ma dc) show a relatively sharp high-energy emission which is band-edge radiation. They also show a broad low-energy emission which is due to transitions to deep impurity levels. This latter emission line may have different shapes and positions depending on the nature of the material. It can sometimes be quite near the band gap giving rise to "bumps" in the high-energy line. The results presented in this paper will be concerned only with the high-energy emission which becomes more intense relative to the low-energy emission when the current is increased and the temperature is lowered.

Jan 1, 1964