Search Documents

By D. T. Peterson, T. Carnahan

The diffusion coejTicients of carbon, nitrogen, and oxyget were determined in $ thorium over the tempernilcre range 1440" io 1715°C. The diffusion coyfiicir?zls are given by: D = 0.022 exp (-27,000/RT) jor carbo)~, D = 0,0032 exp(-l7,00Q/RTj for nitrogen, and D = u.0013 expt(-11,UOU/RT) for oxygen. Cavl~orz was found to increase the hardness of thoriunz nearly linearly with concentration over the range 100 to 1000Ppm carbon. ThORIUM has a fcc structure up to 1365°C and a bcc structure from this temperature to its melting point at 1740°C. Diffusion of carbon, oxygen, and nitrogen in bcc thorium was of interest in connection with the purification of thorium by electrotransport.' In addition, it was possible to measure the diffusion of all three of these interstitial solutes in the same bcc metal. Only in niobium, tantalum, vanadium, and a iron have all three interstitial diffusion coefficients been measured in a given bcc metal. Diffusion coefficients have been measured for carbon and oxygen in a thorium by Peterson2, 3 and for nitrogen by Gerds and Mallett.4 Activation energies for diffusion are reported by the above authors to be 38 kcal per mole for carbon, 22.5 kcal per mole for nitrogen, and 49 kcal per mole for oxygen. Values of the diffusion coefficients of carbon and nitrogen in 3 thorium have been reported by Peterson et al.' However, these were secondary results of their investigation of electrotransport phenomena in thorium and it was hoped that the present study could provide more precise data. EXPERIMENTAL PROCEDURE The specimens used in this study were the well-known pair of semi-infinite bar type. The couple was formed by resistance butt welding two 0.54-cm-diam by 3.0-cm-long bars of thorium together under pure helium, the concentration of the solute being greater in one cylinder than that in the other. The finished couple then contained a concentration step at the weld interface and diffusion proceeded only along the axis of the rod. The thorium used in this study was prepared by the magnesium intermediate alloy method.5 The total impurity content was less than 400 ppm. The major impurities were: carbon, 100 ppm: nitrogen, 50 ppm; and oxygen. 85 ppm. The total metallic impurity content was less than 150 ppm. The high solute concentration portions of the diffusion couples were prepared by adding the solute to the high-purity thorium in a non-consumable electrode arc melting procedure. Carbon and nitrogen were added in the form of spectroscopic graphite and nitrogen gas while a Tho2 layer was dissolved by arc melting to add oxygen. High-purity thorium formed the low concentration portions in the carbon and nitrogen couples. The low oxygen portions were obtained by deoxidizing high-purity thorium with calcium for 3 weeks at 1000°C according to a method reported by Peterson.3 The high C-Th contained 400 ppm C, the high N-Th contained 400 ppm N, the high 0-Th contained 220 ppm 0, and the low 0-Th contained 25 ppm O. The high O-Th was brine-quenched from 1500°C to retain most of the oxygen in solution at room temperature. These concentration levels were all below the solubility limits in 0 thorium at 1400°C. A resistance-heated high-vacuum furnace was used to heat the couples. The samples were mounted horizontally on a tantalum support which had small grooves near each end. Spacer rods of thorium, 0.4 cm in diam, were placed in these grooves to prevent contact between the sample and the tantalum support. This arrangement should have prevented contamination of the sample by contact with the support. In further effort to reduce contamination, the oxygen diffusion couples were sealed inside evacuated outgassed tantalum cylinders lined with thorium foil. Thorium rings around each end of the samples acted as spacers in this case. Pressure during diffusion runs was about 10-6 torr after an initial outgassing stage. Temperature measurements were made by sighting on black body holes in the sample support adjacent to the samples with a Leeds and Northrup disappearing-filament optical pyrometer. Temperatures were constant during a diffusion anneal to ±5C. The observed temperatures were corrected for sight glass absorption after each diffusion run. The pyrometer was checked against a calibrated electronic optical pyrometer and a calibrated tungsten strip lamp with the electronic pyrometer being taken as the standard. All temperature readings agreed to within ±3C over the temperature range 1450" to 1690°C. Time corrections due to diffusion during heating and cooling were necessary because of the short diffusion times. The diffusion times ranged from 6 min for the oxygen sample run at 1690°C to 90 min for the carbon sample run at 1500°C. A series of temperature vs time plots were made for heating and cooling of the samples to the various diffusion temperatures. This data was then used in a method according to shewmon6 to determine the time corrections. The corrections amounted to

Jan 1, 1970

By H. G. Doll

As a continuation of the earlier paper on the general subject of the SP log, a more complete analysis of certain features of the SP log in shaly sands is given. The pseudo-static SP in front of shaly sands is compared, on a theoretical basis, to the static SP in front of clean sands, as a function of the respective amount of shale and sand in the formation, and of the relative resistivities of the shale, of the uncontaminated part of the sand. and of the invaded zone of the sand. As a conclusion, the advantage of using reasonably conduc. tive mud in this case is shown. The discussion is illustrated by field examples. INTRODUCTION The discussion reported in the present paper is based on a theoretical analysis, and not on experiment. The field examples, joined to the text. are shown only as qualitative illustrations of the essential results of this analysis. Although the hypotheses made in the theoretical developments may perhaps be somewhat improved, it seems, nevertheless, that the results obtained account reasonably well for the actual phenomena, and give a fair approximation of their order of magnitude. The paper contains a mathematical analysis of a tri-dimen-sional distribution of potentials and current lines. due to spontaneous electromotive forces arising at the contact of shales and free electrolytes. as a function of the geometry and of the respective resistivities of the different media involved. It is assumed, although this hypothesis is not proven, that the emf's remain the same even if the shale occurs in very thin layers or in dispersed particles. It has already been pointed out 1,2,3 that, all other conditions being the same, the deflection of the SP log in front of a shaly sand is smaller than opposite a clean sand. When the thickness and the conductivity of a clean sand are large enough. the deflection of the SP log reaches a limiting value which is equal to the "static SP" of the clean sand. It is generally convenient to take the static SP of shale as the reference value or "base line." As a consequence, and for the sake of abbreviation, the expression. "static SP of a clean sand," is often used to designate the difference between the static SP of that sand and that of the shales, which difference is a measure of the total electromotive forces involved in the chain mud sand-shale. A similar limiting value is; also observed for the SP deflec-lion opposite a thick shaly sand, but it is smaller. just as if the total electromotive force involved were smaller in that case. This limiting value has been called the "Pseudo-Static SP" of the shaly sand. The static SP of a clean sand depends on the salinity of its connate water with respect to that of the mud, and, to a certain extent. on the differential pressure which controls the electro-filtration potentials, but it does not depend on the resistivity of the sand. On the contrary. the pseudo-static SP of a shaly sand depends not only on the salinity of its connate water and on the differential pressure, but also on the percentage of shale and on the resistivities of the shale, of the uncontaminated part of the sand, and of the zone invaded by the mud filtrate. If the three resistivities above were equal, the pseudo-static SP would be proportional to the percentage of sand in the shaly sand, and its departure from the static SP of a clean sand having the same connate water would simply be proportional to the percentage of shale. In that case, the pseudo-static SP of a shaly sand containing 10 per cent of shale would he 10 per cent less than the static SP of a clean sand. When. however. the sand is. on the average. substantially more resistive than the shale. the percentage of departure of the pseudo-static SP from the static SP of a clean sand is much larger than the percentage of shale. For that reason, the peaks of the SP log opposite shaly sands are systematically of smaller amplitude when the sands are oil-bearing than when they are water-bearing, all other conditions being the same. This feature is observed even when the sand beds are thick. and even when they do not contain a large percentage of shale. All this has already been described in all earlier publication", but mostly in a qualitative way. The present paper will analyze in more detail the action of the local SP currents which are generated inside of the shaly sands, and which are responsible for the abnormally low value of the pseudo-static SP. The quantitative computations have been extended to the general case of thin interbedded layers of sand and shale, where the resistivities of the shale and sand streaks do not have the same value: they are summarized in charts giving values of the pseudo-static SP of a shaly sand as a function of the different parameters involved. DEFINITIONS The static SP of a clean sand has been defined as the potential that would exist in the mud opposite that sand, were the SP current prevented from flowing. Such an ideal condition is represented on Fig. I-A. By analogy, the pseudo-static SP of a shaly sand can be defined as the potential that would exist in the hole, if the circuit shaly sand — surrounding shales — mud column were interrupted by the insulating plugs placed at the boundaries

Jan 1, 1950

By Craig S. Hartley

