Search Documents

By A. S. Yue, A. G. Guy

A study was made of the effects of a prolonged flux of zinc atoms through the a solid solution of zinc in copper. The experimental arrangement consisted essentially of a copper disk about 0.01 in. thick, at one of whose surfaces a gaseous atmosphere containing zinc atoms was maintained, and at the other surface a gaseous atmosphere with a minimum of zinc atoms was maintained. During prolonged theothersurfaceexposure at high temperatures the zinc content of the copper disk gradually built up to the steady-state concentration distribution and then remained at this value. The concentration-distribution curves for various conditions were determined by chemical analyses. The results showed that the condition of steady-state diffusion was achieved. The diffusion coefficients calculated from the experimental data, although not of high precision, agreed with the values obtained by other workers using unsteady-state methods. Relatively slight porosity developed in the specimens in the course of diffusion. A LTHOUGH most diffusion studies have been made under unsteady-state conditions, it is known' that the steady-state method is often superior with respect to the directness and accuracy of interpretation of the data. Steady-state diffusion of gases through metal diaphragms is well known. Also, Harris' and Smith" have used the steady-state method in studying the diffusion of carbon in aus-tenite. The accepted mechanism in this system involves the motion of the interstitial carbon atoms in the rigid framework of the lattice of iron atoms. Thus, there is little difficulty in visualizing the steady flow of the small carbon atoms through the austenite. The situation in substitutional diffusion is quite different. Here the atoms are comparable in size, and it is not evident how a steady flow of one of the atoms through the solid solution might be achieved. At the time the present research was started, it was known that a previous exploratory attempt to produce steady-state diffusion in a substitutional alloy, the Au-Ag system,' had been unsuccessful and had indicated that perhaps there were basic difficulties that could not be ovei-come. Therefore, the present research began as a study of the effect of a prolonged flux of metal atoms through a substitutional solid solution. Eventua.lly, it was possible to produce actual steady-state diffusion in the system chosen for study, the a Cu-Zn alloys. Experimental Procedure The aim in the experiments was to maintain a high zinc content, about 30 pet, at one surface of a copper sheet, and to maintain a low zinc content, near 0 pet, at the opposite surface. The zinc would then diffuse into and through the copper, first building up to the steady-state concentration distribution and then maintaining this distribution. The three types of specimens that were used are shown in Fig. 1. In type A specimens the copper disk through which diffusion occurred was welded to the top of a cylindrical molybdenum tube, the bottom of which also was sealed by welding. At the diffusion temperature the brass chips in the molybdenum container were the source of the zinc vapor which maintained the lower surface of the copper disk at 30 pet Zn. The upper surface was maintained at 0 pet Zn by the vacuum in which type A, and also type B, specimens were diffused. Since the molybdenum container was impervious to zinc vapor, it was intended that the only path of escape for the vapor from the brass chips would be through the thin copper diffusion disk. However, it was found that small leaks often developed at the welded joints during the diffusion treatment, and in most specimens some of the zinc was lost in this manner. Although even small losses of this kind were a serious handicap in attempting to determine the flux through the disk, they did not prevent the maintenance of satisfactory boundary conditions for the attainment of the steady-state condition. Type B specimens differed from type A in having a weight of about 300 g supported on the copper disk by 15 to 20 short quartz rods. This change was made when it was observed that the copper disk was being bowed upward by the difference in the pressures acting on its two surfaces. Since the grain-boundary cracking which occurred in the bowed specimens could be attributed largely to the accompanying creep,3 it was desirable to minimize this effect. The counterweight was effective in significantly decreasing both bowing of the disk and cracking at grain boundaries. Type C specimens differed considerably from the others in that the low-zinc atmosphere at one sur-

Jan 1, 1959

By G. H. Turner, R. C. Bell, E. Peters

Zinc oxide activities in a typical lead blast furnace slag have been calculated from plant operating data. These activities were used to assess the probable effect of fuel composition, oxygen enrichment, and air preheating on the efficiency and capacity of the slag-fuming operation. THE physical chemistry of zinc fuming has been examined with three objectives in mind: 1—to predict conditions favorable to increasing furnace capacity, 2—to predict the changes required to fume zinc more economically, and 3—to explain reported differences in the efficiencies of various slag-fuming plants. This study, made at ail in the plants and laboratories of The Consolidated Mining and Smelting Co. of Canada Ltd., developed from a program undertaken some three years ago on behalf of the AIME Extractive Metallurgy Div. subcommittee on slag fuming. Lead metallurgists first became interested in the recovery of zinc from lead blast furnace slags in 1905 and 1906. An excellent review of the early experimental work has been made by Courtney,' who described blast furnace, reverberatory furnace, and converter methods of fuming zinc from slag. Some of the investigators did not appreciate the importance of reducing the zinc oxide content of the slag to metal in order to fume it, since they tried compressed air blast without fuel in their earliest attempts. However, by 1908, the importance of reducing the zinc was established.' In 1925, the Waelz process for the recovery of zinc oxide from oxidized zinc ores was developed in Germany.' This process was not readily adaptable to lead blast furnace slags because of the difficulty in handling fusible charges in a kiln. What appears to have been the first slag-fuming operation as it is known was commenced by the Anaconda Copper Mining Co. at East Helena, Mont. in 1927." The first Trail furnace was completed in 1930, and this was followed by the construction of several other slag-fuming plants. During the period in which slag fuming has been extensively employed, little development of the chemistry of this process as a whole has taken place. Several good papers on the petrography of lead blast furnace slags have been published,""= but these studies could do little more than establish the forms in which lead and zinc occur in the initial charge and final products of the slag-fuming operation. In recent years, zinc-smelting problems have been ap- proached from a thermodynamic point of view. Maier has published an excellent thermodynamic treatment of zinc smelting." The important thermodynamic properties of zinc and its compounds have been determined and checked by other investigators.' However, to the best of the authors' knowledge, no thermodynamic treatment of the fuming of zinc from slag has been published. A thermodynamic study of any process requires that the essential chemistry of that process be known. In slag fuming there appear to be some differences of opinion as to whether the active reducing agent is elemental carbon or carbon monoxide. Furthermore, some observers have noted that high volatile coals appear to be more efficient than low volatile coals, indicating that hydrogen is also an important factor in the reducing efficiency of a fuel. That both hydrogen and carbon monoxide are effective reducing agents for the zinc oxide content of lead blast furnace slags can be demonstrated readily by introducing these gases into a slag bath held in a neutral vessel at 2100°F (1150°C). Elemental carbon also will reduce zinc oxide, but it is improbable that much free carbon is available for reduction of zinc, as the reaction between the finely powdered coal and air should be largely completed before the solid coal particles reach the slag. Some large-scale fuming experiments using gaseous hydrocarbons have been carried out by other investigators, but, as far as is known, these have not been developed yet into operating processes. The thermodynamic treatment in this paper is based on the following reactions: 1—to supply the thermal requirements C+V2O2- CO [1] C + 0,-CO, [2] H2+ ~z0,-H,O 131 and 2—to reduce ZnO ZnO + CO + Zn + CO, c41 ZnO + H, e Zn + H,O. [51 The furnace-gas composition also is controlled by the equilibrium constant of the familiar water-gas reaction H,O + CO + CO, + H2. C6l In order for the thermodynamic calculations to be quantitatively applicable, it is necessary that the chemical reactions to which they are being applied

Jan 1, 1956

By W. C. Hauber

Techniques previously published to predict the production performance of a water flood when the mobility ratio is not unity have been primarily restricted to a five-spot well pattern. When other types of well patterns are considered, the prediction of performance has required either a model study or a complex numerical solution of a multiple-fluid system with a high-speed digital computer. Simple methods are developed in this report for predicting the areal sweep efficiency and injectivity history of a water flood for arbitrary well patterns and mobility ratios. The effect of an initial gas saturation can be incorporated into this method of prediction. With the layer superposition concepts published by Prats, et al.,1 these methods can be used to predict the production performance of a water flood in a reservoir having a wide range in permeability. A reasonably good agreement is shown by comparisons of the results calculated by this method with those obtained either by model studies or by complex numerical solution of multiple-fluid systems. Therefore, it has been concluded that a reasonably accurate method has been developed for predicting the production performance of a water flood employing an arbitrary type of well pattern for any given mobility ratio. A medium-speed digital computer can be used in applying this method. INTRODUCTION Production performance of multi-fluid five-spot water floods can be predicted by the techniques published by Prats, et al.1 This technique takes into account the effects of initial gas saturation, mobility ratio and vertical variations in horizontal permeability. A basic assumption in applying this technique is that unit displacement efficiency does not vary with time after the passage of the flood front. Prats, et al. demonstrated that this assumption is valid for normal waterflood applications with small mobility ratios. To extend Prats' technique to other types of well patterns has required either an experimental model study or a numerical solution of a multiple-fluid system to describe the performance of an individual homogeneous layer. Sheldon and Dougherty2 presented an exact numerical solution employing any arbitrary well pattern. This technique is complex to apply, requiring (1) numerical solution of the single-fluid system, (2) conformal mapping of the streamlines and iso-potential lines to form a new grid network, and (3) a numerical solution of the multi-fluid system on the new grid network. A method was developed in this paper to approximate areal sweep efficiency and injectivity as functions of time for an individual homogeneous layer for arbitrary well patterns and mobility ratios. Analytic expressions were derived to directly apply this method to five-spot and direct line drive well patterns. For other types of well patterns, this method is similar to the approach of Higgids and Leighton,6 in that it utilizes the stream tubes and iso-potential lines obtained by the solution of differential equations when the mobility ratio is unity. However, this method is easier to apply than the Higgins-Leighton method since it does not require knowledge of the variation between oil and water relative permeabilities and saturation. Using the relationship between areal sweep efficiency, injectivity and the volume of water injected into an individual homogeneous layer as obtained by the method presented in this paper, and Prats' technique for the superposition of individual layers, one can predict the production performance of any water flood. ASSUMPTIONS GENERAL Basic assumptions made in this report are that fluids are incompressible and the layers are horizontal, homogeneous and uniformly thick. It is assumed that as water is injected into the layer, three distinct regions are formed —a water bank, an oil bank, and an unflooded region. Each region is assumed to displace other regions in a pis-tonlike manner, with only one fluid flowing in each region, so that there is no change in unit displacement efficiency with time after the flood front has passed in the reservoir. Therefore, the residual oil saturation in the water bank and the oil saturations in the unflooded region are constant throughout the layer. The ratio (F) of the bulk volume of a layer occupied by the oil and water banks at any time before oil breakthrough, to the corresponding bulk volume of a layer occupied by the water bank, is given' by the equation BREAKTHROUGH AREAL SWEEP EFFICIENCY To predict the behavior of any type of waterflood pattern, one must know the breakthrough areal sweep efficiency Eb(Mc,o, F) for the pattern for any water:oil mo-

