Search Documents

By J. C. Swartz

New measurements of graphite and cementite solubilities in ferrite have been obtained to resolve disagreements among previous data. To measure graphite solubility, specimens of iron were equilibrated with rarious carburizing atmospheres at temperatures be-twee~z 480" and 1435°C. The carbon activity of each atmosphere relative to graphite was taken from published data. After quenching from each equilibration treatment the concentration of carbon (typically 10 pprn by wt) was determined from the height of the Snoek peak in the internal friction spectrum. A few inclividual analyses by combustion support the reszdts of Smith on the solubility of graphite near the eutectoid temperature. Combination of these combustion and internal friction analyses determines In a, (at 0.01 wt pct) = 24,000/RT - 12.72, where a, is the activtity relative to graphite in ferrite at fixed cornposition. The same specimens were used to determine the solubility of self-stressed cementite precipitates. The results indicate solubilities of 6.3, 23.2, and 68 pprn carbon by wt at respectively 400°, 500". and 600°C. The once large discrepancy between the cementite solubility determined by internal frictiori and that calculated from flee energy data has been almost completely resolved. A recent study1 (hereafter designated I) on the solubility of cementite in a iron did not fully resolve the disagreement among previous data. The remaining disagreement is illustrated by the phase diagram in Fig. L The data of I are represented by the line labeled "stress-free" cementite. They determine a heat of solution of cementite of (14.8 * 0.2) kcal per mole which contrasts with the (20 * 1) kcal per mole derived from Smith's graphite solubility2 and the formation energy of cementite.3 Toward resolving this disagreement we have obtained new data on the solubilities of graphite and cementite. METHODS As in I the Snoek peak in the internal friction spectrum was used to measure the concentration of carbon dissolved in ferrite. This method relies on the relation C=kQ-1 [1] in which C is the concentration of dissolved carbon, Q-I is the Snoek peak height, and k is a proportionality factor which is a function of the peak temperature, polycrystalline texture, and other parameters characteristic of the ferrite lattice. To determine true heats of solution via Eq. [1], it is necessary to establish that k is independent of the temperature from which the specimen is quenched. Some evidence on this point is contained in a separate report.* Additional evidence is shown by the data points in Fig. 1. These represent the Snoek peaks of two specimens carburized to about 10 ppm C, copper plated, and quenched first from 888º then from 588°C. The copper plate was removed before each Snoek peak measurement. These data show that there is no detectable dependence of k on quenching temperature. The solubility of graphite was determined by repeatedly carburizing various specimens of iron wire with gas mixtures of known carbon activity relative to graphite. To cover a wide temperature range within the a phase boundaries it is desirable to use as low a carbon concentration as is compatible with reasonable precision. With the present internal friction apparatus the best compromise of these desires is a concentration near 10 wt ppm. The dashed line in Fig. 1 was used to approximate the activity of carbon* for

Jan 1, 1970

By W. W. Whitaker, C. W. Manry, R. H. McLean