An approximate expression is obtained for the energy of a closed dislocation loop in equi1ibriu)n with a constant net stress. The result obtained is valid for loops in isotropic or anisotropzc materials provided that they are suJficiently large that the energy per unit length of a segment of the loop can be approximated by that of an infinite straight dislocation tangent to the loop. It is shown that this approximation leads to very close agreement with a more rigorous calculation of the elastic energy of a circular glide loop. The Gibbs-Wulff Form, GWF, of a dislocation is the closed planar loop which has the smallest elastic self-energy of all possible loops having the same Burgers vector and enclosing a fixed area, A.' The energy of such a loop is related to the net resolved shear stress* required to expand the loop and to the stress required to activate a Frank-Read source.223 In the following sectiorls the problem of determining the form of the GWF is discussed and an approximate method for calculating its elastic self-energy is presented. It is demonstrated that the approximations employed lead to no serious errors when applied to a calculation of the elastic energy of a circular glide loop. This method is then used to obtain a closed form expression for the energy of GWFs in isotropic and anisotropic materials. THEORY Burton, Frank, and cabrera4 have proved that the relationship of the equilibrium shape of a two-dimensional array of atoms under the influence of the Gibbs free energy associated with unit length of its boundary, G(O), is that the polar plot of G(0) vs 0 is proportional to the pedal of the GWF.* The angle 0 is measured "The pedal of the polar graph ofG(0) vs0 is the envelope of tangents to the eraph.relative to some crystallographic reference direction. The difficulty in applying this result to a closed dislocation loop arises from the self-interaction of the loop. For a dislocation the energy analogous to G(0) is a function of the total configurati~n.~ Consequently the relation which determines the GWF is an integro-differential equation rather than the simple differen- tial equation which results when G(8) is a function of 0 alone. Mitchell and smialek3 and Brown~ have used the self-stress concept introduced by ~rown' to calculate the shapes of dislocations in equilibrium with an applied stress. In this approach the glide force on an element of the dislocation loop due to the interaction of the element with the rest of the loop is equated to the glide force exerted by the local applied stress. The shape of the loop is then adjusted so that the two forces above are equal at all points on the loop. It is possible to calculate the energies of such loops by noting that, for equilibrium with an applied stress, the energy is equal to pijbiAj (summation convention) where bi is the Burgers vector, p.. is the local net stress tensor, and Ai is a vector directed perpendicular to the plane of thd loop with magnitude equal to the area of the loop. Also Brown' has calculated the energy of a hypothetical polygonal GWF using the above technique and anisotropic elasticity. However, his indicated solution for the energy in the general case of an arbitrary GWF is only slightly less involved than an iterative solution of the integrodifferential equation referred to earlier. In the present work the approximation employed by DeWit and Koehler' is used to calculate the energy of a closed loop in equilibrium with an applied stress. That is, the energy of a loop segment, ds, is approximated by the product of ds and the energy per unit length of an infinite, straight dislocation in a cylinder coaxial with the tangent to the loop at the angular position of the segment. This is known as the "line tension" approximation. The inner cutoff radius of the elastic solution defines the core radius, while the outer cutoff radius is determined by some characteristic dimension of the loop. Actually, both of these radii vary with the edge-screw character of the segment. The effective core radius changes because of the orientation dependence of the Peierls width of a dislocation,8 and the outer radius should be the radial distance from the circumference of the loop to the center of symmetry of the area enclosed by the loop.g However, since the energy varies logarithmically with the ratio of these radii while depending directly on the effective elastic constants, only the effect of the latter is considered. This approximation also neglects the self-interaction of the loop segments. For small loops this will doubtless be extremely important, but for large glide loops produced by plastic deformation the self-interaction is not nearly so important in determining the energy of the loop. This point is illustrated by the following calculation of the energy of a circular loop. Consider a circular loop of radius R which lies in the XI - x, plane of an infinite isotropic continuum and whose Burgers vector makes an angle $ with xs. The first-order solution for the elastic self-energy is:'

Jan 1, 1970

By P. J. Burkhardt, L. V. Gregor

The oxidation kinetics of single-crystal silicon has been investigated using extremely dry oxygen as the oxidant. Two techniques were used. The first involved a flow system with which incremental thickness measurements were made. This sytem provided an oxygen ambient at atmospheric pressure which had a moisture content of less than 0.1 ppm of water. The second technique involved measuring the change in volume of oxygen with time in a closed system. A temperature range of 900O to 1200 "C and a pressure range of 20 to 760 Tory were covered. Both techniques showed that a simple parabolic rate law holds only at the higher temperatures. At lower temperatures, the data up to 5000Å can be fitted well to a linear-parabolic equation. The two constants in this rate equation are examined in terms of a physical model. The parabolic portion of the low-temperature oxidatLon fits along the high-temperature curve of an Arrhenius plot. The parabolic rate constant for the flow-system technique can be expressed by k = 2 x 10-9 exp (L 31 kcalRT) sq cm per sec. The manostatic technique gave k = 1 x 10-exp (- 23 kcalRT) sq cm per sec. By extending the thickness range to over 104. a rate law of X = At seems to fit all of the data better than the linear-Pavabolic law. Here the exponent n increases from 0.518 at 1150 V to 0.637 at 900°C. PASSIVATION of silicon planar devices through thermal oxidation is of such great importance that considerable effort has been devoted to the study of the process.17 Much of this work involves hydrothermal oxidation, in which H2O is the oxidant: whereby H2 is produced. It has been proposed that the reaction proceeds by the formation of Si-OH groups and hydroxyl diffusion through the film to react with the silicon surface. Thus, the silicon oxide layer probably contains residual hydrogen in some form after oxidation in H2O and it has already been established that hydrogen has a marked effect on the electronic properties of oxide-passivated silicon surfaces. Silicon is oxidized by oxygen at elevated temperatures; but the rate of oxide formation is considerably lower than that for 0." By using dry O2 as the sole oxidant, the presence of residual hydrogen is presumably avoided. At present, the mechanism of thermal oxidation by dry O2 is a matter for conjecture, and the presence of small amounts of H2O was reported to have a marked effect on the rate of oxide formation10 (perhaps on the electronic properties of the surface, as well). Hence, this study of the kinetics of silicon oxide growth in 0, was undertaken: 1) to determine the oxidation rate at atmospheric pressure in a system designed to exclude rigorously all traces of H2O in the ambient 0,; 2) to obtain oxidation-rate data in situ in dry O2 by means of a continuous monitoring of the amount of 0, consumed during oxide growth. These objectives are complementary—the first assesses what effects (if any) are obtained by excluding H2O from the reaction, and the second expands this information while at the same time avoiding the uncertainty inherent in relying upon incremental measurements of oxide growth to furnish kinetic data. EXPERIMENTAL PROCEDURE The conventional silica open tube in a resistance-heated oxidation furnace is subject to sources of trace H2O, even after the input O2 is thoroughly dried and the outlet gas stream is passed through a trap to prevent

Jan 1, 1967

By R. C. Craze

During the past 20 years. research and development in the study of reservoir behavior have dealt principally with flow of oil through sandstones. Many reservoir studies of sand fields have proved valuable in promoting recovery efficiency. This paper discusses fundamental principles governing oil and gas production from sandstone and limestone alike and presents the results of investigations relating to the application of analytical techniques used for sandstone reservoir studies to the study of limestone reservoir performance. The characteristics of limestone porous systems, porosity-permeability relationships, distribution and occurrence of oil, and characteristics of flow through such systems are discussed. Recognition is made of the similarities or differences which these factors exhibit in limestone and sandstone systems. Comparisons between operating data for typical limestone and sand reservoirs are presented. It is indicated that the distribution and movement of fluids in and through porous limestones follow the same fundamental principles underlying such processes in sandstones. This fundamental similarity may readily be discernible in the performance of many limestone reservoirs. The volumetric balance and unsteady state radial flow equation, the fluid displacement equation, use of electrical analogue devices, and other analytical techniques to study the behavior of limestone fields appear fundamentally applicable, but do require thorough understanding of the properties of the formation, of the fluids, their behavior during flow, and adequate production operating data. Need for more complete coring and comprehensive examination of core properties is stressed. The results of "active oil" studies, and of flow and interference tests are presented. Well spacing, well completion, and efficient rates of production in limestone reservoirs are briefly discussed. INTRODUCTION Limestone and dolomite reservoirs constitute the largest source of supply of crude oil in the world. an estimated 60 per cent of present production coming from carbonate reservoirs. In many large geologic provinces such as Mexico. the Middle East, and more recently Canada, almost all the oil is found in this type of rock. In the United States, all of the major oil-producing areas except California and Pennsylvania contain oil-bearing carbonate formations. The discovery in recent years of large oil reserves in the Silurian, Devonian, and Ordo-vician formations in West Texas, in addition to the large reserves in the Permian, has accentuated the interest of operators, geologists, and engineers in limestone formations and in the many problems associated with understanding the performance of these reservoirs. The rapid increase in discovery of oil in limestone formations and the present-day position of prominence held by these fields in the production and reserve picture in all parts of the world emphasize the horizons opened to the reservoir technologist in the field of geological and production research. Pertinent to an interpretation of the behavior of limestone reservoirs are the methods of analysis which may be utilized and a possession of full knowledge of the many factors which influence the analytical procedures. During the past 20 years a well-developed science of reservoir engineering has been built upon the research of many workers who studied the fundamental nature of oil reservoirs, characteristics of the porous media, properties and behavior of the contained fluids, and the mechanics of flow. Application of these studies to production practice has resulted in greater recoveries and more efficient field operation. The major portion of this evolutionary process has been founded upon studies of sand fields. The applicability of these more thoroughly developed techniques for studying sand fields to the study of the behavior of limestone fields becomes a factor of technical and practical significance. In the light of the technological background available to the reservoir analyst, this paper discusses fundamental principles governing oil and gas production from sandstone and limestone alike and presents the results of investigations relating to the application of analytical techniques used for sandstone studies to the study of limestone reservoir performance. LIMESTONE RESERVOIR CHARACTERISTICS The characteristics of limestone and sandstone reservoirs are similar in many respects and they both may occur under similar structural conditions. Fundamentally, the distribution and movement of fluids in and through the porous limestone media follow the same basic principles which dictate such processes in sandstones. Herein lies a fundamental similarity, which may readily be discernible in the performance of many limestone reservoirs. In some limestone fields, as in many sand fields, widely varying formation properties and distribution may reveal themselves in deviations in the performance of the reservoir, in the behavior of wells, and in fluid flow through the rock, and make difficult the delineation of reservoir behavior. Only by reservoir studies, core analyses, and coordinated laboratory and field experimentation can the effects of the many influencing factors upon the nature of limestone production be determined. FORMATION CHARACTERISTICS The chief difference between limestone and sandstone, aside from their chemical compositions, is the difference in the geometry and origin of the porous systems in the two kinds of rock. In sandstone the porous system results entirely from the openings among individual sand grains which occur during deposition. The geometry of the openings between the sand

Jan 1, 1950

By R. F. Dickerson, D. A. Vaughan, A. F. Gerds