Jan 1, 1965

By Ralph Hultgren

IN solid solutions atoms of differing sizes occupy the same crystalline lattice, requiring that some of them be compressed and others expanded. The energy involved has been called misfit strain energy and is an important concept of crystal chemistry. If the atomic sizes and elastic constants of interatomic bonds are known, the misfit energy may be calculated,' provided certain simplifying assumptions are allowable. Usually, isotropic crystals are assumed and interatomic distances are taken to be the statistical average determined from X-ray diffraction. Such calculations yield values of the misfit energy of the order of 1 or 2 kcal per atom in alloys such as Au-Cu at compositions of 50 atomic pct. However, evidence has accumulated in recent times that atoms change their sizes with composition of alloys, implying electronic rearrangement of the bonds. The size changes have been found particularly by application of the X-ray method developed by Warren, Averbach, and Roberts.' Thus, Averbach, Flinn, and Cohen3 determined radii in Au-Cu alloys. Oriani4 showed that these new radii led to a calculated misfit energy in disordered AuCu, which was decreased from the values calculated by the usual theory more than twenty-fold, to only 80 cal per g atom. Thermodynamic calculations from the phase diagram5 also show misfit energy to be no more than a few hundred calories per g atom in this alloy. The question of what electronic rearrangements are possible therefore becomes compelling in estimating misfit energy. In the following pages the results of certain calculations on the AuCu tetragonal superlattice are submitted. Conclusions drawn from these should be applicable in large degree to disordered solid solutions. As in all ordered states, bonding distances in the superlattice are individually known, rather than being merely average distances as found from lattice constants of disordered states. Moreover, only the Au-Au and Cu-Cu distances are strained; the elastic constants of these are known in the elementary state. In the usual calculation it is necessary to assume elastic constants for Au-Cu bonds. Misfit energy has thus been calculable without the need of many simplifying assumptions usually made. It is still assumed that equilibrium bond lengths and elastic properties of the bonds are the same in the alloy as in the pure metals. As previously discussed, this is probably not correct. Also assumed is that the bonds are not affected by strain of neighboring bonds. A calculation of Young's modulus from compressibility data shows this to be far from true; extensive electronic rearrangements take place. It would seem that misfit energy cannot be calculated from elasticity data for the elements. The usual methods may, however, give an upper limit which is often much higher than the true value. The question of electronic rearrangement is, of course, a complex one. Pauling's theory gives a simple, approximate treatment of the relation between type of bond and bond distance. This has been applied with some success to the Au-Cu system, as will be shown in a later section. Misfit Energy in Au-Cu Alloys Hume-Rothery and Raynor6 discuss the Au-CU system as a type example of strain energy. The gold atom is 12.8 pct larger in diameter than the copper atom, near the size factor limit beyond which solid solubility is severely restricted. They therefore consider the misfit energy to be large, a conclusion for which they believe they find evidence in the phase diagram. Gold and copper are completely miscible in the solid state, but the alloy has a minimum melting point at an intermediate composition. From this Hume-Rothery and Raynor conclude that the strain energy is nearly large enough to prevent miscibility; the phase diagram tends toward a eutec-tic type. In Ag-Cu, which has almost identical size relationships, solid miscibility is quite limited; whereas in Au-Ag, where atomic sizes are nearly the same, there is complete miscibility without a minimum in the melting point. From their arguments the heat of formation of Au-Cu would be expected to be endothermic or only slightly exothermic, that of Ag-Cu to be endothermic, and that of Au-Ag to be exothermic. Deviations, from Ve-gard's law of additivity of atomic radii support these conclusions, since Au-Cu and Ag-Cu both have pronounced positive deviations, and Au-Ag has a negative deviation. Nevertheless, Au-Cu alloys form exothermically; indeed, considerably more exothermically than Au-Ag, Table I. Hence, strain energy must be much less important in this case than Hume-Rothery and Raynor have supposed.

Jan 1, 1958

By F. Claisse, R. Angers

Self-diffusion in iron has been measured during rapid a-r transformations using a variant of the Kryukou and Zhukhovitskii diffusion method. The study was performed by thermally cycling iron foils (1 to 6 cpm) through the transformation (=910°C). Some foils have been subjected to over 1000 cycles and some have spent more than 15 pct of their total diffusion time in the process of transformation. The experimental results show that the a-r transformation has no measurable effect on self-diffusion in iron. The study is completed by a quantitative analysis of mechanisms which can affect the diffusion rate during the transformation. The analysis confirms the experimental results. SINCE diffusion is an important factor in many solid-state transformations, it is of interest to study how it is affected by the stresses generated during these transformations. Clinard and Sherby1'2 were the first to make a study along these lines. They measured diffusion coefficients in Fe-FeCoV couples subjected to slow thermal cycling (1.5 cph) through the a-r transformation range. They found an enhancement of diffusion by a factor of about two. The purpose of the present paper is to report measurements of the self-diffusion coefficient of iron during much more rapid thermal cyclings (1 to 6 cpm) through the a-r transformation (-910°C). These more rapid cyclings produce higher strain rates during the transformation and should emphasize any possible influence of transformation upon diffusion. EXPERIMENTAL Iron foils, 25 to 35 µ thick, were cold-rolled from 99.92 pct pure iron and annealed in pure helium for 2 hr at 870°C; the resulting grain diameter was about 150 µ. Specimens 0.5 by 8 cm were cut from the foils and I7e55 was vapor deposited on one of their surfaces. A 38 gage alumel-chrome1 thermocouple was spot welded in the middle of one of the specimen long edges, Fig. 1. Two 38 gage chrome1 wires were also spot welded along the same edge on each side of the thermocouple; they were placed 2.5 cm apart and used for electrical resistance measurements. In order to prevent twisting and crumpling, the specimens were pinched between two quartz plates 0.1 by 1 by 7 cm and the assembly was close fitted into a 1 cm ID quartz tube. Four holes were drilled through the tube to let the 38 gage wires out: these were connected to the recording equipment by means of extension wires. 20-gage nickel wires fixed at both ends of the specimens were used to thermally cycle the foils by Joule heating. The above described device was placed in a 2.7 cm ID quartz tube which in turn was placed in a tubular furnace. Either a pure helium atmosphere or circulating hydrogen was used during the experiments. Specimens were subjected to thermal cycles between a minimum temperature To and a maximum temperature Tm at rates ranging from 1 to 6 cpm. This was obtained by maintaining the furnace at a constant temperature near the minimum temperature To and periodically passing an electric current through the specimen. Cooling was achieved by heat losses to the surroundings. The forms and periods of cycles were varied from one specimen to another; however, each specimen was subjected to one type of cycle only. The temperature and electrical resistance variations of the specimens were recorded as a function of time. The temperature curves were used for diffusion calculations while the electrical resistance curves were used to monitor the transformation and to determine its starting point and its approximate duration. Diffusion was measured by the method developed by Kryukov and zhukhovitskii3 and modified by Angers and Claisse.4,5 In this method a metallic foil is coated on one side with a radioactive isotope and the activity is measured periodically on both sides during the diffusion anneal. The following equation then holds: where: I1 Activity on the surface on which the deposit is made. I, Activity on the opposite surface. t Diffusion time. B Constant. D Diffusion coefficient. d Foil thickness including the deposit. G(t) A function of time; it is a second order correction term which is given graphically in Refs. 4 and 5. The diffusion coefficient D is found by plotting ln[(Il - I2)/(I1 + I,)] -G(t) against t; the resulting slope m leads to an accurate calculation of D through Eq. [2]. The effect of the a-r transformation on diffusion is expressed by the ratio DT/DU where:

