Search Documents

By F. W. Jessen, G. F. Sears

Cavities in massive salt for the purpose of storage of liquid hydrocarbons have assumed a prominent position in recent years. This paper describes a program to facilitate leaching operations for the formation of specifically shaped storage cavities. Various forms and sizes of cavities may be possible through use of the techniques developed. INTRODUCTION The creation of large under ground storage facilities for natural gas and liquified petroleum gases has been practiced for many years. Use of cavities obtained through solution of salt for this purpose is a fairly new and novel approach, but has gained increasing importance due mainly to the economics of this type of mining operation opposed to hard rock mining or surface and pit storage.1 The idea of controlling the shape of any cavity dissolved from massive salt has not been a prime consideration of companies engaged in the formation of storage space. This has been due mainly to an insufficient knowledge of the mechanics of the leaching process and a dearth of published information dealing with both the desirability and ease of control possible for this type of operation. Two distinct advantages are readily ascertainable. From a stability standpoint (i.e., the ability of the completed cavity to withstand stress imposed by overburden pressure and lateral tectonic stresses) a controlled cavity may be generated which will yield the most favorable attitude to these external forces and remain in operable use for longer periods of time.2 A second advantage is that for a given volume, a sphere (which is one controlIed shape possible) represents the minimum surface area exposed. This may become more important in the future when dealing with refrigeration and product losses in underground cavities. The basic solution mining process, without regard for controlling the final shape, is quite simple in that the equipment and materials required to dissolve a cavity in a salt dome or layered salt section are a source of fresh water, a circulating pump, several strings of tubing and a means of disposal of the return brine. To add control measures involves the use of an inert blanket material above the position at which solution proceeds. Initially the annulus of the largest wash pipe is filled with this blanket down to the top of the proposed cavity. The bottom of the proposed cavity is determined by the depth of the original drilled hole or a blanking plug. The blanket material is added incrementally in stages and displaces the water in the enlarging cavity downward. Where the blanket material has displaced the water, no further solution takes place. The rate at which to add this controlling blanket material so as to be able to form specifically shaped cavities of any particular size is of considerable importance. THE PROBLEM It is desirable to have as much information as possible as to the behavior of the mining operation since visual observation obviously is impossible. It would be advantageous, for example, to have a step-by-step program showing how much salt is to be removed, at what rate and the time required. This type of program would facilitate procurement of surface equipment and provide for proper management of manpower requirements, as well as yield a host of intangible benefits. As mentioned earlier, only rule-of-thumb estimates were available until recently and made this type of pre-planned operation haphazardous at best. Recent work3-' has done much to clear up the ambiguities and inaccuracies attendant to solution mining of salt and has helped to place the operation on a more scientific basis. This study attempts to correlate much of the information thus far developed into a complete mining program; to extend the idea of controlled solution mining to include not only spheres, but solid conic sections of the ellipsoidal variety; and to refine the computation of the rate of removal of salt during the various stages of mining. Dommers3 and IIemson4 made extensive laboratory

Jan 1, 1967

By J. E. Hilliard