ALTHOUGH the metallurgy of uranium has been under intensive study since the early 1940's, no systematic effort has been made to identify the non-metallic inclusions in uranium. Uranium carbide (UC), which is probably the most common inclusion found in graphite-melted metal, has been tentatively identified by previous investigators, but the other nonmetallic inclusions have received little attention. Since metallography is a valuable tool in metallurgical studies, the metallographic identification of the nonmetallic inclusions in uranium is important. Such an investigation has been completed and the identification of slag-type inclusions and of uranium monocarbide, uranium hydride, uranium dioxide, uranium monoxide, and uranium mononitride is described. Metallographic Preporation It is often possible to prepare specimens for metal-lographic examination equally well by several methods. The specimens which were examined in this work were prepared by one of two acceptable methods. For the convenience of the reader, both methods will be discussed in detail and will be referred to simply as Method I or Method II in the subsequent sections. For both Methods I and 11, specimens for microscopic examination usually were mounted either in bakelite or in Paraplex room temperature mounting plastic. Method I—Specimens were ground in a spray of water on a revolving disk covered successively with 120-, 240-, and 600-grit silicon carbide papers. It was necessary to perform the final grinding operation carefully on worn 600-grit paper to keep the scratches as fine as possible. After washing and drying, the specimens were polished for 3 to 4 min on a slow speed wheel (250 rpm) covered with a medium nap cloth. Diamet Hyprez Blue diamond polishing paste, Grade 00, 0 to 2 µ, was used as abrasive with kerosene as lubricant on the wheel. Specimens were washed thoroughly in alcohol and final polished electrolytically in an electrolyte composed of 1 part stock solution (118 g CrO, dissolved in 100 cm3 H2O) with 4 parts of glacial acetic acid. A stainless steel cathode was used. At an open circuit potential of 40 v dc, a polishing time of 2 sec retained inclusions well with the bath at room temperature. If additional etching was required to sharpen the interface between the metal and the inclusions, an electrolyte composed of 1 part stock solution (100 g CrO3 and 100 cm8 H20) and 18 parts glacial acetic acid was used at room temperature. Best results were obtained by etching for from 10 to 15 sec at 20 v dc in the open circuit. Surfaces obtained by this method are suitable for microscopic examination. However, if desired, they may be etched further with other chemicals. Method 11—Rough grinding was done on a wet 180- or 240-grit continuous grinding belt. The specimen was then ground by hand successively on 240-, 400-, and 600-grit silicon carbide papers in a stream of water. Final polishing was accomplished on a 4 in. high speed wheel (3400 rpm) covered with Forstmann's cloth. Linde B levigated alumina, suspended in a 1 volume pet chromic acid solution, was the abrasive. Specimens usually were polished in 5 min or less by this technique. Often the inclusions present in the metal were identified in the mechanically polished condition. When etching was required to outline inclusions more sharply, one of the two following methods was used. In the first method, the specimen is etched lightly while electropolishing in the chromic-acetic acid solution described above (1 part of stock solution to 4 parts of acetic acid). The electrolyte was refrigerated in a dry ice-ethyl alcohol bath and specimens were etched at 60 v dc on the open circuit for 2 or 3 cycles of 3 to 4 sec each. The second technique utilizes electrolytical etching at about 10 v dc (open circuit) in a 10 pet citric acid solution at room temperature. X-Ray Diffraction Technique The major problem in the identification of inclusions in metals by X-ray diffraction techniques is the extraction of a sufficient amount of each type of inclusion to obtain an X-ray diffraction pattern. In the present study, X-ray diffraction patterns were obtained from individual inclusions of the order of 10 µ diam. The polished and etched samples shown in the micrographs were examined at a magnification of X54 or XI00 with a binocular microscope. This allowed sufficient working distance to extract the inclusions with a needle probe for powder X-ray diffraction analysis. Friable inclusions such as MgF2, CaF2, UO2, and UH3 could be freed from the metal by probing the as-polished and etched surface. The fine particles then were picked up on the end of a Vistanex-coated glass rod (0.002 in. diam) which was held in a brass adapter made to fit the powder X-ray diffraction camera. The end of the glass rod was centered in the path of the X-ray beam. In the case of the UC, UO, and UN inclusions which are smaller in size, more metallic in appearance, and less friable than the other inclusions, it was necessary to etch the inclusion in relief before extraction. UN inclusions etched sufficiently in relief in the electrolytic polishing solution described in Methods I and II by increasing the polishing time. UN inclusions were relief etched by extending the

Jan 1, 1957

By E. S. Fisher, L. C. R. Alfred

The elastic anisotropy factors, c4,/c6,, c3,/cll, and c12/cl,, for hcp metal crystals vary significantly among the dgferent unalloyed metals. Significant variations with temperature are also found. The effects of elastic anisotropy on the dislocation in an elastic continuum with hexagonal symmetry have been investigated by computing the elasticity factors for the self-energies of dislocations in fourteen different metals at various temperatures where the elastic moduli have been reported. For most of the metals the effects of the orientation of the Burgers vector, dislocation line, and glide plane are small and isotropic conditions can be assumed without significant error. Significant effects of anisotropy are, however, found in Cd, Zn, Co, Tl, Ti, and Zr. The elasticity factors have been applied in the calculations of dislocation line tensions, the repulsive forces between partial dislocations, and the Peierls-Nabarro dislocation widths. It is predicted that the increase in elastic anisotropy with temperature in titanium and zirconium makes edge dislocations with (a), (a + c), and (c) Burgers vectors unstable in basal, pyramidal, and prism planes, respectively. The probability of stacking faults forming by dissociation of Shockley partials in basal planes also decreases with increasing c4,/c6, ratio, when the stacking fault energy is greater than 50 ergs per sq cm. The widths of screw dislocations with b = (a) in titanium and zirconium increase very significantly in prism planes and decrease in basal planes as c4,/c6, increases. The effects of elastic anisotropy on various dislocation properties in cubic crystals have received considerable attention during the past few years. In the case of cubic symmetry the departure from isotropic elasticity depends entirely on the shear modulus ratio, A = 2c4,/(cl, —c12); i.e., the medium is elastically isotropic when A = 1. Foreman1 showed that an increase in the ratio A produces a systematic lowering of the dislocation self-energy for a given orientation and Poisson's ratio. ~eutonico~, has shown that large anisotropy can have a marked effect on the formation of stacking faults by the splitting of glissile dislocations in (111) planes of fcc and (112) planes of bcc crystals. ~iteK' made similar calculations for (110) planes of bcc metals. Both studies of bcc metals showed that the large A values encountered in the alkali metals tend to reduce the repulsive forces between Shockley partial dislocations. In fcc metals, however, A does not vary over the large range encountered in bcc metals; consequently, the effect of A on the forces between Shockley partials is masked somewhat by the differences in Poisson's ratio between metals. The effect of A on the line tension of a bowed out pinned dislocation has also been investigated for cubic crystals, first by dewit and Koehler5 and more recent- ly by Head.6 In both cases the line energy model is applied and the core energy is not taken into account, thus making the conclusions somewhat tenuous with regard to the physical interpretation. Nevertheless, the fact that a large A decreases the effective line tension is clearly evident and the tendency for large A to produce conditions that make a straight dislocation unstable (negative line tensions) also seem evident. Head, in fact, shows visual microscopic evidence that stable V-shaped dislocations occur in 0 brasse6 For hcp metals the definition of elastic anisotropy is more complex and, furthermore, significant deviations from an isotropic continuum are found among a number of real hcp metals, especially at higher temperatures. The present work was carried out to survey the effects of elastic anisotropy on the elasticity factors, K, that enter into the calculations of the stress fields around a dislocation core. Some isolated analytical calculations have previously been carried out for several hcp metals but they are restricted in the dislocation orientations and temperature.8'9 The present computations are based on single-crystal elastic moduli that have appeared in the literature and consider various orientations requiring numerical computations. The results are then applied to survey the effects of temperature on the dislocation line tension and dislocation splitting in hcp metals. PROCEDURE Anisotropy Factors. The degree of elastic anisotropy in hcp crystals cannot be described by a single parameter, such as the A ratio in cubic crystals. The following three ratios must be simultaneously equal to unity in order to have an elastically isotropic hexagonal crystal: The magnitudes of these ratios at several temperatures, as computed from the existing data for the elastic moduli of unalloyed hcp metals, are given in Table I. There are no cases of complete elastic isotropy, but the large anisotropy ratios encountered in the cubic alkali metals are also missing. There are, however, several significant differences among the hcp metals, the most notable being the relatively small A and B ratios in zinc and cadmium and the differences in the magnitudes and temperature dependences of A. It has been noted that the temperature dependence of A has a consistent relationship to the occurrence of the hcp — bcc tran~formation. For cadmium, zinc, magnesium, rhenium, and ruthenium, A is less than unity at 4'~ and, with exception for rhenium, decreases with increasing temperature. In the case of rhenium, A has essentially no temperature dependence between 923' and 1123"~, so that it is clear that A does not approach unity at higher temperatures. Cobalt is similar to the above-mentioned group of metals in that it also does

Jan 1, 1969

By V. Leroy, A. G. Guy