Jan 1, 1970

By W. Hurst, R. E. Leeser, W. C. Goodson

Three aspects of gas deliverability are presented in this paper. The first treats with the gas deliverability or availability of a normal depletion-type dry gas field. Such encompasses not only the period of stabilized constant rate, but more so, the "tailings" when a fixed abandonment pressure is reached and the rate by necessity must decline. A comprehensive work plot is offered, developed from mathematics herein included, that removes the triai-and-errnr computations that attended such undertakings in the past. The second part treats with the discount factor of the open flow potential constant from what is observed initially in testing a gas well to what is evidenced when stabilization is reached. This prevails in tight formations, such as the Kansas Hugoton field which is offered as the example. The means of establishing this factor are pressure build-up curves which, as sustained by analytical deductions, reproduce this entire period of transient flow under conditions of a constant rate influx. Finally, what is offered in this paper is the deliverability performance of an exceedingly rich gas condensate field producing from a tight formation. The example shown is the Knox Bromide field in Oklahoma, producing from the Bromide formations. The results are ominous, showing early reduction in permeability to gas pow, due to the retrograde condensate forming in the pore space, with the attending early logging-up of these wells. The analytics of lowered permeability are incorporated in the gas deliverability formula along with the PVT data that gives the increased condensate liquid saturation as the gas flows to the well bore. This paper would not be complete without a critique oflered at the end. With the many gas wells now in production and those that have completed their life, there has been no factual information collected by any source as to what constitutes that permeability range where a gas well would be unimpaired in its gas deliverability by the presence of rich condensate content, and the lowered range where such would be harmful. This question confronts all producers. INTRODUCTION Various aspects of gas deliverability are presented in this paper that includes depletion-type reservoirs, deteriora- tion factor of the gas deliverability constant, and the performance of a rich gas condensate reservoir producing from a tight sand. With respect to the presentation of gas deliverability and its tailings for depletion-type gas reservoirs, one notes that this is essentially the information offered by every gas transmission company and producer appearing before the Federal Power Commission for Letters of Conveyance in the dedication of reserves. In the ordinary procedure, as many engage upon this study, trial-and-error calculations are included, particularly as apply to the tailings. For many years one of the writers has employed mathematical analyses to perform this step and avoid the complexities so associated. In the preparation of this paper these analyses have been amplified to include any slope n for the open flow potential relationship for which the tailings can be determined from Fig. 1. With reference to the deterioration or discount factor of the open flow potential constant as such occurs in the gas deliverability formula, this for the most part has been an unexplored subject. Although the issue first appeared in the Kansas Hugoton field, where such was surmised but only recently resolved, this situation of a deterioration of the gas deliverability constant can occur wherever dry gas production from a tight sand is encountered. The first concerted attacks upon this problem were the presentations of Hurst' and Goodson' before the Kansas Corporation Commission to show that transient fluid flow and unsteady-state flow formulas prevailed. This was amplified later before the Federal Power Commission3 to show that this deterioration factor could be identified from pressure build-up curves. This has been reported by McMahon.4 Its importance to the industry merits the review of these essential features in completing the program on the aspects of gas deliverability. Finally, as illustrated here, for a low permeability formation such as the Knox Bromide field where the gas is rich, representing some 165 bbl of condensate per MMcf of effluent gas, the gas deliverability can be of limited extent in the life of the field, leaving substantial amounts of condensate and gas unrecovered. In cases such as this, gas cycling is mandatory. This is particularly revealed by the fluid mechanics introduced here, employing factual field as well as laboratory data, to show this-restriction upon gas deliverability. PRESSURE DEPLETION What will now be offered is the study of gas deliver-

By D. L. Feucht, R. L. Longini

The semiconductor heterojunction is considered in terms of simple models which may lead to an understanding of move complex heterojunctions. Metallurgical and electrical properties of hetero-junctions aye discussed including the interface structure, energy -band diagram, and carrier transbovt across the interface. It is found that in a heterojunction all mechanisms such as injection, tunneling, and junction recombination found in simple junctions play modified voles. INTERFACES between materials (grain boundaries, the electrical junction between two differently doped materials in a single crystal, the oxide-metal interface, or metal-metal junctions) are of considerable importance in many situations. These various interfaces all have one very fundamental thing in common. Quantum mechanically speaking, the wave functions of the electrons in one material may penetrate the other material but, in general, only to the extent of angstroms. From an electrical point of view the conduction mechanism changes as a current passes through such junctions. In some cases the change is tremendous, in others almost negligible. The interface, then, is the locus of a change of conduction mechanisms. Some of these, particularly in semiconductors, are well-understood. The ordinary p-n junction in a single crystal can be the locus of an injection mechanism or a tunneling process, depending on conditions. The mechanisms are probably best understood in semiconductors because of the possible simplified view of particlelike conduction. The bands are either nearly filled or nearly empty and band overlap is seldom involved. The same fundamentals are probably important in other situations too but they are very difficult to look at naively. Although the simple look at the semiconductor case only gives us a relatively rough picture which must then be refined, the other systems, which involve a more complex situation, immediately are in many ways too difficult. There are too many initial choices of complex systems and therefore it is not possible to be even reasonably certain of any one model. Because of the relative simplicity of semiconductors, their good and controllable structure, and because of the ability to make many measurements on them not normally available to either metals or insulators! they are probably the best understood materials. It is therefore desirable to use them as a tool to further the understanding of interfaces in general. Semiconductor-heterojunction concepts were first proposed by kroemer1 in 1957. This was followed several years later by reports on the fabrication and experimental characteristics of heterojunction structures by Anderson2 and Diedrich and jotten.3 I) THE HETEROJUNCTION STRUCTURE To get down to hardware, when we refer to a semiconductor heterojunction we imply that there exists an intimate contact between different semiconductor materials. We could put two pieces of material together, complete with oxide layers, we could remove the oxides, or we could even melt the interface and hopefully get wetting and a good "bond" on solidifying. In fact we could by some means grow a crystal of one material using the other as a seed. Essentially we are interested only in the last two because they are the simplest to look at analytically. The degree of perfection of fit varies greatly and is reflected somewhat in the arc welder's joint strength. The lattice match of the two materials, their orientation, and so forth. is obviously necessary for a good bond but so is the continuity of any polar bonds which are involved such as in the III-V semiconductors. The mechanical misfit between two similar lattices can be described in terms of edge dislocations. The edge-type dislocations must be very close together for the usual misfit and there must be dislocations for each of several different Burger's vectors in order to produce a lattice match. The .'dangling bonds'' resulting will be involved in producing interface charge. Order of magnitude estimates of the charge density extrapolated from low densities of dislocations in homogeneous materials give 5 x 1013 cm-2 Ge-Si and 1 X 1012 cm-2 Ge-GaAs electronic charges. Edge dislocations also act as very active recombination centers between holes and electrons. One lattice "matching" difficulty usually exists even if two structures have essentially the same lattice constants as they will have different coefficients of therma1 expansion. Thus, on cooling from the usually high temperature of fabrication to room temperature, dislocations are produced, a good fit not existing at both temperatures. In brittle materials this shrinkage may even result in cracking. For the Ge-Si interface the mismatch is about 2 x 10 -6 per degree whereas it is less than 10"7 per degree between germanium and GaAs. The exact effect of the misfit is dependent on the thickness of the materials involved. For a very

Jan 1, 1965

By H. C. Rogers