AS is well-known, the length per unit volume, Lv, of any lineal feature (such as grain edges in a polycrys-talline specimen) can be estimated from a count of the number of point intersections with a random plane of polish. Specifically:' Lv = 2NA in which NA is the expected number of intersections per unit area of the plane. As with all other measurements in quantitative metallography, the estimation of Lv is subject to a sampling error, the magnitude of which will depend on the number of observations. The analysis is usually performed by counting the number of intersections falling within a closed test figure of known area applied to a succession of different areas on the plane of polish. The test figure can be drawn directly on the screen of a projection microscope or else superimposed on the image of the structure by means of a reticule in a focusing eyepiece of an ordinary microscope. The use of micrographs is very inefficient unless, as in the present study, these are also needed for some other purpose. Consider an analysis in which n applications are made of a test figure within which the expected number of intersections is N. The expected total number of intersections counted is then given by NT = nN If (N) is the standard deviation in N for one application of the figure and if the latter is applied to non-overlapping regions so that the N's are independent in the statistical sense, then the standard deviation in Nt is given by: a(NT) = a(N)vn Eliminating n from Eqs. [2] and [3] we obtain: Thus, k will be a constant for a given specimen and test figure. From Eqs. [ I] and [4] we find As is to be expected, therefore, the fractional error in the estimation of Lv is inversely proportional to the square root of the total number of intersections counted. It is very desirable to have at least a rough estimate of k before starting an analysis, since it will then be possible to decide in advance by use of Eq. [5J how many counts should be made in order to achieve the desired accuracy in Lv. Without such information the investigator, unless lucky, will either do needless (and tedious) work or else obtain a result with too large an error. It is very easy to calculate k for the limiting case of N — 0, i.e., when the test figure is vanishingly small. In this case there is a negligible probability of two or more intersections falling within the figure. There are therefore only two possible outcomes: a count of one with a probability of, say, p and a zero count with a probability of (I-p). Thus and Therefore: We might expect k to depend on N and therefore on the size of the test figure. Such a dependence was observed2 in the analogous case involving the estimation of grain-boundary area per unit volume by a count of the number of intersections of a test circle with traces of the boundaries on the plane of polish. In this case the theoretical value of k at the limit N — 0 was 1.41. It was also predicted that k would display an initial sharp decrease with increasing N. Experimentally, k was observed to decrease to a value of 0.62 at N - 6 and thereafter remain approximately constant. Furthermore, the value of k for N > 6 was found to be insensitive to the form of the structure provided no pronounced anisotropy or grain-size variation was present. In order to determinewhether k would exhibit a similar dependence on N in the estimation of grain-

Jan 1, 1967

By G. R. Countryman, G. H. Thomas, I. Fatt

Procedures are described for scaling up turbulent friction pressure drops of inelastic fluids and laminar friction pressure drops of inelastic/elastic fluids in circular ducts. The laminar scale-up calculation method does not require the use of arbitrarily defined flow models, and can be performed directly from laminar capillary viscometry data, or by means of a numerical solution of the integral equation for Poiseuille flow, inputing true shear rate-shearing stress data from other simple shearing experiments. e.g. Couette and plane Poiseuille flows. The turbulent scale up procedure, based on the Dodge-Metzner correlation for inelastic fluids, requires evaluating the characteristic rheological parameters at the existing wall stress under laminar flow conditions. Since the wall stress is the quantity sought in a pipe flow scale-up problem, a trial-and-err or solution is indicated, and therefore a calculation method has been developed and programmed for a medium-size magnetic-drum-memory digital computer. Laminar flow data from other somces may also be adapted to this turbulent scale-upmetbod for inelastic fluids. Comparisons between predicted and experimentally observed friction losses for generalized Reynolds numbers ranging from 3 — 32.000 indicate an error band width of less than 10 per cent in circular tubes ranging in size from 0.5 to 1.89 in. in diameter. INTRODUCTION Considerable effort has been devoted to the problem of establishing an accurate procedure for predicting the laminar and turbulent flow behavior of non-Newtonian materials in pipes from data obtained with small scale laboratory or pilot plant equipment, as evidenced by the variety of reviews on the subject.1-3 Although there does not appear to be a universally accepted scale-up method for laminar flow, the various methods currently employed can be divided into two principal categories depending upon whether a functional relationship between shear rate and shearing stress is assumed, or whether a generalized correlation is adopted. The methods of Hedstrom4 and of Weltmann are typical of the first category, while the generalized correlation of Metzner and Reed6 illustrates the second category. Under appropriate conditions either approach yields a reasonably satisfactory description of laminar flow behavior. In contrast, the prediction of non-Newtonian transitional and turbulent friction losses is much more complex. For example, certain types of materials have been observed to sustain laminar flow at much higher flow rates than Newtonian liquids. When turbulence does develop, the level of agitation is damped to such an extent that the friction factor is lower than that of the original low viscosity Newtonian medium. Extremely broad transition regions and "diameter effects" have also been noted in the turbulent flow of certain chemically complex materials.' Aside from the latter complications, which are attributed in part to viscoelasticity, the existence of a variety of scale-up procedures for predicting transitional and turbulent friction losses is due principally to the chioce of viscosity parameter used in the definition of the Reynolds number. In contrast, Dodge and Metzner have demonstrated that the friction factorReynolds number relationship for inelastic fluids can be satisfactorily generalized provided the fluid property parameters are evaluated at the existing wall shearing stress.* However, the quantity usually sought in a pipe flow scale-up problem is the wall stress, and hence a trial-and-error procedure is required when dealing with non-power law materials. In this report we discuss two calculation methods: (1) for determining the wall stress, and hence the friction pressure drop, for inelastic fluids under conditions of turbulent flow by means of an iterative procedure, and (2) for predicting laminar friction pressure drops without recourse to arbitrarily defined flow models.

By A. R. Gregory

A shear wave velocity laboratory apparatus and techniques for testing rock samples under simulated subsurface conditions have been developed. In the apparatus, two electromechanical transducers operating in the frequency range 0.5 to 5.0 megahertz (MHz: megacycles per second) are mounted in contact with each end of the sample. Liquid-solid interfaces of Drakeol-aluminum are used as mode converters. In the generator transducer, there is total mode conversion from P-wave energy to plain S-wave energy. S-wave energy is converted back to P-wave energy in the motor transducer. Similar transducers without mode converters are used to measure P-wave velocities. The apparatus is designed for testing rock samples under axial or uniform loading in the pressure range 0 to 12,000 psi. The transducers have certain advantages over those used by King, 1 and the measurement techniques are influenced less by subjective elements than other methods previously reported. An electronic counter-timer having a resolution of 10 nanoseconds measures the transit time of ultrasonic pulses through the sample; elastic wave velocities of most homogeneous materials can be measured with errors of less than 1 percent. 5- and P-wave velocity measurements on Bandera sandstone and Solenholen limestone are reported for the axial pressure range 0 to 6,000 psi and for the uniform pressure range 0 to 10,000 psi. The influence of liquid pore saturants on P- and S-wave velocity is investigated and found to be in broad agreement with Biot's theory. In specific areas, the measurements do not conform to theory. Velocities of samples measured under axial and uniform loading are compared and, in general, velocities measured under uniform stress are higher than those measured under axial stress. Liquid pore fluids cause increases in Poisson's ratio and the bulk modulus but reduce the rigidity modulus, Young's modulus and the bulk compressibility. INTRODUCTION Ultrasonic pulse methods for measuring the shear wave velocity of rock samples in the laboratory have been gradually improved during the last few years. Early experimental pulse techniques reported by Hughes et al.2 and by Gregory 3 were beset by uncertainties in determining the first arrival of the shear wave (S-wave) energy. Much of this ambiguity was caused by the multiple modes propagated by piezoelectric crystals and by boundary conversions in the rock specimens. Shear wave velocity data obtained from the critical angle method, described by Schneider and Burton4 and used later by King and Fate sand by Gregory, 3,6 are of limited accuracy, and interpreting results is too complicated for routine laboratory work. The mode conversion method described by Jamieson and Hoskins7 was recently used by King 1 for measuring the S-wave velocities of dry and liquid-saturated rock samples. Glass-air interfaces acted as mode converters in the apparatus, and much of the compressional (P-wave) energy apparently was eliminated from the desired pure shear mode. A more detailed discussion of the current status of laboratory pulse methods applied to geological specimens is given in a review by Simmons.8 A new shear wave velocity laboratory apparatus has been developed for testing rock samples under elevated pressures. This apparatus employs liquid-solid interfaces as mode converters. The technique for generating pure shear waves at liquid-solid interfaces was developed from studies of critical angle methods.6 Two identical transducers are placed in contact with opposite ends of the rock sample, one transducer acting as an ultrasonic energy emitter or generator and the other acting as an ultrasonic energy receiver or motor. The P-wave energy emitted from the piezoelectric crystal in the generator transducer is converted into pure S-wave energy at the Drakeol-aluminum

By P. K. Weimer

Coinpletely integrated thin-film circuits inco?,porating more than 1000 active and passive elements have been fabricated reproducibly in the laboratory by evaporation of- all components. A 180-stage modified shift register comprising 540 CdSe TFT's, 360 resistovs, and 180 capacitors was deposited upon a glass substrate during one pump-down of the vacuum systeitl. Experimental units have been operated for more than 5000 hr at temperatures uP to 85°C without failure of any slages. The present pupper reviews the charvactevistics of evaporated thin-film transistors (TFT's) and discusses their application to integrnted circuils. INTEGRATED circuits most widely used today are based upon silicon technology, with the circuits formed on the surface of a block of single-crystal silicon1 or on a block of sapphire 2,3 upon which a thin layer of silicon has been deposited. Transistor patterns and interconnections are defined by photographic procedures using photoresist masking. Large groups of transistors are produced in batches which undergo many processing steps, some of which are carried out at elevated temperatures. Passive elements are formed by treatment of the semiconductor surface or are deposited by evaporation alongside the active elements. The size and complexity of the circuits are limited by defects in the devices and by the size of the available single-crystal chips. Capacitative coupling to the semiconducting substrate may also limit useful size and frequency response of the circuits. An alternative approach to integrated circuits is based upon the deposition by evaporation of all circuit components upon a noncrystalline, insulating substrate such as glass.4 Thin-film field-effect transistors,1'6 known as TFT's, can be deposited in large arrays on glass or plastic substrates at the same time the associated passive elements are being deposited. Although still in the research stage the all-evaporated approach is expected to offer certain advantages for many applications. In particular, the use of thin films upon a noncrystalline substrate offers the prospect of producing complex, large-area circuits not now feasible with single-crystal silicon. Important applications of this type include memory arrays for computers and imaging devices for television. The ease with which the evaporated TFT's can be deposited may permit continuous, automated production methods not possible with high-temperature batch processing. Recent advances7'8 in the fabrication of thin-film transistors have led to greatly improved reproduci-bility and life of the TFT's. Experimental circuits having over 1000 active and passive elements deposited upon a common glass substrate are currently being produced in the laboratory. The present paper reviews recent progress in the development of the TFT and in its application to integrated circuits. I) DESCRIPTION OF THE THIN-FILM FIELD-EFFECT TRANSISTOR The TFT is a field-effect transistor of the type having the control gate separated from the semiconductor by a thin insulating spacer. Although its operation resembles that of the silicon MOS transistorg its construction is based entirely on thin-film techniques. Fig. 1 shows a cross-sectional view of an evaporated TFT deposited upon a glass substrate compared with a silicon MOS transistor deposited upon single-crystal sapphire.3 The semiconductor in the TFT is a poly crystal line layer, of a material which can be readily evaporated, such as cadmium sulfide, cadmium selenide,l0 or tellurium." Silicon monoxide layers have usually served as the gate insulator, but other oxides, sulfides, and fluorides may be used. Film thicknesses range from a few hundred to a few thousand angstroms. The source, drain, and gate electrodes in the TFT are of metal. Qpical source-drain gap dimensions are 0.0003 x 0.030 in., with a gate strip whose width slightly overlaps the source and drain. The source and drain electrodes should make low impedance contacts to the semiconductor. With cadmium sulfide films, gold is a satisfactory material for underlying electrodes, while aluminum is preferable for overlying electrodes.12 Considerable variation in design of the TFT is permissible, with the electrodes being

Jan 1, 1967