Serious disagreements are often found between experimentally determined intrinsic diffusion coefficients and those calculated employing the usual form of the vacancy theory. In the new theory it is proposed that the total intrinsic flux, Ji, of component i, is the sum of a part, f, due to the usual random exclzanges of component i with the vacancies, and a second part, Ji, due to exchanges with the uacancies composing the net vacancy flux. The present treatment, while less powerful than that of Manning, has the advantage of easy uisualization and of facilitating the application of the vacancy-flux effect to complex systems. IT is becoming increasingly evident that there are serious deficiencies in the version of the vacancy theory of diffusion that has been widely used for the past 20 years. One type of evidence is the frequent lack of agreement between intrinsic diffusion coefficients and tracer diffusion coefficients, even taking account of the thermodynamic factor. A second kind of evidence is the observation of a Kirkendall shift larger than theoretically possible, that is, larger than can be accounted for without assigning a negative value to one of the two intrinsic diffusion coefficient.'- The thermodynamic factor could conceivably make both coefficients negative, but not just one. It is clear that a cause of these anomalies, apart from any inadequacy of the usual vacancy theory, might lie in an oversimplified treatment of the data. Adequate experimental techniques, including the use of moderate pressure during the diffusion anneal, are now available to insure that porosity, lateral expansion, and so forth, can be kept negligibly small in most cases. The effect of differences in atomic volume can be of major importance, and it is essential that one of the available methods4 be used to account for this factor. In the present treatment this is accomplished by the consistent use of moles per cubic centimeter as the unit of concentration. Of the various possible inadequacies of the vacancy theory, attention will be given here only to effects of the net vacancy flux. annin' has previously considered this question, beginning with an analysis of atomic jumping of tracer atoms. When he added the effect of a concentration gradient, new terms arose that could be associated with the flow of vacancies. The present treatment uses quite a different approach. The usual vacancy flux, J,, is introduced explicitly, and a simple analysis predicts major changes in the intrinsic diffusion coefficients from this cause. The usual assumptions are made that only a vacancy mechanism is operative, that the formation of voids can be neglected, and that changes in the partial molal volumes, vl and v2, are negligible. The significant diffusion coefficients for the present topic are Dl and D,, the intrinsic coefficients, which enter in the equations, where the flux Ji, moles per sq cm per sec, is that crossing the Kirkendall interface. The concentration, ci, is in units of moles per cu cm, and the concentration gradient, aci/ax, is evaluated at the Kirkendall interface. It will be recalled' that the calculation of Dl and D2 involves the measurement of areas on the diffusion curve with respect to the positions of the Kirkendall and Matano interfaces. In the case of the anomalies mentioned earlier, the Kirkendall shift is too large to be accounted for by the diiferetzce in fluxes (J2 -J1), given by Eqs. [I] and [2]. The logical inference is that the flux of the solvent atoms, J1, is actually in the same direction as the flux of the solute atoms, Jz. In terms of Eq. [I] this requires that Dl have a negative value. However, it would be somewhat misleading to state that the solvent atoms are diffusing up their own concentration gradient. The explanation that will be advanced here pictures competing processes producing the net flux of solvent atoms: 1) diffusion of the solvent atoms down their own gradient by random exchanges with vacancies; and 2) diffusion of solvent atoms in the opposite direction by exchanges with the net vacancy flux. ACTION OF THE NET VACANCY FLUX Theories of vacancy diffusion can be formulated with varying degrees of refinement, and the present theory has purposely been kept as simple as appeared adequate to explain the phenomenon in question. In particular the following aspects have been neglected: 1) the gradient of vacancy concentration in comparison to the gradient of the atomic species; 2) departure of vacancy concentration from the local equilibrium value; 3) variation of the jump frequency, LO, with the specific surroundings of the atom-vacancy pair being considered; 4) correlation effects. These and other refinements can be considered once the essential mechanism has been established. The essential idea of the present analysis is to calculate the total intrinsic flux, Ji, of component i as the sum, JlJ?±j{ [3] where J; is attributable to the usual random atomic jumping, and J{ is a contribution arising from the net vacancy flux, J,. The latter quantity, of course,

Jan 1, 1967

By Ralph Hultgren, Prodyot Roy

Manganese vapor pressures from 1250° to 1500°K were measured by conventional Knudsen and torsion-effusion methods in twelve Fe-Mn alloys with compositions from 9 to 80 at. pct Mn. The Knudsen re-sults agveed approximately with previous measurements found in the literature which indicated the solutions were nearly ideal. However, except for the higher manganese compositions, the initial torsion readings indicated much (up to 50 pct) higher vapor presszcres than the Knudsen method. These high readings decreased steadily with time. The results are interpreted as due to depletion of surface concentration of manganese during evaporation. Thus the initial torsion readings are the most nearly correct and Knudsen methods on alloys are to be regarded with suspicion unless it can be demonstrated that diffusion rates are rapid enough to replenish surface concentrations depleted by evaporation. For the lowest manganese contents and highest temperatures, depletion causes initial torsion readings to be too low. FEW thermodynamic data are available for alloys of the high-melting transition metals due primarily to experimental difficulties at the high temperatures involved. Measurement of vapor pressure is one of the most promising techniques to be applied to this problem. From the vapor pressure of a component of an alloy phase, its activity or its partial molar Gibbs energy can be directly calculated. From measurements over a sufficient range of temperatures and compositions the partial and integral Gibbs energies, enthalpies, and entropies may be determined through Gibbs-Duhem integration. For solid phases of variable composition, which commonly are found in alloys, measured vapor pressures are low because of depletion of the surface concentration of the more volatile component as it is selectively vaporized. Diffusion from the interior tends to restore the depleted concentration. The seriousness of the effect depends on the relative rates of diffusion and vaporization. Although depletion has been recognized as a factor,' direct measurements of its rate are lacking and quantitative estimates of its effect on measurements are uncertain. In the present work the vapor pressures above a series of Fe-Mn alloys have been measured by the torsion-effusion method. This method permits a measurement of the vapor pressure of manganese as soon as the sample comes to temperature and continuously throughout the run. The rate of decline of the apparent vapor pressure measures the rate of depletion and the torsion reading at zero time should be correct. At high manganese concentrations (70 and 80 at. pct Mn), the torsion reading remained nearly constant; surface depletion was negligible. There was a slow drop in measured pressure due to bulk loss of manganese from the sample. However, for lower manganese contents, especially at high temperatures, the depletion effect was considerable and the apparent vapor pressure decreased steadily with time. For these alloys the Knudsen method should give pressures which are too low. The measurement of pressure in the Knudsen method cannot be made until enough manganese has been vaporized to be weighable. The Knudsen result is therefore an integrated average between the initial and final vapor pressures. To verify this, conventional Knudsen measurements were also made with results in approximate agreement with previously published Knudsen measurements for these alloys. As expected, for the lower manganese contents the initial torsion readings of pressure were as much as 50 pct higher than the Knudsen. THE APPARATUS The apparatus, shown in Fig. 1, is capable of operating to temperatures up to 2000°K. A vacuum of 5 x 10"6 mm Hg can be maintained. Temperatures were measured to +3° by a W-Re thermocouple placed in the dummy cell, N. The torsion cell consisted of two alumina crucibles with holes in their covers. They were placed in the molybdenum holder shown suspended from the torsion wire. Rectangular torsion wires (1 by 4 mm) were found to have superior sensitivity and less residual distortion than circular (2-mm-diameter) wires. The precision of angular measurement was found to be *0.025 deg; the variation is no doubt due to temperature fluctuation and vibration. For the Knudsen experiments the alumina cell was set on the support, N. SAMPLE PREPARATION 800-g samples were prepared by melting together electrolytic iron and electrolytic manganese and pouring the melt into a water-cooled copper mold. Melting and pouring were done under a helium at-

Jan 1, 1965

By Glenn C. Waterman

THE diamond drill is a very important tool in exploration and development testing and its use is increasing. In almost all cases results of diamond drilling are analyzed on the basis of grade and tons. A proper evaluation of core and sludge assays is important if drilling results are to be acceptable as a basis for geologic and engineering appraisal. The relatively wide variation in assay averages as calculated by various well-known combining methods indicates that the engineering choice of a method may affect the outcome of the drilling in terms of ore and waste. The problem of combining assay results from core and sludge samples has been discussed many times in conference and in the literature.'-' Most writers agree that the field of disagreement in methods is large and that the engineer on the job must consider features unique to his drilling, pick one of several combining methods, and depart from the rules when abnormal results come in. All the discussion to date can be summed up by the admission that as yet there is no fairly simple, generally acceptable combining method that is practicable over a wide range of drilling conditions, ground conditions, and ore occurrence. The combining problem is important in evaluating drilling results at Chuquicamata. Recently a reappraisal has been made of recovery variables and their effect on assays, with the result that a new combining method is offered which fits average drilling conditions and is mathematically reasonable. It is simple in application, fundamentally correct, and an improvement over most combining methods. At Chuquicamata diamond drillholes are used to outline the grade and position of blocks of normal and marginal grade oxide, mixed, and sulphide ore. Most holes penetrate all classes of material (and waste), and it is important for mining programs as well as ore reserves to know almost precisely the soluble and insoluble copper content of mineralized ground. At present three classes of ore are mined and treated differently. For an orderly sequence of mining operations which can provide regular daily tonnages of all three ore types and keep grade at certain levels with minimum variation, ore type and its grade must be predicted. Diamond drilling plus geologic mapping and bench sampling are tools for prediction. And drilling data are largely used to calculate grade of material more than a few meters away from bench faces. The orebody at Chuquicamata" is criss-crossed by millions of barren or mineralized weak to fairly strong slips and fault fissures. Mineralization is diverse and encompasses many quartz and oxide or sulphide-bearing copper veins, as well as seams and disseminated grains. Copper occurs in oxide or sulphide minerals or mixtures of these two mineral types. Rock conditions vary from intensely seri-citized (soft and porous) through clay-altered ground to almost fresh granodiorite. The result is an orebody which offers many obstacles to good and consistent core recovery in diamond drilling. Recovery varies considerably in the several alteration zones, the various types of oxide and sulphide ores, and the position and inclination of the drillhole within the complex fracture pattern. As core recovery drops sludge samples must be used with core samples to calculate grade. Many years of drilling at Chuquicamata indicate that in good grade oxide zones and in the sulphide areas core recovery is good and the ground uniformly mineralized. Moderate loss of core, therefore, does not markedly affect grade calculations based on core assays. An early core-sludge combining method used core assays at face value as indicating grade down to 50 pct core recovery, but below this recovery percentage sludge samples were used and weighted according to the standard Longyear chart. This method apparently did not introduce serious errors, but it abruptly used sludge assays with high weighting factors representing 100 pct return irrespective of actual percentage of sludge recovered. Recent drilling activities have been directed toward outlining the marginal ore areas. The non-uniform mineralization and generally poorer core recovery in such ground indicated that a more exact core-sludge combining method was required to equate wide differences between core and sludge assays and recovery. In fringe ore areas at Chuquicamata core recovery averages about 50 pct, from a minimum of 10 to 15 pct to a maximum of 100 pct. Sludge recovery is likewise variable and averages perhaps 80 pct even though holes are cemented as water return falls off. In a homogeneously mineralized area cut by many slips and faults, with hard and soft ribs, loss of core is loss of ground which has a grade similar to that of core recovered, and core assays approximate true grade. In this case sludge samples need not be used. However, it would be unusual to know beforehand that an area is uniformly mineralized, and in fact this condition is probably uncommon. Generally the distribution of valuable minerals in the ground does not exactly compare with their recoverability in core. Thus, in the usual case, loss of core decreases the