A phenomenological study of the failure of polycry stalline ductile metals at room temperature was carried out using light and electron microscopy. Tensile fractures as well as sections of partially fractured bars of OFHC copper in particular were examined. The initiation and growth of the central crack in the neck of a tensile specimen occurs by void formation. After the formation of the central crack the f'racture may be completed in either of two ways: by further void formation or by an "allernating slip" mechanism. The first leads to a "cup-cone" failure; the second, to a "double-cup" failure. In the past decade or decade and a half there has been a great deal of emphasis on the solution of the problem of the brittle fracture of metals, particularly those which normally exhibit considerable ductility such as steel. Since the problem of the fracture of metals after large plastic strains has less immediate commercial or defense significance, there has been considerably less effort expended in describing the details of the phenomenology and determining the mechanism of this type of fracture. The present research was undertaken to increase our knowledge in this area. The problem of ductile fracture has not been neglected completely, however. Ludwik1 first found by sectioning a necked but unbroken tensile specimen of aluminum that fracture began with a large internal crack which appeared to have started in the center of the neck. Examination of the fracture indicated that the crack had propagated radially with increasing deformation until a point was reached at which the path of the fracture suddenly left this transverse plane and proceeded at approximately 45 deg to the stress axis until the surface was reached. This gives rise to the commonly observed cup-cone tensile fracture. When MacGregor2 was attempting to demonstrate the linearity of the true stress-true strain curve from necking until fracture, he found that copper was anomalous in that the stress dropped off markedly from the straight line value before fracture occurred. Radiography indicated that in the copper an internal crack was formed long before the final fracture, the stress decreasing during the growth of this crack. One of the most significant advances in the understanding of ductile fracture was the result of work by Parker, Flanigan, and Davis.3 By the use of etch-pit orientations they were able to demonstrate conclusively that the fracture surface at the bottom of the cup, although on a gross scale normal to the tensile axis, did not consist of cleavage facets as had been previously supposed by many investigators. Recently, Forscher4 has shown evidence of porosity near the tensile fracture of hydrogenated zirconium which he attributes to hydride decomposition. The workers at the Titanium Metallurgical Laboratory5 have also shown evidence of porosity in a number of the commonly used metals after heavy deformation. Many metals have relatively low ductility during creep tests at high temperature. The fractures are intercrystalline, resulting from the nucleation and growth of grain boundary voids. The work in this area has been recently reviewed by Davies and Dennison.6 It is possible that some of the observations and conclusions may have a bearing on the present study? especially since at least two studies7,' have been extended down to room temperature and below using magnesium alloys. However, since magnesium does exhibit low-temperature cleavage, these results may not be pertinent to the present one. The use of the electron microscope as an aid to the study of fractures has been extensively exploited by Crussard and coworkers.9 The examination of direct carbon replicas of the fractures of a large number of metals and alloys showed that the bulk of the fracture surface was covered with cup-like indentations of the order of 1 to 2 µ in size. These frequently had a directionality by which Crussard claims to be able to tell the direction of the crack propagation. With this rather disconnected background of information, this investigation was undertaken in the hope of presenting a unified picture of the initiation and propagation of a fracture in a ductile metal. To this end all of the techniques previously used were employed simultaneously so that there might be a good correlation of the data obtained by different techniques. EXPERIMENTAL PROCEDURE The metal which was chosen as the starting material for this investigation was OFHC copper. Of the dozen or so materials considered, it best fulfilled the requirements of commercial availability in large sizes, good ductility, relatively high melting point compared with room temperature and

Jan 1, 1961

By W. Allen

FOR the past six years gravity surveys have been used for underground prospecting in the copper mines at Bisbee, Ariz. The primary purpose of the surveys has been to reduce the diamond drilling and crosscutting necessary for exploration. Since many of the orebodies are small, and geologic control is not always apparent, any information that will direct the drilling and crosscutting is highly desirable. Because of extensive development and exploration work in the copper mines at Bisbee, it has been possible to cover more than 630,000 ft of crosscuts on 30 levels with the gravity surveys. In the process the gravity procedures have been refined to a high degree. Density Contrast: For a gravity survey to be successful, a sufficient density contrast must exist between the geologic feature sought and surrounding host rocks. Most mineralized areas will provide this contrast if fairly massive bodies are present. In the Bisbee area the entire sequence of formations, except for alluvium, appears to have specific gravities ranging from 2.65 to 2.70. These values have been determined by means of a large number of cut samples and diamond drill cores. As a further check, vertical gravity differences have been used where nonmineralized sections are known to occur.' The only known major gravity disturbances result from mineralization that has increased the density and the voids that have decreased density. The voids are caused by mining operations and by underground water movement that has developed several areas of caverns. Equipment: While not absolutely essential, a small rugged gravity meter, such as the Worden meter, is highly desirable. A tall tripod, about the height of a transit tripod, permits instrument set-ups in deep water and in locations where fallen timber and muck piles make it impossible to use a short tripod. An additional advantage of a tall tripod is that it places the meter in the center of the crosscut, reducing the error caused by the crosscut void. Size and weight are important, since the only satisfactory means of operating the meter underground is to carry it by hand. A backpack can be used in rare instances but is usually a hindrance because of the close station spacing. The operator's ability to move through tight clearances will improve survey coverage, as it is then possible to move through raises and caved areas and to pass mine cars and machinery with a minimum of trouble. Station Control: Gravity stations are normally located every 100 ft along the crosscuts, at each intersection, and in the face of all stub crosscuts. In areas of high gravity relief, or where small anomalies might be expected, stations may be located at 25 or 50-ft intervals. When possible, the stations should be offset to avoid effects of raises or other voids. The gravity stations on a level are tied to one or more base stations, which are usually located at the shaft or near the portal of an adit. The base stations may be part of a gravity control net that extends to each level in the mine as well as to the surface. Such a net extending throughout the potential area of the surveys is highly desirable, as it is then possible to compare all gravity stations on a uniform basis. The stations that are part of the base net should be carefully established by multiple readings and, if necessary, by a least squares adjustment of the loops. In some instances where levels do not have a shaft station, or where access may be blocked by caving, it may be necessary to establish secondary bases at the top and bottom of the raises that are between levels. Under fair conditions 70 to 90 gravity stations can be located and run in 6 hr by a two-man crew. The best field procedures depend on conditions. Reduction of Field Data: Most of the time required to produce a final gravity map is consumed in processing the data. Each meter reading must be corrected for a minimum of five factors that affect the gravity value in addition to the density contrast being sought. These factors are 1) instrumental drift, 2) station elevation, 3) topography, 4) latitude, and 5) regional gravity gradient. Mine openings, such as stopes and raises, will affect the value. However, it is seldom practical to make corrections for these voids. Usually a rotation is made on the field note on the station, and any

Jan 1, 1957

By G. F. Leaming

The paper by R.V. Ramani and R.J. Sweigard is a wonderful description of the labyrinthine web that has been spun about the mining industry by energetic bureaucrats and politicians over the past 50 years. The remedy for the problem, however, is not more of the same, but less. That may be difficult for the industry to achieve, for it is not a technical solution but a political one. And the current fervor for more detailed planning at all levels of government and private enterprise has become deeply ingrained. The authors recommend the provision of more information about mining and mineral resources to "macro" (i.e., government) land use planners. They apparently overlook, however, the already strong tendency on the part of most government land use planners to consider themselves omniscient. Thus, giving them more information about the technical problems of mining will only make them want to get more and more involved in the "micro" (private, site specific) mine development and production plans of the individual mining firm. In fact, this has already happened at all levels of jurisdiction from municipal to federal government. Examples are legion. The most effective way to ameliorate the adverse impacts of government land use planning on existing and potential mining operations is to: (1) introduce greater flexibility in the definition of land use zones by local and state governments; (2) adopt realistic and relevant ambient environmental performance standards in governing relationships between mineral land uses and concurrent or subsequent nonmining land uses; (3) allow greater leeway for economic considerations in land use decisions in contrast to the explicit legalistic approach now in vogue; (4) recognize that all minerals are not the same and that sand and gravel mining should not be treated the same as underground metal mining, coal stripping, oil field production, or in situ leaching; and (5) eliminate the notion that mining operators should be responsible for determining in detail the use of land by subsequent owners of mined land. This last bit of conventional ethic really makes no more sense than requiring the builders of every shopping center or government office complex to provide detailed plans for the use of that land when its use for shopping or government is ended. Did the builder of Ebbetts Field plan for Brooklyn after the Dodgers went to Los Angeles? Should the developer of the Bingham Pit plan for suburban Salt Lake City after the copper mining goes to Chile? The nation's mining industry must address these questions before further bankrupting itself to provide more data to planners and spending thousands of dollars per acre to create land that when reclaimed is worth only a few hundred dollars per acre. ? Reply by R.V. Ramani and R.J. Sweigard We thank Mr. Learning for his valuable contribution. His views on the problems of land use planning and mineral resources are most welcome additions to our paper. As the title indicates, our paper was more concerned with the impacts of land use planning on mineral resource conservation than with the details of the planning process. On the whole, his five recommendations would be helpful for mineral resource conservation. However, we would suggest that the argument he presents for his final recommendation does not address the differences between mining as a land use and commercial or institutional uses. We believe that this difference is the crux of the issue. We share Mr. Learning's desire to ameliorate the adverse impacts of land use planning. Possibly the most detrimental impact is the loss of mineral resources. Any development, whether mineral or community, that does not give proper consideration to other resources can result in permanent loss or sterilization of resources. With proper planning, some of these losses can be avoided. As our paper indicated, one factor that limits the consideration of mineral resources, and ultimately leads to their sterilization, is the generally inadequate levels of resource characterization and understanding of the unique nature of mineral resources and mining operations. The last point raised by Mr. Learning is also important. In terms of reclamation and land use planning in mining districts, we certainly do not advocate spending more than what the results are worth. The main thrust of the paper was to explore the avenues for conserving the mineral resources so that, at some appropriate time, the issue of mining and reclamation can still be addressed. ?

Jan 1, 1986

By E. Teghtsoonian, E. M. Schulson