By B. B. Wood, M. J. Rzasa, E. W. Hough

An apparatus for the determination of fluid-fluid interfacial tension by the pendant drop method has been constructed. The apparatus is refined beyond those previously described in that the samples are introduced into an evacuated windowed cell without contacting air or mercury. Interfacial tension data are reported for the water-methane system between 74 and 280°F and 15 and 15,000 psi. The variation of wettability of the stainless-steel dropper tip with pressure and temperature was observed incidentally to the measurements. INTRODUCTION The most refined previous apparatus --for the determination of fluid-fluid interfacial tension at reservoir conditions was built by Hauser and Michaels.' These authors used the pendant drop method, previously investigated by Andreas, Hauser, and Tucker.' Recently, Fordham3 has made the pendant drop method an absolute one by numerical integration of the differential equation of drop shape. It is necessary to know the difference in density of the two fluids and the acceleration due to gravity in order to evaluate the interfacial tension from measurements of drop dimensions. Contact angle is not involved in the measurements by this method. Some precise tests of the method have been made at atmospheric pressure and 77°F by Douglas," who found agreement with established values to within about 0.1 per cent for water and air, cyclohexane and air, benzene and air, and benzene and water. No values of interfacial tension from the recent pendant drop apparatusJ have been published by Hauser and Michaels. However, values for the water-benzene and water-decane systems are given in Michaels' thesis5 or Dressures from 15 .to 10,000 psi and temperatures from 73 to 268ºF. Some determinations on gas-oil and gas-water systems in the range of 15 to 4,000 psi and 78 to 178°F have been made by Hozott,' by the drop-volume method. The capillary rise method has been used to investigate the interfacial tension in gas-oil systems at 88°F to several thousand psi by Swartz.7 Other determinations in gas-oil systems at pressures under 1,000 psia have been made by Jones: and Beecher and Parkhurst.9 The scarcity of data for pressures above 5.000 psi is apparent. APPARATUS An assembly drawing of the working section of the interfacial tenqion apparatus is given in Fig. 2. 'The endanant drop (A) is housed in a pressure chamber having glass windows (B). An agitator (6) is provided. The metal parts of the apparatus coming into contact with the sample are made of AISI: Type 316 steel. Drain grooves are provided around the whdows so that the sample in contact with the O-rings is not returned to the working section when the pressure is reduced. A schematic arrangement of the apparatus is given in Fig. 3. The working section containing the pendant drop (A) is filled by samples from bombs (B) with the aid of the vacuum system (C). The O-ring sealed pistons (D) used to separate the sample from the hydraulic fluid are shown together with the pressure production and measurement system (E). A front view of the apparatus is shown in Fig. 1. Most of the hydraulic system is at the top of the panel. The pump and displacement piston are at the lower right. The Variacs for the air bath heaters. the camera. and the vacuum gauge are a1so visible on the panel. The working section together with a window closure and a drop-forming tip closure are shown in Fig. 9. The carbon arc light source with filter is visible in the upper part of the picture. Two sides of the air bath with the impeller and a part of the convection shield can also be seen. The assembled air bath is shown in Fig. 10.

Jan 1, 1951

By E. W. Hough, B. B. Wood, M. J. Rzasa

An apparatus for the determination of fluid-fluid interfacial tension by the pendant drop method has been constructed. The apparatus is refined beyond those previously described in that the samples are introduced into an evacuated windowed cell without contacting air or mercury. Interfacial tension data are reported for the water-methane system between 74 and 280°F and 15 and 15,000 psi. The variation of wettability of the stainless-steel dropper tip with pressure and temperature was observed incidentally to the measurements. INTRODUCTION The most refined previous apparatus --for the determination of fluid-fluid interfacial tension at reservoir conditions was built by Hauser and Michaels.' These authors used the pendant drop method, previously investigated by Andreas, Hauser, and Tucker.' Recently, Fordham3 has made the pendant drop method an absolute one by numerical integration of the differential equation of drop shape. It is necessary to know the difference in density of the two fluids and the acceleration due to gravity in order to evaluate the interfacial tension from measurements of drop dimensions. Contact angle is not involved in the measurements by this method. Some precise tests of the method have been made at atmospheric pressure and 77°F by Douglas," who found agreement with established values to within about 0.1 per cent for water and air, cyclohexane and air, benzene and air, and benzene and water. No values of interfacial tension from the recent pendant drop apparatusJ have been published by Hauser and Michaels. However, values for the water-benzene and water-decane systems are given in Michaels' thesis5 or Dressures from 15 .to 10,000 psi and temperatures from 73 to 268ºF. Some determinations on gas-oil and gas-water systems in the range of 15 to 4,000 psi and 78 to 178°F have been made by Hozott,' by the drop-volume method. The capillary rise method has been used to investigate the interfacial tension in gas-oil systems at 88°F to several thousand psi by Swartz.7 Other determinations in gas-oil systems at pressures under 1,000 psia have been made by Jones: and Beecher and Parkhurst.9 The scarcity of data for pressures above 5.000 psi is apparent. APPARATUS An assembly drawing of the working section of the interfacial tenqion apparatus is given in Fig. 2. 'The endanant drop (A) is housed in a pressure chamber having glass windows (B). An agitator (6) is provided. The metal parts of the apparatus coming into contact with the sample are made of AISI: Type 316 steel. Drain grooves are provided around the whdows so that the sample in contact with the O-rings is not returned to the working section when the pressure is reduced. A schematic arrangement of the apparatus is given in Fig. 3. The working section containing the pendant drop (A) is filled by samples from bombs (B) with the aid of the vacuum system (C). The O-ring sealed pistons (D) used to separate the sample from the hydraulic fluid are shown together with the pressure production and measurement system (E). A front view of the apparatus is shown in Fig. 1. Most of the hydraulic system is at the top of the panel. The pump and displacement piston are at the lower right. The Variacs for the air bath heaters. the camera. and the vacuum gauge are a1so visible on the panel. The working section together with a window closure and a drop-forming tip closure are shown in Fig. 9. The carbon arc light source with filter is visible in the upper part of the picture. Two sides of the air bath with the impeller and a part of the convection shield can also be seen. The assembled air bath is shown in Fig. 10.

Jan 1, 1951

By P. J. Root, J. E. Warren