Jan 1, 1956

By G. Bruggeman, K. G. Kreider

Grain boundary dij]usion coefficienls were measured in tungsten between 1400° and 2200° C and can be expressed by the equation sq cm per sec This activation energy confirms some eavlier estimates made .from tungsten sintering experiments. Grain boundary diffusion was found to occur in sub-bozrndavies having -misorientations of less than 10 deg. The actiuation energy for this subboundavy diffusion is equal to that for dijjusion in incoherent grain boundaries with in the limits of error. This is shown to be consistent with the dislocation model of Low-angle boundaries wheve diffusion occlcvs along- the dislocation 'YPipes" comprising -tile boundary. RECENT investigations of the sintering of tungsten powders all report activation energies which are considerably less than the activation energy for tungsten volume diffusion. Kothari' reports a value of 100 * 5 kcal per mole, Hayden and Brophy' obtained 90 kcal per mole, and Vasilos and smith3 found 110.7 kcal per mole from their sintering studies. Since most determinations of the activation energy for volume diffu-sion4-' fall between 120 and 160 kcal per mole (the true value seems most likely to be nearer 150 kcal per mole), the conclusion is drawn that the mass transport leading to densification during sintering is accomplished by grain boundary diffusion. This interpretation is consistent with various diffusion models of the sintering process. 10-12 Vasilos and Smith calculate diffusion coefficients from their data which fit the equation D * 1.36 x 10* exp(-llO,700/HD However, no direct measurements of tungsten grain boundary diffusion have been made. Furthermore, considerable disagreement exists between the directly measured values of tungsten volume diffusion.'-' In order to corroborate the inferred results of the sintering experiments concerning grain boundary diffusion and to provide accurate diffusion data essential to the analysis of the kinetics of creep, oxidation, precipitation, and so forth, the present work was undertaken to measure self-diffusion in single-crystal and polycrystalline tungsten between 1400" and 2200°C. It is within this temperature range that tungsten sintering is done, the re crystallization of tungsten occurs, and the widest application of tungsten as a high-temperature material will probably be made. EXPERIMENTAL PROCEDURE Radioactive WlE5 was produced by irradiating tungstic acid in a neutron flux of 1.2 x 1012 neutrons per sq cm per sec for 36 hr. A 2-week waiting period was allowed for the decay of w"~ also produced by the irradiation. (w"~ has a half-life of 24 hr.) The half-life of the remaining isotope was determined to be 75 days confirming the presence of w lE5 and the absence of any undesired radionuclide. Specimens 4 in. in diam and $ in. thick were cut from polycrystalline swaged tungsten rods (recrystal-lized) and from Linde single-crystal rods. Chemical analyses of these materials appear in Table I. Actually upon closer examination, the single-crystal specimens were found to consist of several subgrains separated primarily by tilt boundaries in which the misor-ientation ranged from 3 to 10 deg. Thus, it was possible to measure boundary diffusion coefficients in these low-angle subboundaries as well as in the incoherent boundaries of the polycrystalline specimens. The two faces of each specimen were ground flat and parallel within 0.0001 in. The radioactive tungstic acid was dissolved in concentrated ammonium hydroxide, placed on the ground flat of the specimen, and evaporated to dryness. The oxide was then reduced in hydrogen at 1000°C resulting in a layer of wlE5 approximately 1 p thick. The diffusion anneals were performed in vacuum in a tantalum resistance furnace. Time at temperature ranged from 10 hr at 1400°C to 2 hr at 2200°C. The penetration profile was determined by measuring the residual activity after successive removal of surface layers by grinding on metallographic polishing paper. Extreme care was exercised to insure that sections were always taken normal to the diffusion direction; this was verified repeatedly by checking that front and back surfaces of the specimen remained parallel. The activity was measured with an end-window Geiger-Mueller counter. The sides and edges of the specimen were well-shielded to eliminate possible effects due to surface diffusion. The weight of the

Jan 1, 1968

By R. S. Wagner, W. C. Ellis

A new mechanism of crystal growth involving oapor, liquid, crnd solid phases explains many observations of the effect of implurities in crystal growth from the vapor. The role of the impuuitq is to form a liquid Solution with the crystalline tnalerial to be grown from the vapor. Since the solution is n prefevred site for deposition firorti the uapor, the liquid becorrles supersaturated. Crystal growth occurs by precipitatzon from the supersaturated liquid crt tlie solid-liquid zntevfnce. A crystalline defect, such as a screw dislocation, is not essetztial for VLS (vapor -liquid-solid) growth. The concept of the VLS mechanism is discussed in detail with reference to tire controlled growth of silicon crystals using gold, platinum, palladium, nickel, silver, or copper as an implurity agent. RECENTLY a short communication' described a new concept of crystal growth from the vapor, the VLS mechanism. In this paper we present a detailed description of the process and its application to the growth of silicon crystals and we discuss its relevance to existing concepts of .'whisker" crystal growth. Crystal growth from the vapor is usually explained by a theory proposed by Frank2 and developed in detail by Burton, Cabrera, and Frank.3 In this theory a screw dislocation terminating at the growth surface provides a self-perpetuating step. Accommodation of atoms at the step is energetically favorable, and is possible of much lower supersatu-ration than required for two-dimensional nucleation. Crystals of a unique form resulting from aniso-tropic growth from the vapor are "whisker" or filamentary ones. Such crystals have a lengthwise dimension orders of magnitude larger than those of the cross section. For most filamentary crystals both the fast-growth direction and directions of lateral growth have small Miller indices. The special growth form for a whisker crystal implies that the tip surface of the crystal must be a preferred growth site. sears4 proposed that, according to the Frank theory. a whisker contains a screw dislocation emergent at the growing tip. Such an axial defect provides a preferred growth site and accounts for unidirectional growth. The hypothesis was extended by Price. Vermilyea. and Webb," still implying the presence of a dislocation at the whisker tip. They postulated that impurities arriving at the fast-growing tip face become buried while those arriving on the surface of slow-growing lateral faces accumulate and thereby hinder growth. These considerations led to a whisker morphology. There is increasing evidence that most whisker crystals grown from the vapor are dislocation-free. Webb and his coworkers6 searched for an Eshelby twist7 in zinc? cadmium, iron. copper, silver, and palladium whisker crystals. They found unequivocal evidence for an axial screw dislocation in only one element, palladium. However, not every palladium crystal examined contained a dislocation. Observations with the electron microscope have failed to show dislocations in whisker crystals of zinc, silicon.9 and one morphology of AlN.10 Since many whiskers are completely free of dislocations, an axial dislocation does not appear to be required for whisker growth of many substances. A significant advance in understanding whisker growth has been a recognition of the need for impurities. This requirement has been clearly demonstrated for copper,11 iron,13 and silicon9-1 whiskers. For silicon, detailed studies proved conclusively that certain impurities, for example, nickel or gold, are essential. Another pertinent phenomenon which has received little attention is the presence of a liquid layer or droplets on the surface of some crystals growing from the vapor. Crystals in which this has been observed include p-toluidine,14 MoO3,15 ferrites,16 and silicon carbide.'" The liquid layers or globules were considered to be metastable phases, molecular complexes, or intermediate polymers originating from condensation of the vapor phase. The possibility has been suggested that the halide being reduced is condensed at the tip18 or adsorbed on the surface11 of a growing metal whisker, for example copper. The literature on whiskers discloses illustrations of rounded terminations at the tips. These appear. for example, on crystals of A12O3,19,20 sic,21 and BeO.22 For BeO, Edwards and Happel suggested that during growth of the whisker the rounded termination consisted of molten beryllium enclosed in a solid shell of BeO. A recent paper9 on the growth of silicon whiskers contains many observations pertinent to an understanding of the mechanisnl of whisker growth. These observations are summarized as follows. 1) Silicon whiskers are dislocation-free. 2) Certain impurities are essential for whisker growth. Without such impurities the silicon deposit is in the form of a film or consists of discrete polyhedral crystals.

Jan 1, 1965

By D. F. Stein

In spite of the apparent usefulness of autoradiography in demonstrating segregation, it has had very limited success in demonstrating grain boundary segregation. Because of this limited success, a model system amenable to mathematical analysis was devised to determine which variables in the experiment are most important. As a result of these calculations, it was concluded that autoradiograPhy is a rather insensitive technique to measure grain boundary segregation. The range (energies) of the emitted particles (ß and a) must be low, and the concentration of the radioactive speczes at the grain boundary must be (in general) two or three orders of magnitude greater than the concentration within the grain. Because of these very restricted conditions , the limited success of the technzque is not surprising. In spite of the apparent usefulness of autoradiography in demonstrating segregation, it has had very limited success in demonstrating grain boundary segregation. Sulfur segregation in iron1 and possibly polonium segregation in Pb-5 pct B1 2 alloys are the only autoradi-ography experiments that have demonstrated segregation to grain boundaries in metals without the formation of a second phase. Segregation to a boundary without the formation of a second phase is often called Gibbs' absorption and is discussed by McLean in Chapter 5 of Ref. 10. There are several394 experiments showing grain boundary diffusion of a radioisotope, but this type experiment is not representative of equilibrium between the concentration of an element at a grain boundary and that in bulk, so it will not be discussed in this paper. In an attempt to determine if temper embrittle ment of low-alloy steels was associated with segregation of antimony to grain boundaries, a program to use auto-radiography (using Sb-125) was initiated. Even though other measurements strongly suggest that segregation is occurring during embrittlement, no evidence of grain boundary segregation was observed in the auto-radiography experiments. An attempt was also made to detect segregation of carbon (using C-14) to grain boundaries in iron during slow cooling. There is again strong indirect evidence7 that segregation occurs during such a treatment, but the autoradiography experiment gave no evidence of such segregation. Because of the failure of these experiments and the general lack of success by our metallographic unit in measuring grain boundary segregation using autoradi-ography, a model system amenable to mathematical analysis was devised to determine which variables in the experiment are most important. MODEL SYSTEM The model system is illustrated in Figs. 1 and 2. It is a semi-infinite bicrystal with a single grain boundary perpendicular to the top face. The width of the grain boundary, W, is defined as the region to which segregation has occurred and it is assumed that this region is of constant composition. The assumption of an infinite dimension is for mathematical reasons but the results of such an analysis (as will be demonstrated later) are valid for even very small specimens. It also is assumed that the measurement is not limited by means of detecting the radiation (photographic emulsions, counters, and so forth) but that the means of detecting radiation is linearly sensitive to the radiation and has infinite resolving power. The distance y is the perpendicular distance from the center of the grain boundary to the point at which the background radiation is measured and x is the integration variable. CALCULATION FOR BETA PARTICLES Two values of the intensity of radiation will be calculated, the intensity at the center of the grain boundary and the intensity far from the grain boundary. The radiation at the center of the grain boundary can be calculated in the following way. Assume a grain boundary of width W having a uniform concentration equal to Cgb + Cb where CGb is the excess concentration of the element under consideration at the grain boundary and Cb is the bulk concentration. The intensity of radiation, IgB, at the center of the grain boundary, is a consequence of the amount of radioactive material, decay rate, absorption, and geometrical factors which can be represented mathematically by the following expression:

Jan 1, 1968

By A. J. Griest, P. D. Frost, H. A. Robinson

The flow characteristics and uniform elongation of alpha-beta titanium alloys in the solution-treated condition were shown to be markedly affected by the solution temperature. Two classes of alpha-beta alloys were distinguished, based on the effect of solution temperature on the uniform elongation. Martensitic alpha-beta alloys, which are relatively weakly beta-stabilized, showed decreasing values of uniform elongation as the solution temperature was initially increased. However, after solution-treatment at still higher temperatures in the alpha-beta field, the uniform elongation of these alloys was again high, and often attained its maximum value. Nonmartensitic alpha-beta alloys, which are more heavily beta-stabilized, showed a continuous decrease of uniform elongation with increasing solution temperature. The behavior was rationalized in terms of the alpha-beta ratio and the alloy content of the beta existing at the solution temperature. The data indicate that in commercial practice close temperature control may be required to get maximum uniform elongotion, and hence, maximum potential formability, in solution-treated sheet. Bend-test data also were obtained for the solution-treated conditions. A MOST important objective of current titanium research is the development of age-hardenable alloys with superior formability in the solution -treated condition. Alpha-beta titanium alloys have been shown to be amenable to high-strength aging treatments after solution-treatment in the alpha-beta field.' However, heat-treatment data on the properties of solution-treated sheet, such as uniform elongation, minimum bend radius, and flow properties, which have a bearing on the potential formability, have not been widely available. This paper presents only the most significant results of a research program conducted at Battelle for the Air Force in an effort to develop sheet alloys that can be formed readily as solution-treated and subsequently aged to very high strengths. Data for 10 compositions representative of the more than 25 alloys studied in this program are presented. A primary objective of the paper is to demonstrate the pronounced dependence of the flow characteristics and uniform elongation on solution-treatment temperature. This effect is believed to be vitally important in the commercial heat-treatment of age-hardenable titanium sheet alloys. Experimental Procedures Melting and Fabrication—The nominal and actual compositions, beta transi temperatures, and rolling temperatures of 10 of the 30 alloys studied are presented in Table I. Stock :€or two of the alloys, E (Ti- 8Mn) and J (Ti-6A1-4V) was obtained from commercial heats. The remaining heats, except D and I, were double melted as 10-lb heats in a consumable-electrode furnace. Melting stock was 120 Brine11 hardness sponge. The oxygen and hydrogen analyses spotted throughout the heats indicate the range of interstitial levels resulting from the melting practice used. Heats D and I were melted as 20-lb ingots with the same grade of sponge. These ingots were cut in half and one portion used in the present evaluation. The ingots, which were about 4 in. diam and 4 in. in height, were upset forged and drawn at 1750°F to sheet bar 5/8 in. thick by 2 1/2 in. wide in such a manner that the top of the ingot corresponded to the top of the sheet bar. Slabs of the sheet bar, 5 1/2 in. in length were scalped to 1/2 in. thickness and rolled at the temperatures indicated in Table I to approximately 0.055 in. sheet. Elongation during rolling took place at 90" to the original forged length. The fabrication procedure used is described in some detail since the experimental heats showed considerable directionality in mechanical properties. The rolling procedures for the commercial heats were unknown, but were such as to produce essentially isotropic properties. All heats, as rolled, showed microstructures typical of alpha-beta alloys worked in the alpha-beta field; that is, the structures consisted of equiaxed alpha particles in a beta matrix. Heat Treatment— The effects of solution and aging heat-treatments on the properties of the alloys were evaluated uSing longitudinal and transverse tensile and bend coupons sheared from the sheet. Two or three solution-treated conditions were evaluated for each alloy. In addition, the alloys were also tested after annealing at 1200°F for 4 hr and furnace cooling. The solution-treated and annealed conditions were those in which the alloys might be expected to

Jan 1, 1959

By R. R. Jones, R. W. Kraft

Zn-Mg2Znn eutectic alloys nzay freeze willr either rodlike or lanzellar rnorphology. Alloys with slighlly more than /he eutectic arrzount of rnagnesillrn usually contain three-cnned dendrjles of MgzZnll in a eutec-lic ttlulris. All three morphologies haue the same cryslallographic orientution relationship: (0UOl) zn - 11 (111) Mg2Znll and (2310)Zn 11(101) Mg2Znll, but u3ith different prej-erred groulth direclions. The lurnellae lo rods transifion in con/rolled ingols qf euleclic cotnposition occurs because lhe large kinelic undercooling due to MgzZnll minirrzizes /he ejj-ecl of the solid-solid inlerface energy. The eutectic morphology is influenced by the presence of lhree-nned dendrites 0-f MgzZn11 which may conlrol /he rricroslrccture by acting as nuclealion sites. In recent years there has been much interest in eutectic solidification and several theories have been proposed. One of the confusing factors is the existence of various morphologies in which the solidified phases may form. The lamellar microstructure seems to be most common in metal eutectics, and it has been claimed' that all regular eutectics should be lamellar if sufficiently pure. However, there still remain eutectic alloys which are not lamellar or which change their morphology as a function of growth conditions. The eutectic between zinc and the intermetallic phase Mg2Znll was chosen for this investigation because it has been found to solidify in more than one morphology. The diagram in anssen' locates the eutectic point at 3.0 wt pct Mg and 367°C. lliott gives 364°C as the eutectic temperature, leaving the phase compositions unaltered. Since the growth conditions determine the micro-structure of the solidified alloy, the factors controlling the transition from one morphology to another could be studied. The lamellae to rods transition is of particular interest. PROCEDURE Alloys were prepared from carefully weighed portions of 99.999 pct Zn and 99.97 pct Mg by melting in Pyrex containers under argon and casting into graphite boats. The resulting ingots were remelted under argon and solidified unidirectionally in a horizontal tube furnace at growth rates ranging from 2.0 to more than 50 cm per hr under a temperature gradient, measured over a 5-cm length, of 9" to 14°C per cm. The solid-liquid interface appeared to be planar at all growth rates although no attempt was made to confirm this by decantation or quenching. A few ingots were allowed to freeze uncontrolled. Most alloys were of the nominal eutectic composition, 3.0 wt pct Mg according to Hansen2 and lliott, but some contained as much as 3.35 wt pct Mg. Chemical analyses were not run since metallographic examination confirmed that the desired composition was achieved. Specimens were cut from the middle portion of the ingot normal to the growth axis, polished mechanically, and etched with 2 pct Nital. Suitable areas were selected for the determination of crystallographic orientation relationships by a tiontechniqueof described previously by one of the authors.4 The (2310) planes of zinc and the (8701, {944}? (1032) planes of Mg2Znll were found suitable for orientation determination; experimental error was on the order of 2 or 3 deg. RESULTS Three different morphologies were found in the unidirectionally solidified alloys: lamellar eutectic, rod-like eutectic, and a structure whose most predominant characteristic was the presence of three-vaned (cellular) dendrites of Mg2Znll. These dendrites were only found in alloys with more than the eutectic amount of magnesium. In some ingots fine hexagonal needles of Mg2Znll surrounding a core of MgZn2 were observed. They were probably due to incomplete alloying and seemed to have no effect on the eutectic morphology. In addition hexagonal spirals like those discussed by Fullman and wood5 and Hunt and acksonh ere observed in some ingots frozen without directional control. Both MgZZn,, and MgZnz were detected by X-ray diffraction in these alloys. Since the morphology could not be grown unidirectionally and no characteristic orientation relationship between the phases was found, further study was limited to the lamellar: rodlike, and three-vaned dendrite morphologies. Alloys of Eutectic Composition, No Dendrites. The mcrostructures of allovs with no three-vaned dendrites were either lamellar or rodlike depending on the growth rate. At rates below 10 cm per hr the morphology was lamellar, consisting of two sets of parallel plates intersecting at about 54 deg like the Mg-MgzSn eutectic described by raft.7 At growth rates faster than 14 cm per hr the microstructure showed rods of zinc in a matrix of MgnZnll, while intermediate rates yielded mixtures of rods and lamellae in small groups. The lamellar "grains" were often several millimeters in cross section, but contained small irregular areas which divided each grain into perfect islands 100 or 200 p in diam. Lamellae were parallel to each other throughout the grain in spite of these defects in the structure, Fig. 1. Rods, on the other hand, could only be produced in small groups arranged like fish scales and separated by irregular areas of appreciable thickness, Fig. 2. Alloys Not of Eutectic Composition, With Dendrites. In alloys with 3.1 to 3.35 wt pct ME,-. three-vaned dendrites bf MgzZnll were usually found surrounded by eutectic. At growth rates slower than about 10 cm per hr the dendrites were separated from each other by small areas of both lamellar and rod eutectic, Fig. 3.

Jan 1, 1969

By E. Mille, R. Geffken

The Mg-Si and Mg-Ge phase diagrams were rede-levtnined by thermal analysis, and the existence of a single congruent melting compound in each system was confirmed. The melting points of the two compounds Mg2Ge and ,Wg2Si were found to be 1117.4° and 1085.0°C respectively. The euteclics for the Mg-Ge system occur at 635.6°C (1.15 at. pcl Ge) and 696. 7°C (64.3 at. pct Ge); for the Mg-Si system the eutectics are at 6376°C (1.16 at. pct Si) and 945.6°C (53.0 al. pcl Si). The phase diagrams and known thermodynamic data were used to calculate activity values for both systems. The activities calculated for the Mg-Ge system agreed very well with those previously published. Partial molar enthalpy values for the Mg-Si systetn were calculated from the phase diagram for the composition region where no experimental values have been reported. THE phase diagram for any system is an important source of thermodynamic data. Steiner, Miller, and Komarek1 have derived equations which permit calculation of the activity in binary systems with an inter-metallic compound! if the liquidus and enthalpy data are known. The thermodynamic properties of the Mg-Ge and Mg-Si systems have recently been determined in this by by an isopiestic method, and it was considered that it would be interesting to compare these directly determined values with those computed from the phase diagram. The basic features of the Mg-Ge and Mg-Si systems are essentially similar. The one intermediate compound present in each system. Mg2X, crystallizes in the antifluorite structure and melts congruently. Raynor4 has accurately determined the temperature and composition of the magnesium-rich eutectic in both the Mg-Ge and Mg-Si systems. Klemm and West-linning5 investigated the entire Mg-Ge liquidus, employing sintered alumina crucibles; the purity of the magnesium and germanium starting materials was not reported. The melt was not stirred, and the temperature was automatically recorded to an accuracy of ±3°C. The authors reported large weight changes due to magnesium evaporation between 50 and 67 at. pct Mg. The Mg-Si system has been studied by a number of investigators, and the results have been compiled by Hansen and Anderko.6 Significant discrepancies exist between the two principle investigations of voge17 and Wohler and Schliephake.8 Two different grades of silicon were used by Vogel, one of 99.2 pct purity and the other quite impure, containing 6 pct Fe and 1.7 pct Al. The magnesium purity was not specified. The melts were contained in graphite crucibles with porcelain thermocouple protection tubes under an atmosphere of hydrogen. Samples weighing 10 g were rapidly heated to 50° to 100°C above the liquidus: held, and then rapidly cooled without stirring. Accuracy was ±1 at. pct which is equivalent to a maximum error in temperature of ±18°C. Wohler and Schliephake used 97.9 pct Mg and 99.48 pct Si. The graphite crucibles contained a stirrer and the 15-g samples were melted under an atmosphere of streaming hydrogen. The samples were chemically analyzed after each run. Because of the scarcity of the data, the impurity of the starting materials, and the resultant uncertainty and inconsistency in the published liquidus values, it was decided to undertake a reevaluation of the Mg-Ge and Mg-Si phase diagrams by thermal analysis. EXPERIMENTAL PROCEDURE Alloys were prepared from 99.99+ pct Mg (Dominion Magnesium Ltd.) with impurities in ppm: 20 Al, 30 Zn, 10 Si, <1 Ni, <1 Cu. <10 Fe; 99.999 pct Ge (United Mineral and Chemical Corp.), and 99.999 pct Si (Wacker Chemie Corp.). All graphite parts were machined from high-density (1.89 g per cu cm) G-grade graphite obtained from Basic Carbon Corp. with a total ash content of 0.04 pct. Boron nitride parts were machined from rods of National-grade H.B.N. boron nitride. All graphite and boron nitride pieces were baked out under running vacuum at 1100°C for 24 hr before us Cylindrical graphite crucibles (1; in. OD, 23/4 in. long, l3/8 in. ID) were tightly closed with threaded graphite covers which had 21/4-in.-long thermocouple wells and 1/4-in.-diam off-center holes for stirrers. The cover and thermocouple well were machined from a single piece of graphite. A stirrer was made from a flat cylindrical graphite plate perforated with five 3/16-in.-diam holes and a 1/2-in.-diam central hole, and was held parallel to the crucible bottom by a 1/4-in.-diam. 4-in.-long graphite rod which screwed into the plate and extended up through a tightly fitting hole in the crucible cover. An iron core enclosed in a glass capsule was attached to the stirrer with an 18-in.-long molybdenum wire, so that the stirrer could be magnetically raised and lowered from outside the system. One crucible and stirrer with essentially the same dimensions given above was made entirely of boron nitride. Chunks of magnesium were premelted, cast into 11/2-in.-diam. rods, and then cut into lengths varying from a to 1 in. A 5/16-in. hole was drilled through the center of each piece to accommodate the thermocouple well and the individual pieces were then cleaned and rinsed with acetone. The total weight of an alloy was 50 to 70 g in the Mg-Ge system and 40 to 60 g in the Mg-Si system. The pure components were weighed to an accuracy

Jan 1, 1969

By H. Conrad, V. V. Damiano, G. J. London

The slip mode activated during the c axis compression of single crystals of commercial-purity ingot SR beryllium, high-purity (twelve-zone-pass) beryllium, and Be-4.4 wt pct Cu and Be-5.2 wt pct Ni alloys in the temperature range of 25° to 364°C was determined using two-surface slip trace analysis, slip-step height analysis, and electron transmission microscopy. All three techniques indicated the occurrence of copious pyramidal {1 122) (1123) slip in the alloys over the entire temperature range, the amount increasing with temperature. Pyramidal slip was also indicated in the high-purity beryllium by slip trace analysis and electron transmission microscopy, but the amount was somewhat less than in the alloys. For the commercial-purity ingot crystals, only a very small number of pyramidal slip lines were observed, and these were in the immediate vicinity of the fracture surface. No pyramidal dislocations could be detected by electron transmission microscopy in this material. Dislocatransmissiontions with Burgers vectors [0001] and +(ll20) were identified by electron transmission microscopy inthe (1122) slip bands, as well as those with the j (1123) vector. This was interpreted to indicate that the edge components of the 3(1123) vector dislocations activated during c axis compression dissociate upon unloading according to the reaction i (1123) — [0001] + 3(1120) THE microstrain c axis compression of single crystals of commercial-purity ingot SR beryllium (99.6 pct), high-purity twelve-zone-pass beryllium (99.98 pct), Be-5.24 pct Ni and Be-4.37 pct Cu alloys was described in a previous paper.1 This paper covers in detail the analysis of slip traces observed on two mutually perpendicular lateral surfaces of these specimens, and a detailed description of transmission electron microscopy studies performed on foils cut from the bulk crystals after they had been deformed to fracture in the c axis compression. Observation of slip traces on single surfaces of deformed single crystals are generally insufficient to positively identify slip or twinning modes. The use of two carefully cut and oriented perpendicular surfaces can greatly aid in the positive identification and index- ing of slip traces, although even this technique may be quite inadequate if more than one type of slip system operates and if an insufficient number of traces are observed on the surfaces. The problem is greatly simplified for symmetric cases like that for c axis compression of an hep crystal such as beryllium, in which the operating slip systems are all equally inclined to the direction of the applied stress, and each slip system of a given slip mode has an equal chance of operating. For such cases, the traces of any given slip mode observed on the surfaces cut parallel to the c axis are symmetrically tilted about the c axis. It is therefore possible to quickly determine whether one or more slip modes are operating. Confirmatory evidence in support of the observations made on the external surfaces can be obtained from foils cut from the deformed crystals and examined by transmission electron microscopy. This latter technique serves to identify not only the operating slip plane but also the Burgers vector of the dislocations which participate in the slip. For this purpose, a simplified technique based upon a double tetrahedron notation is used in the present paper. The planes and directions in the hep lattice are all designated by letters rather than indices and extinction conditions are easily determined if the Burgers vector lies in the plane contributing to the diffraction. RESULTS 1) Slip Trace Analysis. The standard (0001) stereo-graphic projection of beryllium is shown in Fig. 1. The two mutually perpendicular, lateral surfaces of the compression specimen are represented by the diametrical planes AA&apos; and BB&apos;, also referred to as surface A and surface B. For the specific case represented (a Be-5.24 pct Ni specimen deformed by c axis compression at room temperature), the A surface is tilted 5 deg to the (10i0&apos;) plane and the B surface is tilted 5 deg to the (1120) plane. Two surface trace analyses may be facilitated by examining in turn the intersection of various great circle traces of specific pyramidal planes with two surfaces and comparing the angles made with the (0001) plane with those actually observed on the two surfaces. One then identifies the slip traces by trial and error on a best-fit basis. The (1122) type planes (it was found that slip occurred on these planes) are shown plotted on the stereographic projection in Fig. 1. One obtains directly the angles between the (0001) plane and the {1122) traces by measuring the angle from the periphery to the point of intersection along the lines

Jan 1, 1969

By H. L. Hartman

The estimation of the magnitude of pressure changes which occur in mine ventilation circuits is of primary importance to the mining engineer in making changes in an existing mine or in projecting the ventilation requirements of a new mine. It is a formidable task in many cases, because no adequate theory is available and empirical data are unreliable or even lacking. The common practice of increasing one&apos;s estimate of the mine static pressure by an arbitrary "safety factor" to compensate for the unknown is a poor substitute for accuracy. Generally little difficulty is encountered in calculating friction losses in mine airways. As long as the friction factor can be reliably estimated, the well-known Atkinson formula gives accurate results. However, it is in the determination of shock losses that the ventilation engineer encounters difficulty and is apt to resort to an arbitrary allowance. Shock losses may constitute as little as five or as much as twenty-five percent or more of the overall pressure drop in the mine, and such allowances are generally completely unsatisfactory. Since other shock losses have been covered by McElroy (1), this paper deals with a long-neglected source of shock and pressure change in ventilation circuits, divergence and convergence of airflows, here referred to by their more common mining terms, splits and junctions. PREVIOUS WORK There is practically no fundamental information available in ventilation literature on the splitting and junction of airflows, either in mining or mechanical engineering periodicals or handbooks. Fluid mechanics textbooks seemingly avoid the subject of divergence and convergence, except for brief empirical treatment. The only investigation reported for mine ventilation systems was by Weeks, et a1 (21), in 1933. Although applicable to the case of line brattices in coal mines, it cannot be applied to the general case of mine splitting and junctions involving change of direction with little or no area change. In dealing briefly with the subject of shock losses at splits and junctions, McElroy (1) suggests only that they be considered comparable to bends, with losses at junctions being increased 50 pct to allow for interference of merging streams. He recognizes the importance of quantity ratio but concludes that exact quantitative data are lacking to compute losses at splits and junctions. A very capable summary and analysis of the state of knowledge in industrial ventilation and fluid mechanics regarding losses due to divergence has recently been reported by Gilman (3). In considering all earlier work dealing with both air and water as fluids, he compares results on a common basis and obtains good empirical agreement. For divergent flow, the controlling variable in each branch is the quantity ratio, although the deflection angle is of secondary importance. Empirical equations are presented for straight and 90-deg branches with varying velocity ratios, assuming branches smaller than the main duct. These formulas are modifications of the Borda equation for abrupt contraction and expansion. The only information available on shock losses in convergent flow are approximate experimental data presented by Alden (4) and the Manual of Recommended Practice for Industrial Ventilation of the American Conference of Governmental Industrial Hygienists. Alden considers the effects of both quantity ratio and deflection angle. The data indicate an effect opposite to that observed in divergence, the shock loss varying directly with quantity ratio. SHOCK LOSS THEORY Because pressure losses due to shock have been found to bear a constant relation to the mean velocity of flow in a given conduit, they are frequently expressed as a dimensionless function of the velocity head, termed the shock loss factor. For airflow in ducts and mine airways, the shock loss may be represented by the formula (1): Hx - XHv, (1)

Jan 1, 1961

By J. H. Polhems, R. B. Brackin, D. B. Grove

THE heavy media separation process plays an outstanding role in the concentration of 4000 tons of zinc ore per day at the Mascot mill of the American Zinc Co. of Tennessee. Of the total tonnage, 72 pct is treated in the heavy media separation plant to reject 56 pct of the ore as a coarse tailing, which has a ready market. Concentrates from this separation are beneficiated further by jigging and flotation. Approximately 25 pct of the total zinc concentrate production is made in the jig mill. Jig tailings are ground and pumped to the flotation circuit where the balance of the production is made. Fig. 1 shows a generalized flowsheet of the mill. The Mascot ore is a lead-free, honey-colored sphalerite in dolomitic limestone, with lesser amounts of chert and some pyrite. A mineralogical analysis is given in Table I. After 10 years of successful operation with galena medium and treatment of nearly 10,000,000 tons of ore, a decision to convert to ferrosilicon was made early in 1948 because of the increasing price of galena and consequent high operating costs. The conversion was made on Nov. 6, 1948, and the results obtained since that time have shown remarkable improvement over those made with galena. The Table I. Mineralogical Analysis of Mill Feed, Pct Calcium carbonate 49.5 Magnesium carbonate 35.2 Iron oxide and aluminum oxide 1.5 Zinc sulphide 4.5 Insoluble 9.3 100.0 Table II. Comparative Data, Galena and Ferrosilicon Ferro- Diner-Gelenaa siliconb ence Operating costs per ton milled, ct. 21.21 9.12 12.09 Medium consumption per ton milled, lb 0.80c 0.15 0.65 Reagent consumption per ton milled, lb 0.45 0.02 0.43 Tailing assay, pct Zn 0.310 0.297 0.013 Concentrate. oct Zn 12.08 10.33 1.75 Heavy medla ieparatlon recovery. pct 89.38 90.22 0.84 Mill feed rate, tons per hr 153 166 13 Heavy mesa separation feed rate. tons per hr 100 10 0 Tons milled per heavy media separation man shift 350 620 270 Mill feed to coarse tailings, pct 51.0 56.7 5.7 Lost mill time, pct 5.6 5.0 0.6 Power consumption, kw-hr per ton 2.06 1.92 0.14 a 1947. " First 6 months of 1950. c Net consumption after deducting credit for reclaimed waste galena. Consumption of new galena was 1.320 lb per ton milled. For entire life of galena operation, a credit of 40 pct of the value of the new galena added was realized from the sale of waste galena. comparisons given in this report cover the first 6 months of 1950 as representing the ferrosilicon operation, and the year 1947 as representing the galena operation. This was the last full year in which galena was used exclusively and is representative of the best work done during the 10 years of operation with this medium. After only 2 years&apos; operating experience, with ferrosilicon and treatment of 1,807,585 tons many advantages have been revealed and are summarized in Table 11. Development Prior to the introduction of the heavy media process, all the mill feed was crushed through 5/8 in. and treated by jigging. A finished tailing assaying 0.66 pct Zn was made on rougher bull jigs, and cleaner jig tailings were ground for treatment by flotation. The first test work on the sink-and-float method of mineral beneficiation was carried out at Mascot in 1935, using a 3-ft cone and galena medium for batch tests. The following year a 6-ft cone was installed for pilot-plant work. This unit became a part of the mill circuit on March 1, 1936, and handled a gradually increasing tonnage in the next 2 years as the process developed to the point where it could treat all the + 3/8-in. material in the mill feed. Coarse jigging was then discontinued on March 1, 1939, and all coarse tailings have been made by the heavy media separation plant since that time. Feed Preparation: The original feed preparation plant consisted of a drag washer followed by two 4x10-ft Allis-Chalmers washing screens. A surge bin and two additional 5x12-ft AC washing screens were added in 1943. Use of primary and secondary washing screens was found essential to provide the cleanest possible feed for the cone and thereby avoid excessive contamination of the galena medium. Improved washing was obtained by replacing the drag washer with a 7x20-ft Allis-Chalmers scrubber, shown in Fig. 2, which has been in service since May 1944. Throughout the life of the galena operation, delivery of extremely muddy ore to the mill overloaded the medium cleaning system, and it frequently was necessary to cut off the feed and clean the medium for several hours until its normal viscosity had been re-established. The cleaning circuit

Jan 1, 1952

By J. K. Elliott, P. H. Kelly

Many engineering functions such as surface metering work and laboratory compressibility check points involve the use of liquid densities of light hydrocarbon mixtures at various pressures and temperatures. However, at the present time, no simple reliable method exists for determining density variation, particularly if the composition of the liquid is unknown. Consequently, a study was undertaken to develop and present a simple and accurate method of predicting density variation of a light hydrocarbon liquid with pressure and temperature, knowing only the density of the liquid at some condition. The experimental liquid compressibility data from API Project 37 by Sage and Lacey&apos; have been considered to be accurate within 0.5 per cent and cover a wide range of pressure (14.7 to 10,000 psia), temperature (100" to 400°F) and molecular weight (up to 150). From these data, a set of liquid density curves, which relate density to pressure, temperature and molecular weight, was developed. These curves make it possible to predict density variation with pressure and temperature. Compared to extensive laboratory compressibility data on a complex, light hydrocarbon liquid, the use of the liquid density curves resulted in an average error of less than 0.5 per cent. Based on the results of this analysis, it is concluded that the set of liquid density curves developed from the data of Sage and Lacey provides an accurate and simple method for predicting the density variation of light hydrocarbon liquids when the density at some condition is known. These curves should be very helpful in many engineering calculations, particularly in the surface metering of light hydrocarbon liquids. INTRODUCTION Many situations arise in field and engineering laboratory work, such as reservoir engineering studies, check of experimentally determined laboratory data and orifice flow-meter formulas, where liquid density factors at various pressure-temperature conditions are required. Also, the need for accurate light hydrocarbon liquid information has become more important with the advent of miscible-type displacements for secondary recovery purposes in oilfield operations. Several reliable methods are available1 - "or determining the density of liquid hydrocarbons if the composition of the liquid is known. However, there is a definite lack of methods for accurately determining the variation of density when the composition of the liquid is unknown. The purpose of this study is to review the various methods for determining hydrocarbon liquid densities and to develop a simple and reliable method of determining variation in density of light hydrocarbon liquids with pressure and temperature when the compositio~n of the liquid is unknown. METHODS FOR DETERMINING DENSITY OF LIQUIDS OF KNOWN COMPOSITION Sage, Lacey and Hicks&apos; have proposed a method to predict the density of light liquid hydrocarbons by using partial molal volumes. Data are available on experimentally developed partial liquid volumes of hydrocarbons over a rather limited range of temperature, pressure and composition. The partial mold volume method has proved satisfactory for determining the density of some hydrocarbon liquids when the composition is known. Within the range covered in the Sage, Lacey and Hicks1 data, the results agree within about 3 per cent of the experimental values. Hanson mentions the limitation of this method to a composition range of approximately 10 per cent by weight of methane, which will not allow this correction to cover most low molecular weight-light hydrocarbon liquids. Standing and Katz2 studied data on light hydrocarbon-liquid systems containing methane and ethane at high temperature and pressure and have presented a method for determining liquid densities, assuming additive volumes for all components less volatile than ethane and using apparent densities for methane and ethane. The compressibility and thermal-expansion curves used by Standing are based on assumptions that compressibility of a hydrocarbon liquid at temperatures below 300°F is a function of the liquid density at 60°F and that thermal expansion of the liquid is affected little by pressure. The information required to use this technique with an example problem is furnished by Standing.&apos; Hanson eports an average error of - 0.5 per cent using the method of apparent densities in calculating liquid densities of several volatile hydrocarbon mixtures. However, as implied, the apparent density method is not applicable for liquid density calculations when the composition of the liquid is unknown. Watson- as presented a method