The tensile behavior of bcc ordered P' AuZn single crystals (CsCl structure) has been investigated under varying conditions of temperature, composition, and orientation. Between -0.2 and 0.4 T, multi-stage hardening occurs fm stoichiometric and nonstoichio-metric crystals oriented near the middle of the primary stereographic triangle. At higher and lower temperatures, parabolic type hardening occurs, followed by work - softening at the higher temperatwes. Deviations from stoichiometry give rise to increased flow stresses. Multi-stage hardening was observed for most orientations, except along the [loll-[lll] boundary and near the [001] corner of the stereo -graphic triangle, where parabolic type hardening occurs. Along two slip systems, (hk0)[001] and (, operate simultaneously while in the [001] comer, slip occurs mainly on the system. Electron microscopy of deformed crystals revealed bundles of edge dislocations forming walls approximately Perpendicular to the glide plane. In general the plasticity of 4' AuZn closely resembles the plasticity of bcc crystals. In recent years, considerable interest has arisen concerning the mechanical properties of the CsCl type intermetallic compounds Ag Mg,'- Fe co,' and Ni Al.'-' The compound P'AuZn is structurally similar. It has a low and congruent melting point of 725"~,'" remains ordered up to the melting point,16 and pos-esses a range of solid solubility from 47.5 to 52.0 at. pct Au at room temperature.15 The present paper reports the results of an investigation on the general tensile behavior of material in single crystal form. Some dislocation configurations characteristic of the deformed state are also reported. The results of a detailed study of the slip geometry in AuZn are presented in a separate paper.17 PROCEDURE Alloy preparation, crystal growing techniques, and the procedure followed in selecting specimens of minimum composition variation are reported elsewhere.17 Dumb-bell shaped tensile specimens were prepared by carefully machining single crystals in a jewellers' lathe to a gage length of 0.80 in. and diam of 0.090 in. Back-reflection Laue X-ray patterns and room temperature tensile tests revealed that machining damage could be eliminated by electrochemically polishing 0.005 in. from the machined surface followed by annealing at 300°C for 1 hr. Specimens were polished in fresh 5 pct KCN solution (40°C, 12 v). Experiments were performed by gripping specimens in a self-aligning pin-chuck and threaded collet system, then straining in a floor model Instron tensile machine. All tests were performed in duplicate. Experimental variables included temperature, composition, and orientation. Unless otherwise stated the strain rate was 2.5 x 10"3 per sec. Liquid testing environments included nitrogen (WOK), nitrogen cooled petroleum ether (133" to 293"K), and silicone oil (293" to 488°K). Resolved shear stress-shear strain curves were electronically computed from autographically recorded load-elongation curves. Stress and strain were resolved on the macroscopic noncrystallographic (hkO) [001] system operative under the specific test conditions of temperature, strain rate, and orientation reported earlier.17 RESULTS The temperature dependence of the work-hardening curves is shown in Fig. 1 for gold-rich crystals of 51.0 at. pct Au oriented near the center of the stereo-graphic triangle. Over the range of intermediate temperatures from -200" to 400°K, they are very similar to those classically observed for fcc metals (reviewed by Nabarro et al.).'' The beginning of deformation is characterized by a region of decreasing hardening rate, stage 0, which is followed by a region of low linear hardening, stage I, and then a region of higher linear hardening, stage 11. At the higher temperatures, stage 111 is observed, a region of decreasing hardening rate. Over the intermediate temperature range, the extent of stage 0 and of the slow transition between stages I and I1 decreases with increasing temperature. Total ductility is large, often greater than 300 pct shear. As the temperature is either increased or decreased, the extent of stage I is decreased, giving rise to parabolic type flow and reduced ductility. Similar temperature effects have been reported for bcc ~r~stals.~~-~~ Below -14O°K, hardening is terminated in brittle fracture while above -400°K. initial hardening is followed first by work-softening and then by chisel-edge type ductile fracture. Stoichiometric (50.0 at. pct Au) and Zn-rich (51.0 at. pct Zn) crystals were also tested from 77" to -500°K. The effect of composition on the flow behavior is illus-

Jan 1, 1970

By E. Eichen, G. S. Ansell. L. R. Sefton

A Technique for the Preparation of Thin Films of Two-Phase Alloys Suitable for Use in Transmission Electron Microscopy In order to obtain foil sufficiently thin to permit transmission electron microscopy of two-phase alloys, particularly alloys in which the second phase is present as a very finely dispersed second phase, a technique has been developed to permit the thinning of these materials without either the large degree of pitting present in the electrochemical or chemical thinning process or the slicing of the second phase particles accompanying the mi-crotoming process. The metal or alloy to be thinned is first ground to initial sheet thickness of approximately 0.002 in. by hand grinding on 4/0 emery paper. This starting sheet of material is then reduced in cross section by an air-abrasion technique using an S. S. White dental abrasive unit. This unit consists of a supply of compressed CO, which provides a stream of high velocity gas impinging on the metal surface. Incorporated in the machine is a hopper-type of device in which fine abrasive particles of alumina are introduced into this high velocity gas stream. The abrasive particles are irregular in shape, approximately 30 across. If any particles are introduced into the sample during thinning, they should be readily observable. Impingement of this gas abrasive mixture on the metal surface causes rapid removal of the metal in a rather uniform manner. The rates of metal removal may be adjusted by both regulating the gas pressure and amount of abrasive material introduced into the gas stream. By moving the gas stream abrasive mixture slowly about the specimen or sheet surface, a large area of the sheet may be uniformly reduced to a thickness which is then suitable for use in transmission electron microscopy. Attendant with this metal removal is a result in cold working of the metal surface due to the impingement of the abrasive particles. The resultant foils, therefore, if made from materials which are desirable to view in the annealed condition, have to be subsequently heat treated. If the alloy from which these foils have been made has been initially heavily worked or deformed during manufacture, then no subsequent heat treatment is, necessary as in the case of the SAP-type alloys. Following this initial abrasive thinning technique, the material is then given a final electropolish which is just sufficient to remove any oxide which has formed on the metal surface during the abrasion technique. Care has to be exercised in this technique to make sure that none of the attendant problems associated with the chemical or electrochemical removal process are encountered. An electrolyte solution which has been found suitable for aluminum-aluminum oxide SAP materials is one containing 78 ml of perchloric acid, 120 ml of distilled water, 700 ml of ethanol, and 100 ml of butylcellosolve. The current density used is 1 amp per sq cm and a polishing time of 3 sec on each side of the foil is used. Extreme care must be used in choosing the current density and time to ensure that just the oxide film on the specimen is removed. Times or current density which are greater than this will lead to perforation of the film due to the removal of the A1,0, particles or a complete loss of the film due to complete solution of the matrix. Once the oxide film is removed the specimens must be kept in alcohol until placed in the microscope to prevent any further oxidation. Fig. 1 shows -an electron micrograph at X28,000 of a thin foil of an aluminum-aluminum oxide SAP-type alloy prepared using this technique. The alloy consists of a matrix of commercial purity aluminum containing a very finely dispersed second phase of aluminum oxide flakes. In the electron micrograph, the dark patches are the aluminum oxide particles. The lighter portions of the micrograph are the aluminum matrix. Within the aluminum matrix, the dislocation structure of the alloy is clearly visible. This structure consists of intersecting twin boundaries and single dislocations. No abrasive particles were observed in the samples. In addition, thin films of commercial purity aluminum were prepared using this method and observed in the electron microscope. Here also, there was no evidence of abrasive particles introduced in the samples. The technique presented provides a satisfactory method for the production of thin films of aluminum-aluminum oxide SAP-type alloys for use in transmission electron microscopy where none of the usual thinning techniques appears to be satisfactory.

Jan 1, 1962

By W. Rostoker, R. F. Domagala, R. P. Elliott

The phase equilibria of the Hg-Th system over the composition range 0-100 pct Th and temperatures up to 1000°C have been studied for a small-volume, closed system. The solubility of Th in liquid Hg is about 5 pct at 300°C and decreases sharply with decreasing temperature. Two intermediate phases occur, Hg3 Th and HgTh. The structures of these are hexagonal (nonideally close-packed) and face-centered cubic, respectively. The HgTh phase decomposes eutectoidally at 400°-500°C. The solubility of Hg in solid thorium seems to be negligible. AFULL-phase diagram for this system would have to be defined on temperature-composition-pressure co-ordinates. This paper describes the pseudo phase diagram of a closed system, that is, where the alloy enclosed in a small volume equilibrates with a vapor pressure of mercury dictated by composition and temperature. Because of the experimental difficulties in studying a system of this nature, many of the phase relations can only be sketched. Alloy Preparation Alloys over the full range of composition were made from triple distilled mercury and one of two grades of thorium. For the bulk of the work, a calcium-reduced metal in sintered pellet form of reported 99+ pct total thorium content was used. Arc-melted specimens of this thorium gave a hardness of 135 VPN. The microstructure showed small primary dendrites of ThO2. A number of alloy compositions were made with a high-purity, iodide-decomposition thorium metal. The are-melted hardness of a button of this material was 35 VPN. Although the microstructure of the arc-melted specimens showed no dendrites of ThO2, there was definite evidence of an unidentified phase enveloping the grain bound-aries. There were no distinguishable differences between the constitution of alloys made with the two grades of thorium metal. Under normal conditions thorium is not wetted by liquid mercury. The film of ThO2 on all thorium metal cannot be penetrated by either liquid or vaporous mercury. It was therefore necessary to comminute thorium in the presence of mercury under such conditions that oxide films could not reform on the newly exposed metal surfaces. This was accomplished by the use of a high-speed, carbide-tipped rotary cutter incorporated in a chamber purged with argon and connected at the bottom to a demountable Vycor bulb containing a weighed amount of mercury. This experimental device is fully described in a separate paper.1 Alloy compositions were calculated by weighing the empty bulb, the bulb containing the mercury, and the bulb containing the mercury and the thorium chips. Many alloys were analyzed chemically for thorium and/or mercury after subsequent homogenization; the agreement between analyzed and calculated compositions was invariably very close. Bulbs containing the requisite amounts of mercury and fine thorium chips were clamped off, removed to a sealing unit, evacuated and sealed. Amalgamation under these conditions proceeded rapidly even at room temperature. To insure homogeneity, the specimens were annealed to 300-400°C. Alloys containing less than 30 pct Th remained pasty after all treatments, indicating an equilibrium condition of liquid plus solid. Alloys with more than 30 pct Th were transformed into a dark powdery product. These latter specimens were annealed for times of up to 1 week to complete interdiffusion. Many of the alloy compositions are pyrophoric. On exposure to air they oxidize with considerable evolution of heat to a mixture of ThO2 and free mercury. It was mandatory that alloy specimens be handled in a "dry box" purged thoroughly with argon. All X-ray diffraction specimens were powdered, screened, and sealed in capillary tubes within the dry box. Experimental Procedures Thermal analysis experiments, useful only in the mercury-rich region of the system, were conducted with the alloys in their original containers. A reentrant thermocouple well formed an integral part of the bulb. These bulbs were heated in a silicone oil bath and cooled in a dry ice-acetone mixture. The rates of heating and cooling were slowed by immersing the specimen bulb in a larger tube containing silicone oil. This provided a suitable thermal lag. In all tests, pure mercury was run as a basic standard. While the invariant reaction at about the melting point of mercury was detected by thermal analysis, the heat effect at the liquidus was not sufficient to produce an inflection in the cooling curve. It was necessary to determine the liquidus temperatures at the mercury-rich end of the system by "breaks" in electrical reslstivity versus temperature curves for individual alloys. The apparatus for this purpose consisted of a pyrex tube about 2 in. diam and 12 in

Jan 1, 1959

By ED. E. Furman, R. H. Atkinson

HITHERTO the practical creep testing of precious metals has received little or no attention. The only previous creep tests of precious metals have been made with wires under conditions such as to yield much more rapid rates of creep than in engineering tests.', ' Up to the present time the value of creep bars of adequate size, in the absence of real need for engineering data, has deterred investigators. However, the increasing use of platinum at high temperatures has demonstrated the need for reliable creep data for the guidance of engineers, especially those engaged in designing certain specialized chemical plant equipment. In order to supply this need, creep tests were conducted at 1382°F (750°C) on 0.290 in. diam specimens of platinum, 90 pct Pt, 10 pct Rh and palladium. The platinum was high purity, nominally 99.95 pct Pt. The 90 pct Pt, 10 pct Rh was of the same high quality as is used for making gauzes for the catalytic oxidation of ammonia. The palladium was also of high purity; two batches of palladium bars were tested, one deoxidized with calcium boride and the other with aluminum. Spectrographic examination of the palladium confirmed its good quality; the only significant impurities apart from the residual deoxidizers were traces of silicon and lead. Procedure The creep bars, which were furnished by Baker and Co. to our specification, were 6 ¾ in. in overall length with a 4½ in. (4 in. gage length) reduced section 0.290 in. in diam and had the ends threaded (?-NC16). It may be of interest that the bars were valued at up to $600 each. The specimens were supplied in a 50 pct cold-worked condition to facilitate attachment of the creep extensometer, which was of the push rod type. Because of the softness of the platinum and palladium, the extensometer rings were secured to the test section by means of circular knife edges instead of the usual pointed set screws. The extensometer rods extended through the bottom of the furnace and readings were taken with a 0.0001 in. "Last Word" dial gage fastened to the rods for the duration of the test. The bars were directly loaded by hanging weights from the lower specimen grip. All tests were conducted at 1382°F ± 2°F, and an effort was made to maintain the temperature gradient over the test section within 2°F. The ends of the furnace tube were packed with asbestos wool, which allowed a very slow circulation of air through the tube. Annealing was accomplished in the creep furnace before the load was applied. The platinum and palladium specimens were annealed at the test tem- perature for about 17 and 24 hr respectively; in the case of the rhodioplatinum it was found expedient to anneal for 1 hr at 1922°F (1050°C). Pilot samples cut from the same stock as the bars were used to check annealing procedures. Pertinent measurements of grain size and hardness were recorded. Results and Discussion The creep data obtained are given in Table I and the creep curves are plotted in Figs. 1, 2, and 3. Two platinum specimens, tested under a stress of 250 psi, had almost identical creep rates at 2000 hr, namely 0.000008 and 0.000009 pct per hr. A third platinum specimen, stressed at 400 psi, had a creep rate at 2000 hr of 0.000026 pct per hr; the reason for a rather sharp decrease in creep rate during the period from 1200 to 1600 hr is unknown. As it was thought that 90 pct Pt, 10 pct Rh would have a lower creep rate than platinum, the first sample was tested at 400 psi; however, the creep rate was approximately 50 pct greater. Microex-amination revealed that differences in grain size might be responsible for the unexpected result, as annealing at 1382°F developed an average grain diameter of 0.0021 in. in the rhodioplatinum specimen compared with 0.004 in. in the platinum bar. Annealing the alloy for 1 hr at 1922°F (1050°C) increased the average grain diameter to 0.0032 in. and materially improved the creep resistance, making it much better than platinum. A second specimen annealed at 1922°F (1050°C) and tested under a stress of 550 psi had a creep rate of 0.000022 pct per hr at 2000 hr, which was still substantially lower than that shown by the specimen annealed at 1382°F (750°C) and stressed at only 400 psi. In contrast to the creep behavior of the platinum and rhodioplatinum specimens, the palladium bars, whether deoxidized with calcium boride or aluminum, were characterized by high first stages of creep. However, after about 1200 hr of test, the creep

Jan 1, 1952

By H. R. Cooper, T. S. Mika

T. S. Mika (Department of Mineral Technology, University of California, Berkeley, Calif.) - Dr. Cooper's attempt to establish a correlation between process behavior and operational variables on the basis of a statistical analysis after imposing a reasonable process model is a very commendable improvement on the use of standard regression techniques. However, it must be recognized that the imposition of a model has the potential of yielding a poorer representation if its basic assumptions or mathematical formulation are invalid. It appears that at least two aspects of his treatment require some comment. First, the limitations on the kinetic law where xta represents a hypothetical terminal floatable solids concentration (cf. Bushell1), should be mentioned. Most current investigations2-9 appear to utilize the concept of a distribution of rate constants rather than a single unique value, k, to describe flotation kinetics. A distributed rate constant is certainly a more physically meaningful concept than that of a terminal concentration. The study of Jowett and safvi10 strongly indicates that xta is merely an empirical parameter, whose actual behavior does not correspond to that expected from a true terminal concentration. Rather than being a strictly mineralogical variable, as Dr. Cooper's treatment implies, it apparently represents the hydromechanical nature of the test cell as well as the flotation chemistry. The extension of batch cell kinetic results to full-scale continuous cell operation is a suspect procedure if the effect of such nonmineralogical influences on x,, remain unevaluated. There is evidence that introduction of a terminal concentration is necessitated by the inherent errors which arise in batch testing and are eliminated by continuous testing methods.' Possible lack of validity of the author's use of Eq. 1 is indicated by two unexpected results of the statistical analysis of his batch data. The first is the apparent corroboration of the assumption that the rate constant, k, is independent of particle size, i.e., of changes in the size distribution of floatable material. This assumption directly contradicts numerous results 2,4,11-l8 for cases where first order kinetics prevailed and ignores the phenomenological basis for the analysis of flotation in terms of a distribution of k's. It must be recognized that, if the rate constant is size dependent, the lumped over-all k would be time dependent; Eq. 1 would then no longer be valid. Cooper's x,, is determined by batch flotation of a distribution of sizes for an arbitrary period of time. If the size dependence of k is artificially suppressed, x,, will become a function of the experimental flotation time used in its determination. Upon reviewing the rather extensive literature concerning batch flotation kinetics, there appear to be few instances where constant k and x,, adequately adsorb variations in floatability due to particle size. The second surprising result is the low values of the distribution modulus, n, determined. Contrary to Cooper's assertion, most batch grinding (ball or rod mill) products yield values of n > 0.6, which increase as the material becomes harder.'' It is likely that the values of n = 0.25 and n = 0.42 for Trials 1 and 2, respectively, are completely unreasonable, and even the value n = 0.54 obtained for Trial 3 is unexpectedly low. Possibly, this indicates inherent flaws in the three trial models considered, in particular the assumed particle size independence of the rate constant, k. The above does not necessitate that Eq. 1 (and the terminal concentration concept) is invalid; it could constitute a good first approximation. However, the qualitative arguments used by Dr. Cooper in its justification are somewhat frail and require verification, particularly since much of the flotation kinetics literature is in opposition. Apparently, no effort was made to test these hypotheses on the actual data; in fact, since they pertain to a single batch test time, his data cannot be utilized to evaluate the kinetics of flotation. To evolve a control algorithm on the basis of this infirm foundation seems a questionable procedure. Another difficulty in his analysis arises in consideration of the froth concentrating process. As Bushel1 ' notes, for Eq. 1 to be valid it is necessary that the rate of recycle from the froth be directly proportional (independent of particle size) to the rate of flotation transport from the pulp to the froth, a restrictive condition." Harris suggests that it is more realistic to assume that depletion occurs in proportion to the amount of floatable material in the pertinent froth phase volume (treating that volume as perfectly mixed).12,21,22 The physical implications of

Jan 1, 1968

By David Turnbull

IN isothermal recrystallization processes, new crystals generally grow into the matrix until they impinge upon other new crystals or an external surface, at constant linear rates G. Before impingement the perceptible course of growth can be described by the equation: 1 = G(t-7) C1I where, G = dl/dt, 1 is a crystal dimension measured in a constant direction, t is the time, and 7, the nucleation period, is a positive intercept on the time axis. Fig. 1 is a schematic representation of I as a function of time for a recrystallizing grain. G is dependent upon temperature, driving energy (strain or surface energy), relative grain and boundary orientations, but is generally independent of time. The frequency of nucleation, fi, (time" volume") can be defined by the equation: N = 1/fV [2] where ? is the mean nucleation period and V is the volume of the specimen that has not recrystallized. The kinetics of primary and secondary recrystallization generally can be described satisfactorily in terms of the parameters N and G only.'-" After recrystallization is complete the average grain size 7 increases with time by "normal grain growth;" didt, the average rate of grain growth, is strongly time dependent and has not yet been precisely related to G for the motion of the individual grain boundaries constituting the system. It has been suggested4* " that the elementary act in grain boundary migration is closely related to the elementary act in grain boundary self-diffusion. Although the distance of atom movement in the two processes may be somewhat different, there is reason to expect that the activated states may be very similar, so that the free energy of activation for grain boundary migration should be of the same order of magnitude as for grain boundary self-diffusion. Therefore, it is desirable to develop a satisfactory basis for comparing data on self-diffusion and grain boundary migration and to make such comparisons where possible. Theory The formal theory of grain boundary migration rates is analogous to the theory for the rate of growth of crystals into supercooled liquids reviewed elsewhere 6-8. Boreliuss has shown that the latter theory describes, within the theoretical uncertainty, the growth of selenium crystals into supercooled liquid selenium. Motto and more recently Smolu-chowski" have derived expressions for the rate of boundary migration in recrystallization. The treatment to be presented is similar to Mott's excepting that the formalism of the absolute reaction rate theory will be used. The atomic mobility, M, in grain boundary migration is defined by: G = -M6p/6x where p is the chemical potential per atom and x is the coordinate measured in the direction of grain boundary movement. Let AF be the free energy difference per gram atom on the two sides of the boundary and k the thickness of the boundary. For RT>>AF the potential gradient across the boundary (6p/6x) is essentially linear and it follows that: SF/8x = - aF/N\ [4] where N is Avogadro's number. According to the Nernst-Einstein equation, M is related to a diffusion coefficient, Do, for matter transport in grain boundary migration by the equation: M = Da/kT [5] Substituting eqs 4 and 5 into eq 3 gives the basic relation between Do and G: G = DoaF/\RT [6] Do values may be calculated from experimental values of G from eq 6 and directly compared with the coefficient of self-diffusion within the crystal, DL, or the grain boundary self-diffusion coefficient D,. However, a more convenient, though equivalent, basis for comparing atomic mobility in grain boundary migration and self-diffusion is through the constants of the absolute reaction rate theory. According to this theory diffusion coefficients may be written:" D = k2(kT/h) exp [-AF,/RT] 171 aFa, the free energy of activation, is related to the measured energy of activation, Q, by the equation: AFA = Q - T aSx - RT [8] where aSa is the entropy of activation. Substituting eqs 8 and 7 into eq 6 gives: G = ek(kT/h) (aF/RT) exp [(AS,,)C/R] exp C-Qc/RTI C91 where the subscript G refers to boundary migration. The relationship between the driving free energy and the free energy of activation in boundary migration is indicated schematically in Fig. 2. Experience indicates that the variation of G with temperature can be described by: G= Go exp [- Qc/RT] [10] where Go and Qc are generally temperature independent over wide ranges of temperature. Comparison of eq 9 with eq 10 gives:

Jan 1, 1952

By A. Bauer, P. Calder, N. H. Carr, G. R. Harris

Since both rotary and jet piercing drills are used by the Iron Ore Co. at Smallwood, it is often desirable in planning to know in which regions of the orebody or new orebodies a particular drill will be the most economic. This makes it necessary to establish a correlation between drillability and pierceability and some physical rock properties. For rotary drills a good correlation was found with penetration rate and grinding factor index. The jet piercers were found to have a reciprocal relationship in the sense that the best rotary ground was the worst jet ground and vice versa. It is also indicated how an economic comparison could be made using these penetration rate versus grinding factor index curves, the hole size distribution curves for single pass and chambered holes and the mine distribution curve for grinding factor index. A discussion is presented on the fuel oxygen ratios to be used in jet piercing and on the site gas sampling and analysis which has been used to set up the drills. The fuel has been cut back so that stoichio-metric conditions exist, carbon monoxide is drastically reduced and pop-up or exploding holes eliminated. No decrease in penetration rate has been observed contrary to the published results of previous workers. The blasting procedure and results at Smallwood are discussed and the operation of Iron Ore Co.'s slurry pump-mix truck is also described briefly. Smallwood mine is part of the Iron Ore Co.'s Carol Lake operation and is situated in Labrador, 240 miles north of Sept-Iles, Quebec. Last year 15 million tons of crude ore were crushed to yield 6.3 million tons of concentrate and pellets. This year the figures will be 17 million tons of crude and 7% million tons of concentrate and pellets which is the full plant capacity. Carol Lake ores consist primarily of specularite and magnetite mixed with quartz. For convenience the ore has been split-into the following classifications depending on the percentage of magnetics in the sample, shown in brackets: specularite (0 to 10%), specularite-magnetite (10 to 20%), magnetite- specularite (20 to 30%), magnetite (>30%). The order of classification also represents the order of increasing grinding difficulty - the specularite generally being the easiest and the magnetite the hardest. The orebody also contains a small percentage of waste materials consisting of limonite carbonate, quartz carbonate and quartz magnetite. The first two materials are among the softest in the mine, generally softer than the specularite, and the quartz magnetite is amongst the hardest. The bulk of the material in the mine is of the specularite-magnetite and magnetite-specularite classifications. As a result of test drilling at Smallwood in 1960 with rotary, jet and percussion drills, the Iron Ore Co. purchased four JPM-4 jet piercers for the bulk of production drilling and set up an oxygen plant to supply 20 tons of oxygen per day. This oxygen is sufficient for two machines operating full time and one part time. In addition, there are two 50-R, one 60-R and one 40-R machines in use. The benches are 45 ft high and 50 ft holes are generally drilled. JET DRILLING At the onset of jet drilling in the late fall of 1962, two major problems were encountered: 1) freezing due to winter operations; experience and the use of heat at more places, such as the rotary head, has eliminated this,'" and 2) exploding or "popping" drilled holes; this happened frequently (several holes "popping" each day) and was the cause of two lost time accidents. In one instance a hole was being measured with a tape which fell down the hole causing it to "pop." Safety glasses though pulverized saved the wearer's eyesight. Various methods were then employed to detonate the holes before measuring or loading (dropping lighted rags of fusees down, or sparking across a spark gap). These methods were time consuming and far from completely successful. Consideration was given to the fuel oxygen ratio on the machines and what this would produce in the way of product gases. A fuel oxygen weight ratio of 0.35 which was quite oxygen negative was being used. Theoretically appreciable carbon monoxide would be produced at this fuel oxygen ratio. On the close down procedure of the jet which calls for low oxygen after flame out, oxygen would be left in the hole along with this carbon monoxide. This is an explosive mixture. The fuel oxygen ratio was cut back to stoichiometric