An idealized model has been developed for the purpose of studying the characteristic behavior of a permeable medium which contains regions which contribute significantly to the pore volume of the system but contribute negligibly to the flow capacity; e.g., a naturally fractured or vugular reservoir. Unsteady-state flow in this model reservoir has been investigated analytically. The pressure build-up performance has been examined in some detail; and, a technique for analyzing the build-up data to evaluate the desired parameters has been suggested. The use of this approach in the interpretation of field data has been discussed. As a result of this study, the following general conclusions can be drawn: 1. Two parameters are sufficient to characterize the deviation of the behavior of a medium with "double porositym from that of a homogeneously porous medium. 2. These parameters can be evaluated by the proper analysis of pressure build-up data obtained from adequately designed tests. 3. Since the build-up curve associated with this type of porous system is similar to that obtained from a stratified reservoir, an unambiguous interpretation is not possible without additional information. 4. Differencing methods which utilize pressure data from the final stages of a build-up test should be used with extreme caution. INTRODUCTION In order to plan a sound exploitation program or a successful secondary-recovery project, sufficient reliable information concerning the nature of the reservoir-fluid system must be available. Since it is evident that an adequate description of the reservoir rock is necessary if this condition is to be fulfilled, the present investigation was undertaken for the purpose of improving the fluid-flow characterization, based on normally available data, of a particular porous medium. DISCUSSION OF THE PROBLEM For many years it was widely assumed that, for the purpose of making engineering studies, two param- eters were sufficient to describe the single-phase flow properties of a producing formation, i.e., the absolute permeability and the effective porosity. It later became evident that the concept of directional permeability was of more than academic interest; consequently, the degree of permeability anisotropy and the orientation of the principal axes of permeability were accepted as basic parameters governing reservoir performance. 1,2 More recently, 3-6 it was recognized that at least one additional parameter was required to depict the behavior of a porous system containing regions which contributed significantly to the pore volume but contributed negligibly to the flow capacity. Microscopically, these regions could be "dead-end" or "storage" pores or, macroscopi-cally, they could be discrete volumes of low-permeability matrix rock combined with natural fissures in a reservoir. It is obvious that some provision for the inclusion of all the indicated parameters, as well as their spatial variations, must be made if a truly useful, conceptual model of a reservoir is to be developed. A dichotomy of the internal voids of reservoir rocks has been suggested.7,8 These two classes of porosity can be described as follows: a. Primary porosity is intergranular and controlled by deposition and lithification. It is highly interconnected and usually can be correlated with permeability since it is largely dependent on the geometry, size distribution and spatial distribution of the grains. The void systems of sands, sandstones and oolitic limestones are typical of this type. b. Secondary porosity is foramenular and is controlled by fracturing, jointing and/or solution in circulating water although it may be modified by infilling as a result of precipitation. It is not highly interconnected and usually cannot be correlated with permeability. Solution channels or vugular voids developed during weathering or burial in sedimentary basins are indigenous to carbonate rocks such as limestones or dolomites. Joints or fissures which occur in massive, extensive formations composed of shale, siltstone, schist, limestone or dolomite are generally vertical, and they are ascribed to tensional failure during mechanical deformation (the permeability associated with this type of void system is often anisotropic). Shrinkage cracks are the result

By C. J. Rosa, W. C. Hagel

An attempt is made to determine the diffusion coefficients of oxygen ions in hypostoichiometric zirconium oxide. A phenomenological theory is developed for three-phase, unidirectional volume diffusion of a gaseous species from a planar source. Knowledge about oxygen diffusion in a and ß zirconium and experimentally measured rates of movements of interfaces permit the evaluation of oxygen ion diffusivi-ties in ZrO2-x. The results conform to the relationship: sq cm per sec, for the temperature range of 875" to 1050° C. DIFFUSION coefficients of oxygen ions in stoichio-metric oxides have been studied by many investigators. The obtained results not only supplement the knowledge about the lattice-defect structures of oxides and conduction mechanisms but also provide information of technological importance. Rate processes such as oxide sintering, high-temperature creep deformation, and oxidation of metals and their alloys are directly dependent on diffusion phenomena. Unfortunately, diffusion coefficients obtained from studies on stoichio-metric oxides cannot be treated without restraints when applied to oxidation kinetics, if precise evaluations of oxidation rates are desired. Oxides formed on metals are usually nonstoichiometric and exhibit nonuniform distribution of either cation or anion vacancies across their thicknesses. Debuigne and Lehr1 extended Danckwerts'2 mathematical analyses and furnished a semitheoretical treatment for evaluating oxygen diffusion coefficients in zirconium oxide in terms of weight gain measurements during oxidation of zirconium in the temperature range 400° to 850°C. Their calculations are based on theoretical partitioning3 of the total amount of oxygen between a zirconium and its oxide phase. An inherent discrepancy appears to continue throughout their calculations. The amount of oxygen confined to the oxide phase should be corrected by multiplying the calculated oxide thickness by the inverse of Pilling-Bedworth ratio, if expansion of zirconium dioxide is taken into account. The purpose of the present investigation is to provide a new semitheoretical treatment for evaluating diffusion coefficients of moving species in hypostoichiometric oxide. The analyses are relevant to unidirectional volume diffusion. Although the presented method is quite general, its applicability will be discussed in terms of kinetic measurements from oxidation of zirconium at temperatures above the a/ß transition, for reasons stated in the later part of this paper. EXPERIMENTAL Pure, 99.93 wt pct, zirconium samples of approximately 2 by 1 by 0.35 cm were exposed to oxygen for various times at four temperature levels: 875°, 950°, 1000°, and 1050°C. High-purity, better than 99.9 vol pct oxygen was passed through a cold trap and columns containing anhydrous CaSO4 to remove residual water vapor. Before exposing to oxygen at a pressure of 500 Torr, specimens were subjected to wet abrasion on 200- through 600-grit Sic papers and polished with 1—p diamond dust. Prior to oxidizing, all zirconium samples were vacuum-annealed at 10-6 mm Hg for about 15 min at the test temperature. The temperature variation across the uniform heat zone was ±2°C from the preset value. Final polishing, for phase thicknesses measurements, was performed by using 8- and 1-p diamond dust and followed by vibrator polishing with 0.05-µ ? Al2O3 solution. Measurements of phase thicknesses were made by projection on a metallographic viewing screen. Conventional powder X-ray analyses for the oxide were employed. RESULTS Fig. 1 is an example of the cross section of oxidized zirconium. Designations of phases given below this figure are the same as those used in the three-phase diffusion system represented by Fig. 5 in the following section. Zirconium oxidized for longer exposure times showed two types of oxides: white oxide at the surface and an underlaying layer of gray oxide. X-ray powder patterns of these oxides at room tern-

Jan 1, 1969

By A. Gittens

In two recent papers stevens13,14 has compared the methods used by various investigators to calculate the contribution of grain boundary sliding to the total creep strain. In previous work Gittins15 has shown that the strain due to grain boundary sliding Egb may be estimated from the lateral, d, and vertical, v, displacements of a longitudinal marker line as it crosses a grain boundary. He obtained the formula, where nL is the number of grains per unit length measured on the specimens before creep, and \p and ?2 are, respectively, the angle between the grain boundary and the stress axis in the plane of the surface and normal to the surface (see Ref. 13, Fig. 1). The subscript L indicates that observations are made in a longitudinal direction parallel to the stress axis. It was generally found more convenient to measure the grain size of a specimen before creep; this was necessary to characterize the specimen. Stevens13,14 however prefers to define Egb in terms of NL, the number of grains per unit length after creep. Since nL - NL(l + €1) where cr is the total strain, Eq. 1 11 becomes The quantities d and ?, can be determined directly from the marker lines, while v and ?2 are found by sectioning normal to the surface. In the interior when symmetry considerations suggest that v = d and ?1 - ?2 the strain due to grain boundary sliding may be estimated from a suitable internal longitudinal marker line by the equation To determine cgb at the surface using either Eqs. |1| or [2] requires a large number of observations, particularly if many estimates of cgb are required, and also it is impossible to determine the variation of cgh with creep strain on the same specimen because the specimen must be sectioned to determine ?2. Gittins15 found that for a number of materials the sliding vectors d and v were related to each other so that provided one was measured the other could be estimated. Also during creep, grain boundaries tend to migrate to meet the surface at angles close to 90 deg,16 so that the term nL(v/tan?2)L was comparatively small when 4,- 90 deg. By determining the frequency distributions of ?2 on control specimens this term can be estimated without sectioning the creep specimen. In this way an empirical equation based on Eqs. |1| and |2| enables cgb to be determined with reasonable accuracy with fewer observations of the angles and sliding components. For B brass, the relationship between the sliding vectors indicated that ceb = 0.85nLVL . In Table I of Stevens paper this equation is stated to be incorrectly derived. However it should be pointed out that this is an empirical equation derived from Eq. |1|, based on a knowledge of the relationship between d and v and on the frequency distributions of ?1 and ?2 for a particular material under given experimental conditions. It was used only when these relationships were known, to reduce the number of observations required to estimate cb. Therefore the comment "derivation incorrect does not apply. Gittins did in fact make use of Eqs. |1| and 131. Author's Reply I must apologize to Dr. Gittins for the inadvertent misrepresentation of his work. Other slight errors in Table I have been pointed out to me. Langdon (No. 12 Table I), Gifkins (No. 7, Table I), and possibly other workers did not measure the number of boundaries per unit length after creep (A') but the corresponding quantity bejove creep (n). As Dr. Gittins points out, a completely analogous set of equations to those given in the paper can be derived by the same method: using the number of boundaries per unit length before instead of after creep. The basic equations corresponding to Eqs. [2] and 171 of the paper, from which all the others can be derived, are

