Search Documents

By AIME

BORN in 1899, in Bridgeport, Conn., E. A. Anderson grew up in a center of the nonferrous metal industry. Perhaps that had something to do with his selection of mining as a career while an undergraduate in Yale. After securing his Ph.B. in mining he worked for a year at Handy & Harman. By this time physical metallurgy had supplanted his interest in mining and he entered the Yale graduate school where he secured his M.S. degree in metallography under C. H. Mathewson. After completing his graduate work he entered the employ of The New Jersey Zinc Co. in their research division where he has remained ever since. The transition from graduate study to industry involved no disruption of technical interests. Academic research relating to the zinc-aluminum alloys, then coming into use for die-casting, was followed by extensive and productive industrial research in the perfection and commercialization of these alloys. Research on rolled zinc and rolled zinc alloys, corrosion and stress corrosion, dry cells, electroplating on zinc, and electrodeposition of zinc followed in overlapping succession as time passed.

Jan 1, 1947

By J. D. Sudbury, J. E. Landers, D. A. Shock

The external corrosion of casing is one of the most important problems facing the present day production man. A 1953 NACE report' estimated the annual cost of casing corrosion at $2,200,000 in 22 fields studied. The fields represented only 10 per cent of the wells in the United States. It was significant that 77 per cent of the failures reported were due to external corrosion. A variety of means for combatting the casing corrosion problem have been discussed in papers such as the NACE survey' mentioned previously. No attempt to revicw these will be made in this paper. Thc scope of this paper has been limited to a discussion of the following major items: (I) an instrument for studying external electrolytic corrosion; (2) use of this instrument to determine the quantity of protective current required to prevent corrosion of casing; and (3) problem of accelerated external corrosion of the oil string at the base of the surface string. ELECTROCHEMICAL MECHANISM OF CORROSION It IS gcncrally belleved that one major cause of external corrosion is the existence of electrolytic cells operating between various formation waters. A theoretical discussion of the role of earth potentials in casing corrosion has been presented by L. de Witte and F. J. Radd.2 A brief picture of the electrolytic corrosion action between two differcnt formations (A and B) is presented in Fig. I. This picture is a review to many but is given to insure clarity. In the anode area, oxidation takes place as iron ions go into solution in formation A as a product of corrosion. In the cathode area. reduction takes place and hydrogen atoms are formed at the boundary between the pipe and formation B. The flow of electrons. as indicated. is into the pipe at the anode or corroding area and out of the pipe into the formation at the cathode. In order to conform with regulations of this journal, the flow of a conventional current (I! is also indicated in this sketch. All subsequent references will be to this so-called conventional current flow and will be indicated by the symbol (I). As this action takes place the layer of hydrogen developed on the cathode offers a resistance to the flow of current in the cell. It is said (hat the cathode has been polarized. This polarizing film can be rcmoved by addition of a so-called dcpolarizing agent which reacts with and removes this hydrogen film. Such an agent must be present in ordcr for this corrosion reaction to go at a rate sufficient to cause a severe corrosion problem. The two common depolarizing agents are: (I) oxy-sen. which readily reacts with hydro gen to form water, and (2) certain bacteria which use atomic hydrogen as part of their life process. It is not difficult to explain the presence of these agents. Oxygen may be present in some very shallow formations. Sulfate reducing bacteria are found to exist in nearly any ground where mud pits might be dug. THE CASING POTENTIAl-PROFILE Instrumentation A sketch of the equipment* uscd in determining casing potential profiles is shown in Fig. 2. The tool is run in the hole using a hoist truck with a three-conductor cable. This tool must be run with tubing and rods out of the hole and should be run with the flow line disconnected. By proper choice of electrode size it can be run in any size casing. The electrodes must be run in an open hole or a hole filled with oil. Measurements cannot be made while immersed in water. The tungsten carbide-tipped knives serve as the electrical robes and arc spaced 25 ft apart. The potential difference along a 25-ft section of casing is measured with a potentinmeter to an accuracy of about 5 mv. Current flowing in the pipe can then

Jan 1, 1956

By E. L. Claridge, K. S. Lee

Areal sweep efficiency of oil displacement by enhanced-viscosity water exhibiting pseudoplastic behavior was measured in a Hele-Shaw model representing one-quarter of a five-spot pattern. The pseudoplasticity of polymer solutions and the velocity distribution in the five-spot pattern produced a condition under which the mobility ratio between the displacing and the displaced fluid could not be assigned a single value. Instead, the movement of the displacement front is governed by local mobility ratios which are also time dependent. The areal sweep at breakthrough with polymer solutions was poorer than the sweep obtained with Newtonian fluids of comparable viscosity. However, the areal sweep at 1 PV throughput was greatly improved as compared to flood water without polymer. It was also demonstrated that, even after the oil-cut had declined to a low value during a regular waterflood, switching to polymer flood efficiently swept out the oil remaining in the model. INTRODUCTION The behavior of fluid displacements in isotropic porous media for various patterns of injection and production wells has been extensively investigated.1-7 These investigations all concerned Newtonian fluids, i.e., the viscosity of each fluid was constant regardless of flow rate. The generally unfavorable influence on areal sweep efficiency of higher mobility of the displacing fluid as compared to the mobility of the displaced fluid has been established for both miscible and immiscible fluids. The principle was also established that a close correspondence exists between miscible and immiscible flood front behavior,2, 3 although oil recovery in a waterflood at unfavorable mobility ratio may be less than that observed in a miscible displacement at the same mobility ratio. This is true even when oil recovery is expressed on the basis of movable oil. The reason is that oil saturation only slowly achieves its final value behind the water-flood front in accordance with the Buckley-Leverett8 simultaneous flow relations. It is convenient to use miscible displacements for laboratory simulation of waterflood frontal advance since the interfacial tension forces which are negligible in proportion to viscous forces on a reservoir scale9 are thus made nonoperative in the laboratory model. For miscible displacements, the Hele-Shaw type of model adequately represents a porous medium so long as the appropriate scaling rules5-10 are observed in its design and operation. During simulation of waterflood front behavior in the laboratory by using miscible displacements, the behavior of connate water may ordinarily be disregarded since it is usually indistinguishable from flood water in this process. However, when the flood water is deliberately thickened to improve the mobility ratio between water and oil. the effect on the sweep efficiency due to generation of a connate water bank during the process must be considered. In a uniform porous medium, such a bank is generated11 and efficiently displaced by injection of thickened water.l2 The oil originally in-place at the start of the waterflood is then displaced by connate water followed by thickened water. If the flood water must be thickened to obtain a favorable mobility ratio, the mobility of the oil phase is appreciably less than that of the connate water. Hence, the oil phase is inefficiently displaced by the connate water bank, and a considerable proportion of the oil comes in contact with and is displaced by the thickened waterflood front' A recent study 12 of five-spot displacements using miscible Newtonian fluids in which connate water banks of various sizes were simulated shows that, for connate water banks as large as 40 percent of 1 PV, major improvements in sweep efficiency resulted from use of a thickened

Jan 1, 1969

By I. Javandel, P. A. Witherspoon

The finite element method was originally developed in the aircraft industry to handle problems of stress distribution in complex airframe configurations. This paper describes how the method can be extended to problems of transient flow in porous media. In this approach, the continuum is replaced by a system of finite elements. By employing the variational principle, one can obtain time dependent solutions for the potential at every point in the system by minimizing a potential energy functional. The theory of the method is reviewed. To demonstrate its validity, nonsteady-state results obtained by the finite clement method are compared with those of typical boundary value problems for which rigorous analytical solutions are available. To demonstrate the power of this approach, solutions for the more complex problem of transient flow in layered systems with crossflow are also presented. The generality of this approach with respect to arbitrary boundary conditions and changes in rock properties provides a new method of handling problems of fluid flow in complex systems. INTRODUCTION Problems of transient flow in porous media often can be handled by the methods of analytical mathematics as long as the geometry or properties of the flow system do not become too complex. When the analytical approach becomes intractable, it is customary to resort to numerical methods, and a great variety of problems have been handled in this manner. One such method relies on the finite difference approach wherein the system is divided into a network of elements, and a finite difference equation for the flow into and out of each element is developed. The solution of the resulting set of equations usually requires a high speed computer. When heterogeneous systems of arbitrary geometry must be considered, however, this approach is sometimes difficult to apply and may require large amounts of computer time. The finite element method is a new approach that avoids these difficulties. It was developed originally in the aircraft industry to provide a refined solution for stress distributions in extremely complex airframe configurations. 27 Clough has recently reviewed the application of the finite element method in the field of structural mechanics.6 The technique has been applied successfully in the stress analysis of many complex structures.l, 27, 28 Recognition that this procedure can be interpreted in terms of variational procedures involving minimizing a potential energy functional7 leads naturally to its extension to other boundary value problems. In the field of heat flow, there recently have been introduced several approximate methods of solution that are based on variational principles.2-4, 17 By employing the variational principle in conjunction with the finite element idealization, a powerful solution technique is now available for determining the potential distribution within complex bodies of arbitrary geometry. In the finite element approximation of solids, the continuum is replaced by a system of elements. An approximate solution for the potential field within each element is assumed, and flux equilibrium equations are developed at a discrete number of points within the network of finite elements. For the case of steady-state heat flow, the technique is completely described by Zienkiewicz and Cheung.33 Since the flow of fluids in porous media is analogous to the flow of heat, Zienkiewicz et al. have employed the finite element method in obtaining steady-state solutions to heterogeneous and anisotropic seepage problems. 34 Taylor and Brown have used this method to investigate steady-state flow problems involving a free surface.25 The work of Gurtin has been instrumental in laying the groundwork for the application of finite element methods to linear initial-value problems.12 As a result, Wilson and Nickell have recently

Jan 1, 1969

By R. Lagneborg, R. Gullberg

An X-ray method for determinig volume fractions of phases in textured materials has been developed. The method involves measurements of the integrated intensities of reflections of each phase for one single orientation of the sample. The method was checked by evaluating the phase contents in samples from two stainless Cr-Ni steels containing austenite and ferrite in known amounts. In the present case measurements of about five reflections from each phase were required. The agreement between calculated and known volume fractions was found to be very close. AN X-ray diffraction method is often employed for determining the amounts of two or more phases in a material. The details of the method have been described at several instances.'" In short, the method is based upon the fact that the integrated intensity of a reflection is directly proportional to the volume fraction of the phase considered. In principle, the method only requires the measurement of the intensity of one reflection for each phase. The most severe shortcoming with this method is that it disregards influences from preferred orientation in the material. In a diffractometer, where this type of analysis is most conveniently carried out, the geometry is such that the diffracting planes for a specific Bragg peak are all parallel to the surface of the specimen. This implies that it must be required that the occurrence of the diffracting planes parallel to the specimen surface must be as frequent as for any arbitrarily oriented physical plane of the sample. This requirement will generally only be fulfilled for a material free of texture. Since all materials in reality always exhibit a preferred orientation to some extent. the deficiency of the method is quite serious. In cold-deformed material the texture is appreciable, but also annealed material will in many cases display a considerable texture. The problem of quantitative phase analysis in textured material is one of obtaining the integrated intensity averaged over all orientations of the specimen with respect to the X-ray beam. Lopata and ~ula~ have shown that this can be accomplished by integration of intensities over the whole pole figure, or by random rotation of the specimen in the X-ray beam, or by a combination of both. Hence, in this method the average integrated intensity for a specific reflection is obtained by directly determining the intensity for a sufficient number of physical orientations of the sample. However, since an excess of diffracting planes in a certain surface plane necessarily must be accompanied by a deficiency of another type of diffracting plane, it appears to be possible to obtain the correct average intensity for a specific diffracting plane by performing an appropriate averaging of intensities belonging to a large number of diffraction planes of different hkl when all intensities have been produced from one single orientation of the specimen surface. In the present work such a method has been attempted and tested experimentally for two stainless Cr-Ni steels containing austenite and ferrite. 1) THEORY The expression for the integrated intensity of a Bragg reflection diffracted from a phase in a polycrys-talline sample in a diffractometer is: where C is a constant only dependent on the incident X-ray power and the diffractometer geometry, v the volume of the unit cell, F the structure factor, m multiplicity of the diffracting plane, L Lorentz factor, P polarization factor, e-2 temperature factor, p average linear absorption coefficient of the irradiated vO1umej and V volume fraction of the phase 'Onsidered. Hence, the factor R in the abbreviated form of Eq. [I] can be calculated and depends only on 0 and the diffracting plane {hkl}. For an alloy containing two phases, a and y, we can write

Jan 1, 1967

By F. L. Toy, C. E. Sims

PNEUMATIC processes for converting molten pig iron to steel were the major producers of steel during the latter half of the 19th Century and until shortly after the turn of the century, when these processes began losing ground to their rival—the basic open hearth process. The first pneumatic process for converting molten pig iron to steel was introduced almost simultaneously by Sir Henry Bessemer in England and by William Kelley in the United States during the decade 1850-1860, but it was Bessemer who pushed it to a successful commercial operation. The process itself is so basically simple that it has not lent itself to many extensive or revolutionary improvements since its inception. During the first 30 yr, numerous' designs of vessel and many modifications of practice were tried before the industry standardized on something resembling current practice. The evolution in construction has resulted in proportions and arrangements based more on mechanical expedience than on metallurgical requirements. Many of the early converters were stationary and had to be tapped similarly to the cupola or open hearth furnace. Such types were used for many years in Sweden and Germany, but the rotating or tilting types were favored in England and the United States. There was evident trouble with bottom tuyeres, becAuse about half of the early designs were constructed for side blowing with tuyeres submerged below the metal-bath surface. Even portable tuyeres, introduced from above, were tried. The principal advance, which made bottom blowing practical, was the invention, in 1863, by Henry Bessemer, and subsequent improvement by A. L. Holley, of a detachable bottom so that the tuyere section could be renewed. In the meantime, it was the demonstration, by Robert Mushet, of the value of manganese in overcoming hot-shortness in converter steel, that made the process metallurgically feasible. At first, only the acid practice was used, and Swedish pig iron, because of its low sulphur and phosphorus, was in great demand. The development of a suitable basic lining, by Sidney G. Thomas and Percy G. Gilchrist in 1876-1879, made the process applicable to the high-phosphorus ores of Europe. From the beginning of the use of the Thomas-Gilchrist process, trouble was experienced in obtaining sufficiently high temperatures in the blown metal, particularly in small heats. It is necessary, of course, to have the low-carbon blown metal at a considerably higher temperature for handling from vessel to ladle to molds than is needed for the molten pig iron. All of the increase in temperature in a pneumatic process must be obtained from the heat of combustion of the elements contained in the pig iron plus that from silicon when ferrosilicon is added. In bottom blowing practice, carbon is almost completely burned, but inasmuch as it cannot be oxidized to a higher state of oxidation than to CO when in contact with liquid iron, only a minor portion of the heat-producing capacity is utilized. Under such conditions, it is not considered to be an efficient source of heat. In the acid practice, therefore, silicon and manganese are the large heat producers, are quickly oxidized, and are efficient fuels. Too much manganese, however, thins the slag greatly and tends to promote excessive slopping or throwing of slag and metal from the converter during the blow. Most operators, therefore, like to keep the manganese content of the iron below 0.60 pct, and thus the main source of heat is the silicon, which must be at least 1 pct and sometimes 2 to 3 pct for small blows such as for castings. Ferrosilicon is sometimes added during the blow to obtain an adequate temperature level. In the basic (Thomas) practice, even more heat is required, because cold lime is added to flux the phosphorus after it is oxidized, and the basic slag must be maintained at a proper temperature. Silicon

Jan 1, 1951

By A. Long, A. Silverman, J. L. Kulp

Garth Crosby has written an excellent geologic description of the deposits surrounding the Gem stocks in the Coeur d'Alene district. The ore deposits in the area of the Gem stocks may hold the key to the age and mineralization history of the entire Coeur d'Alene district. A recent isotopic study10 of the Coeur d'Alene district in this laboratory illustrates that a classical geologic interpretation of these deposits, with its implications for the whole district, may lead to erroneous conclusions. For this reason it seems pertinent to submit a summary of these results and an interpretation of the mineralizing history of the Coeur d'Alene districts. Pitchblende from the Sunshine mine has been dated as 1100-1200 million years (m.y.) by Kerr and Kulp,11 and Eckelmann and Kulp.12 This pitchblende is known to predate the major sulphide introduction in the Sunshine mine, but it clearly follows the metamorphism and folding of the Belt Series. Recently, Goldich (personal communication) dated sericite from the Belt rocks as 1200 m.y. by the potassium-argon method. This date establishes a minimum Belt age if the sericite is of metamorphic origin and a maximum age if detrital. The Belt Series, therefore, is not latest Precambrian. Lead isotope ratios of galenas have been measured and model ages calculated by Bate, et all3 for the Sunshine mine, Farquhar and Cummingl4 or the Bunker Hill and Sullivan mine, and Cannonl6,18 for a few other deposits in the district. Galenas from similar deposits in Southern British Columbia have been measured by Bate, et al,23 Farquhar and Gumming,11 and Wanless and Leech." The new work in this laboratory includes leads from ten widely scattered mines north and south of the Osburn fault, and through a vertical distance of over 5000 ft (Star mine). Included in this group are the Hercules, Sunset, and Mountain Goat mines, which are adjacent to the Gem stocks. In addition, galenas from stringers which cut the Gem stocks in the Success and Sunrise mines and the Hercules hanging wall dike in the lower western portion of the Hercules mine, were included in the study. All of the above leads, except one, give model ages of 1300-1400 m.y. for the time of segregation from the parent material. The concordance of the common lead and U-Pb ages, the great horizontal and vertical uniformity of the lead ratios, the sulphur ratios of two Sunshine mine galenas and three galenas from the Edith, Page, and Hercules mines indicating average crustal sulphur, and the U23/Pb204 ratio inferred for the original lead environment, indicate that the entire mineralization of the Coeur d'Alene-British Columbia base metal belt was derived from a deep source and emplaced in the Precambrian. One sample of galena from sulphide stringers cutting the North Gem stock at the Sunrise Prospect mine gave a modern (i.e., Tertiary) lead ratio. This mine is located at the northern tip of the stock. The presence of modern lead suggests mobilization of some rock lead during the intrusion of the Gem stocks and their outliers. Larsen et al.l8 have dated the Gem stocks as Tertiary by the lead alpha method. If Coeur d'Alene mineralization was genetically related to Tertiary magma formation and intrusion

Jan 1, 1961

By D. M. Boulin, E. S. Greiner

Single crystals of vanadium were bent at 293" and 77°K to a maximum fiber strain of about 1 pct. Analysis of slip line and etch markings of- crystals deformed at either temperature show the slip planes are {112}. The slip directions at 298° and 77°K, deterwined from the asterism axes of rotation of crystals compressed in a {110} direclion, were (111). Etch pits, associated with dislocations, were observed in the plastically deformed vanadium. There was no evidence of mechanical twinning. THIS paper presents results on single crystals of vanadium plastically deformed by bending and compressing at 298° and 77°K. Bcc metals have been observed to slip on {110}, (1121, and (123) planes while the slip direction has invariably been found to be ( 111) 1 Lindley and smallman2 have studied the deformation intension of polycrystalline vanadium between 298° and 20°K and reported that at and above 77°K the deformation is predominantly by slip, but at 20°K the deformation is associated with twinning. Later, Edington and smallman3 found that single crystals of vanadium compressed in a (111) direction at 77°K deformed by twinning on { 112}( 111) systems. Clough and pavlovic4 have reported that mechanical twins were observed in vanadium loaded by impact at 77°, 173°, and 195°K, respectively. PREPARATION OF SINGLE CRYSTALS The vanadium single crystals were made from commercial vanadium reported to contain 0.03 pct C, 0.03 pct 0, 0.02 pct N, 0.0X* pct Ca, 0.0X pct Si, proximately 0.8 cm diam, 4 cm long and were grown in (100) and (110) directions. A crystal for studies of etch pits was gound to have a face parallel to a {l00} plane by the Bond crystal-orienting technique.7 In order to avoid straining, the crystals were placed on the steel mount with a low-melting organic compound (glycol phthalate) before grinding. Crystals for the bending experiments were ground to have the four side faces parallel to ( 100) planes. The crystals were electropolished8 in 10 pct H2So4 using a current density of about 3 amp per sq cm to remove the strained surface material. The crystals were bent slowly by three-point loading in a stainless-steel jig designed by Sylwestrowicz.9 The end rollers of the jig were 3.8 cm apart. DEFORMATION AT 298°K The bent vanadium crystal, after removal from the jig, had a radius of curvature of 20.4 cm over the 2-cm length close to the center loading roller, where the maximum bending moment occurred. This radius of curvature corresponded to a maximum fiber strain of 0.9 pet. The side faces contained slip bands as shown in Fig. 1. These bands are straight where the strain is high, top and bottom of micrograph, but are curved on either side of the neutral axis in the center. This curvature of the slip bands in the low-stress area may be due to cross slip. The plastically deformed crystal was subsequently electropolished and then electroetched in 10 pct H2SO4 solution, using a current density of 0.25 amp per sq cm. The structure is shown in Fig. 2. The directions of the traces on two adjoining faces were plotted stereo graphically, Fig. 3.* Since there were two sets of

Jan 1, 1968

By W. A. Backofen, B. B. Hundy

THE pole figure interpretation in a recent paper' on the torsion texture of copper was questioned in discussion,2 nd a simpler alternative interpretation was proposed. In the hope of reaching a positive decision in favor of one or the other, further studies were undertaken with 70-30 brass and Armco iron. This additional work has not answered all the questions about torsion textures, but the new pole figures do offer support for the original interpretation and show that the texture is more complicated than the alternative suggests. From the previous paper,' the appropriate method of describing the texture in a twisted bar is known to be a pole figure plotted on a plane tangent to the surface of the bar. The axes of the pole figure are chosen to coincide with directions in the tangent plane which are parallel and at right angles to the bar axis (the longitudinal direction, L.D., and the transverse direction, T.D.), for these directions define symmetry axes of the texture. Therefore any ideal orientation that might be selected to represent the texture gives the crystallographic plane parallel to the tangent plane and the crystallographic direction parallel to the axis of the torsion specimen. A unique texture was assumed in the original interpretation, and the four ideal orientations, (112) [111], (111) [112], (110) [001], and (112) [131], were used in its description. Of the four orientations, the first three provide possible slip directions along the direction labeled T.D. in the pole figure, which is one of maximum shear stress in a specimen being twisted. The alternative, however, simply considered the texture to be made up of two components, one a [Ill] alignment along the specimen axis with a random orientation about the axis as in the wire texture, and the other a superimposed orientation (110) [0011. Experimental evidence for this simpler interpretation was derived from two observations: 1—twisting did not alter the orientation of a single-crystal torsion specimen of a brass in which an octahedral plane was perpendicular (and a [Ill] direction parallel) to the axis, and 2—deformation bands with approximately the (110) [001] orientation were seen to develop in the crystal during the test. The new pole figures were obtained with a brass and an iron torsion specimen of 3/4 and 1 in. gage length diameter, respectively. After annealing, polishing, and twisting, seven rods, about 0.030 in. in diameter and with different orientations, were prepared from material near the surface of each specimen; these were required for the method of constructing pole figures by using a Geiger counter as developed by Norton.5 Since the shear strain gradient in plastic torsion is linear, it was readily computed that the material supplying the brass rods had been subjected to a shear strain of about 3.2, and that supplying the iron rods had received a shear strain of about 2. All details concerning specimen design, testing procedure, rod preparation, and use of the Norton technique can be found in ref. 1. Only this procedure was followed in working out the (111) and (100) pole figures for brass, Figs. 1 and 2. However, X-ray photograms were also involved in the construction of the (110) pole figure for Armco iron, Fig. 3. The photograms were made according to the conventional transmission diffraction technique using a small sheet taken parallel to the specimen surface and etched to a thickness of about 0.0025 in. They were especially important in locating the boundaries

Jan 1, 1954

By H. C. Chang, N. J. Grant

REEP studies and measurements in most in-V> stances are based on a relatively gross gage length. Even in some recent theoretical studies on the mechanism of creep, changes were followed by means of X-rays and metallography without regard for the behavior of the individual grains and grain boundaries.'" Analytical studies of creep have frequently indicated that more than one competing process is operative and that the ultimate creep performance depends on the effects of temperature, time, and strain rate as well as grain size.'-" he grains and grain boundaries are known to respond differently to these variables but the extent of the response is still little known, especially in poly-crystalline materials where the relative orientations of grains and grain boundaries are important. This paper reports some experimental observations of the behavior of the grain boundaries and their effect on the deformation in the grains. High purity aluminum was used to avoid complications of aging reactions, second phase, inclusions, etc. Experimental Procedure A small vertical nichrome wound furnace was built to permit direct loading on the test specimen. This furnace is capable of operating up to about 1300°F with reasonable temperature distribution. A ground rectangular quartz window, 1.25x0.5 in., is imbedded in the external stainless steel tube to permit observations of the creep process with a microscope. An argon atmosphere ,was provided to prevent oxidation. The quartz window caused a temperature gradient between the back and front surfaces of the specimen of 15" at 700°F. The temperature gradient could be cut down to 3°F along the length of the specimen by adjusting shunts across the furnace windings. An ordinary metallurgical microscope with a working distance of 14 mm was rebuilt so that it could be moved in three mutual perpendicular directions. Visual examination could be made with either a 30X eyepiece or a 15X filar eyepiece. Using an objective of focal length 24.3 mm and a numerical aperture 0.20 a total magnification of 120X or 240X could be obtained. High purity aluminum rated at 99.995 pct was supplied by the Aluminum Co. of America as % in. then the gage portion was further machined to 3/16 in. diam. The ends were Y4 in. diam and threaded to fit available specimen holders. Two opposite flat surfaces of the gage portion were obtained by milling, giving final dimensions of the gage portion of 3/4x3/16x3/32 in. Machined specimens were first electropolished in Jacquet's solution.' Reference marks were produced by pressing a thin sewing needle into the surface of the specimen. The specimen was then annealed at 900°F for 24 hr and at 1150°F for 12 hr. Typical macrostructures of specimens treated in this manner are shown in Fig. 1. The range of the grain size was 1 to 5.5 mm (see Fig. l), but most of the grains were of 2.5 mm and went through the whole thickness of the specimen. The specimen was re-electropolished and etched electrolytically to reveal the grain boundaries. In order to reveal the orientation difference of the grains, etch pits were produced on the surface of the grains by etching in a mixture of hydrochloric and nitric acids (HC1, 45 pct; HNO,, 10 pct; H,O, 45 , by volume). The specimens were very carefully put into the holders and furnace; yet even with this care, fine slip lines were developed on the specimen surface.

Jan 1, 1953

By E. P. Miesch, P. B. Crawford

A study was made of the effect of permeable and impermeable lenses in a reservoir on the production capacity of a well. Both steady-state and unsteady-state data were obtained. An electrical resistivity model was used to obtain the steady-state data and thermal models. were constructed to obtain the unsteady-state data. The productivity of a well is affected very greatly only when the lenses are close to the well. The effect of circular lenses on the productivity ratio can be correlated with the distance from the center of the lens to the center of the well divided by the radius of the lens. When this dimensionless distance is equal to six or greater, the effect of the lenses on production capacity will be negligible. The pseudo steady-state productivity of a heterogeneous reservoir can be predicted using steady-state data. INTRODUCTION Many analytical solutions of reservoir behavior assume that reservoir rock is uniform and homogeneous. Although this assumption is used, all of the data from core analyses and well logging indicate that the reservoirs are heterogeneous. Very little work has been done on the performance of heterogeneous reservoirs. The work of Landrum,l et al. showed that transient phenomena in oil reservoirs could be studied with thermal models. Pickering and Cotman2 used thermal models to study flow in stratified reservoirs and investigated the effect of inhomo-geneities in oil reservoirs on transient flow performance. Loucks3 made a mathematical study of the pressure build-up in a system composed of two concentric regions of different permeability. Root, Silberberg and Pirson4 studied the effect of the growth of the flooded region on water influx predictions using a thermal model consisting of three concentric cylindrical regions of different thermal properties which simulated the aquifer, the flooded region and the unflooded portion of the original hydrocarbon region. Tomme, et a1.5 made a mathematical study of vertical fractures. The object of this investigation was to study the effect of highly permeable and impermeable lenses in the vicinity of the wellbore on the pressure depletion history of the well. Steady-state data were obtained for both conductive and nonconductive lenses that completely penetrated the formation. The lenses were symmetrically located at various distances from the wellbore. The unsteady-state data were obtained on seven thermal models. EXPERIMENTAL EQUIPMENT AND PROCEDURE STEADY-STATE DATA The steady-state data were obtained from an electrical resistivity model 30 in. in diameter and approximately 11/2 in. deep. The outside of the model was lined with a 30-in. diameter copper strip, which served as the outer boundary of the reservoir. The bottom was covered with a sheet of plexiglass so that it would be nonconductive. The model was filled with a slightly saline solution. The well size was varied from an 0.064-in. diameter copper wire to a 10-in. diameter copper cylinder. Readings were taken with an impedance bridge using AC current to prevent polarization at the contacts. Copper and wax lenses were used to represent infinitely conductive and nonconductive lenses, respectively. The resistance was first measured for each well diameter with no lenses in the reservoir. Then the conductive and nonconductive lenses were spaced symmetrically at various distances from the well and the resistance read from each lens location. The diameters of the conductive lenses were 3, 1.022 and 0.624 in., and those of the nonconductive lenses were 3, 2.25 and 1.563 in. The ratio of each resistance reading to the resistance reading with no lenses in the reservoir

Jan 1, 1965

By E. L. Dougherty, J. W. Sheldon

Using the numerical techniques shown in this paper it is possible to compute the simultaneous dynamical behavior of multiple fluid-fluid interfaces in two dimensions. Hence, fluid-fluid inter face models of several oil recovery processes can be constructed which allow prediction of the effect of fluid flow in two dimensions on the behavior of these processes while taking into account various physical effects, such as saturation gradients, phase changes and thermal stimulation. Techniques for constructing fluid-fluid interface approximations for several oil recovery processes are reuiewed. The validity of these techniques is established by comparison to experimental results. The availability of a computer program for computing the behavior of multiple fluid-fluid interface problems makes the use of potentiometric models uneconomic. The results indicate the areal sweep efficiencies obtained from potentiometric model studies at high mobility ratios may be somewhat optimistic. The results also indicate that the areal sweep efficiency curves obtained at high viscosity ratios 2—4 using miscible fluids in porous media are misleading and that considerable care must be taken in defining the mobility ratio for a displacement process, especially for miscible fluids. INTRODUCTION Many of the methods used to estimate areal sweep efficiency assume that the oil recovery process can be represented by one or more fluid-fluid interfaces. When the mobilities of the displaced and displacing phases are assumed equal, this quantity can be computed using analytical mathematical methods.l.2 For non-unit mobility ratio, experimental methods based upon the fluid-fluid interface concept include potentiometric and electrolytic blotter models. Scaled models of porous media have also been used to study this problem. Fay and Prats3 and Aronofsky4 used numerical techniques to compute breakthrough sweep efficiency, but neither of these studies were extended beyond breakthrough. Numerical techniques which have been developed for solving a general fluid-fluid interface problem are presented in a companion paper.5 These techniques have been incorporated into a program for a large-scale digital computer which can treat simultaneously the dynamical behavior of as many as six fluid-fluid interfaces. Thus it is possible to simulate the behavior of several oil recovery processes taking into account such effects as saturation gradients. Techniques for representing several oil recovery processes by fluid-fluid interfaces are presented. Processes considered are conventional water flooding with and without a mobile gas phase, hot water flooding, miscible flooding and enriched gas drive. Other processes, whose effects could be studied by this technique, but are not considered in detail here, are underground combustion and steam injection processes and various slug processes such as gas-water injection and alcohol-solvent processes. Theoretical results obtained with single interfaces for mobility ratios of 0.25, 1, 2, 4, 8, 16, 32 and 64 in a repeated five-spot are compared to experimental results for this geometrical configuration obtained using miscible fluids and conductive analogs. Predicted recovery curves for water flooding in a repeated five-spot are compared to experimental results reported in the literature. Example calculations are shown for a conventional water flood in the presence of a free gas phase and for a hot-water flood. MATHEMATICAL PRELIMINARIES We assume that the rock matrix is homogeneous with constant thickness, absolute permeability and porosity. The fluids are incompressible. We also neglect the effects of gravity.

Jan 1, 1965

By A. E. Dwight, A. F. Berndt

ThE published U-Ru equilibrium diagram indicates the existence of six intermetallic compounds in this system.' The structure of only two of the compounds, u2Ru2,3 and URU3,4 have been identified. This communication reports on the crystal structure and density of Uau3, which is formed by a peritectoid transformation at 1025° ± 5°C: U3Ru4 + U3Ru5 = U2Ru3. An alloy with the stoichiometric formula U2Ru3 was prepared by arc melting electrolytic uranium and ruthenium of 99.9 pct purity on a water-cooled copper hearth under an argon atmosphere. After homo-genization at 1000°C for 3 days, a metallographic examination revealed that the peritectoid transformation was complete and less than 5 pct of a second phase was present. A powdered sample, obtained by crushing the brittle compound, was annealed for 1 hr at 1000°C and air-cooled. Debye-Scherrer patterns were obtained with filtered CrKa, FeKa, CoKa, or CuKa radiation. The X-ray patterns of U2Ru3 weFe indexed as an fcc structure with a. = 12.895 0.00lA. The unit-cell constant and its standard deviation were calculated by a computer program described by Mueller, Heaton, and Miller.5 The density, measured by an immersion method, was 14.45 g cm-3. A calculated density of 14.48 g cm-3 was obtained on the assumption that twenty-four formula units (120 atoms) comprise a unit cell. The space group is either oh7 Fd3m, or the less symmetric, T4h, Fd3, since reflections indexed as hk0 with h + k + 4n are absent. All values of h2 + k2 + 12 permitted by both space groups are listed in Table I together with dc, do, and I. values. The calculated d spacings were obtained by the computer program of Mueller, Meyer, and simonsen.6 The observed d spacings above 0.77A were obtained from X-ray" patterns prepared with CuKa radiation, above 0.90A with CoKa radiaaion, above 0.98A with FeKa radiation, and above 1.16A with CrKa radiation. An attempt was made to identify the crystal structure of U2Ru3 by comparison with published data on other A2B3 compounds. A recent compilation7 of intermetallic compounds lists thirty-two structure types with the A2B3 stoichiometry. Only two of the structure types are fcc and they are Sb2O3(r), senarmontite, with 16 Sb2O3 per unit cell, and In2Te3 for which several structures have been proposed.8 However, the two cases cited do not have twenty-four formula units in the unit cell. We therefore, cannot identify UPRU3 with any of the other known A2B3 compounds. Ternary alloys were made in which 5 at. pct Os, Rh, or A1 were substituted for 5 pct of the Ru in U2Ru3. The osmium substitution increased the unit cell constant to 12.899 * .001A, which is to be expected since the atomic volume of osmium is slightly larger than

Jan 1, 1969

By L. G. Austin

My first contact with industrial milling was during the time I worked in the electricity generating industry in the United Kingdom. In visits to power stations to investigate either deposits in the boiler furnaces or polluting deposits settling around the stacks. I had to check the performance of the vertical coal pulverizers, since poor pulverization aggravated both problems. Naturally, then, when I came to the USA in 1957 to take a PhD in fuel technology at Penn State, I was put to work to review the science of coal pulverization. After this reviewing, I was completely confused. On one hand, there was a well-developed understanding of stress-strain equations, and a rapidly developing knowledge of how stressed, brittle solids fractured, based on the Grifith crack theory. On the other hand, reading in the grinding literature gave me: -Kick's Law, which was clearly not correct in the light of modem fracture theory; -Rittinger's Law, which was also clearly not correct; -Bond's Third Law of Comminution, which was claimed to have something to do with the Grifith crack theory, but where the connection between the two was made by intuitive pseudo-scientific reasoning I could not accept; -the choice of mill motor power for the most common type of coal mill, the Raymond pulverizer, was calculated from the fan power required to move air through the mill. Although I could accept the empirical connection between the two, it made no sense from the point of view of fracture energy. Even today, most books or review chapters on size reduction start from these laws. Incorrect statements abound in the literature, such as "the Hardgrove Index is based on Rittinger's Law," which it is not, "The Bond theory states that work input is proportional to new crack tip length produced in particle breakage," which is not true, etc. My own test work showed that these "laws" did not fit the data for grinding of coal. At about this time, Epstein (1948) and Broadbent and Callcott (1956), following the original work by R.L. Brown (1941) at the British Coal Utilization Research Association, proposed describing breakage as a series of fracture stages. I took their concepts and developed the basic differential equation for a batch grinding process continuous in time, analogous to a batch chemical reactor. Robin Gardner then joined the project and did his PhD on treating batch grinding in the same way as a batch chemical reactor. He found that the basic equation had already been partially derived by Sedlatschek and Bass (1953) in Ger- many. We confirmed experimentally the validity of the equations for describing batch grinding (1962) and formulated the equation describing steady-state continuous grinding in a fully-mixed mill. At about the time this work was published, Gaudin and Meloy (1962) and Filippov (1961) independently published essentially the same equations, but without experimental proof of the validity of the concepts. I will give a brief overview of what these beginnings had led to in the design of mills for size and power, and show some of the results of this more detailed understanding of grinding processes. Concepts of Fracture Mills such as tumbling ball, rod, pebble and autogenous mills and vertical mills such as the Raymond, and E-type, apply compressive stress to lumps or particles relatively slowly. Compressive stress applied to a particle of an elastic brittle solid imparts overall strain energy to the solid and produces local regions of tensile stress, [(Fig. 1)] (Berenbaum and Brodie, 1959).

Jan 1, 1985

By J. G. Early, F. V. Lenel, G. S. Ansell

The isothermal shrinkage rates of copper-powder compacts were determined in the temperature range from 760o to 1060oC. The rates for compacts fabricated from a pure spherical copper powder were compared, on the one hand, with those fabricated from a spherical powder containing a dispersion of aluminum oxide in the copper matrix and, on the other hand, with those fabricated from a pure copper powder having irregularly shaped particles. The irregular copper-powder compacts shrank much faster than the spherical copper-powder compacts and these, in turn, shrank faster than dispersion -strengthened copper -powder compacts. The experiments were run under a small external load. By varying this load, the stress dependence of the isothermal shrinkage rate could be measured for compacts from spherical-copper powders. This stress dependence was found to be approximately proportional to the 4th (dispersion-strengthened alloys) or 4.5th (pure metals) power of the stress. It is concluded that the material-trans port mechanism is due to dislocation climb rather than to Herring-Naharro creep. In this investigation the isothermal shrinkage rates of compacts made of copper powder with an average particle size of 70 g in the temperature range from 780o to 1060°C were studied. The aim of the investigation was to determine the mechanism of material transport responsible for the shrinkage of these compacts. Since the phenomenon observed in this investigation was shrinkage of compacts at high temperatures, it will be necessary to consider all the material-transport mechanisms possible in sintering. It is unlikely that evaporation and condensation, surface diffusion, or even volume diffusion of vacancies, in which surfaces of small radius of curvature are vacancy sources and those of larger radius of curvature vacancy sinks, will materially contribute to a sintering process involving shrinkage, i.e., one in which the centers of adjacent particles move towards each other. The two principal mechanisms which were considered were the so-called Nabarro-Herring mechanism and a plastic deformation mechanism by dislocation movement. Both of these mechanisms were originally conceived to explain creep, i.e., the plastic deformation of metals at elevated temperatures under a constant applied stress. The shrinking of metal-powder compacts which was studied in this investigation can be readily treated as a creep process. This means that the shrinkage rates observed are considered as creep rates under the influence of an applied stress. It is true that the stresses in sintering experiments are generally not constant but decrease continually with time, but, by considering creep rates or shrinkage rates over a relatively short time interval, the small changes in stress during this time interval can generally be neglected. Nabarro-Herring microcreep is a high-temperature deformation mechanism which is caused by a diffusional flow of vacancies from boundaries under a tensile stress towards those under a compres-sive stress. Herring1 has shown that in this type of microcreep the strain rate e is proportional to the applied stress and follows an equation of the form ma, aexp(-Q/RT) [l| where <j is applied stress, Q is the activation energy for self-diffusion, R is the gas constant, and T is the absolute temperature This type of creep can explain the shrinkage of a porous compact. The surface of the pores are subjected to the tensile stress of surface tension, while the grain boundaries are under a compressive stress also due to surface tension. Vacancies diffuse from the pore surface through the lattice to the grain boundaries, where they are eliminated. In order to maintain vacancy equilibrium at the pore surface the vacancies leaving the region of the pore surface are replenished from the volume of the pore, which means that the pores shrink.

Jan 1, 1964

By H. J. Morel-Seytoux

Methods of predicting the influence of pattern geometry and mobility ratio on water flooding recovery predictions are discussed. Two methods of calculation are used separately or concurrently. The analytical method yields exact solutions in a convenient form for a unit mobility ratio piston-like displacement. A few typical pressure distributions, sweep efficiencies and oil recoveries are presented for various patterns. For non-unit mobility ratio, one may resort to a numerical method, such as that of Sheldon and Douherty. 1,2 Because the domains of applicability of the analytical and numerical techniques overlap, the exact solutions provide estimates of the errors in the numerical procedures. The advantages of the analytical and numerical methods can be combined. To develop a numerical technique as independent of geometry as possible, the physical space is transformed into a standard rectangle. The entire effect of geometry is rendered through one term, the "scale-factor", derived from mapping relations. The scale factor can be calcu-lated from the exact unit-mobility ratio solution for the particular pattern of interest. By this means recovery performances for arbitrary mobility ratio con be obtained for many patterns. A sample of tesults obtained in this manner is presented. INTRODUCTION Pattern geometry and mobility ratio are two major factors in making a waterflood recovery prediction. Because assisted recovery has become increasingly important to the oil industry, pattern configuration and mobility ratio also assume a greater significance in the assessment of the economic value of recovery projects. The influence of pattern geometry and mobility ratio in shaping a recovery curve and on the other quantities of interest to the reservoir engineer is the main subject of this paper. Much effort has already been spent on estimating quantitatively the influence of either pattern or mobility ratio or both on oil recovery. The literature reports many investigations of this nature. 3-9 However, many results or methods of recovery prediction presented in the literature cannot be considered fully satisfactory. Even for unit mobility ratio and piston-like displacement, where analytical solutions are available, the literature shows discrepancies. For non-unit mobility ratio, the divergence in the results is extreme. For infinite mobility ratio in a repeated five-spot, depending on the investigator, the sweep efficiency ranges from 0 per cent to 60 per cent. With respect to the influence of pattern on recovery, only the repeated five-spot has received much attention. Other confined patterns and pilot configurations have received very little attention. Two calculation methods are presented in this paper, either separately or concurrently: the analytical method of potential theory and the numerical method of finite-difference approximation. The analytical method is more restricted in scope than the finite-difference method, but it has the definite advantage of providing exact solutions within its range of applicability. If a unit-mobility ratio piston-like displacement is assumed, the analytical approach is possible. A few typical results are reported in this paper; the detailed description of the general method and of a great variety of results will be the subject of other articles. Fornon-unity mobility ratio, we must resort to a numerical scheme. The numerical technique is that which was described by Sheldon and Dougherty.l,2 It is not limited to piston-like displacement. However, mainly single interface results will be presented here. Because the respective domains of applicability of the analytical and the numerical method overlap, useful comparisons of exact and numerical solutions can be made for a variety of patterns. The advantages of the analytical and numerical approaches can be combined. The reason for the success of this analytical-numerical approach can be summarized in the following two points: 1. The numerical solution for arbitrary mobility ratio can be programmed most efficiently when the physical space in which the displacement actually takes place is transformed into a standard shape; and 2. This can be done with remarkable simplicity whenever an analytical expression for the pressure

Jan 1, 1966

By FREDERICK J. ALCOCK

THE high prices which lead and zinc have commanded during recent months have given a great impetus to search for workable deposits of these metals, and there has accordingly risen a demand for information about deposits which have been worked in the past and which might possibly be opened up anew under these more favorable conditions. Production of zinc and lead in Canada has been on the increase during the last few years. In 1924, 175,- 485,499 lb. of lead with a computed value of $14,221,- 345 were mined. Of this, 96 per cent. was produced in British Columbia, 3 per cent. in Ontario, while Quebec and the Mayo district of the Yukon accounted for the remaining one per cent. During the same year the production of zinc amounted to 98,909,077 lb.. valued at $6,274,791. Of this 96,000,069 lb. were produced in British Columbia and the remaining 2,909,008 in Quebec. The large producer of zinc and also of lead is the Sullivan Mine at Kimberley, B. C. In comparison with the Sullivan, the deposits which are known to occur in eastern Canada are relatively small. Two, however, the Kingdoll lead mine at Galetta, Ontario, and the Tetreault zinc property at Notre Dame des Anges, Quebec, are steady producers. The object of this paper is to describe briefly these and other deposits of eastern Canada, to show their relation to geological events and to point out what are the chances of any of those which are idle at present developing into active producers or what are the possibilities of new deposits being discovered.

Jan 1, 1926

By D. I. Bardos, F. X. Spiegel, R. K. Malik, Paul A. Beck

Study of the occurrence of elongated ß-Mn phase fields and of their direction of extension in fifteen ternary systems indicates that the ß-,mn phases formed by manganese with other transition elements and with either silicon, germanium, or tin may be considered as electron phases. In these phases the alloying behavior of silicon on the one hand and of germanium and tin on the other shows a significant difference . Formally all observations can be satisfactorily described by attributing to germanium and tin an apparent electron concentration of nine electrons per atom. This behavior may be interpreted by assuming that germanium and tin exclude itinerant electvons because of a strong repulsive potential arising fromthe requirement that their wave junctions be orthogonal to those of- the closed d shells of these atoms. The absence of a d shell in silicon allows this atom to admit itinerant d electrons. The following new ternary phases of unknown structure were identified by means of their X-vay difjraction patterns: (Mn,Co,Si)I , (Mn,Co ,Si)Y, (Mn.Co,Si)Q , (Mn ,Ni ,Si)W , and (Mn, Ti ,Si)A. The X-ray data suggest that these phases have large unit cells and that to some extent they may he related structurally to the 0 phase. The a phase and certain other phases formed by transition elements are characterized by their preferred occurrence at a constant electron concentration. AS a result, these "electron phases" have relation to electron concentration. The tendency of the B phase to form elongated phase fields in ternary systems along the directions of constant electron concentration has been previously demonstrated. The occurrence of an elongated 0-Mn phase field in the ternary system Mn-Co-Si was schematically indicated by Cherkashin et al.l0 The present study includes fifteen ternary systems, each containing in addition to manganese a second transition element from the group Ti, V, Cr, Fe, Co, Ni and one of the nontransition elements silicon, germanium, or tin. In each system the boundaries of the B-Mn phase were determined in the 900°, 950°, or 1000°C isothermal section, the temperature being dictated on the one hand by the requirement of sufficient rates of diffusion to approach equilibrium conditions in reasonably short annealing periods and on the other hand by the need to avoid as much as possible the partial melting of alloys during annealing.

Jan 1, 1967

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

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

Jan 1, 1967

By T. S. Lundy, R. A. Padgett

The fact that cobalt diffuses in silver at a much slower rate in a region near a free surface than in bulk material has been demonstrated in a variety of experiments. Various possible mechanisms of this near-surface effect (NSE) are examined in light of experimental evidences. The tendency for surface segregation of some solute atoms appears to be a likely mechanism to explain the effect. UNIDIRECTIONAL and isothermal diffusion from a plane source of radioactive tracer into an isotropic, homogeneous medium follows the Gaussian solution to Fick&apos;s second law, or: where A(x,t) is the position and time-dependent concentration of the diffusing species, M is the quantity of diffusant per unit area initially at x = 0, and D is the diffusion coefficient. For a given time, say t = tl, one can easily see that: d In Aixji)_____1_ r2] d(x2) 4Dti Thus, if diffusion data obtained by serial sectioning are plotted as In A(x,tl) vs x2, the diffusion coefficient is obtained from Eq. [2]. The ideal behavior seldom is observed in real systems. Deviations from ideality generally occur either at small or large values of penetration distance. The overall picture is shown schematically in Fig. 1 where Regions I, 11, and 111 are due, respectively, to some near-surface effect or effects (NSE), volume diffusion, and short-circuiting effects. For different systems and conditions of specimen preparation and analysis, the relative importance of Regions I, 11, and I11 may be drastically different. Region I is completely absent in many systems and there is no NSE. In others, however, there may be considerable difficulty in determination of volume diffusion coefficients due to the overshadowing influences of Region I. Such is particularly true when indirect techniques such as the surface concentration decrease or the Gruzin method&apos; are used to find D. In this paper several experiments which had as their goal a better understanding of near-surface effects are reported. Specifically, it is well-known that unusual penetration behavior results for the diffusion of several impurities in ilver- and that considerable controversy exists with regard to interpretation of the published diffusion data. We chose to examine those unusual results which occur for cobalt diffusion in monocrystalline silver395 primarily because knowledge of experimental procedures and materials were readily available. Both radioactive co60 and co5&apos; could be utilized to offer a wide range of available concentration of the diffusing species. Experiments were also performed on the diffusion of AgnO in silver. They gave base-line data for the system—both for technique evaluation and for examination of the vacancy distribution near the surface. I) PROCEDURES Cylindrical-shaped specimens 1.25 cm diam and 1 to 1.25 cm long were cut from 99.999 pct Ag single crystals obtained from Monocrystals Co., Cleveland. One end surface of each specimen was polished flat and etched with a mixture of 20 ml saturated Cr03 solu-tion-20 ml H20-3 ml concentrated HC1. Radioactive tracer (co60 or co5&apos; or Ag&apos;lO, each having a radio-chemical purity of greater than 99 pct) was deposited onto the prepared surface either dropwise with evaporation to dryness of the solution containing the tracer (0.5 N HC1 solution for the co60 and co5&apos; and 0.6 N HNO3 solution for the A"&apos;), by electroplating, or by evaporation from a tantalum or tungsten filament in an evacuated bell jar. The amount of deposited cobalt was varied from many atomic layers of cobalt containing the co60 to approximately one co5&apos; atom per lo6 surface atoms of silver. This large variation in amount of deposited cobalt yielded information about both the influence of a thick-source layer of tracer and solubility effects on the NSE. Such variations were possible because of the large differences in half-lives,

Jan 1, 1969