Jan 1, 1967

By E. P. Burtchaell

A method is presented for evaluating the effect of gravity drive upon the reservoir performance of a high relief pool. Conventional forms of reservoir analysis do not consider the alterations in the basic material balance data caused by gravity segregation of reservoir fluids. A procedure is outlined for structurally weighting physical and chemical data for use in the material balance equation. It is demonstrated how actual pool performance data can be utilized to evaluate the future reservoir performance of a gravity drive pool. INTRODUCTION Conventional reservoir engineering. procedure is inadequate for the analysis of an oil pool which has considerable structural relief, steep dips, and good permeability development. In, pools of this type, gravity drainage has an important part in the movement of oil to the wells and the effects of gravity on the overall pool performance should be included in any analysis of reservoir behavior. Many engineers have the opinion that the force of gravity in the movement of oil is not important until the later life of a pool.' Probably the basis for this belief is that gravitational effects may not be readily discernible until a pool is nearing depletion. This would be especially true for pools not having a high degree of structural relief and permeability development. Actually the effects of gravitational forces are at a maximum when the pool pressure is high, for during this period the hydrostatic head of the oil column is at a maximum and the viscosity of the oil is at a minimum. Oil recoveries from pools having favorable gravity drive characteristics may equal or even exceed recoveries which might be expected from water displacement. Field evidence indicates that in some reservoirs gravity drive has resulted in recoveries greater than that which could have been expected from gas expansion or water drive.'.3 Unfortunately, the possible effects of gravity drive on pool performance have been underestimated and other reasons have been sought to explain the high recoveries obtained. There are unquestionably many reservoirs to which the principles of gravity drainage can be effectively applied. It is the purpose of this paper to illustrate one method whereby gravity drive is included in the reservoir analysis of an oil pool. A hypothetical pool, typical of many California reservoirs, is used as an example. As used in this paper, "gravity drive" is defined as the overall effect of gravitational influences on the recovery of petroleum from the reservoir; "gravitational segregation" as the gravity separation of oil and gas within the reservoir; and "gravity drainage" as the downward movement of oil as caused by the force of gravity. SAND VOLUME DATA Fig. 1 presents a structural contour map of the pool under study. Maximum closure is 1950 feet with dips on the south flank approaching 45". The original gas-oil interface was set at -5200 feet. Average thickness of the producing sand was 200 feet. For use in subsequent calculations ill this paper, the pool was subdivided into 100-foot vertical increments and the sand-volume content of each increment was obtained. If the gross sand thickness is small, under 100 feet, the sand-volume content can be obtained by superimposing an isopachous map upon a structural contour map and planimetering the average thickness of each 100-foot increment. For sand thicknesses over 100 feet, one approacli would be to construct a sufficient number of cross-sections of the pool from which the weighted sand-volume of each 100-foot increment could be obtained. Variations in the sand body with depth, as determined by core data, can also be included in the above process. Table I presents a summary of sand-volume calculations, core data, and the original distribution of reservoir hydrocarbons in the pool. Fig. 2 illustrates the structural distribution of the sand-volume content. A total of 171,398 acre-feet is contained within the productive limits of the pool. Assuming an average porosity of 25% and an interstitial water content of 20%, the original hydrocarbon content was computed to be 227,075,000 barrels. DEPTH-PRESSURE DATA The determination of the initial vertical pressure arrangement in the pool is necessary for PVT and material balance calculations. Whenever sufficient data are available, a plot of pressure versus subsea depth of measurement should be made. From this plot a representative fluid pressure gradient can be established. Lacking sufficient initial pressure data, an initial pressure gradient may he estimated or calculated from avail-