Jan 1, 1968

By A. Schild

The effect of heating a producing well on the rate of oil recovery has been analyzed in the simple case of a well producing oil by a radial drive and in the steady state. Differences of production rates with and without borehole heating have been calculated and plotted; these measure some of the benefits of the thermal recovery process. A striking feature is that the production increase tends rapidly to a constant value for large production rates. Some concrete examples, including a rough comparison with field data. are discussed. Heating requirements have been calculated and plotted on the assumption that the heaters used are 100 per cent efficient. Heating has been thought of as applying only to low volume wells producing heavy oils. However, the calculations suggest that the thermal process considered might be profitable even for some wells producing large quantities of light oil, especially in unconsolidated formations where fracturing is not possible. In the case of wells producing other fluids with the oil or of wells producing partly or wholly by gravity drainage, the theory developed does not apply in detail. It can be expected, however, that any general qualitative be-havior will remain valid. Under normal temperature conditions, high viscosity emulsions sometimes form around the wellbore, lowering the permeability in the vicinity of the well and impeding the production of oil. Borehole heating often prevents such emulsification. This important effect is not easily amenable to mathematical analysis and has not been included in the theory presented here. It should therefore be kept in mind that the benefits due to borehole heating may in fact exceed the values given by the following calculations. INTRODUCTION The recovery of viscous, low gravity oils can often be stimulated by applying heat to the production zone. The higher temperature lowers the viscosity of the oil, increasing its mobility, and therefore its rate of production. One thermal method of recovery 1,2,3,4. consists in operating a heater inside a producing well and, at the same time, producing oil from the well. In this paper a simple theoretical model is analyzed to study the effect of such heating on the recovery of oil. The reservoir mechanics of the process is the only aspect of the problem which is considered here. The details of the heating system are ignored, and it is simply assumed that the heater raises the temperature in the immediate vicinity of the producing well to a pre-assigned temperature T w. Although heating requirements are calculated, these refer only to the heat supplied to the formation and to the oil being produced; they do not include heat losses due to the inefficiency of the heater. A short outline of the theory follows: The heat supplied to the producing well flows out into the formation by means of thermal conduction. At the same time, the flow of oil carries heat by convection from the formation towards the producing well. After sufficient time has elapsed, a steady-state is reached where the temperature, the velocity of the oil, etc., no longer vary with time but are functions of position only. In the steady state the temperature decays like an inverse power from the high value Tw at the wellbore to the original reservoir temperature T. at some distance away from the well. The effect of the high temperature and low oil viscosity near the producing well is the same as if the oil were produced at reservoir temperature To but through a producing well of a larger radius. From this point of view, the energy and money being expended to provide heat

Jan 1, 1958

By T. B. Massalski, C. S. Barett

THAT metals and alloys of the body-centered-cubic structure tend to become unstable at low temperatures is so nearly universal that any exceptions are worthy of special attention. Studies of the exceptions may disclose unrecognized examples of metastability or may shed light on the nature of some of the factors that stabilize the body-centered-cubic phase. Among the structurally analogous body-centered-cubic 3/2 electron compounds, the ß-phase in the Cu-Zn system deserves further attention for, although this phase is known to undergo partial mar-tensitic transformation on cooling below room temperature to the temperature of liquid nitrogen when the composition is less than about 42 wt pct Zn, it does not transform on cooling .when there is more than this amount of zinc present.1-4 The present authors recently found that an alloy containing 48.35 wt pct Zn could be transformed, however, if deformed severely at room temperature or lower temperatures.V his discovery led to the suggestion that the jerky flow of ß-brass single crystals6 and some of the unusual temperature variations of the yield strength and hardness are due to strain-induced transformation during the tensile tests. The present investigation is an extension of that work to cover the entire available range of compositions of Cu-Zn ß-phase alloys and certain interesting compositions of the ß-phase in the Cu-Zn-Ga system. Particular attention is given to the structure of the transformation product resulting from cold work, which is found to vary with alloy composition, and to the temperature range in which transformation can be induced in the various alloys. The composition limits of the body-centered-cubic phase in the Cu-Zn system, based on a recent review by Raynor7 and the precision determinations of Beck and Smith: are shown in Fig. la; in the ß region the alloys are disordered and in the ß region they are ordered. The rate of ordering is so rapid that a high degree of order is found in any of these alloys immediately after quenching. Experimental The nine binary and four ternary alloys of this investigation were prepared from cathode copper (99.99 pct purity), electrolytic zinc (99.99+ pct), and gallium (99.98+ pct). All alloys were cast in transparent Vycor tubing, homogenized at temperatures between 800° and 890°C, and quenched in water. Ternary alloys were cast in small amounts (-5 g) to predetermined compositions to fall along a line of constant electron concentration in the ternary system. Changes of weight during the casting operation of the ternary alloys were not greater than 1 part in 2500 from the intended nominal weights and therefore their compositions were accepted without chemical analyses. All binary alloys were cast in larger amounts (100 to 150 g), homogenized, and hot rolled after preheating to 800°C. The central portion of the rolled strip of each alloy was annealed for eight days and quenched in brine. All binary alloys were analyzed chemically. The compositions, heat treatments, and metallographic appearances of all alloys are tabulated in Table I. Various samples were subjected to heat treatment or cold work at several temperatures and were then examined metallographically or by means of X-rays. Specimens about 6x6x3 mm in size were examined metallographic ally to confirm the presence of the ß-phase after quenching. Cold work was introduced by fastening each sample between two pieces of steel, cooling the assembly to the desired temperature, and transfering it rapidly to a vise where moderate or severe compression was immediately applied. Cold working at liquid helium temperature was done by a chisel-shaped tool mounted on an X-ray spectrometer, described elsewhere.% A!! metal-

Jan 1, 1958

By F. A. Crossley