In an eccentric annulus, cement may favor the widest side and bypass slower-moving mud in the narrowest side. Tendency of the cement to bypass mud is a function of the geometry of the annulus, the density and flow properties of the mud and cement and the rate of flow. Bypassing can be prevented if the pressure gradient protluced from circulation of the cement and buoyant forces exceeds the pressure gradient necessary to drive the mud through the narrowest side of the annulus at the same velocity as the cement. In the absence of buoyant forces, one requirement for this balance is maintenance of the yield strength of the cement greater than the yield strength of the mud multiplied by the maximum distance from the casing to the wall of the borehole and divided by the minimum distance. If the yield strength of the cement is below this value, bypassing of mud cannot be prevented unless buoyant forces or motion of the casing significantly aid the displacement. INTRODUCTION Successful primary cementing leaves no continuous channels of mud capable of flow during well treatment and production. Prevention of channels requires care. Tep-litz and Hassebroek provide evidence of channels of mud after primary cementing in the field.' Channeling of cement through mud in laboratory experiments has also been reported.'-' Recommendations for improving the displacement of mud include (1) centralizing the casing in the borehole,'-" 2) attaching centralizers and scratchers to the casing and moving it during displacement,18 "3) thinning the isolating the cement by plugs while it is circulated down the casing,%( (5 establishing turbulence in the cement," and (6) holding the cement slurry at least 2 lb/gal heavier than the mud and circulating the cement slurry at a very low rate of flow.' Although much has been written about the above parameters, the relative importance of each has not been well defined. In this investigation, the mechanics of mud displacement are described through results from analytical models and experiments. The model chosen — a single string of casing eccentric in a round, smooth-walled, impermeable borehole — is analagous to casing centralized in a borehole which is not round and to placing more than one string of casing in a borehole. In each, some paths for flow are more restricted than others. A fluid flowing in the borehole may seek the least restricted, or most open, path. This tendency for uneven flow can lead to channeling of cement through mud unless preventive measures are taken. The analytical models describe channeling and give means of balancing the flow. Experimental data test the analytical models and illustrate effects of motion of the casing, differences in density and mud's tendency to gel. Results are encouraging. Piston-like displacement of mud by an equal density cement slurry is possible through proper balance of the flow properties of the mud and cement slurries to the eccentricity of the annulus. The more eccentric the annulus, the thicker must be the cement relative to the mud. If proper balance is not achieved. bypassing of mud by cement cannot be prevented without assistance from motion of the casing or buoyant forces. Increasing the rate of flow can help to start all mud flowing but cannot prevent channeling of cement through slower moving mud in an eccentric annulus. Thinning the cement slurry tends to increase channeling although the extent of turbulence in the annulus may be increased. Description of flow in an eccentric annulus begins in the next section. It is assumed that (1) the casing is eccentric and is stationary, (2) the mud and cement slurries have the same density and (3) the gel structure of the mud has been broken and the mud and cement follow the Bingham flow model. Effects related to these restrictions will be discussed. FLOW PATTERNS SlNGLE FLUID IN ANNULUS Flow of a single fluid through an eccentric annulus is illustrated in Fig. 1. Part A shows laminar flow of a Newtonian fluid. This distribution of flow was calculated by Piercy, Hooper and Winney.' In fully developed turbulent flow, the velocity distribution around the annulus is less distorted, but the flow still favors the widest part of the annulus Parts B, C and D of Fig. 1 are a qualitative representation of the flow of a Bingham fluid. The yield strength of the fluid increases the severity of bypassing compared to Newtonian flow. At a very low rate of flow, all flow is confined to that portion of the annulus which has the minimum perimeter-to-area ratio. The fluid shears on the perimeter of that area when the pressure gradient multiplied by the area just exceeds the yield stress of the fluid multiplied by the perimeter. Whether or not the minimum perimeter-to-area region encompasses all of the annulus or only a part (as shown in Part B) depends on the geometry of the annulus. If only a part begins to flow, increasing the rate of flow increases the area flowing until finally there is flow throughout the annulus.

By Joseph W. Leonard

The purpose of this paper is to review and discuss some of the less evident and sometimes neglected opportunities for progressive developments in coal research. While a great deal of both promotional and technical information flows from some areas of coal research, output deficiencies in other areas of activity have reached a magnitude where important developments have been, and will increasingly be, unfavorably affected. These areas mainly involve coal mining and preparation. Some recommendations for the intensification of effort in these areas follow: Coal Mining While a huge tonnage of in-the-ground coal is assured, the location and distribution of these tonnages are becoming less favorable. The easy-to-mine coal which is located in or near population centers has been, or is being, mined. The vigor with which the less accessible reserves are recovered by the mining industry depends largely on the condition of the coal market at the time of mining. Hence, during a buyer's market, the commercially oriented mining industry is compelled to mine the easier and less costly reserves. Conversely, during a seller's market, the need to rapidly expand production results in more difficult mining and higher cost coal as few obstacles are encountered in finding markets. Hence, a seller's market tends to enhance the recovery of reserves while a buyer's market does not. One reason for today's fuel supply problems is that the Nation has recently emerged from a long-term coal buyer's market which lasted from about 1950 to 1968. During that period, national policy caused severe production cutbacks which regretably drove the industry to mining only the more accessible and better quality reserves. Often in order to remain in business, many hundreds of millions of tons of more difficult to mine reserves were abandoned and lost behind caved areas. Many of these reserves are close to population areas and would not have been lost in a more stable economic climate. It is difficult to fully account for all the impacts that were caused by the great buyer's market of the 1950s and 1960s. Besides the obvious loss of reserves that were once considered national wealth, the mining of better reserves tended to produce a generation of technically optimistic mining people. Mining people frequently became accustomed to looking at nothing less than outstanding mining conditions as a result of the declining market. Many are now and have long since received a re-education in the other half of mining. Going from many years of mining accessible, select and easy-to-win reserves, to the crash-driving of development entries in reserves that were considered unworthy of mining during 50s and 60s, frequently results in a much higher rate of encounter with in-seam and out-of-seam rock as well as with coal-deficient areas or "washouts." Intensive entry driving activity and compulsory non-selective mining in sometimes lean reserves were brought on by the need to rapidly open up new supplies of coal. Working under these requirements presents a continuing reminder that much more needs to be known about the relatively esoteric art of planning the best direction for driving entries in order to insure that a more consistent and greater supply of coal is available during early mine development. All of the preceding discussion tends to point to a need for a better estimate of those reserves of coal that are likely to be mined in the future. Such estimates should not be limited to the compilation of the amount of coal in the ground; but, where possible, should also include information concerning the capability for producing this coal. After all, a coal seam of ample thickness may have a degree of thickness variability, undulation, bad roof or floor, so as to make what would otherwise appear to be an attractive mining condition untenable. Underlying the problem involving the feasibility of producing known reserves is the need to develop better methods for the characterization of coal seams and associated lithotypes, based on drill core data, once at area is selected for mining. Reserves and their characterization involve aspects of exploration technology that are frequently considered mature. The resulting technological deficiencies may be the main reason why coal exploration frequently does not end with core drilling of a property, as it should, but extends into the mining operation during the driving of development entries. When exploration is extended to the driving of development entries, the near absence of integrated decision-making theory involving mining, geology, mathematics, and economics becomes, once again, all too painfully apparent and frequently results in very costly rationalizations. Hence, by the formal initiation of a concentrated program to combine the cyclical effects of economics with geology and mining, more relevant estimates of reserve distribution, tonnages, and production capability should be forthcoming. Moreover, a similar formal effort is needed to develop a combination of the most advanced concepts of mathematics, geology, and mining to better "see" coal seams as a means to favorably implement many long-range decisions involving mine safety and productivity. Much more applied research needs to be done on coal mining systems for mining in thin seams and/or under bad roof. Current difficulties in both of these areas at recently opened coal mines should provide a sobering glimpse into the future. Full-scale applied research, sponsored by appropriate federal agencies, is urgently needed on a scheme involving a new combination of established mining and preparation elements. The scheme may include: (1) a continuous mining machine remotely operated by a miner stationed at some distance behind the machine using a cord attached control box; (2) hydraulic transport of coal through pipes from the mining machine to a coarse refuse removal grid, crusher, and then on to portable concentrating equipment; (3) the hydraulic transport of clean coal out of the mine in pipes to the surface for thermal dewatering, if neces-

Jan 1, 1974

By W. Vickers

Single-point rotating cantilever fatigue tests have been carried out at 550" to 650°C on beryllium produced by a variety of fabrication routes. All the specimens gave similar plots of stress against number of cycles to failure but failure occurred at a lower stress range in the hot-pressed material than in the extruded material. A discontinuity or "kink" appears to exist in the curves and this is considered. A true endurance limit was not evident up to 10 cycles. When Brush hot-pressed warm-extruded beryllium was tested at 600oC a relationship where Nf = number of cycles to failure, ?ep = plastic strain range, and c = constant, was shown to be applicable for values of ?ep > 6.3 x 10-4. Further for this moterial and temperature and for frequencies between 150 and 4500 cycles per min, the variation of number of cycles to failure with frequency obeyed the expression N = Af0.44 where N =number of cycles to failure, f = frequency, and A = constant. Examination of the fracture surfaces indicated that they were generally intergranular and in all cases the higher the stress the more irregular was the fracture. Microexamination has shown that grain boundary cavities and fine intergranular cracks were associated with the fractures. REVERSING or cyclic stresses set up in fuel elements can lead to failure by fatigue. Fatigue stresses may be developed by instability of the fuel element in the coolant flow leading to rattling against the channel wall or fin vibration. Thermal cycling also sets up reversing stresses due to differential expansion of fuel and can. Fatigue resistance is of some importance therefore in the evaluation of canning materials. At the time the work described in this paper was initiated, beryllium was being considered for the can and endcaps of an advanced gas-cooled reactor fuel element. With this in mind, beryllium obtained from different manufacturing routes was tested on single-point rotating cantilever fatigue machines at temperatures the same as those envisaged in the reactor to obtain appropriate fatigue data. The results obtained indicated a discontinuity in the S-N curves coinciding with the yield stress in the specimens. Although in some cases the data is insufficient to warrant a detailed quantitative analysis, a qualitative analysis of this effect was thought to be justified. 1) EXPERIMENTAL DETAILS The materials used are given in Table I together with the production route and available analyses. The specimen shape is shown in Fig. 1. The fatigue specimens were machined slightly oversize, after which they were given a stress-relieving anneal of 1 hr at '750° C in vacuum. This was followed by an etch in sulfuric-phosphoric-chromic acid mixture down to the finished dimensions given in Fig. 1. Etching served to remove the interference colored oxide film formed on the surface during annealing and very small recrystallized surface grains too. The specimen was mounted in the grips of a cantilever machine and rotated inside a furnace controlled at the testing temperature. The specimen was loaded in increments as the speed was increased until the required load was applied at 5000 cycles per min. A revolution counter was used to record the total number of cycles.

Jan 1, 1964

By John C. Bilello

The brittle fracture behavior of cold-worked sintered tungsten was studied over the temperature range 4.2° to 298°K using a high-sensitivity strain measuring system and electronfractography. Similar observations were made on a swaged electron beam zone-refined monocrystal. In sintered tungsten irreversible plastic deformation was observed during cyclic load-unload tests at stress levels well below the fracture stress for all temperatures, but general microyielding could be detected only down to 202°K. For the zone-refined samples macroyielding occurred at all test temperatures with evidence for twinning below -202°K. The fracture stress of the sintered tmgsten was virtually independent of temperature, while the zone-refined crystal showed a 2.3 times increase over the same temperature range. Electronfractography confirmed the presence of numerous rod-shaped and spherical submicroscopic voids which ranged in diameter from 1400 to 4300A in the sintered tungsten; no voids could be found in the zone-refined tungsten. Contrast effects observed on the replicas in the vicinity of certain voids indicated that plastic deformation could be induced by the local stress concentration. It has been suggested that the presence of these voids may be responsible for the low-temperature brittle failure of sintered tungsten. Based m this suggestim und on the evidence obtained here, a dislocatim model is presented to account for the brittle behavior of sintered tungsten. In this model slip, which is induced by the local high stress concentration in the region at the edge of a favorably oriented void, could cause the void to grow to a microcrack of critical size. STUDIES of brittle fracture in bcc metals have led to the well-known experimental relationships between grain size, yield stress, fracture stress, and temperature which have formed the basis for the various dislocation pile-up1-3 or interaction4'= models for slip-induced microcrack nucleation. While microcracks can be nucleated by deformation twins,6,7 there has been no direct evidence furnished by transmission electron microscopy to support conclusively either the Zener pile-up or Cottrell dislocation reaction models for producing micro-cracks in all "brittle" materials. In addition to the "inverse" grain size relationship for yield and fracture stresses the cottrel14 theory predicts that the fracture stress below the transition temperature should behave in a fashion similar to that of the yield stress above this temperature. Such behavior has been verified for several bcc metals.8-10 With reference to both grain size effects and the tem- perature dependence of the fracture stress below the transition temperature, the behavior of sintered tungsten appears anomalous. Early work by Bechtold and Shewrnon 11 showed no apparent temperature dependence of the fracture stress below the ductile-brittle transition temperature (DBTT). They attributed this result to the intergranular nature of the fractures observed. More recent work by Wronski and Four-deux12'13 on considerably purer material did not show any systematic relationship between the fracture stress and temperature below DBTT. The dependence of flow and fracture stresses on grain size is also not clearly established for sintered tungsten. Koo, for example, has shown that the DBTT for sintered tungsten depended chiefly on the annealing temperature and was relatively insensitive to the actual grain size achieved. Using electrofractography and transmission electron microscopy, Wronski and Fourdeuxl3 showed that numerous spherical and rod-shaped submicroscopic voids could be found in sintered tungsten but not in melted tungsten of nominally the same purity. They suggested that these voids could be responsible for the temperature insensitivity of the fracture stress below the DBTT. In the present work the temperature dependence of the fracture stress for high-purity commercially sintered tungsten has been determined. The presence of submicroscopic voids in sintered materials was confirmed, and these were studied in detail to examine the role they could play in nucleating brittle fracture. A dislocation model is suggested which could cause an inherent spherical void to lengthen into a Griffith crack of critical size. EXPERIMENTAL PROCEDURE Commercially sintered tungsten rod was obtained in the as-swaged condition from Sylvania. A zone-refined crystal was obtained from the same source. This crystal was grown by giving three zone passes (at 25.4 cm per hr) to a sintered rod of high-purity tungsten. The rod axis prior to cold working was -15 deg from the [110] direction. Originally the zone-refined rod was -6 mm in diam; it was reduced to -3 mm by eight swaging passes, at high temperatures, with each step having about the same reduction of area. The final swaging step gave a 7.5 pct reduction of area at 1050°C. All swaging operations were performed in a hydrogen atmosphere. For the sintered rod a similar working schedule was employed. Metal-lographic examination of the sintered material revealed that the cold-worked structure had an apparent grain diameter of -25 u transverse to the swaging direction (obtained by the intercept method). In the longitudinal direction cold-worked grains were approximately 1.5 to 2 times their diameter. No distinct fiber structure could be observed optically for the zone-refined rod. The cold-worked structure in the

Jan 1, 1969

By R. W. Parsons, P. R. Chaney

Oil recovery by the imbibition mechanism can be important in fractured carbonate reservoirs with a bottom water drive. Laboratory experiments were performed on water-wet carbonate rocks to model this process. The results are interpreted in light of the applicable scaling laws. The imbibition production behavior of a pillar of rock subjected to a slowly rising water table can be synthesized from data on total immersion experiments on small rock samples. If the rate of water table rise is slow enough and if the physical properties are such that the imbibition zone is of a small vertical extent, then only the final oil recovery on the small samples need be known. Applying laboratory imbibition data to a real reservoir should be done with caution because of possible wettability problems and the unknown behavior of two-phase flow in fractures. INTRODUCTION The spontaneous taking up of a preferential wetting fluid into a porous medium with the simultaneous expulsion of the contained fluid is termed imbibition. One of the many physical imbibition processes is the production of oil by the taking up of water into a water-wet reservoir rock. In certain types of heterogeneous reservoirs, e.g., fractured, it is realized that imbibition alone may be the dominant mechanism. Theoretical studies of imbibition oil recovery have followed several lines of attack. The behavior of certain mathematically tractable and conceptually simple models have been analyzed. Birks1 considered a bundle of vertical capillaries model, while Perotti et al.2 and Kelemen3 used a model of a series of vertically oriented fractures. If the oil production behavior from a block of reservoir rock is known, then the cumulative behavior of a vertical stack of these blocks subjected to a rising water table can be determined by a method outlined by Aronofsky et al. 4 Knowing the capillary pressure and relative permeability functions should uniquely set the imbibition behavior for a given sample. Solutions of the partial differential equations describing the process, however, present some formidable problems. Digital computer techniques have been applied to displacement processes in which the capillary pressure has been included.5"8 More recently, similar techniques have been used by Blair 9 to calculate countercurrent imbibition behavior. The fact that imbibition behavior is described by the usual fluid flow equations (Darcy's Law and continuity equation) implies that the derived scaling laws should apply equally to this situation. Use of model experiments to study imbibition has been employed by Graham and Richardson 10 and Mattax and Kyte,11 who actually tested the scaling law applicability. This study is concerned with scaled model imbibition experiments on water-wet carbonate rocks. Of particular interest is the slowly rising water table experiment. This simulates a highly fractured reservoir with the water table rising uniformly through the vertical fracture system. Experiments were performed on different sized samples and with different rates of rise of the water-oil interface surrounding the rock pillar. Total immersion imbibition experiments on smaller rock samples were performed to test the basic linear scaling laws and to see if these simpler tests could be used to synthesize the rising water table experimental results. THEORY IMBIBITION SCALING LAWS Scaling laws for the flow of two incompressible immiscible fluids in a porous medium have been presented by Rapoport.12 These are derived directly from Darcy's Law for the individual phases and the continuity equation. Results from a laboratory model experiment will duplicate those from some reservoir prototype if a constant proportion is maintained among the three forces acting on the fluids — the capillary pressure gradient, the gravitational gradient and the flowing pressure gradient. Operation of the model (i.e., fluid injection or withdrawal) must be conducted in accordance with specific equations. If the

Jan 1, 1967

By P. T. Chiang

Chemical etching methods for the simultaneous revealing of the tellurium or selenium Phase and the chalcogenide grain boundaries of the alloy systems are given. A tellurium eutectic was found Present in zone-melted ingots. Similarly, a selenium monotectic was present in ingots. In general, the second phase (tellurium or seleniumn) occubies three different sites; viz., along the chalcogenide grain boundaries, as inclusions within the chalcogenide grain, and on the undersurface of the ingot. The detection limit for the tellurium phase is about 1 u in width. THERMOELECTRIC materials based on Group V (bismuth, antimony) and Group VI (selenium, tellurium) elements have aroused considerable interest in recent years in the practical application of thermoelectric cooling. In many cases, a small amount of excess tellurium (or selenium) was added to the material to optimize its thermoelectric properties. Then the question immediately arises as to the number of phases present in the resultant alloy. In the binary systems of Bi-Te, Sb-Te, and Bi-Se, the congruent melting compositions have been reported to be non-stoichiometric and are represented by Bi~Te respectively. It is to beexpected and known that Bi2Te3 and SbzTe3 crystallize from the melt with an excess of bismuth and antimony in the lattice and that tellurium forms a eutectic.~' The same could be assumed to take place in the pseudo binary systems of (Bi,Sb)zTe3 and Bi2(Se,Te)3 as well as in the system studiedby puotinen5 and other workers. Likewise, BiaSe3 crystallizes from the melt with an excess of bismuth in the lattice and selenium forms a monotectic.~ Therefore, in practice, alloys solidified from the melt often contain a second phase (tellurium or selenium) in one region or another of the solid mass even without the addition of excess tellurium (or selenium). ~u~~recht' studied the thermoelectric properties of (Bi,Sb)2Te3 alloys with excess tellurium and simultaneous additions of selenium. He mentioned that the materials show two phases because of the considerable excess of tellurium or selenium. However, he did not report as to how the tellurium or selenium phase was identified. It is generally believed that the presence of an excessive amount of tellurium or selenium phase in the alloy would adversely affect its thermoelectric properties and its uniformity. Consequently, there is a need for a simple method for the identification of the tellurium and selenium phase. The quantity of the second phase present is usually too small to be detected either by chemical analysis or by normal X-ray techniques. This investigation was therefore carried out, first, to devise a simple metallographic method for the identification of the tellurium or selenium phase coexisting with the chalcogenides and, second, to determine the distribution and specific location of the tellurium or selenium phase in the ingots. EXPERIMENTAL PROCEDURE The starting materials used for the alloy preparations were 99.999 pct pure bismuth, antimony, and tellurium and 99.997 pct pure selenium. The bismuth and antimony were obtained from Consolidated Mining and Smelting Co. of Canada Ltd., while the selenium and tellurium were obtained from Canadian Copper Refiners Ltd. The tellurium was purified further in the laboratory by zone refining. The elements were pulverized in a stainless-steel pestle and mortar. The amounts for the desired composition were weighed out each time on an analytical balance to make up a 100-g sample. Then the sample was introduced into a Vycor ampule (19 by 150 mm), pumped down to a vacuum of 10"5 Torr for 15 min, and sealed off. The ampule was then heated in a horizontal resistance furnace at 800" to 900°C for about 20 hr. During this period the assembly was rocked back and forth several times to ensure good mixing. At the end of the heating period, the ampule was quenched in cold water and then transferred to the zone-melting apparatus described in a previous publications to grow large-size aligned polycrystals. The background and ring-heater temperatures were adjusted to make the freezing solid-liquid interface slightly convex to the liquid. The recorded temperature gradient in the vicinity of the freezing solid-liquid interface was around 15°C per cm. The ampule was moved horizontally at a speed varying from 0.4 to 2 cm per hr so that the ring heater would cover the whole ingot length from end to end. A single zone-melting pass was used for the Bi-Te, Sb-Te, and Bi-Sb-Te ingots. Two passes in the forward and reverse directions were carried out for the Bi-Se and Bi-Se-Te ingots. Six passes in the forward and reverse directions were performed for the Bi-Sb-Se-Te ingot. The zone-melted ingots were found to contain several large crystals, with their basal planes (0001) approximately parallel to the growth axis. Samples of bismuth and antimony tellurides coated with a layer of tellurium, and bismuth selenide coated with a layer of selenium, were prepared for comparison in phase identification. These coatings were made by dropping a piece of the zone-melted ingot into some molten tellurium or selenium under argon atmosphere and allowing them to cool slowly to room temperature. The metallographic specimens were prepared by

Jan 1, 1969

By Peter Soo

Pvestrained Ferrovac E iron has been neutron-irradiated at approximately 90°C to an integrated flux of 1020 nut (E > 0.82 mev]. The irradiation was found to produce an incveased temperature dependence of the flow stress in addition to a greatly increased athemal stress. Measurements of the flow stress and stress relaxation, from which the activation volume and activation energy for slip were deduced, show that neutron irradiation changes the rate -controlling slip process to one based on dislocation interactions with tetragonal distortions which are Produced around submicroscopic interstitial loops in the lattice. The study indicates that without prestraining prior to irradiation the chances of detecting a change in the rate -controlling slip process are greatly reduced because in the initial stages of slip a substantial fraction of the radiation defects are swept out of the slip plane by gliding dislocations. Thus, activation parameters which are subsequently measured are representative of a greatly reduced defect density and would not differ appreciably from those for unirradi-ated material. The large increase in the athermal component of the flow stress is probably connected with the presence of depleted zones in the lattice which are introduced by irradiation. ALTHOUGH fast neutron-irradiation has not been observed to markedly alter the activation parameters for slip in bcc metals,' small but significant changes do occur. Most experimenters agree that irradiation predominantly increases the athermal component of the yield stress.'-= In addition to this, Laidler and smidt7 have shown that in iron irradiated to 5 X 10" nvt and molybdenum irradiated to 10" nvt, changes occur in the activation volumes for slip. A similar conclusion has been reached by Milasin and Malkin8 for irradiated iron. Work by Ohr et a1.5 shows that for Ferrovac E iron, irradiated to 1.2 X 1016 nvt, small increases in the activation energy for slip also occur. So far these changes in the activation parameters have not been explained on a firm theoretical basis. One important factor which would minimize the chances of detecting a change in the slip mechanism upon irradiation is the presence of "channeling" which has been observed in molybdenum,9 niobium,10 and iron.11These channels are formed by gliding dislocations which sweep irradiation defects out of the active slip planes and thereby create zones in which continued dislocation motion is encouraged. The activation parameters for the dislocations gliding in the defect-free channels would, therefore, be similar to those for unirradiated iron and a change in the rate-controlling slip process would be difficult to detect. In the present work, an attempt has been made to reduce the effect of uneven deformation on the measured activation parameters for slip in neutron-irradiated Ferrovac E iron polycrystals, so that a more realistic assessment of the effects of neutron-irradiation could be made. Primarily, the experiments involve the irradiation of specimens which had been prestrained to 9 pct elongation at room temperature prior to insertion into the reactor. It was hoped that the introduction of a large number of evenly distributed dislocations would substantially decrease any channeling effect which might otherwise occur. MATERIAL AND EXPERIMENTAL PROCEDURE The starting material was vacuum-melted Ferrovac E iron, an analysis of which is given in Table I. The standard tensile specimen had a gage length of 1.125 in., a cross-sectional diameter of 0.120 in., and a re-crystallized grain size of 1.2 x 10-3 in. All tensile tests were conducted on a floor model "Instron" tensile machine at a strain rate of 3 x 10-4 per sec. The irradiation of the prestrained specimens was performed in the Brookhaven High Flux Beam Reactor to an integrated flux of 1020 nvt (E > 0.82 mev) at a temperature of about 90°C. All specimens were excap-sulated in high-purity aluminum sheaths which were lightly swaged around the samples to ensure good thermal contact. Subsequent measurements on the irradiated specimens showed that within experimental accuracy the swaging had not deformed them. EXPERIMENTAL RESULTS Fig. 1 shows the flow stresses for a series of unirradiated control samples. In order to produce a comparable dislocation substructure throughout the test sm range, all specimens were prestrained

Jan 1, 1970

By O. J. Hassel, W. T. Maidens, J. H. Calbeck

The rotary gas fired reheating furnace used by the American Zinc Oxide Co. at Columbus, Ohio for Therotarygasfiredreheatingfurnacerefining lead-free zinc oxide is described. The outstanding features of this operation are that the color of the zinc oxide is greatly improved, sulphur is eliminated, and cadmium arethatrecovered without densifying the product to an objectionable degree. IN 1919 Leland S. Wemple obtained a patent for a process of reheating zinc oxide wherein the "coarsening of grain due to excessive heating was avoided." He taught in his specification that if solid carbonaceous material, such as lamp black, was added to the zinc oxide in proper amounts prior to reheating, objectionable sulphur compounds could be removed and the color would accordingly be improved and no objectionable densification would occur because of the relatively low temperature required. The situation that made this invention imperative was the newly opened zinc oxide plant of the American Zinc, Lead & Smelting Co. in Hills-boro, Ill. This was one of the early Western Type American Process zinc oxide operations. Characteristic of all of these early Western operations using Tri-State and Western ores was the great difficulty encountered in obtaining a product low enough in sulphur to compete with the Eastern Type American Process zinc oxides which were made from ores containing very low sulphur percentages. Wemple demonstrated that the refining process of his invention produced a superior color and although this was true and a most welcome feature, the primary purpose of the early refining operations at Hillsboro was to reduce substantially the high sulphur content of the crude zinc oxide. Although many and varied attempts had been made for refining zinc oxide none of the processes had a commercial history of any consequence until Wemple's invention became standard practice for the American Zinc, Lead & Smelting Co. in 1919 and their operations have been unique in that substantially all of their lead-free zinc oxide has been reheated since the first installation at Hillsboro. This process has become known in the industry as refining. The furnace developed by Wemple and continued in use by the company from 1919 until 1943 was unusual and merits some consideration by way of review in this paper. The furnace was essentially a double hearth coal-fired muffle furnace with a mechanical raking system consisting of a central shaft supporting six rabble arms in each muffle. The untreated or "crude" zinc oxide was fed onto the outer rim of the top muffle, moved to the center where it dropped to the lower muffle and progressed to the outer rim where it was discharged into an alloy screw conveyor. The retention in this furnace was extremely short, about 5 min, and the shallow zinc oxide bed on the hearths of the muffles was being continuously turned by the fast moving rabbles. Soft coal was burned on the grates below the lower muffle and the long yellow flame necessary to carry the heat around both muffles resulted in a very inefficient combustion of the fuel. The temperature of the top of the lower muffle seldom exceeded 65 °C although the oxide itself often reached 700°C before discharge. The capacity of this furnace was approximately 1/2 ton per hr. In our plant at Columbus it was necessary to keep four of these furnaces running in parallel to take care of the production because, as mentioned above, every pound of zinc oxide produced during these 24 yr passed through one of these refining furnaces. An essential part of this refining operation was the use of carbonaceous material admixed with the zinc oxide fed to the furnaces. Between 1 and 2 pct of a bran produced in the processing of cotton seed was added to all zinc oxide charged to the furnaces. The bran ignited on the top hearth and was still burning when the charge fell from the top hearth to the bottom hearth making a cascade of sparks. The rapid turning of the zinc oxide caused these particles of bran to flash on the hearths behind each rabble; but the combustion, of necessity, had to be complete by the time the charge reached the outer rim of the bottom hearth, otherwise the finished product would be contaminated with the charred particles of bran which would give the zinc oxide an unsatisfactory color. Although this operation was initiated to reduce objectionable sulphur percentages, as time went on new properties of the product were appreciated which made advisable continuing the refining process long after other methods of sulphur reduction became known in the industry. The particle size and particle size distribution, the absence of colloidal fines and perhaps a unique surface condition gave this product an outstanding performance when used in paints. The Wemple furnaces installed in Columbus in 1919 had to be rebuilt frequently and were extravagant in the use of fuel. The raking mechanism and the muffles required excessive maintenance expense and as the furnaces wore out the problem arose whether to continue along this line or to explore the possibilities of obtaining similar or better results in the simpler and more commonly used rotary furnace. To this end special research was initiated in 1941 on a small laboratory rotary

Jan 1, 1951

By D. T. Peterson, T. Carnahan

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

Jan 1, 1970

By T. E. Seidel, A. U. MacRae

The dependence of the electrical and crystalline properties of silicon containing ion implanted boron atoms have been studied as a function of the incident dose, substrate temperature, and annealing tempera-ture. Hall effect studies show that the compensation due to defects is negligible after annealing at and above 500°C for a dose of 1013 ions per sq cm. Steps occur in the isochronal annealing curve for the sheet resistance at 300°C and 800°C and the mechanisms responsible for these steps are discussed. The range of the ions at 150 kev is 3500 * 5001. Electron diffraction studies on samples implanted with 150 ker ions show the existence of damaged (polycrystalline) silicon closer to the surface than the range of the ions. Etching and subsequent annealing of etched and unetched specimens suggest that the poly crystalline silicon is closer to the surface than the bulk of the distribution of ions. THE purpose of this study is to investigate some electrical and crystalline properties of ion implanted boron in silicon. Sheet resistance, Hall measurement, depth distribution by chemical staining, thermal probe and electron diffraction have been used to study the properties of samples implanted in a random direction with ions having energy up to 150 kev. The boron-silicon system has been the subject of considerable study'-7 due to its possible application in the fabrication of electronic devices. The implanted ions produce considerable crystalline damage to the substrate crystal before they lose most of their primary energy and come to rest. Fortunately, most, if not all, of this damage can be removed by subsequent annealing or by holding the substrate at an elevated temperature during the implantation process. Sheet resistivity measurements of the p-type layers produced by boron implantation into room-temperature silicon, reveal that the annealing occurs in a two-step process: the first occurring at -300°C and the second at -800°C.' Electron-microscope and low-angle electron diffraction studies of this damaged material indicate that the bombarding ions produce an amorphous silicon region about the trajectory of each ion.' At high doses, these damaged regions overlap and an amorphous surface layer results. Furthermore, the 800°C annealing step, appears to be associated with the elimination of this gross damage. The depth distribution of the as-implanted and also the annealed samples have been studied by an angle lap and staining technique as well as by Hall measurements.3-7 We find the experimentally determined ranges of the boron are less than those calculated by the Lindhard-Scharff-Schiott theorys (hereafter referred to as the LSS theory). Our results are in essential agreement with pre- vious results. In addition, we have made Hall measurements on samples bombarded with a low dose (2 x 1013 ions per sq cm) to minimize impurity banding effects, and find that deep ionized levels influence the carrier concentration and room-temperature hole mobility until the sample is annealed at temperatures above 500°C. Above this temperature, the data corresponds to that obtained when silicon is doped with boron by conventional techniques. We have also made a depth distribution determination of the crystalline damage due to the implanted boron by using glancing angle electron diffraction techniques. It was found that the gross damage depth distribution of heavily implanted samples is closer to the surface than the depth distribution of the electrically active boron. By removing this gross damage, without affecting the distribution of the implanted boron, we have determined that the 800°C annealing stage is not exclusively due to the improvement in the crystalline perfection of the amorphous material but also appears to be associated with the annealing of some point defects. These results and their relation to previous results are discussed in the following sections. IMPLANTATIONS Boron ions were implanted into silicon samples using an accelerator capable of operating from 10 to 150

Jan 1, 1970

By Charles Smeal, Robert J. Schafer, Max Quantinetz

Normally metal powders cannot be ground to sub-micron sizes because of welding and agglomeration phenomena. Through the use of selected grinding aids and grinding fluids, nickel and other metal powders have been ball-milled as fine as 0.1 It was found that certain inorganic salts are more effective grinding aids for metal powders than conventionally used surfactants. METAL and alloy powders are used to produce reagents, pigments, coatings, solders, brazes, and parts for industry by powder metallurgy techniques. They are also combined with refractory compounds to produce cermets and dispersion-hardened products. One of the interests at the Lewis Research Center has been to explore the potentialities of the dispersion strengthening process. Since the work of Ir-mann, many investigators have shown that the strength of dispersion-hardened products may increase with decreasing interparticle spacing.24 One approach to achieving small interparticle spacing is to combine fine refractory compounds and metal powders, preferably below 1 in size. In attempting to obtain fine metal powders, it was found that until very recently the best that could be obtained from commercial suppliers, particularly of ductile metals, was about 1.0 . Interest was therefore developed in providing finer metal powders for dispersion-hardening studies. Information obtained from the literature, from others working in the field, and from prior experimental work performed at NASA, led to a consideration of ball-milling as a technique to produce the desired materials. Some of the variables associated with ball-milling are the size, material, and nature of construction of the grinding container; the nature and amount of the grinding material and material to be ground; and the nature and amount of the grinding liquid and grinding aid, if employed, and the grinding time. In all ball-milling, welding and agglomeration can oc- cur as well as grinding. Because of the tendency of ductile metals to weld together, they are difficult to grind.5 The ultimate particle size obtained on grinding is generally the one at which the rate of grinding becomes equal to the rate of welding. To help delay welding and thus obtain smaller particle sizes, grinding aids are often employed. The principal objective of this investigation was to produce submicron metal powders by ball-milling the powders with selected grinding aids and grinding fluids. A secondary objective has been to attempt to explain the variations observed in grinding behavior by considering possible grinding mechanisms and correlating various parameters with grinding effectiveness. Three groups of ball-milling experiments were run; one in which the grinding aid was varied, a second in which the grinding fluid was varied, and a third in which the material being ground was varied MATERIALS, APPARATUS, AND PROCEDURE In the first group of experiments in which the grinding aid was varied Inco Carbonyl Grade B nickel powder initially 2.5 (all sizes refer to average particle size as measured by Fisher Sub -sieve Sizer) was used as the material being ground and 200-proof ethyl-alcohol as the grinding fluid. Surfactants, representative of typical organic structures, were selected as grinding aids. Inorganic salts used as grinding aids were chosen on the basis of the size and valence of their ions. Water soluble salts were used in order to facilitate their removal from the slurry after grinding. In the second set of experiments, grinding of the 2.5- Ni powder was tried in four different grinding liquids; water, cyclohexane, n -heptane, and methy-lene chloride. In this study seven grinding aids selected from those tried in the first group of experiments were employed. In the third group of experiments 200-proof ethyl alcohol was again used as the grinding fluid to mill Cu, Cr, Fe, Ag, and Ni powders of various initial particle sizes' All mill charges contained 300 ml of grinding liquid and 3000 g of 1/2 in. stainless steel balls. When inorganic salts were used as the grinding aid, 70 g of salt and 210 g of metal powder were employed, and with surfactants 6 g of grinding aid and 300 g of

Jan 1, 1962

By F. H. Vitovec, D. W. Hoeppner

Pusk-pull fatigue tests were conducted on copper tricrystals of 99.988 pct purity to ascertain the role of grain boundaries in the initiation and propagation of fatigue cracks. Significant differences in behavior were found for specimens which possessed different transverse-boundary misorientation. In speciwens with low boundary angles cracks initiated within the transverse boundary, while higher angles led to transcrystalline fatigue failure. It is suggested that at low angular misorientation moving dislocations may interact with dislocations of the boundary or dislocations present in adjacent pains on favorably oriented glide planes, thus initiating a fatigue crack. MANY fatigue studies have been concerned with fatigue-crack initiation within grains and the mechanism causing initiation and propagation1-3 Although the initiation of fatigue cracks in or near grain boundaries of pure metals has been observed and reported in the literature, the mechanism of this phenomenon has received little attention.4"11 EXPERIMENTAL PROGRAM Testing Procedure. To investigate the role of grain boundaries regarding initiation and propagation of fatigue cracks, copper tricrystals were tested in push-pull. Axial-stress tests were used to avoid the stress gradients introduced by stressing of some other nature. Copper was selected as the most suitable test material since extensive work has been done on the formation of fatigue-induced slip in copper. Tricrystal specimens were used to provide a grain boundary geometry with one boundary transverse to the principal stress axis and one or more boundaries nearly coincident with the direction of maximum resolved shear stress. The boundary energy and the slip characteristics in the vicinity of the transverse boundary depend on the relative orientations of grains across the boundary. Testing was done at room temperature, at a frequency of 700 cpm, in an atmosphere of either air or argon. It is known that the incidence of grain boundary cracking increases with increasing temperature and decreasing frequency.4 The fatigue machine applied a uniaxial tension-compression load to the specimen by a flat spring which was actuated by an adjustable eccentric. Use was made of an adjustable head and a load cell which had strain gages mounted at 120-deg intervals around its periphery, thus providing a means for eliminating any detectable superimposed bending moment. A clip gage was mounted between the gripping heads to record the cyclic strain amplitude applied to the specimen. Each specimen's hysteresis characteristics were recorded by supplying the load and strain signals to an oscilloscope. Microstructural changes were observed and recorded with an optical microscope which was mounted on the fatigue machine. Immediately prior to insertion in the machine, each specimen was chemically polished. Extreme care was exercised while inserting the specimen in the test machine to avoid either bending the specimen or introducing a mean load. Each specimen was stressed in the tensile direction first and subsequently the load was reversed. Specimen Preparation. Tricrystal fatigue specimens of 99.988 pct Cu were grown from the melt using a modified Bridgman technique. Ingots about 1 in. high were grown to the configuration shown in Fig. 1. Upon sectioning the resultant ingot, several specimens were provided with an identical shape and grain boundary orientation. Except for plane sectioning and final polishing, this method eliminated machining the specimens. The apparatus used for growing the tricrystals is shown in Fig. 2. A spectrographically pure graphite mold, Fig. 3, was inserted in a vycor tube which was mounted on a vertical zone-refining table. Prior to insertion in the tube the mold was assembled as follows. The insert and graphite rods were fixed in position in the mold. The copper was then placed in the mold and the entire mold assembly was positioned in the vycor tube, Fig. 2. At this point the tube was sealed, evacuated, purged with helium gas, and finally placed under a slightly positive helium pressure. The heating coil was then adjusted to melt the copper charge, after which it was raised at a speed of 4 in. per hr. It is also possible to use

Jan 1, 1964

By V. Leroy, A. G. Guy

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

Jan 1, 1967

By C. H. Stewart, R. A. Fagin

A two-dimensional, three-phasereservoir simulator was programmed for a large memory digital computer. It was designed to provide a practical solution to describing the complex physical relation between the natural forces and the physical properties of a heterogeneous reservoir when subjected to a specific set of conditions. A reservoir study is briefly described to illustrate application of the model. A full volumetric account of three phases (oil, gas and water) is performed simultaneously throughout an integration net representing the reservoir. Absence of one or two of the phases is treated as a special case of the more general situation. Expansion (or contraction) of all phases, including rock expansion, is performed so that the pressure calculation is the general unsteady-state case. To account for the large variations of subsurface elevation encountered in some reseruoirs, and to allow for segregation of the various phases, a gravity head term is included in the basic drive potential. Appropriate fluid and rock properties are used in polynominal surface form (functions of pressure and/or depth) or they can be entered as space variables at each position of the integration net. An unsteady-state water influx calculation, based on the method of van Everdingen and Hurst, l was connected to the boundary of the matrix to simulate aquifers of various sizes. In addition to reservoir calculations, three-phase fluid flow from the producing depth to the wellhead, including provisions for gas lift, was incorporated in the simulator. A workover routine was also built which can automatically switch to a different set of production relations when a gas-oil ratio or water fraction reaches a limit; or it can shut-in the well if prescribed. INTRODUCTION This paper describes a reservoir engineering mathematical simulator used to represent the complex interaction of natural forces and physical properties of a reservoir during natural depletion or with various injection schemes. The simulator, which was programmed for a large memory digital computer, is a two-dimensional calculation which handles three mobile fluid phases simultaneously (oil, gas and water). Basic requisites for the method are individual well production and pressure data, hydrocarbon fluid properties, geological data (producing depth and net sand), capillary pressure data, relative permeability data and permeability and porosity information. Matching the past performance of a combination drive reservoir often has yielded information concerning continuity and the validity of basic data. Detailed predictions of future performance can be made for continuation of current depletion methods (natural depletion) as well as for various types of recovery by gas or water injection. Combination injection cases and pattern studies can also be performed. Workover programs, gas lift and different types of artificial lift programs have been investigated using a technique similar to that described by Kern and Nicholson2 except that conditions of pressure and saturation at the block within which the well is located are used rather than average reservoir conditions. Drilling additional wells to optimize profit was explored, both as to number and location, by placing wells at different spots within the reservoir matrix. Special depletion processes can be examined, such as upstructure drainage and lateral (or strike) waterfloods in thin oil columns. In one case the mathematics of the simulator were modified to calculate the displacement in the vertical plane rather than in the horizontal plane. In this manner specific reservoir problems can be studied, such as coning of gas and/or water around production points, fingering along permeable stringers or, more generally, frontal advances in a heterogeneous section.

Jan 1, 1967

By E. F. Reed

CONSIDERABLE deep hole prospect drilling has been done in the last few years in the Globe-Miami mining district about 70 miles east of Phoenix, Arizona, and in the San Manuel-Tiger area about 50 miles south of the Globe-Miami region. More than 205,000 ft of churn drilling have been completed by the San Manuel Copper Corp. at their property in the Old Hat Mining District in southern Pinal County. The deepest hole on this property is 2850 ft; there are 49 holes deeper than 2000 ft. At the adjoining Houghton property of the Anaconda Copper Mining Co., where only one hole reached 2000-ft depth, there were 27,472 ft of churn drilling and 3436 ft of diamond drilling. Three churn drill holes were deepened by diamond drilling methods. Near Miami in the Globe-Miami district the Amico Mining Corp. drilled four holes by combined churn and rotary drilling methods, the total amounting to 13,879 ft, of which 2256 ft were drilled with a portable rotary rig. In the same district, besides doing a large amount of shallow prospect drilling, the Miami Copper Co. drilled two holes of 2560 and 3787 ft, respectively, which were completed by churn drilling methods. The rocks encountered in drilling at San Manuel and Tiger are described by Steele and Rubly in their paper on the San Manuel Prospect' and by Chapman in a report on the San Manuel Copper Deposit? The rocks are well-consolidated Gila conglomerate, quartz , monzonite, and monzonite porphyry. In some places these formations stand very well while being drilled, and three holes were drilled without casing, the deepest of which was 2200 ft. In other holes faulted and fractured ground made drilling difficult. In the Globe-Miami district the deep drilling was done in the down-faulted block of Gila conglomerate east of the Miami fault and in the underlying Pinal schist. The geology of this area is described by Ransome.3 In the Amico holes the conglomerate varied from material consisting entirely of granite boulders and fragments to a rock made up of schist fragments in a sandy matrix; in the Miami Copper Co. holes there were more granite boulders and the material was poorly consolidated. Drilling was much more difficult and expensive in the Miami area than in the San Manuel district, mainly because of the depth of the holes and the formations drilled. All the deep hole prospecting described in this paper was done with portable rigs. The churn drill rigs were of several types, of which the Bucyrus-Erie were the most popular. Bucyrus-Erie 28L, 29W, and 36L rigs were used on some of the deeper holes on the San Manuel property. A few Fort Worth spudder types were tried, and the deepest hole at San Manuel was drilled with a Fort Worth Jumbo H. The spudder type is considerably larger than most other rigs used on this work and required a larger location site. The spudders were belt-driven machines with separate power units, and time required for setting up and moving was much longer than with the more portable drills. All the churn drilling was done by contractors or with machinery leased from them. A few of the contractors had complete equipment, including most of the necessary fishing tools. Unusual and special, fishing tools were obtainable from the supply companies in the oil fields of New Mexico or in the Los Angeles area. Most of the contractors used equipment with standard API tool joints, so that much of it was interchangeable. Failure of tool joints is one of the principal causes of fishing jobs. It can be minimized if the joints are kept to the API specifications and the proper sized joints are used in the various holes. The minimum sizes that should be used with various bits are as follows: 12-in. and larger bits, 4x5-in. tool joints; 10-in. bits, 3 1/4x4 1/4-in. tool joints; 8-in. bits, 2 3/4x 3 3/4-in. tool joints; 6-in. bits, 2 1/4 x3 l/4 -in. tool joints; 4-in. bits, 1 5/8 x2 5/8-in. tool joints. Two rotary drill rigs were tried at San Manuel on the same hole, and a portable rotary drill rig was used on the Amico drilling for test coring the formation and for drilling in holes 3 and 4. Rotary drilling differs from churn drilling or cable tool drilling in that the bit is revolved by a string of drill pipe and the cuttings are removed from the hole by a thin solution of mud pumped through the drill pipe. The principal parts of a rotary rig are the power unit, a rotating table to revolve the drill pipe, hoists to raise and lower the pipe and to handle casing, and a pumping system to circulate the drilling liquid. The rig used on the Amico property at Miami was mounted on a truck. The larger rig used on the San Manuel property was hauled by several trucks and had separate turntable and pumping units. Diamond drill coring equipment was used successfully with the rotary rig in the holes on the Amico property, To allow for 2 3/8-in. drill pipe with tool joints, 3 1/2-in. core barrels and bits were used. With the standard 3 1/2-in. core barrel there was considerable difficulty in maintaining circulation with mud, so a barrel was designed with a smaller inner tube and a broad-faced bit. This allowed coarser material to circulate between the barrels. Rock bits of 5 5/8 to 3 7/8 in. were used with the rotary rig for drilling between core runs. Diamond drill equipment is much lighter than churn drill tools, so that fishing tools can usually be obtained from supply houses by air express when needed. Three churn drill holes on the Houghton property at Tiger were deepened by diamond drilling with Longyear UG Straitline gasoline-driven-machines. The open churn drill hole was cased with 2 1/2-in. black pipe. In deep hole churn drilling, casing is one of the most important items, especially in drilling in unconsolidated material like the formations drilled by

Jan 1, 1952

By E. A. Gulbransen, K. F. Andrew

Weight loss measurements were made using a sensitive microbalance operating in a high vacuum system. The Langmuir equation was used to Calculate the vapor pressures of the several metals and alloys. Aluminum lowered the vapor pressure of i9,on in a 5.4 Al-94.6 Fe alloy, and of iron and chromium in a 4.8 Al-21.5 Cr-73.7 Fe alloy. The influence of oxide films formed on the alloys on the effective vapor pressures of the alloys was also studied. These studies were used to aid in the inte9,pretation of the high temperature oxidation of alloys based on iron, chromium, and aluminum. RECENT studies on the oxidation of chromium1 and the heat resistant alloy 5 Al-22 Cr-73 Fe showed transitions in the rate of oxidation at 900o and 1050oC respectively. The transition for chromium occurred at a temperature where the rate of evaporation of chromium from oxide-free surfaces equalled the rate of chromium atoms reacting with oxygen. This paper presents new vapor pressure studies on iron, chromium, and the specially prepared alloys 5.4 Al-94.6 Fe, 21.9 Cr-78.1 Fe, and 4.8 Al-21.5 Cr-73.7 Fe. The ternary alloy is similar in basic composition to the patented heat-resistant commercial alloy known as Kanthal. The binary alloys were studied to show the individual effects of aluminum on the vapor pressure of iron and of iron on the vapor pressure of chromium. The purpose of the proposed studies was to test the influence of metal volatility on the oxidation behavior of heat-resistant alloys based on iron, chromium, and aluminum. Since oxide films formed on the metal may limit metal transfer to the surface, the influence of oxide films on metal volatility was also studied. LITERATURE A) Vapor Pressure. The vapor pressure of iron has been studied by Jones, Langmuir, and Mackay,' Dornte and Nrton. Edwards. Johnston. and Ditmar. and McCabe, Hudson, and Paton. Stull and SinkeT have calculated a heat of sublimation at 298° K of 99.83 kcalper g atom. The vapor pressure of chromium has been reported by Bauer and Brunner,' Speiser, Johnston, and Blakburn, Gulbransen and Andrew,'' Vintaikin,o McCabe, Hudson, and Paton, and Kubaschewski and Heymer.12 Good agreement was obtained except for the early work of Bauer and Brunner.' Stull and Sinke7 calculate a heat of sublimation at 298°K of 95.0 kcal per g atom. Vapor pressure measurements on aluminum have been made by Brewer and Searcy,13 Bauer and Brunner,aand Farkas.I4 Stull and Sinke,7 giving Brewer and Seary's' data the most weight, derive a heat of sublimation of 77.5 kcal per g atom at 298 oK. B) Method. Two methods15 are used for measuring the vapor pressure of metals: 1) The Langhuir free evaporation method, and 2) the Knudsen effusion method, In the Langmuir method, used in the present work, the vapor pressure of the metal P is given by the equation: Here M is the molecular weight of the vapor species. T is the absolute temperature. (dw/dt) is the rate of sublimation in g per sq cm per sec and a is the condensation coefficient. a is a measure of the efficiency of condensation of molecules striking the metal surface from the vapor. If condensation results from each collision a is unity. This coefficient has been found to be unity,a, le for metals. These results establish the general validity of the Langmuir free evaporation method as applied to metals. EXPERIMENTAL A) Vacuum Microbalance Method. A modification of the Langmuir free evaporation method was used. Strip specimens were suspended from a sensitive microbalance operating inside of a high vacuum system. For alloys, microbalance methods17 must be used on large area samples at relatively low temperatures to minimize composition changes. The

Jan 1, 1962

By P. L. Rittenhouse, M. L. Picklesimer

A rapid and seniquantitative method of determining prefered orientation on large numbers of. Zircaloy-2 specimens was desired. knoop microhardness measurerrzetzls were irvestigated as a solldtion to this pro6lem. The variation of Knoop microhardness measurerlerzts on selected planes as a function of 'irzdenle axis relutilje lo cryslallograplzic or fabrication divectals were used to corzstrrcl a polar coordinate hardness contour map. With use of an empirical relationslip between the single-crystal hardnesses and those of the polycrystalline material conentional pole figures could be constructed which compare favorably also obtained. To determine preferred oriention qualitatively from hardness data requires a minimum of twelve measurements per plane on three, preferattention to grain size, specimen prepartion, and in-rreinetzt, and analjlsis is of the order of- 45 to 60 rrin. THE design of structures from Zr alloys requires consideration of preferred orientation and the resulting anisotropy of mechanical properties. Rapid and semi-quantitative methods of evaluating anisotropy and of determining preferred orientation are needed for quality control and for examining the large number of test specimens required in development programs. The variation of Knoop microhardness with indenter orientation has been studied in single crystals of several hcp metals including titanium.' beryllium,' magneium, and zinc.4 The maximum hardness observed on any crystallographic plane other than the basal plane invariably occurred when the long axis of the Knoop indenter was either parallel or perpendicular to the projection of the [ OOOI.] direction on the plane of examination. When the major slip mode was (0001)(i2i0j the maximum hardness occurred at the parallel positio! but when the operating slip system was {10i0) (i210) the hardness was a maximum perpendicular to the [0001] projection. A minimum hardness was always observed at a rotation of 90 deg from the maximum. It follows that the orientation of the projection of the [0001] can be determined on any non-basal plane by making hardness measurements at a number of indenter orientations. The determination of the orientation of the [0001] projection on two non-parallel, preferably orthogonal, surfaces of the specimen will allow location of the [0001] direction in the specimen. It seems possible that such measurements could be used to examine, at least qualitatively, the preferred orientation existing in polycrystalline hcp materials. An investigation of the microhardness anisotropy in Zircaloy-2 was undertaken to ascertain whether these measurements could be used for this task. EXPERIMENTAL PROCEDURE Single crystals of Zircaloy-2 were grown by an a-8-a annealing sequence using electron-beam heating.= The orientation of the crystals was determined by a conventional back-reflection Laue technique. The crystals were then mounted on a goniometer head and the desired crystallographic faces were milled and chemically polished. Polycrystalline Zircaloy-2 specimens were prepared from fabricated sheets or plates. Inverse pole figures for these materials were obtained using the X-ray diffraction techniaue described by Jetter, McHargue, and illiams.' A Wolpert-Greis Micro-Reflex hardness-testing machine was used to make the Knoop microhardness measurements. The single-crystal and polycrystalline specimens were loaded to 0.5 and 2.0 kg, respectively. Determination of the scatter of hardness numbers as a function of applied load for several indenter orientations on several specimen surfaces showed that the loads selected were the lightest consistent with minimum scatter. Heavier loads did not appreciably decrease the scatter and produced hardness impressions too large to be conveniently measured with the equipment used. Seven crystallographic planes of the single crystals were examined, while six planes of examination were used in studying Zircaloy-2 polycrystals. The specimens were rotated 10 to 15 deg after each measurement and from four to eight impressions were made at each angle of rotation. RESULTS The two angles which were used to relate crystallographic directions and planes in the hcp cell to the planes of examination of the single crystals are shown in Fig. 1. 8 is the angle between the c axis. [0001]. and the normal, N, to the plane of examination. a is the angle between the long diagonal of the Knoop indenter and the projection of the c axis on the plane of examination. For exampIe, if P = 90 deg and a = 0 deg, the plane of examination is of the family

Jan 1, 1967

By R. A. Vandermeer, J. C. Ogle

The preferred crystallographic orientations (texture) developed in randomly oriented, poly crystalline niobium during rolling were studied by means of X-ray diflraction techniques. The evolution of texture at both the surface and center regions of the rolled strip was carefully examined as a function of increasing defamation throughout the range 43 to 99.5 pct reduction in thickness. Certain aspects of the center texture development in niobium are in agreement with the predictions of a theory by Dillamore and Roberts, but others cannot be explained by the theory in its present form. Above 87 pct reduction by rolling, a distinctly different texture appeared in the surface layers which was unlike the center texture. The present results are compared with previous results obtained from other bcc metals and alloys. RANDOMLY oriented, poly crystalline metal aggregates when plastically deformed to a sufficiently large extent develop preferred orientations or textures. In a recent review article, Dillamore and Roberts1 pointed out that the nature of the developed texture may be influenced by a large number of variables. These include both material variables such as crystal structure and composition and treatment variables such as stress system, amount of deformation, deformation temperature, strain rate, prior thermal-mechanical history, and so forth. From a practical point of view, the control of preferred orientation may often be important for the successful fabrication of metals into usable components. During the past few decades many experiments have been devoted to the study of deformation textures. This work, however, has been confined in large part to metals and alloys that have an fcc crystal lattice. By comparison, bcc metals and alloys have received much less attention, and consequently our understanding of preferred orientations in these materials is only shallow. This state of affairs worsens when it is realized that almost all of our present howledge about this class of materials derives from studies on irons and steels.' The bcc refractory metals, which are relative newcomers to the industrial world, have, on the other hand, been given at best only passing glances in the area of texture development. Our understanding of the evolution of preferred orientations in bcc metals can only remain fairly limited until systematic studies of metals and alloys other than the irons and steels have been carried out and the influence of the many variables has been determined. To that end a program was initiated to investigate in detail texture development in niobium. The present paper reports some of the results of this study. Textures were determined at both the center and surface of strips rolled variously to as much as 99.5 pct reduction in thickness at subzero temperatures. Emphasis in this paper is on texture description and on texture evolution during rolling to progressively heavier deformation. EXPERIMENTAL PROCEDURE The niobium was purchased from the Wah Chang Corp. as a 3-in.-diam electron-beam-melted billet. Chemical analysis indicated the impurities to be less than 300 ppm Ta, 40 ppm C, 10 ppm H, 170 ppm 0, and 110 ppm N. All other impurities were below the limits of detection by spectrochemical analysis. This large-grained billet was fabricated into specimen stock so that a fine-grained randomly oriented grain structure resulted. This was accomplished in three deformation steps alternated with recrystalli-zation anneals of 1 hr at 1200°C in a vacuum of low 10"6 Torr range after each deformation step. The first step was to alternately compress the billet 10 to 20 pct in each of three orthogonal directions. The second step was to compress in only two directions 90 deg apart to produce a 2-in.-sq bar. The final step was to roll this bar 50 pct to give a 1-in. by 2-in. cross section. After the final anneal, metallo-graphic examination showed the material to have an average grain size equivalent to ASTM No. 5 at 100 times (i.e., 0.065 in. diam). Specimens cut from the center and edges of this bar gave no indication of detectable preferred orientation when examined by X-ray diffraction. Samples 1.5 in. long, either 0.625 or 0.750 in. wide, and approximately 0.400 in. thick were machined from this fabricated ingot. The surfaces corresponding to the rolling planes were ground so as to be parallel. The samples were chemically polished in a solution of 60 pct nitric acid and 40 pct hydrofluoric acid (48 pct solution) prior to rolling to remove any cold work introduced in the machining operations. Rolling was accomplished with a 2-high hand-operated laboratory rolling mill that had 2.72-in.-diam rolls. Prior to operation, the rolls were polished with 600 grit paper, cleaned with acetone, and then soaked in a container of liquid nitrogen for several hours. The samples were also soaked in liquid nitrogen prior to rolling and were recooled between each pass. While some slight heating of the samples occurred during rolling, this procedure maintained the sample temperature well below 0°C at all times. The samples were rolled unidirectionally, and the rolling plane surfaces were not inverted during any phase of the operation. The draft per pass averaged between 0.010 to 0.012 in. After 96 or 97 pct reduction the draft was reduced to 0.001 to 0.002 in. per pass. Samples were rolled to various reductions in thickness between 43 and 99.5 pct.

Jan 1, 1969

By O. D. Gaither

This paper is an analytical study of the flow of fluids through small vertical conduits. Small conduits are defined as 11/4-in. nominal diameter tubing size and smaller, and approximately twice this area for annular conduits (i.e., 1- X 21/2-in. annulus and smaller). Experimental data are presented for the 1-X2-in. and 11/4- X 2%-in. annuli, and the I-in. and 11/4-in. tubing, since these represent the small conduit sizes and configurations generally encountered in oilfield applications. Data have been gathered for these conduits for single-phase water, single-phase gas and two-phase water-gas mixtures, with particular emphasis on high gas-liquid ratios. Water rates in excess of 2,000 BID and gas rates in excess of 2.5 MMcf/D, and two-phase flow ratios in between these two, represent the scope of the data gathered. Existing equations have been applied to predict flowing pressures and compared with experimental data. New correlations have been developed. INTRODUCTION The increased economic pressure on the domestic oil industry in the United States has constantly required the use of new techniques and equipment designed to reduce the cost of finding and producing oil and gas. Since tangible items are most readily apparent in economic analysis, the advent of lower-cost well completions was inevitable. One of the methods used to reduce costs which has received widespread attention is the slim-hole completion technique where tubing is used as the well casing and in which small conduits are used for tubing if necessary. Small conduits, defined by Kirkpatrick1 as "11/4-in. diameter nominal tubing and smaller for tubing flow and less than twice the 11/4-in. diameter nominal tubing internal flow area for annulus flow", have also found widespread usage as siphon strings for de-watering gas wells and as "kill" strings in deep high-pressure oil and gas wells. The growing use of small-diameter tubing has resulted in an increased need for development of improved methods to measure or predict flowing bottom-hole pressures since the physical dimensions generally preclude the use of subsurface-recording pressure gauges. Even in the cases where small bombs are available, the relatively high velocities encountered at nominal flow rates make it necessary to use excessive weight bars or special hold-down devices. Attempts to use recognized correlations to accurately predict flowing or gas-lift performance in wells equipped with small conduits have been generally unsuccessful. Insufficient field data were available to allow the development of a correlation on this basis, and an experimental approach was applied in an attempt to obtain a workable relation. The experimental approach used to obtain the data presented in this paper was actually a compromise between a field installation and a laboratory study. A test well 1,000 ft in length was used to obtain flow data on single-phase liquid, single-phase gas and two-phase water-gas flowing mixtures. Liquid rates up to 2,200 B/D and gas rates up to 3 MMcf/D were used in the single-phase flow studies. Two-phase flow rates from 100 to 600 B/D with gas-liquid ratios from 500 to 8,000 cu ft/bbl were recorded. Experimental data were obtained for single- and two-phase flow through 1-in and 11/4-in. nominal tubing, and through the annuli between 1- and 2-in. and 11/4- and 2%-in. nominal tubing strings. Experimental results for the two-phase flow are compared to the Poettmann-Carpenter correlation2 which is widely used as a comparative standard for development of multiphase flow predictions in flowing and gas-lift wells. Correlations developed by Tek,3 Baxendell and Thomas" were also investigated. The experimental data recorded herein fell in between the two flow regimes as defined by Ros," and this correlation also failed to yield satisfactory results. The fact that existing correlations failed to confirm the experimental data led to the need for development of a new correlation. Although a two-phase flow study was the primary objective of this investigation, data were also recorded for single-phase flow of water and gas, and constants were developed relating to pipe roughness and equivalent diameters for annular flow. These single-phase studies assisted materially in the development of certain of the two-phase flow results. Considerable previous work has been published which presented relationship of surface measurements to bottom-hole condition. The works of Buthod and Whiteley,6 Jones,' Poettmannb and the Texas Railroad Commission" are classic examples of the successful use of mathematical relationships which allow acceptable predictions of subsurface pressures, when gas is the flowing fluid. Darcy and others have derived relationships which may be used with minor modifications to predict subsurface flowing conditions in injection and water-supply wells. As previously stated, the application of the single-phase flow relationships