Jan 1, 1949

By E. P. Burtchaell

A method is presented for evaluating the effect of gravity drive upon the reservoir performance of a high relief pool. Conventional forms of reservoir analysis do not consider the alterations in the basic material balance data caused by gravity segregation of reservoir fluids. A procedure is outlined for structurally weighting physical and chemical data for use in the material balance equation. It is demonstrated how actual pool performance data can be utilized to evaluate the future reservoir performance of a gravity drive pool. INTRODUCTION Conventional reservoir engineering. procedure is inadequate for the analysis of an oil pool which has considerable structural relief, steep dips, and good permeability development. In, pools of this type, gravity drainage has an important part in the movement of oil to the wells and the effects of gravity on the overall pool performance should be included in any analysis of reservoir behavior. Many engineers have the opinion that the force of gravity in the movement of oil is not important until the later life of a pool.' Probably the basis for this belief is that gravitational effects may not be readily discernible until a pool is nearing depletion. This would be especially true for pools not having a high degree of structural relief and permeability development. Actually the effects of gravitational forces are at a maximum when the pool pressure is high, for during this period the hydrostatic head of the oil column is at a maximum and the viscosity of the oil is at a minimum. Oil recoveries from pools having favorable gravity drive characteristics may equal or even exceed recoveries which might be expected from water displacement. Field evidence indicates that in some reservoirs gravity drive has resulted in recoveries greater than that which could have been expected from gas expansion or water drive.'.3 Unfortunately, the possible effects of gravity drive on pool performance have been underestimated and other reasons have been sought to explain the high recoveries obtained. There are unquestionably many reservoirs to which the principles of gravity drainage can be effectively applied. It is the purpose of this paper to illustrate one method whereby gravity drive is included in the reservoir analysis of an oil pool. A hypothetical pool, typical of many California reservoirs, is used as an example. As used in this paper, "gravity drive" is defined as the overall effect of gravitational influences on the recovery of petroleum from the reservoir; "gravitational segregation" as the gravity separation of oil and gas within the reservoir; and "gravity drainage" as the downward movement of oil as caused by the force of gravity. SAND VOLUME DATA Fig. 1 presents a structural contour map of the pool under study. Maximum closure is 1950 feet with dips on the south flank approaching 45". The original gas-oil interface was set at -5200 feet. Average thickness of the producing sand was 200 feet. For use in subsequent calculations ill this paper, the pool was subdivided into 100-foot vertical increments and the sand-volume content of each increment was obtained. If the gross sand thickness is small, under 100 feet, the sand-volume content can be obtained by superimposing an isopachous map upon a structural contour map and planimetering the average thickness of each 100-foot increment. For sand thicknesses over 100 feet, one approacli would be to construct a sufficient number of cross-sections of the pool from which the weighted sand-volume of each 100-foot increment could be obtained. Variations in the sand body with depth, as determined by core data, can also be included in the above process. Table I presents a summary of sand-volume calculations, core data, and the original distribution of reservoir hydrocarbons in the pool. Fig. 2 illustrates the structural distribution of the sand-volume content. A total of 171,398 acre-feet is contained within the productive limits of the pool. Assuming an average porosity of 25% and an interstitial water content of 20%, the original hydrocarbon content was computed to be 227,075,000 barrels. DEPTH-PRESSURE DATA The determination of the initial vertical pressure arrangement in the pool is necessary for PVT and material balance calculations. Whenever sufficient data are available, a plot of pressure versus subsea depth of measurement should be made. From this plot a representative fluid pressure gradient can be established. Lacking sufficient initial pressure data, an initial pressure gradient may he estimated or calculated from avail-