The additions: o, Sn, Zr, Cb, Mo, and V were studied for their effects on the a/a + Ti3A1 boundary of Ti-A1 base ternary alloys. These additions were chosen because of their importance to commercial titanium alloy practice—oxygen as the major impurity (in terms of influence on mechanical properties) and the others as major alloying additions. The composition and temperature ranges of investigation varied according to the system as shown in Table I. The principal means of investigation was light-microscopy. The basis of the ternary studies was the recently determined titanium-rich end of the Ti-A1 system.' MATERIALS, METHODS, AND PROCEDURES Ingots were prepared in 15 g size for all alloys except those of the Ti-Al-Mo and Ti-A1-V systems which were prepared in 20 g size. These alloys were prepared from high-purity, electrolytically refined titanium obtained from the Bureau of Mines, Boulder City Metallurgical Laboratory, Boulder City, Col., and Chicago Development Corp., Riverdale, Md., for the 15 g and 20 g ingots, respectively. The alloying additions were obtained from commercial sources, and were at least of 99.9 pct purity, except the V-15 wt pct A1 master alloy for vanadium additions. This alloy was of 98.2 pct purity; the major impurities were 0.37 wt pct Fe and 0.38 wt pct Si. The ingots were prepared by multiple melting in a nonconsum-able-electrode arc furnace. The ingots were covered with Markal coating "CR-T22" (Markal Co., Chicago) to protect them from scaling. Markal coating is applied as a slurry of solids in a volatile carrier liquid. The coated article is air dried. Upon heating to the hot working temperature the coating fuses to form a glassy, protective film. The ingots were heated to 1000°C and rolled to 0.16 in. thickness, reheating after each pass. The glassy film was removed by sand blasting. Ten mils were removed from all surfaces by grinding and/or pickling (3 pct HF, 30 pct HN03, balance water). Grinding was always followed by pickling to remove traces of the grindmg medium. Hydrogen removal was accomplished by vacuum annealing at 1000°C for 2 hr to a final pressure of about 0.02 µ Hg. This was followed by annealing in evacuated Vycor bulbs (800°C, 2 hr, ice brine quenched (IBQ)) in order to maximize ductility. The materials were then lightly pickled and cold rolled to 0.12 in. thickness or cracking, whichever occurred first. After rolling, 1/2 by 1/4 in. specimens were cut for the equilibration anneals given in Table 11 and carried out as previously described.1 Final polishing of the specimens for metallographic examination was done by electrolytic means. Specimens were first examined without etching under bright field and polarized light illumination, for this is the only means of detecting the "two-phase" phenomenon.1,2 The term "two-phase syndrome" refers to Ti-A1 al-loys which appear on the bases of their chemical, phys-ical, and mechanical properties, to be two-phase, but are actually single phase and inhomogeneous in composition. Since processing and heat treatment of the alloys of the first group were completed before the two-phase syndrome was discovered, a number of the specimens were found to exhibit this phenomenon—80 pct of those of the Ti-Al-Zr system, 5 pct of those of the Ti-Al-Sn system, 30 pct of those of the Ti-Al-Zr system. The affected specimens were homogenized by heating above the ß transus temperature and quenching. The heat treatment and working schedules to accomplish this are given in Table 111. After rolling, specimens were pickled (3 pctHF, 30 pct HNO3 in water) to remove 1 to 2 mils from the

Jan 1, 1970

By M. A. Grossman

The harden ability of most steels can be predicted within 10 to I5 per cent provided the complete chemical composition is known, including "incidental" elements; and provided the as-quenched grain size is known; and provided, finally, that the composition and heating temperatures for hardening are such as to result in austenite free from carbide particles. In the method proposed herein, the steel is considered as having a base hardenability due to its carbon content alone (the hardenability of a "pure steel" of the given carbon content, without any other elements), and this base hardenability is multiplied by a multiplying factor for each chemical element present. After multiplying all these together, the final product is the hardenability. Grain size may be taken into account either in the base hardenability or after the multiplication. Hardenability is stated in terms of "ideal critical diameter", namely, the diameter of bar, in inches, that will just harden all the way through (absence of un-hardened core) in an "ideal" quench, and the calculation may also be related to the Jominy test The data bring out certain features rather clearly. For example: i. It is quite useless to attempt to predict hardenability unless all elements, including "incidentals,,) are known; thus an "incidental,, chromium content of 0.20 per cent increases the hardenability by about 50 per cent. 2. The relative effectiveness of different alloys is sometimes unexpected; thus molybdenum when calculated in this way appears to be of the same order of effectiveness as manganese, rather than much more powerful, as would appear to be common experience; observe that an increase from no manaanese to 0.20 per cent Mn provides a multiplying factor of 1.67, and an increase from no molybdenum to 0.20 per cent Mo provides a factor of 1.63, an increase of over 60 per cent in each case. However, if to a steel containing 0.50 per cent Mn there is added "20 points of manganese," the factor is raised from 2.65 (for 0.50 per cent Mn) to 3.35 (for 0.70 per cent Mn), or an increase of only 26 per cent. Thus the fist small addition of an element has a much more powerful percentage effect than an equal further addition when some is already present, and in most cases the effect of molybdenum is considered in relation to a steel in which molybdenum is absent. 3. If two elements are equally effective, a greater hardenability will be obtained by using, for example, 0.5 per cent of each than by using I.o per cent of either of them alone. 4. A knowledge of the as-quenched grain size is essential for precise work, since a difference of only one grain size number (say No. 7 instead of No. 6) makes a difference of almost 10 per cent in hardenability (in these units); this, however, does not apply to certain It should be emphasized that in chromium steels (over 0.30 per cent Cr) and chrome-molybdenum and chrome-vanadium steels, un-dissolved carbides are likely to be present in the steel as quenched, and that in such cases the charts can indicate only a maximum Passible hardenability, whereas the extent of hardening actually obtained may be much less. Thus tests on a number of chrome-molybdenum steels have indicated a degree of hardening amounting to only 50 to 65 per cent of the maximum possible, and in chromium steels from full hardening (100 per cent) down to as low as 70 per cent. On the other hand, when the amount of such elements is small (Cr under 0.30 per cent, Mo up to 0.25 per cent in the

Jan 1, 1942

By M. A. Grossman

The harden ability of most steels can be predicted within 10 to I5 per cent provided the complete chemical composition is known, including "incidental" elements; and provided the as-quenched grain size is known; and provided, finally, that the composition and heating temperatures for hardening are such as to result in austenite free from carbide particles. In the method proposed herein, the steel is considered as having a base hardenability due to its carbon content alone (the hardenability of a "pure steel" of the given carbon content, without any other elements), and this base hardenability is multiplied by a multiplying factor for each chemical element present. After multiplying all these together, the final product is the hardenability. Grain size may be taken into account either in the base hardenability or after the multiplication. Hardenability is stated in terms of "ideal critical diameter", namely, the diameter of bar, in inches, that will just harden all the way through (absence of un-hardened core) in an "ideal" quench, and the calculation may also be related to the Jominy test The data bring out certain features rather clearly. For example: i. It is quite useless to attempt to predict hardenability unless all elements, including "incidentals,,) are known; thus an "incidental,, chromium content of 0.20 per cent increases the hardenability by about 50 per cent. 2. The relative effectiveness of different alloys is sometimes unexpected; thus molybdenum when calculated in this way appears to be of the same order of effectiveness as manganese, rather than much more powerful, as would appear to be common experience; observe that an increase from no manaanese to 0.20 per cent Mn provides a multiplying factor of 1.67, and an increase from no molybdenum to 0.20 per cent Mo provides a factor of 1.63, an increase of over 60 per cent in each case. However, if to a steel containing 0.50 per cent Mn there is added "20 points of manganese," the factor is raised from 2.65 (for 0.50 per cent Mn) to 3.35 (for 0.70 per cent Mn), or an increase of only 26 per cent. Thus the fist small addition of an element has a much more powerful percentage effect than an equal further addition when some is already present, and in most cases the effect of molybdenum is considered in relation to a steel in which molybdenum is absent. 3. If two elements are equally effective, a greater hardenability will be obtained by using, for example, 0.5 per cent of each than by using I.o per cent of either of them alone. 4. A knowledge of the as-quenched grain size is essential for precise work, since a difference of only one grain size number (say No. 7 instead of No. 6) makes a difference of almost 10 per cent in hardenability (in these units); this, however, does not apply to certain It should be emphasized that in chromium steels (over 0.30 per cent Cr) and chrome-molybdenum and chrome-vanadium steels, un-dissolved carbides are likely to be present in the steel as quenched, and that in such cases the charts can indicate only a maximum Passible hardenability, whereas the extent of hardening actually obtained may be much less. Thus tests on a number of chrome-molybdenum steels have indicated a degree of hardening amounting to only 50 to 65 per cent of the maximum possible, and in chromium steels from full hardening (100 per cent) down to as low as 70 per cent. On the other hand, when the amount of such elements is small (Cr under 0.30 per cent, Mo up to 0.25 per cent in the

Jan 1, 1942

By G. H. Enzian

THE effects of phosphorus and nitrogen on the properties and behavior of low carbon steels are important considerations to both the manufacturer and the user of such material. For one thing, these two elements are the principal points of distinction between the chemical compositions of Bessemer and open-hearth grades. Further, open-hearth steels themselves may vary considerably in nitrogen content depending on the raw materials and ..practice used. In addition, appreciable amounts of rephosphorized open-hearth steels are being used at the present time in some low alloy, high strength grades and in certain other cases. For these and other reasons, it was felt desirable to obtain more complete knowledge of the effects of simultaneous small variations in phosphorus and nitrogen content on such properties of low carbon steel as tensile strength, impact toughness, sensitivity to cold work, and strain aging. A long range investigational program designed to furnish such information, was carried out, and the present paper summarizes some of the data currently obtained. General Discussion The strengthening effect of relatively small amounts of phosphorus in steel has been recognized, and to some extent utilized, for many years. A number of investigators have pointed out, from time to time, that under proper conditions of control, phosphorus can be a useful alloying element in low carbon steels.1-3 Nevertheless, in spite of the volume of data on the subject, the use of phosphorus for this purpose has always been handicapped by widespread feeling that its presence makes the steel brittle, difficult to cold work, and unreliable at low temperatures. This belief on the part of some steel users is not entirely justified in many cases, and often stems from early experience with acid Bessemer steels which have appreciably higher phosphorus (and nitrogen) contents than basic open-hearth grades. Nitrogen in steel, and its effect on properties, has also received considerable attention in the literature.4 In 1906, Braune5 ascribed the brittleness of Bessemer rails to the nitrogen content, and proposed a schedule of nitrogen limits for steels for various uses. This was one of the first papers to recommend directly the use of nitrogen content as a quality control measure, but it did not meet with much favor or support; as a matter of fact, the Jernkon-toret was rather sharply criticized for publishing the paper. In the intervening years a number of papers have been published on the subject of nitrogen in steel. As a result, it is now generally recognized that nitrogen can have an important effect on the properties and behavior of steel and that the presence of this element should be considered in any systematic study of factors affecting steel quality. In 1932, Graham6 introduced the concept that the quality of various steels could conveniently be evaluated and compared in terms of "sensitivity", which was defined as "the characteristic reaction rate of a piece of steel toward brittleness or hardness." This concept formed the basis for a broad research program into factors causing variations in the quality of commercial steels. At an early stage of this work, it became apparent that the observed effect of either phosphorus or nitrogen on the sensitivity of steel depended, in no small measure, on how much of the other of these two elements was also present. Indeed, the combined effect of phosphorus and nitrogen in steel may often be more important than the specific effect of either. With a few exceptions, it is not possible to predict accurately the effect on steel properties of a change in the amount of either phosphorus or nitrogen in the steel without also specifying the approximate amount of the other element. It seems desirable, at this point, to discuss briefly the distinction between strain-sensitivity and strain aging as used in this paper. In the early stages of its development, the work-brittleness test7 (which will be described in more detail later) was used

Jan 1, 1951

By V. A. Josendal, J. W. Wilson, B. B. Sandiford

Two radioactive tracers have been tested as a means of determining core saturation in multiphase flow studies. Cesium chloride was tried as a water-phase tracer, but complications in its use in low permeability cores resulted from sorption of cesium by the core or water-wet pads. IodoI3'ben-zene proved very satisfactory as an oil-phase tracer. The synthesis of iodobenzene from the sodium iodide as received from Oak Ridge is simple and direct. The tracer is insoluble in water and there was no evidence of sorption by any of the core materials used. use of the method to determine saturation profiles during capillary and dynamic desaturations and relative permeability measurements on oil-water and oil-gas systems is described. Comparisons of the dynamic and capillary methods of relative permeability determination were made using the tracer to check core saturation and saturation distribution. Other experiments are also described in which mobility of the oil phase at various saturations was measured by displacing labeled oil by flowing inactive oil. Similar experiments were made using water labeled with cesium 134. INTRODUCTION To be fundamentally sound. any laboratory method of determining relative permeability must meet, among others, the following requirements: (1) the core saturation between the pressure taps must be uniform, and (2) there must be uniform pressure difference between the phases in the region between the pressure taps. In practice it may eventually be shown that appreciable departure from either or both of these conditions may be permitted without introducing sensible error. but the magnitude of such error remains to be established. There is need then for methods of saturation determina- tion which will determine the saturation profile in a core as well as the overall average saturation given by the gravimetric or material balance methods usually employed. Methods which have received consideration include resistivity, X-ray absorption, gamma ray absorption, neutron diffraction, and radioactive tracers. The resistivity method has use in the determination of profile, but needs checking by independent means.' The X-ray method has been used successfully, but requires elaborate equipment and calibration.23 4 At present gamma ray absorption is practical only for very large cores. The neutron methodo does not appear adaptable to measurement of saturation variations in short cores since it is a scattering method. Russell, Morgan and Muskat' employed radiovanadium in a study of the mobility of interstitial water. Coomber and Tiratsook used radioiodine as an oil-phase tracer and measured profiles in unconsolidated sand packs. ADVANTAGES OF RADIOACTIVITY METHOD A thorough discussion of the theory of radioactive tracers and the technique of using them to study the movement of fluids in sands is given by Coomber and Tiratsoos and need not be detailed here. The method requires less equipment than the X-ray method; the counting problem is similar, but the generation and regulation of the radiation, which is a major problem in the X-ray method, is not required in the radioactivity method. In the X-ray method it is customary to add up to 20 per cent of absorber, usually an iodine compound, to the phase to be labeled. while in the radioactivity method a mere trace is sufficient. A possible disadvantage of the radioactivity method is that a disproportionate amount of the activity comes from the portion of the core nearest the counter, whereas the X-ray gives a true average saturation in the region traversed by the beam. 011 the other hand, this characteristic of the radioactivity method renders it peculiarly adaptable to study of certain types of radial variations in saturation. As with the X-ray method, a tracer may be added to either the oil or water-phase. An oil-phase tracer is more generally useful. however. since it permits work on oil and gas in the

Jan 1, 1952

By F. A. Schaufelberger

METAL can be recovered from a leach solution either indirectly by precipitation as a compound that is later reduced or directly by electrolysis, cementation, or chemical reduction, for example, with hydrogen. The widely used electrowinning processes suffer from the inefficiency of converting fuel to low voltage dc and from low current efficiency when complete recovery is attempted. With the batch hydrogen reduction methods outlined here, a unit of fuel converted to hydrogen in modern gas reform plants will make two to three times as much metal as can be obtained with electrolysis. With continuous hydrogen reduction, four to six times as much metal is obtained. Precipitation by reducing a metal salt solution with H, has been known for almost 100 years.' Commercial use of the process, however, awaited the research and development program initiated by Chemical Construction Corp. Within the last few years this program has brought about construction of commercial plants, listed on page 545. Ideal hydrogen reduction will precipitate a pure metal from a solution obtained by commercial leaching methods at a rapid rate without excessive temperatures and pressures. The metal precipitate will be of desired size and density, less than a percent left in the vessel being deposited on the wetted parts. The present discussion will outline Chemical Construction Corp.'s early development program and will discuss the chemistry and mechanics of reducing copper, nickel, cobalt, and cadmium from solution by H2. Work on selective reduction of nickel from cobalt has been described earlier.' The article is not concerned with precipitation of copper by disproportionation of cuprous solutions where yield is limited3 to 50 pet.* It does not discuss use of gaseous reducing agents other than hydrogen, such as SO2,4 which may lead to contamination with sulfur, or C0,5 which is considerably more expensive than hydrogen and produces a gaseous reaction product, CO2.† The article does not include use of other nongaseous reducing agents,7 which are far more expensive than CO. Also, notes on the history, engineering design, and performance of commercial plants have been given elsewhere.' It will be shown that a pure metal can be precipitated by reduction with hydrogen from solutions obtained by commercial leaching methods when the solution composition is controlled and proper acidity, proper metal ion concentration through controlled complex formation and hydrolysis, and adequate agitation for suspension of metal and for transfer of reducing gas are maintained. Reasonable temperatures and pressures are used which are less excessive than those used by previous workers, although they are maintained at above equilibrium values because of process kinetics.' It has also been found necessary to induce nuclea-tion and to control growth and agglomeration in order to get a powder product of suitable physical properties. Review of the Literature Muller, Schlecht, and Schubardt of I. G. Farben10 claimed a successive reduction of silver, copper, nickel, cobalt, and zinc from ammoniacal solution by applying progressively higher temperatures and hydrogen pressures. Metals produced by this technique were not pure, since no attempt was made to adjust the solution composition between the different reductions. It is doubtful that the zinc product contained any metallic zinc. The major contribution to the literature" was made by the Ipatiews,12 who prepared platinum, iridium, copper, nickel, cobalt, lead, tin, arsenic, antimony, and bismuth at often extreme conditions, up to the critical temperature of water (373ºC), at 1500 to 7500 psi, for as long as several days. Even cadmium and zinc were claimed to have been observed in trace quantities precipitated together with basic salts. In the early experiments solutions were not stirred. Extensive formation of oxides and basic salts occurred, probably long before sufficient hydrogen could be supplied in the unstirred liquid to decrease the metal concentration in solution by reduction. The importance of hydrogen pressure was

Jan 1, 1957

By M. van Swaay, C. E. Birchenall

A. S. lBrling (Johnson, Matthey & Co. Ltd., Laboratories)— Because of its initial emphasis upon the production of membranes by powder metallurgy and by normal casting and rolling techniques, this paper is of considerable interest to those concerned with the design and application of industrial diffusion units. As the subsequent performances of the three grades of palladium were not compared in detail, it can perhaps be concluded that specific permeabilities were very similar. A short table in which the absolute diffusion rates obtained by the Authors with their three grades of palladium were compared with those of other investigators18-25 would, however, have been of particular value in view of the very variable results obtained with this membrane material. The Authors state early in the paper that the permeability of the thoriated samples decreased in much the same way as that of the nonthoriated samples. Later, they emphasize that sintered palladium without thoria was very susceptible to poisoning and did not respond as well as the other materials to the decontamination treatment. These results, which suggest that sintered membranes are less effective than those produced by normal methods, have considerable industrial implications, and amplified experimental details would be very welcome. The results of the Authors differ in one important respect from those of other investigators. Thus, no "bleeding" was required to maintain a constant rate of diffusion, although Davis, Darling, and De Rosset found this to be necessary. Davisz3 found that "bleeding" was necessary even when the hydrogen had already been purified by diffusion through a palladium tube, and it would be of value to know the purity of the gas stream emerging from the purification train shown on Fig. 1 of the Authors' paper. The Authors show that diffusion rates varied according to the square root of the pressure difference providing that the palladium membrane was not contaminated in any way. This finding is in agreement with the results of those workers2,18,19,23,24 who diffused their hydrogen through the membrane into a partly or completely evacuated system.Those workers who did not find this square-root law dependence did not have a partial vacuum on the downstream-side of the membrane.'' The difference between the two experimental con- ditions appears to be significant, as the absolute diffusion rates per sq cm of diffusion area are lower when the square-root law dependence is not observed. While the effect could be attributable to poisoning of the membrane it is conceivable that a positive pressure on the downstream side of the membrane might influence the mechanism of diffusion. In view of the capricious behavior of palladium, factors of this type, which are likely to influence its industrial significance, should be carefully considered. M. van Swaay and C. E. Birchenall (authors' 9-eply) — Our Eq. 191 may be compared with several others for hydrogen diffusion in palladium: D = 5.95 x exp (-5770/1)' D = 3.55 x 10"3 exp (-5550/)' D = 4.3 x 10"3 exp (-5620/1')' Other values are given in the literature, but in many cases surface poisoning may have affected the results. No analysis was performed on the hydrogen gas stream after purification. It should be noted that the presence of liquid nitrogen traps near the heated membranes may have set up convection currents and provided the equivalent of "bleeding." The question of whether a partial vacuum existed on the downstream side of the membrane is not clear. Does it imply that there is something unique about the permeability measured near one atmosphere pressure? We do not believe that there is. When square root pressure dependence is not observed, diffusion rates are not measured. Is "absolute diffusion rate per sq. cm. of diffusion area" equivalent to permeability? The terminology is new to us. In any case we doubt very much that pressures of the order of several atmospheres can affect the mechanism of diffusion.

Jan 1, 1962

By C. G. Dunn, F. W. Daniels, M. J. Bolton

IN recent investigations1. a data on relative grain boundary energies in silicon iron have been obtained. The present investigation is a continuation of this work along similar lines for the purpose of obtaining additional information on boundary energies. One aim has been to extend the data for a (110) series, which was formerly studied,' to small differences in orientation, since a recent theory3,4 based on a dislocation model of grain boundaries predicts a E. ? (A-ln ? )* form of energy curve for small differences in orientation. Another aim has been to provide data for a (100) series—all grains having (100) planes parallel with the surface of sheet speci-mens—and to compare the results with the predicted form of energy curve. Finally, information on the nature of approach to equilibrium conditions was to be obtained through observations on the movements of grain boundaries. Recently published papers2,3,5 have treated the mathematical problem of expressing equilibrium conditions for grain boundaries for the general case when boundary energy depends upon boundary orientation. The equilibrium equations, which relate equilibrium angles and grain boundary -free energies, contain additional terms that express the variation of energy with boundary rotation. It has been necessary in the present investigation, however, to neglect these additional terms since no data for their evaluation were available. Dropping the additional terms leads to the following approximate equations which were used: E12/sin ?3 = E18/sin?2 =E28/sin?1 or those of the form: E28 + E13 cos ?3 + E12 cos ?2 = 0 Experimental Procedure Except for sample F1, all specimens were made from two lots of silicon iron called C and L respectively. The compositions of these are listed in table I. In the preparation of 12 three-grain specimens for the present investigation, a controlled grain growth technique, which has been described previously,' was used. After preparation, the specimens were notched at the boundaries to shorten them as a means of providing more rapid approach to equilibrium in the anneals. Initial grain boundary angles were determined from micrographs taken at X500. Seven anneals, totaling about two days at 1300°C and two to four days at 1400°C, were run in an atmosphere of pure dry argon as described in a previous paper.' Boundary changes generally could be followed without metallographic surface preparation, because thermal etching occurred during the anneals. These observations indicated angle changes of as much as 35" in some instances, with the major changes occurring during the first anneal. Annealing was discontinued when the angles reached stationary values. Fig. 1 gives a schematic diagram for purposes of defining crystallographic and boundary directions for three grains with common plane parallel with the sheet. Differences in orientations (expressed by A) are the angles between [00l] directions. They can be calculated with the aid of the ?'s. The grain boundary angles can be calculated from the w'S. Results Results are given in table 11. Types of crystal-lographic planes refer to those parallel with the sheet. Orientations are given in terms of [00l] direc-tions as shown in fig. 1. Except for grain 3 of specimen H6, which was unintentionally tilted 7" out of a (100) plane, all grains generally were within 1" or 2" of the plane specified. The directions of boundaries were measured generally to within 1" on both upper and lower surfaces, and the average direc- tions were determined for use in calculating the grain boundary angles. Variations in grain boundary directions refer to deviations from average directions. Boundaries near the junction point for all specimens generally were very close to perpendicu-

Jan 1, 1951