Jan 1, 1949

By Jr. F. M. Perkins.

Capillary forces play a controlling role in water-drive displacement processes both in laboratory experiments and in actual reservoirs, but their quantitative importance may be quite different in the two cases. Because of the importance of conducting laboratory experiments which are representative of field conditions, it is necessary to understand exactly the role of capillary forces in the displacement process. Though a number of experimental investigations related to this subject are contained in the literature, there appears to be a lack of information pertaining to unsteady-state experiments in water-wet media. This experimental study was conducted to obtain additional laboratory data to clarify further the role of capillary forces in both the macroscopic and microscopic flow of oil and water in porous materials. THEORETICAL CONSlDERATlONS The capillary pressure is defined as the difference in pressure between a continuous oil phase and a continuous water phase in a porous material.' The magnitude of this pressure difference depends on the interfacial curvature and the interfacial tension. The interfacial curvature is determined by the geometry of the pore spaces, the wettability of the rock surfaces, and the quantity of cach phase present. Capillary forces are involved in a water-drive displacement process in that they exert a controlling influence on the microscopic fluid distribution which in turn is reflected in the saturation or macroscopic flow behavior. Microscopic Fluid Distribution Because of the microscopic nature of the displacement of oil by water, it is necessary to consider the flow and the fluid distribution in individual pores. On this microscopic scale the capillary Forces, which act over a distance of one or two sand grain diameters, control the distribution of oil and water under static equilibrium conditions. When an external force is applied to the fluids, such as in a water-injection experiment, the applied forces tend to distort the oil-water interfaces. However, in most fine-grained, water-wet sands, the applied pressure difference across one or two grain diameters is usually several orders of magnitude less than the capillary pressure difference. These con-siderations lead to the theory' that even during flow the capillary forces continue to control the microscopic distribution of oil and water within the pores of a porous material for all practical reservoir and laboratory flow rates. This concept of capillary forces controlling the microscopic distribution of fluids has been substantially verified by other investigators3-7 who have found a lack of dependence of relative permeability and residual oil saturation on rate of fluid injection. Macroscopic Distribution The microscopic influence of capillary forces cannot be observed easily and only the effect on the macroscopic or average saturation can be detected. The saturation, of course, is really the point of interest. During a water flood, large differences in saturation at the flood front cause large capillary pressure gradients. This, in turn, causes water to advance ahead of the flood front thereby reducing the capillary pressure gradient in this region. The result is that in homogeneous porous media capillary pressure gradients tend to cause a diffuse displacement front. At low rates in laboratory columns, the front may extend over the entire column length. When the advancing water first reaches the outflow face of a core, the water, which is the wetting phase, cannot be produced because the pressure in the water just inside the core is lower than the pressure in the oil-filled space around the outflow face. This difference in pressure is equal to the capillary pressure for the water saturation existing at the outflow face. Water, therefore, accumulates at the outflow end of the core which causes a reduction in the capillary pressure. Because the capillary pressure does not vanish except at the residual oil saturation,7,8,9 water will not be produced until the residual oil saturation exists at the outflow face. This entire effect,"' which is called the "boundary effect", results in a region of relatively high water saturation near the outflow face. At low rates of injection in a short column, this region of high water saturation may extend over a considerable portion of the column. The influence of capillary forces on the macroscopic flow of oil and water have been described by Leverett.10 For unidirectional, viscous flow in the absence of gravity segregation, the expression in dimensionless form for the fraction of water in the flowing stream, f, is

Jan 1, 1958