Search Documents

By Irving Cadoff, Richard Nolder

An analysis oj a series of samples of single -crystal silicon grown on sapphire shows that ,four distinct orientation relations exist. There are at least thirteen crystallographic planes which serve as substrates ,for epitaxial growth The model which assumes that silicon substitutes for aluminium and bonds with the oxygen atoms of the sapphire to form the first later for subsequent growth can account for the epitaxial geometry in all four cases. Multiplici-ties 01 specific orientations are found which are explained by the symmetry of both the sapphire and silicon crystal structures. It was reported in the previous paper' that single-crystal silicon films can be deposited on single-crystal a A1,03 (sapphire) with a wide variety of orientations. The silicon films, however, were found to deposit in only four different orientation relationships with respect to the sapphire substrate. The X-ray data have been analyzed to determine the plane and direction of atomic fit and an atomic model to account for the epitaxial relationships at the silicon-sapphire interface has been developed. These results are described in the present paper. EXPERIMENTAL The Laue X-ray data described in the previous paper were reduced to stereographic projections in a standard manner. These data were supplemented with X-ray reflection data obtained with -& E and A Full-Circle Goniometer (Electronics and Alloys, Inc.) mounted on a Phillips diffractometer unit. The silicon films were thin enough to result in characteristic reflections from the substrate as well as the silicon. The superimposed stereographic projections were analyzed to determine the orientation relationships. Ambiguities were resolved by referring to constructed models of various lattice planes for both silicon and sapphire. These lattice models were also used to develop the model for atom fit at the interface. DATA AND ANALYSIS 1) Crystal Structures. The crystallography of silico-n, which is of the diamond cubic class with a. = 5.4~, is straightforward. The sapphire structure is more complex and generally difficult to represent in pictorial form. It is of the rhombohedra1 class but usually referred to hexagonal axes. The transformation of indices between these two systems can lead to ambiguities which will have bearing on the analysis presented herein. Therefore, the basis used in this paper is described in the appendix. Referred to the hexagonal system, the parameters of sapphire are co = 12.95.&, a, = 4.75~, c/a = 2.73. All parameters are for room temperature. Figs. 1 and 2 illustrate schematically the hexagonal unit cell model of sapphire used for this analysis. The oxygen atoms are assumed to occupy positions in an idealized hexagonal array on basal planes. Their actual positions are in a slightly distorted hexagonal array. The aluminum-atom positions are scaled to their actual positions in the lattice with respect to the basal layers of oxygen atoms. It was found that an idealized model where the aluminum atoms are assumed to lie exactly at the octahedral interstices of the oxygen lattice, similar to that used by ~ronberg,' was not adequate for this study.

Jan 1, 1965

By C. A. Bruch, J. F. Duffey

On the basis of metallographic and cooling-rate studies of 15 U-Zr alloys, a new version of the 6 phase region of the U-Zr phase diagram is proposed. IN several published diagrams of U-Zr alloy system:' there is good agreement that at elevated temperatures complete solid solubility exists between y uranium and ß zirconium. Similarly, there is reasonably good agreement concerning the uranium-rich and zirconium-rich ends of the diagram. However, there is no agreement concerning the nature of the diagram at intermediate compositions and at temperatures below about 650°C. The early work of Saller et al.8 hows an intermediate phase, 6. In a later diagram by Saller and Rough' shown in Fig. 1, the d phase field is drawn with dotted lines, indicating doubt that this phase exists. The diagram of Summers-Smith2 does not show the d phase. More recently Holden and Seymour' have shown that in the composition range 40 to 60 wt pet Zr, the high-temperature ? phase, which is body-cen-tered-cubic, transforms to an ordered structure, 8. They found that the order-disorder transformation occurs between 595° and 610°C, and presented a diagram based on the study of four different alloy compositions. Most recently Saller et a1.5 presented a new diagram which differs with respect to the composition limits of the 8 field from the others proposed by these authors. The purpose of the present investigation was to determine the boundaries of the d phase field. Experimental Procedures and Results Alloys—Buttons, weighing 80 g of 15 different alloy compositions were prepared by arc melting Mal-linckrodt uranium and iodide-process zirconium in a water-cooled crucible. A controlled atmosphere of argon and helium was maintained in the furnace. Prior to making each button, a sacrificial button of zirconium was melted to increase the purity of the furnace atmosphere. Each alloy was then melted six times to improve its homogeneity. Subsequently the alloy buttons were jacketed in copper and hot rolled to slabs at 775°C. From these slabs, specimens were machined for chemical, thermal, and metallographic analyses, Table I gives both the estimated chemical composition of the alloys, and the values which were obtained by analysis. The oxygen contents varied from 0.013 to 0.032 pet, nitrogen from 0.0014 to 0.014 pet, and carbon from 0.0013 to 0.011 pet. The purity of the alloys ranged from 99.947 to 99.975 pet. Thermal Analysis—The specimens for thermal analysis were blocks 11/4 x 1/4 x % in. The bead of a calibrated thermocouple was placed between two of these blocks and molybdenum wire was wrapped around the assembly. Each assembly was sealed in a quartz capsule having an argon atmosphere. The capsule was placed inside a refractory crucible which was heated in a furnace. Using a method described by Smith,' the emf of a differential thermocouple, with junctions inside and outside of the refractory crucible, was used to activate the

Jan 1, 1959

By D. Havlena, A. S. Odeh

The use of the straight-line method of solving the material balance equation is illustrated by means of six field cases. Also, the application of statistical criteria to arrive at the most probable answer is shown. The theory underlying the straight-line method of solution and the applicability of the statistical criteria was presented in a previous paper.' The field cases include saturated and undersaturated oil reservoirs with and without water drive. The aquifers discussed are: limited radial, infinite radial, very small aquifer and infinite linear. The field cases also include a gas reservoir producing under water drive. INTRODUCTION In a previous paper,' the authors presented the theory underlying the solution of the material balance as an equation of a straight line. The appropriate equations for various material balance cases as well as the methods of analysis and interpretation with comments and discussion were also included. To illustrate the various theoretical cases treated previously, five field cases are analyzed in this paper by employing the straight-line method of solving the material balance equation (MBE) and one example previously published is referred to. The use of statistical criteria to arrive at the most probable answer is also shown. All the field examples presented, except Case 2, are excerpts from complete reservoir studies. To illustrate the method, only sections specifically dealing with the material balance principles are included. Additional geologic information and basic data are reported to better acquire an understanding of the cases and thus to better follow the reasoning that suggested the successful application of the straight-line method of solving the MBE. The six cases are: (1) saturated reservoir, small gas cap, limited aquifer; (2) saturated reservoir, very small gas cap, infinite aquifer; (3) undersaturated-saturated reservoir, very small aquifer; (4) highly undersaturated reservoir, no water drive; (5) high undersaturated one-well reservoir, limited aquifer; and (6) gas reservoir, infinite linear aquifer. THE D4 SAND, GUICO FIELD, VENEZUELA The D, sand, which was discovered in 1943, is presently in a depleted state. Since its discovery it has produced under water drive, gas-cap-gas expansion, and solution gas drive. In Nov., 1947, water injection was initiated to arrest further pressure decline. When discovered, the D, sand was a saturated reservoir with a gas cap /oil zone volume ratio m estimated volu-metrically at 0.0731, an average permeability of 500 md, a porosity value of 25 per cent, and an oil viscosity at reservoir conditions of 0.3 cp. The volumetrically determined stock-tank oil initially in place was 23.1 million bbl. The volumetrically weighted physical data and production data available until Nov., 1953 are reported in Table 1. In Ref. 1, the effects on the straight-line plot of various values of ro/r for a constant ?to, or of various dimen-sionless times for a constant r8/r10 were theorized and were illustrated in Fig. 3A of that reference. In this field case, the previously theoretically predicted effects are established. Thus, the MBE calculations using Eq. 3c of Ref. 1 were performed for various r6/r10 and dimension-less time values. Eq. 3c of Ref. 1 is:

Jan 1, 1965

By B. S. Gottfried

A generalized mathematical model is presented, which describes the thermal recovery of oil in linear systems with convective external heat loss. Three-phase fluid flow, conduction-convection heat transfer, chemical reaction between oxygen and oil and aqueous phase change are included, making the model applicable to a variety of thermal recovery processes. The complexity of the model requires that numerical solutions be obtained with a digital computer. The numerical approximations of the partial differential equations are presented and the algorithm is discussed in detail. The analysis is applied to the detailed simulation of a forward-combustion laboratory tube experiment. The results of the simulation exhibit with reasonable accuracy all of the principal thermal and bydrodynamic characteristics observed in the laboratory (i.e., the propagation of a combustion zone, a steam plateau, and oil and water banks). Numerical variations in the parameters appearing in the external heat-loss term, the aqueous phase-change term, and the Arrhenius combustion term are examined, and their effects on the predicted temperature profiles, oil and water production histories, and the air-oil ratio are presented. In spite of an oversimplification in the chemical reaction mechanism and the considerable computer time required to obtain solutions, the model represents a major advance in the ability to simulate the phenomena observed in thermal recovery experiments. INTRODUCTION Thermal oil recovery refers to a class of recovery processes where heat is supplied to a reservoir to provide the necessary expulsive energy. This thermal energy can be supplied externally as steam or hot water, or it can be generated in situ by forward or reverse combustion. In either case, however, thermal recovery processes are charac- terized by the simultaneous flow of two or three fluid phases in a variable - temperature field, accompanied by possible chemical reaction or phase-change effects.2-5, 7- 10, 14,16-18,20-22 Although a physical understanding of the thermal recovery processes is far from complete, it is possible to construct mathematical models which describe approximately all of the principal physical and chemical phenomena. However, attempts to solve such models, even with high-speed computers, involve formidable mathematical difficulties. Consequently, theoretical solutions have been obtained only for idealized cases in which important physical phenomena are neglected. For example, consider the process of forward in situ combustion. All such theories which have been developed consider only certain aspects of the process, such as heat transfer, 2-4,8,9,l6,l8 heat transfer with phase change,10 heat transfer with chemical reaction,7 or the hydrodynamics of three-phase flow.21 A general theory including all of the above phenomena has not been developed to date. This paper presents a unified theory of thermal recovery processes in linear systems. A mathematical model is developed which explicitly includes conduction-convection heat transfer with convective external heat loss, chemical reaction between air and oil, aqueous phase change, and the hydrodynamics of three-phase flow. A system of equations is developed which can be solved numerically on a high-speed digital computer, resulting in predicted temperature, pressure, and saturation histories in space and time. The model allows a more detailed simulation of thermal recovery tube experiments than had previously been possible. THEORETICAL DEVELOPMENT Consider the linear flow of gas, water and oil in a homogeneous porous medium. Assume that the oil will react with gaseous oxygen, and that mass is transferred between the water and gas phase by evaporation or condensation. In the absence of hydrocarbon phase change, intraphase diffusion,

Jan 1, 1966

By M. Semchyshen, A. P. Coldren, W. G. Scholz

A survey of the Fe-Mo-B system was made to determine the extent to which boron might affect the microstructure and strength properties of iron-rich Fe-Mo alloys. Seventeen vacuum-induc tion melted ingots were prepared with approximate mo-lybednum contents 01 2, 4, 9, and 19 wt pct and approximate boron contents of 5, 14,50, and 630 ppm by weight. Phases extracted from annealed specimens and identified by X-ray diffraction analysis included E Fe3MO2, Mo2FeB2, and an oxide hazling a diffraction pattern practically identical to that of MnFez04. Three types of alloy strengthening were observed: a) Fe-2 pct Mo-B alloys were hardenable by a bainitic-type transformtion of austenite to ferrile; b) Fe-4 pct Mo-B alloys were hardened by the solid-solution mechanism; and C) Fe-9 pct Mo-B and Fe-19 pct Mo-B alloys were hardened by precipitation of the E Fe3Moz phase. Boron was observed to have a very strong hardening effect in the Fe-2 pct Mo alloys and a mild strengthening effect on the 1200°F properties of the Fe-4 pct Mo alloys. In the Fe-9 pct Mo alloys there was no consistent effect of boron on strength at either room temperature or 1200° F. The Fe-19 pct Mo alloys were so brittle that meaningful tensile or creep-rupture data could not he obtained. No dispersion hardening from borides was recognized in any of the alloys. THE high-temperature properties of iron-rich Fe-Mo alloys were studied by Reiter and Hibbard.' They found that the strength of the alloys increased markedly as the molybdenum content was raised from 0 to 15.9 pct by weight. The authors explained the various degrees of strengthening in terms of a) solid-solution strengthening, b) strengthening from a martensitic-type transformation, and C) precipitation strengthening. The objective of the present investigation was to conduct a survey of the iron-rich portion of the Fe-Mo-B system to determine whether the beneficial effects of molybdenum on the strength of iron can be enhanced by the presence of boron. It was thought that perhaps fine dispersions of molybdenum- rich borides would be found which could improve the high-temperature strength of ferritic alloys beyond the strengths attainable by molybdenum additions alone. Also, it was considered worthwhile to study the effects of boron in alloys which are essentially free of carbon. Most, if not all, of our present knowledge about boron in high-temperature alloys pertains to alloys containing carbon and carbides. EXPERIMENTAL PROCEDURES Alloy Preparation. Seventeen 8-lb ingots were prepared from electrolytic Plastiron, unalloyed molybdenum chips, and boron powder by a vacuum-induction melting procedure employing hydrogen as the main deoxidizer. The alloys were melted in alumina crucibles using 35-lb melts, with each melt being split four ways to achieve varying boron contents within four Fe-Mo base compositions. The results of chemical analyses conducted on chips machined from the chilled ends of the ingots are presented in Table I. The boron analyses were performed by a spectrographic method utilizing primary standards prepared from boric acid. They indicated that in ingots to which no boron was

Jan 1, 1964

By Frederic L. Kadey

Diatomite is a siliceous, sedimentary rock consisting principally of the fossilized skeletal remains of the diatom, a unicellular aquatic plant related to the algae. Thus, it has been formed by the induration of diatomaceous ooze, and consists mainly of diatomaceous silica, a form or variety of opal which is first formed in the cell walls of the living diatom. Diatomaceous silica is not generally regarded as a synonym or the equivalent for diatomite, although it has been so used at various times. Accurately, diatomaceous silica is the preferred name for the principal mineral component of which the rock, diatomite, is composed. The terms diatomaceous earth and kieselguhr are used as synonymous with diatomite. The designations tripoli, tripolite, infusorial earth, etc., were used at one time but are now obsolete. With the changing nomenclature, these terms that were at one time correct when proposed and used for generations would be considered incorrect if used today in the light of current knowledge. The designation diatomite is reserved for those accumulations of diatomaceous silica that are of sufficient quality, size, and minability to be considered of potential commercial value. Processed diatomite possesses an unusual particulate structure and chemical stability that lends itself to applications not filled by any other form of silica. Foremost among these applications is its use as a filter aid, which accounts for over half of its current consumption. Its unique diatom structure, low bulk density, high absorptive capacity, high surface area, and relatively low abrasion are attributes responsible for its utility as a functional filler and as an extender in paint, paper, rubber, and in plastics; and as an anti-caking agent; thermal insulating material; catalyst carrier; and chromatographic support; polish, abrasive, and pesticide extender to name a few representative applications. The United States is the principal producing country, although diatomite is found in numerous other locations. Geology Composition and Morphology Diatomaceous silica qualifies as a. mineral of organic origin in much the same way that aragonite and collophane do. The silica of the fossilized diatom skeleton closely resembles opal or hydrous silica in composition (SiO2 • nH2O) .4, 12, 203 The silica is of acute biological significance, not only for the cell wall component, but also for the basic life process.45, 51 Without silica, cell development ceases .31 In addition to bound water, varying between 3.5 and 8%, the siliceous skeleton may also contain, in solid solution, or as part of the Si02 complex, small amounts of associated inorganic components-alumina, principally-and lesser amounts of iron, 46, 49, 55 alkaline earths, alkali metals, and other minor constituents.59, 63, 67 Boron is reported to be an essential element for diatom growth .52, 58 Since diatomaceous silica is not pure hydrous silica but contains other intimately associated elements, there is good reason to consider it a distinct type or variety.4 Associated with the diatomaceous silica, and integrated as part of the diatomite, may be variable amounts of organic matter, soluble salts, and particles of rock-forming minerals that were syngenetically deposited or precipitated with the diatom frustules. Sand, clay, carbonate, and volcanic ash are typical common contaminants. Other contaminating minerals may be present, such as feldspar, mica, amphiboles, pyroxenes, rutile, zircon-the result of weathering, then transporting, and subsequent redeposition of surrounding land masses. Commercial diatomite may also contain fragments and particles of

Jan 1, 1983

By John H. Schemel

Large Zircaloy-2 hammer or press forged bars did not exhibit the uniform excellent corrosion resistance to steam normally expected of the alloy in wrought form. Weight gains of coupons cut from forged bars were 40 to 60 mg per sq dm compared to the 28 mg per sq dm obtained on hot rolled sheet and strip. The mechanism of this corrosion and its relation to microstructure is discussed. After these observations, a solution heat treatment followed by rapid cooling was tried on coupons from eight heats of forged bars. The basket-weave structure resulting from the quench did not show agglomeration of intermetallic compounds. These coupons showed uniformly good corrosion resistance with weight gains very close to the 28 mg per sq dm expected of Zircaloy-2 after the 14-day steam test. A large forging was solution heat-treated and tested for structure, mechanical properties, and corrosion resistance. Corrosion resistance was excellent and uniform throughout the section. The normally anisotropic mechanical properties were changed to completely iso-tropic. Strength levels were raised while there was a small loss in ductility at the service temperature. Tempering of the quenched structure below the ß transus did not improve the ductility at the test temperature of 600°F. ZIRCALOY-2 end caps for nuclear fuel elements Zircaloy-2 is an alloy of 1.5 pct Sn, 0.1 pct Fe, 0.1 pct Cr, and 0.05 pct Ni, balance hafnium-free zirconium patented by Westing-house Electric Co. used in pressurized water reactors are machined from large hammer-finished forged bars. The bars forged for this application range from 3 1/2 in. to 7 1/2 in. squares. One of the tests used to evaluate Zircaloy-2 for this service is a 14-day exposure to very high-purity steam at 75oF and 1500 psi pressure. Normally wrought products will exhibit a black lustrous film of zirconium oxide after the test and show a weight gain of approximately 28 mg per sq dm. Coupons representing the forged bars exhibited a wide variety of corrosion results ranging from acceptable black coupons to some covered with white crystals of zirconium oxide that had weight gains of nearly 100 mg per sq dm. One of the most common effects was a white corrosion product that looked like a stain to some observers and the outline of massive metal grains to others. Very careful specimen preparation prior to the corrosion test did not affect the result and metallographic examination did not reveal a structural feature of a size and shape that would correlate with the corrosion product. In addition to the relatively poor corrosion resistance, the forged product was not as ductile as desired. The problem of obtaining uniform and more acceptable properties in these heavy-forged sections seemed to be a function of the microstructure of the metal and, in particular, the distribution of the intermetallic compounds. None of the four alloying elements however are appreciably soluble in a zirconium and are present as compounds which usually appear scattered randomly throughout the structure. Moudryl showed that the distribution of these intermetallic compounds in Zircaloy could be related to a ostringero corrosion failure. Grozier et al. discussed another structural defect that causes a similar corrosion effect. In this case, an elongated gas-void was determined to be the cause. M. L. Picklesimer held in his discussion of the Grozier paper that at least a part of the observed =stringersn were caused by the distribution of the intermetallic compounds. He proposed a

Jan 1, 1962

By L. R. Duncan, O. M. Wicken

Magnesium, the eighth most abundant element in the earth's crust, is found widely distributed in a variety of minerals. Among the more commercially important ones are magnesite (MgCO3), brucite (Mg[OH]2), dolomite (CaMg[C03]2), and the salts magnesium sulfate and chloride often found in natural brines. Among other commercially important magnesium-containing minerals are olivine ([MgFe]2Si04), talc (H2Mg3[Si03]4), and serpentine (H4Mg3Si2O9). These are valuable because of desirable characteristics given them by the magnesium content and its placement in the particular crystal structure. These minerals are the raw materials for a host of products including magnesium metal and several grades of magnesia or magnesia-containing materials for refractories, fluxes, fillers, insulation, cements, decolorants, fertilizers, and chemicals. Considering only those minerals valued for their magnesia or magnesium content, over 85% of the material tonnages are processed to be used as refractories. Dolomite, which is used in enormous tonnages for its physical properties in construction, has an important use as a refractory raw material in its own right and, after calcination, as a precipitant for magnesium hydroxide, which is an intermediate for the production of magnesia or magnesium metal. Magnesia from the reaction of calcined dolomite (or limestone) with seawater or magnesium-containing brines has supplanted much of the production of magnesia from mined magnesite, particularly in the United States. The best known of the minerals directly and widely exploited for magnesia content is magnesite, one of the calcite group of rhombohedral carbonates which includes calcite (CaC03), siderite (FeCO3), and rhodochrosite (MnCO3) among others. The members of this group enter into a wide range of substitutional solid solutions when the positive ions have similar radii. The radii of magnesium and iron ions are within 6% of each other; hence, magnesite and siderite form a complete series of which breunnerite (ferroan magnesite) is a well-known end member. However, the radius of calcium ion is 36% larger than that of magnesium ion, and only limited substitution exists at each end of the MgCO3-CaCO3 series. Dolomite is not a member of the calcite group, but results when calcium and magnesium ions alternate in equal number in an ordered structure among carbonate ions. The result of these relationships is that calcite and dolomite, often found intermixed with magnesite, occur commonly as identifiable crystal entities which can often be separated to a varying degree from the magnesite by various beneficiation techniques. Magnesite, when pure, contains 47.8% MgO and 52.2% CO.. The pure mineral is sometimes, but rarely, found as transparent crystals resembling calcite, but, preponderantly, magnesite contains variable amounts of the carbonates, oxides, and silicates of iron, calcium, manganese, and aluminum. Magnesite may be either crystalline or amorphous (cryptocrystalline). The crystalline form has a hardness of 3.5 to 4.0. The color may range from white to black with shades of yellow, blue, red, or gray. The color is not a fundamentally significant indicator of purity, but in a given deposit, an experienced person can often roughly grade magnesite by assessing color and crystallinity. The crystalline form of magnesite occurs in

Jan 1, 1983

By R. F. Domagala, D. J. McPherson

Iodide zirconium was combined with calculated amounts of ZrO2 or master alloys and arc-melted. Annealing treatments were carried out at 21 temperature levels. Metallographic examination of the heat treated specimens permitted construction of the binary phase diagram from zirconium to ZrO2. Oxygen additions to zirconium raise the transformation temperature as well as the melting point. Features of the diagram include the peritectic formation of beta, the formation of alpha directly from the melt, an intermediate phase ZrOB with a' range of homogeneity, and a eutectic between alpha and Zr02. PHASE relationships in the Zr-0 system were determined in the range 0 to 66. 7 atomic pct 0 (ZrO,). Arc-melted alloys were annealed at temperature levels between 600" and 2000°C. Determination of the phase boundaries was accomplished by metallographic evaluation of specimens quenched from the various temperatures. Incipient melting techniques were used to determine solidus curves, and X-ray diffraction work was employed to study the lattice parameters of the a solid solution and the phase ZrO,. Experimental Procedures Westinghouse "Grade 1" iodide zirconium crystal bar (approximately 99.8 pct pure) was employed for the investigation. This material is substantially free of the impurity inclusions characteristic of most zirconium metal. Previous work has shown it to be the most satisfactory grade for phase diagram studies. The zirconium crystal bar, as-received, was coated with corrosion product from autoclave tests by which its grade designation is determined. The bars were sand blasted lightly, pickled for 1 min in a 20 pct HN0,-5 pct HF aqueous solution, rinsed in water and acetone, and dried. The bars were rolled to about 1/32 in. strip, cut into 10 in. lengths, and pickled again. The material was sheared to approximately Y4 in. squares, cleaned with acetone, and stored for use. Two pounds of specially prepared "chemically pure," and 100 grams of spectrographic grade ZrO, were obtained for the preparation of alloys. The "Specpure" ZrO, was obtained from Johnson, Mat-they Co.. Ltd.. and Titanium Alloy Manufacturing Div. of National Lead Co. supplied the other grade of oxide. The "chemically pure" material (henceforth referred to as C.P. ZrO,) was used for the preparation of all alloys employed in the determination of phase relationships. A limited number of alloys were prepared with the "Specpure" oxide to recheck certain phase boundaries. Analyses of the two grades of oxide are given in Table I. The C.P. ZrOI was received in powder form, not suitable for arc melting. Therefore, this material was compressed to wafers 1 in. in diameter and about Vn in. thick. These were segmented and stored for use. The "Specpure" material, in the form of small granules, was suitable for use in the as-received condition. A nonconsumable electrode arc-melting furnace was used to prepare alloys. A water-cooled copper crucible and tungsten-tipped electrode were employed. The material was placed in the crucible and rapidly melted under a protective atmosphere of helium gas by striking the arc on a tungsten stud and directing it upon the charge material. Details of operation, including drawings of this type furnace, have been published.'

Jan 1, 1955

By R. W. Heckel

The growth and shrinkage rates of second-phase particles of various size distributions are analyzed in terms of five mathematical models. These models apply to various forms of diffusion-controlled growth (one model) and interjace-controlled growth (four models). The expression for diffusion-controlled growth assumes local equilibrium at the interface between the second Phase and the matrix and through the choice of the proper flux area perrnits the definition of minimum diffusion-controlled growth kinetics. The interface-controlled kinetics models consider the situations where the rate-controlling step occurs at either the interlace undergoing solution or the interface undergoing deposition. In addition, each of these interface-controlled kinetics models is analyzed for interface reaction rates which are independent of the thermodynamic activity gradient across the interface or proportional to the gradient. THE variation of the size distribution of second-phase particles in a solute-saturated matrix has been a stimulating topic for many years. Research has been carried out in this area for situations where the matrix was either liquid or solid. It is well-understood that the over-all process results in an increase in the average particle size, a decrease in the total number of particles per unit volume, and a net decrease in surface area. Greenwood1 has presented a mathematical model based upon the mass balance existing between growing and shrinking particles. His analysis is stated in terms of Fick's First Law using concentration gradients set up between particles of different size as a result of the effect of surface curvature on solubility. Greenwood's analysis indicates that the mean particle size should increase approximately as the cube root of time. Asimow2 has shown that several of the simplifying assumptions employed by Greenwood may be eliminated from this diffusion-controlled model with little increase in complexity. Lifshits and Slyozov3-5 have also considered the growth-shrinkage problem involving diffusion-controlled kinetics. Their analysis was extended to include the effects of interacting particles ("kinetic approximation"), particle shape, and elastic strains. These authors also show that the mean particle size varies as the cube root of the reaction time. wagner6 has formulated expressions for both diffusion-controlled kinetics and interface-controlled kinetics in the growth-shrinkage problem. Wagner's diffusion-controlled kinetics is based upon the same form of mass balance using Fick's First Law as used by others and is shown to be applicable to situations where kr >> D (k is the interface reaction constant, r is the particle radius, and D is the diffusion coefficient). This analysis also yields a mean particle size which varies as the cube root of time, in agreement with other workers. Wagner also developed interface-controlled kinetics to handle kr « D. This model was based upon the fact that the interface reaction rate is proportional to the interface area. Wagner assumed that the rate was also proportional to the thermodynamic activity difference, ?a, existing across the interface. Thus, the specific reaction rate at the surface of radius r was defined as k • 4nr2 . ?a. It is the purpose of the present paper to reconsider the processes of diffusion-controlled growth and interface-controlled growth in the variation of particle-size distributions as a function of time at constant second-phase volume. The analyses are formulated in a manner which permits data obtained

Jan 1, 1965

By Gareth Thomas

The analysis of Kikuchi pattersns of exct ovientalions from single cryslals and paired Kikuchi lines from single and overlapping crystals is shown to be useful and quanlitalve and is applied to Phase transfovnzcitions including ordering, spinodals, mavten-silic, and nuclealion and growth pyocesses. 112 pinciple, the analysis of exacl orientations enables the crystal system and the Bravais lattice of a crystal to he determined. The advantages of the davk-field imaging technigure for detecting vevy small precipitates are also described. ALTHOUGH Kikuchi electron-diffraction patterns were first observed nearly 40 years ago,' little systematic application seems to have been made of them, until fairly recently during electron microscopy, when their usefulness in contrast experimen and for determining exact orientations8'9 has been pointed out. With the availability of gonio-metric specimen tilting stages it has now become possible to make much wider use of diffraction patterns, particularly the Kikuchi pattern, which is the main subject of this paper. The treatment presented here is not exhaustive but it is hoped that it will stimulate more use of Kikuchi patterns during electron-microscopy investigations. The Kikuchi pattern is formed as a result of Bragg diffraction of the inelastically scattered electrons produced during the interaction of the beam and thick specimens, The important feature of these patterns is that they give an accurate representation of the symmetry of the crystal being investigated, so that it is possible to identify crystal systems and even the Bravais lattice. This means that new structures, e.g., formed in phase transformations, may be identified during normal electron microscopy, so that Kikuchi-pattern analyses considerably extend the uses of the electron microscope. Recent work has also shown that dark-field images are more informative than bright-field images, particularly, for observing small precipitates. The second part of this paper discusses some applications and advantages of dark-field imaging in studies of two-phase systems. 1) DIFFRACTION PATTERNS 1.1) Spot Patterns. Electron-diffraction spot patterns have their limitations because of the importance of the form factor on the intensities and shape of the reciprocal lattice points. Because of the extension of these points into relrods, reflections are possible over a large angular range (+5 deg) and the patterns from thin regions can be complicated because second-layer relrods intersect the reflecting sphere. Spot patterns can thus give only an approximate idea of the crystallography unless the foil is tilted into exact orientation. In this case the spot pattern is symmetrical, with equal numbers of spots on the positive and negative zone directions about the origin. Such cases are necessary for structural analyses and have been used recently to determine the crystal structure of the ordered Ta64C phase." Exact orientation means that the plane of the reciprocal lattice lies exactly normal to the incident beam as shown in Fig. 1(b). It should be noted that in exact orientations the angle of reflection is less than the exact Bragg angle 9 so that the reciprocal lattice points lie to the outside of the sphere, Fig. 1(b). This deviation is denoted by the parameter s and in the usual convention4 s is negative for exact orientations and, of course, zero at the exact reflecting position shown in Fig. l(a). In order to avoid secondary reflections from thickness relrods it is advantageous to work in thicker regions of foils. Diffraction from thick regions may also produce Kikuchi patterns and as dis-

Jan 1, 1965

By Langtry E. Lynd

The predominant tonnage use of titanium is as a white pigment, in the form of rutile or anatase, which are allotropic modifications of titanium dioxide. Because of its whiteness, high refractive index and resulting light scattering ability, titanium dioxide is unequaled as a white pigment for paints, paper, rubber, and a wide variety of other materials. Consumption of titanium dioxide pigments for these uses has been rising steadily for many years, United States consumption in 1956 amounting to about 486,000 tons, worth $250,000,000. Production was off a little in 1957 and 1958. An outstanding development in the titanium field in recent years is the expansion of titanium metal production. Virtually nonexistent before 1948, United States production of pure titanium metal sponge reached a high in 1957 of about 17,200 tons, worth about $77,000,000. United States titanium sponge metal production in 1958 was cut back to about 4,600 tons because of reductions in military requirements. Another major use of titanium, in the form of the mineral, rutile, is for welding rod coatings. United States consumption of rutile for this purpose in 1958 amounted to about 11,000 tons, worth $1,100,000. Titanium-Bearing Minerals Although several minerals contain titanium as a major constituent, only a few of these occur in sufficient quantities and concentrations to be of commercial importance. These are ilmenite (FeTiO3), rutile (TiO2), and alteration products of ilmenite which have been referred to as weathered ilmenite, leucoxene, or arizonite. The other polymorphs of titanium dioxide, brookite and anatase, occur in a few deposits, of which only one has been of even minor importance. Considerable concentrations of titanite or sphene (CaTiSiO5) and perovskite (CaTiO3) are known, but do not compete with ilmenite and rutile deposits in grade or size. Rutile is very widely distributed as an accessory mineral in most rocks, but large deposits of commercial importance are rare. A few rock deposits have been worked for rutile but only one, in Virginia, is in operation (1959). Development of one in Mexico was discontinued recently. At present, the main sources of rutile are the beach and dune sands of Australia and Florida. A major use of rutile is for the manufacture of titanium tetrachloride, which is the intermediate compound universally used in titanium metal manufacture. Another main use of rutile is for welding rod coatings. Although the whiteness and brightness of titanium dioxide pigment is due to the high refractive index of rutile or anatase, natural rutile is too impure to be used directly as a white pigment, is too expensive to use for manufacturing pure oxide pigment, and is too scarce to provide the required tonnages. Instead, the cheaper, more abundant, mixed oxide of iron and titanium, ilmenite, which is readily decomposed by sulfuric acid, is used for manufacture of rutile and anatase pigments. Ilmenite is mined from both rock and sand deposits. Rock deposits contain ilmenite in massive or heavily disseminated form along with magnetite or hematite and are almost invariably associated with anorthositic and gabbroic rocks. Ilmenite also occurs as an accessory mineral in a great many igneous and metamorphic rock types. Concentrations of

Jan 1, 1960

By J. F. Breedis

The morphology and crystallography of marten -site formed during quenclzing were examined by transmission electron microscopy in alloys whose compositions lie between Fe-19 wt pct Cr-11 wt pet Ni and Fe-33 wt pet Ni. Depending upon composition, four types of transformation products were observed in this portion of the Fe-Cr-Ni system. Three of these products are denoted as lath, plate, and surface martensite. The hcp structure found in association with martensite it? certain higlz-clzromium stairzless steels does not influence the transformation crystallogyaphy and cannot he Considered to act as a transitional structure in the transformation from austenite to martensite. DISTINCT changes in the characteristics of mar-tensitic transformations are produced by the partial replacement of nickel by chromium in austenitic iron alloys. Large, discrete martensite plates form in Fe-30 pct Ni, while aggregates of small crystals form in sheets parallel with (111)A austenite planes in stainless steels containing in the vicinity of 18 pct Cr and 8 pct Ni. These transformation structures are termed plate and lath martensite, respectively. In addition, an hcp structure (E) has been observed in certain high-chromium stainless steels in association with martensite. It has been suggested that the E structure acts as an intermediate step in the nucleation of martensitel-' and influences the morphology and crystallography of the final transformation product. However, this structure may be the consequence of the deformation of austenite arising from the shape distortion accompanying the transformation to bcc martensite.' An hcp stacking of atoms is produced in an fcc lattice by the introduction of an intrinsic stacking fault on alternate (111)A planes. At present, the arguments behind each explanation for the E structure are based on circumstantial evidence. This investigation seeks to determine the importance of the hcp structure to the martensitic transformation, and to examine the crystallography and morphology of the structures produced in bulk samples by quenching in alloys whose compositions lie between Fe-19 wt pct Cr-11 wt pct Ni and Fe-33 wt pct Ni. These aspects of the martensitic transformations occurring in these alloys have been studied by electron metallography and selected-area diffraction of foils examined in transmission in the electron microscope. EXPERIMENTAL METHODS The alloys have been prepared from thoroughly mixed powders of pure iron (99.6+ pct, 0.06 max pct C, 0.15 max pct N), chromium (99.9+ pct), and nickel (99.9+ pct). After sintering in an Ar-H2 atmosphere, the mixtures were melted by induction heating in recrystallized alumina boats under purified argon. The ingots were cold-rolled to 0.4-mm sheets (100 pct reduction in thickness) after the removal of at least 1.5 mm from their surfaces by grinding. Three-millimeter discs were punched from the cold-rolled sheet and annealed at 1150°C for approximately 24 hr in evacuated quartz capsules. The average grain size was 40 µ Thin foils for examination in transmission in the electron microscope were prepared from discs quenched to liquid-nitrogen temperature. A depression was first jet-machined into the center of a disc and a thin area (less than approximately 3000A in thickness) was produced at the base of the depression by electropolishing. Electron-microscope studies utilized the RCA EMU-3 and the Philips EM-100B microscopes operating at 100 kv. X-ray examination of the polycrystalline alloys was performed either in a diffractometer or in a

Jan 1, 1964

By F. D. Wright

Production of a coal preparation plant is governed by many restrictions, such as the tonnage of different products and blends that can be sold within a given period, capacities and output proportions of the cleaning and sizing units, blending proportions, quality specifications, and costs and prices for the various products. Determination of the tonnage of each product and blend that should be made in order to obtain the maximum profit can be difficult unless a systematic method such as linear programming is used. In this paper the basic method of linear programming is described briefly. The Old Ben No. 9 preparation plant is used as an example to illustrate in detail how eqwations can be written to form a linear programming model of a coal preparation plant. Three sample problems, each requiring 56 or more equations and 63 structural variables, were solved with an IBM 650 computer. Linear programming is one of several mathematical tools used for operations research. It has been applied to many fields and has been used by a number of industries either to maximize the profits of certain operations or to minimize costs. P. B. Nalle and L. W. weeks' have described the use of the method by the Riverside Cement Co. to minimize the cost of blending raw materials to make portland cement. Their problem is to obtain at minimum cost a mix with certain specifications from a number of possible materials which have various costs and various amounts of CaO, SiO2, Fe2O3, and other constituents. In a paper on the use of linear programming by the National Coal Board in England, K. B. Williams and K. B. Halley2 describe how the transportation method, a variation of linear programming, is used to minimize the cost of sending 37 grades of coal from 28 mines to seven central washing plants which produce coal for furnace coke and foundry coke. The various coals have different percentages of volatile matter, moisture, sulfur, ash, and phosphorous, so there is considerable choice in how they can be blended to meet the specifications of the two products. The purpose of this paper is to show how linear programming can be used to maximize the profits of a coal washing plant which produces individual final products as well as blends. The Old Ben Corp. furnished assumed sample data from their Old Ben Mine No. 9 preparation plant for this investigation. However, the data that have been used are entirely the author's responsibility. METHOD OF LINEAR PROGRAMMING Numerous articles and books have been written on the theory and applications of linear programming.3-5 However, since the method has not been widely applied by the mining industry, a brief, nontheo-retical discussion of its basic method seems to be in order. Linear programming (Table I) is used to determine the best possible solution to a number of interdependent activities. It is essentially a method for making systematic selections from a number of possible solutions to determine positively the optimum solution. There may be other equally good solutions but no better ones. Each activity must have linear coefficients and there must be a criterion to judge how good the solution actually is. Linear programming is different from the solution

Jan 1, 1961

By Mark J. Klein

An internal-friction profile induced by nitrogen in Cr-35 at. pet Re was studied as a function of nitrogen concentration and heat treatment. From these studies, the solubility of nitrogen in this alloy was deduced. The equilibrium solubility of nitrogen in chromium does not appear to be substantially changed by the addition of 35 at. pet Re to chromium. However the precipitation kinetics of nitrogen from solution are significantly decreased. IT has been found that rhenium has a beneficial effect on the low-temperature ductility of Group VI-A metals. As part of a study to investigate the mechanism of this increased ductility, it was of interest to compare the concentration of interstitials retained in solution in one of these ductile bee alloys with the concentration retained in solution in its base metal. In this regard, a study of the solubility of nitrogen in ductile Cr-35 at. pet Re was selected for this phase of the investigation because the solubility of nitrogen in the base metal, chromium, has been determined over a fairly large temperature span.1-3 In addition, nitrogen is known to have an embrittling effect on chromium,4 but does not appear to have a comparable embrittling effect on Cr-35 at. pet Re.5 The experimental method selected for this study was the measurement of internal friction. This technique allows interstitial-concentration measurements to be made in bee metals that are sensitive to small changes in the concentration of a particular interstitial, even when other interstitials are present. When subjected to a cyclic stress, interstitial atoms in solid solution give rise to an internal-friction peak whose origin is the stress-induced ordering of interstitials into preferred lat- tice sites. These internal-friction peaks signify a dissipation of elastic energy which is greatest when the period of oscillation is comparable with the relaxation time for the attainment of an equilibrium interstitial distribution. For dilute interstitial binary alloys, the height of the internal-friction peak is proportional to the concentration of interstitial atoms in solid solutions. snoek8 has shown that the internal friction, Q-1, may be expressed by the equation where A is the relaxation strength, ? is the angular frequency, and 7 is the relaxation time. The relaxation strength, A, is proportional to the concentration of interstitials taking part in the relaxation if no interaction occurs between solute atoms. Thus values of peak height, which can be determined as a function of temperature, give a measure of the relative concentration of interstitial atoms in solution at various temperatures. If equilibrium has been attained, the peak heights will be proportional to the interstitial solubilities. The constant that relates A to the interstitial concentration can be determined experimentally from peak-height measurements as a function of concentration for specimens of constant crystallographic orientation. When this has been done, Q-1 measurements can be used to determine solubilities according to the equation C - KQ-1 where C is the soluble interstitial concentration and K is a constant that can be experimentally determined. It is possible, however, that lattice imperfections may bind interstitials to certain lattice sites. Internal friction does not measure these interstitial~, but does measure interstitials which are free to move in the normal lattice sites when subjected to an applied oscillating stress. On alloying the pure metal with a substitutional

Jan 1, 1965

By C. D. Stahl, S. M. Farouq Ali

This investigation attempts to describe and simulate the alcohol displacement process by means of a cell model, as employed in chemical engineering practice. The proposed model is more simple than previously proposed models. and utilizes parameters chosen on a theoretical basis. The model successfully reproduced the formation of the stabilized bank and the breakthrough of alcohol, the latter depending on one of the model parameters, which may be correlated with the length of the porous medium. Moreover, the effects of the phase behavior of the liquid system involved, as observed in experimental studies, were reproduced. Several variations of the basic model were devised and tested on a digital computer. These included the cases in which: (1) the actual value of fractional flow was used in cell-to-cell computations; (2) the number of cells was varied within the same run; and (3) incomplete rather than complete phase equilibrium was assumed within each cell. The proposed cell model clarifies the basic mechanism of the process. Detailed concentration profiles obtained for each cell, for instance, showed the mechanism of bank formation in relation to the phase behavior characteristics. The results obtained indicated a varying degree of phase equilibrium concommitant with changes in the velocities of the phases in an actual alcohol displacement. This condition was approximated by changing the number of cells during the simulation. Interesting information was obtained on the influence of path length on the efficiency of alcohol displacement, which has been the subject of some controversy. Certain limitations preclude the use of the proposed model as a substitute for experimental studies. The results obtained were, nevertheless, of value in interpreting the experimental results. INTRODUCTION During recent years considerable effort has been directed toward an understanding of alcohol displacement, the process whereby oil and water are recovered from a porous medium by the continuous injection of a solvent. The complex nature of the physical process involved has so far defied a complete mathematical treatment. Other methods of approach, amounting to an overall material balance, have been proposed, yielding useful information on certain aspects of the process.1-3 Taber et al, in particular, defined the displacement mechanism in terms of the phase behavior of the alcohol-oil-brine system invoIved.2 Wachmann reported a mathematical treatment of alcohol displacement subject to certain simplifying assumptions? Donohue proposed the use of a "cell model" for simulating alcohol disPlacement.5 The nature of the assumptions involved limited the utility of the model. The present work attempts to examine the variables involved in the simulation of alcohol displacement and discusses several possible versions of the basic cell model. Under certain conditions the model results are similar to the experimental results. In particular, the spontaneous formation of of the stabilized bank and the effects of the system phase behavior were successfully reproduced. PREVIOUS WORK ON CELL MODELS Cell models and the theoretical plate concept are often used in solving chemical engineering problems in which an explicit mathematical solution may be difficult or impossible to obtain. Examples of such applications occur in distillation, gas-liquid chroma tography,6 reactor technology, absorption, etc. In petroleum engineering, such a model was used by Attra7 to describe non-equilibrium gas drive, and by Higgins and Leigh ton8 to calculate sweep efficiency in water flooding. Aris and Amundsen pointed out the equivalence between the diffusion model and perfectly mixed cells connected in series.9 Deans proposed a three-parameter cell model to simulate two-component

Jan 1, 1966

By G. R. Speich, Morris Cohen

The growth rate of bainite has been determined by hot-stage metallography in five hypereutectoid high-purity iron-carbon, iron-carbon-chromium, and iron-carbon-nickel alloys. The studies have been confined essentially to lower bainite. The growth rate is independent of time at constant temperature, but both the edgewise and sidewise growth rates increase with increasing temperature. Increasing the carbon content, or adding chromium or nickel, decreases the growth rate without greatly affecting the temperature dependence. The results lead to separate diffusion-controlled models for the edgewise and sidewise growth rates. ALTHOUGH bainite has been subject of many detailed investigations,1,2 the growth mechanism of this intermediate transformation product of steel remains unclear. One difficulty has been that practically all the kinetic measurements have been restricted to studies of the overall reaction rate. Only a limited amount of information3, 4-6 is available on the growth process itself. The purpose of the present work was to investigate the effect of time, temperature and composition on the growth rate of bainite using hot-stage metallography. This method makes it possible to observe continuously the isothermal growth of individual bainitic units because of the surface upheavals accompanying their formation.3 This information is needed to test the existing growth models and to obtain a better understanding of the factors controlling the bainitic transformation. EXPERIMENTAL PROCEDURE Alloys—The compositions of the alloys used in the present work are shown in Table I. The five alloys were vacuum-melted and cast, and were prepared from high-purity materials. The three iron-carbon alloys were available from the work of Roberts et al.7 in the form of 3/8-in. diam rods.* To facilitate machining, the rods were ing. The rods were then tempered for 3 hr at 650°C to allow machining. The iron-carbon alloys in Table I were selected for investigating the effect of carbon on the bainitic growth rate. Due to the rapidity of the transformations involved, it was not feasible to measure the growth rates in iron-carbon alloys having carbon contents much less than 0.9 pct C. The iron-carbon-chromium and iron-carbon-nickel alloys were selected to study the influence of chromium and nickel on the growth rate at a carbon level of about 1.1 pct C, since this carbon content was intermediate among the iron-carbon alloys. The concentrations of chromium and nickel were chosen to match isothermal trans formation diagrams in the literature.8'9 Hot-Stage—Since the iron-carbon alloys had such low hardenabilities, it was necessary to quench the hot-stage specimens drastically from the austentiz-ing temperature to the isothermal transformation temperature to prevent the formation of pearlite or proeutectoid constituents. A rapid-quenching metal-lographic hot-stage was built for this purpose. A maximum temperature of about 1100°C and quenching rates up to 500°C per sec could be attained with this equipment. The details of the apparatus are described in a recent publication. After being polished and mounted on the hot-stage, the specimen (0.020 in. thick) was protected from oxidation and decarburization during austenitizing by a vacuum of 10-4 to 10-5 mm. The quench to the reaction temperature was accomplished with a helium blast, and the isothermal run at this temperature was conducted in a helium atmosphere. A fine-wire thermocouple spot-welded to the sample and a high-speed

Jan 1, 1961

By C. K. Eilerts

Finite difference equations were programmed and used to integrate the second-order, second-degree, partial differential equation with variable coefficients that represents the transient linear flow of gas-condensate fluids. Effect was given to the change with pressure of the compressibility factor, the viscosity, and the effective permeability and to change of the absolute permeability with distance. Integrations used as illustrations include recovery of fluid from a reservoir at a constant production rate followed by recovery at a declining rate calculated to maintain a constant pressure at the producing boundary. The time required to attain such a limiting pressure and the fraction of the reserve recovered in that time vary markedly with properties of the fluid represented by the coefficients. Fluid also is returned to the reservoir at a constant rate, until initial formation pressure is attained at the input boundary, and then at a calculated rate that will maintain but not exceed the limiting pressure. The computing programs were used to calculate the results that would be obtained in a series of back-pressure tests made at selected intervals of reservoir depletion. If effect is given to the variations in properties of the fluid that actually occur, then a series of back- pressure curves — one for each stage of reserve depletion — is required to indicate open-flow capacity and related flow characteristics dependably. The isochronal performnce method for determining flow characteristics of a well was simulated by computation. INTRODUCTION The back-pressure test procedurel is based on a derivation of the equation for steady-state radial flow of a gas, the properties of which are of necessity assumed to remain unchanged in applying the test results. The properties of most natural gases being recovered from reservoirs change as the reserve is depleted and pressures decline, and the results of an early back-pressure test may not be dependable for estimating the future delivery capacity of a well. The compressibility factor of a fluid under an initial pressure of 10,000 psia can change 45 per cent and the viscosity can change 70 per cent during the productive life of the reservoir. There are indications that the effective permeability to the flowing fluid can become 50 per cent of the original absolute permeability before enough liquid collects in the structure about a well as pressure declines to permit flow of liquid into the well. The advent of programmed electronic computing made practicable the integration of nonlinear, second-order, partial differentia1 equations2.3 pertaining to flow in reservoirs. Aronofsky and Porter4 represented the compressibility factor and the viscosity by a linear relationship, and integrated the equation for radial flow of gas for pressures up to 1,200 psi. Bruce, Peaceman, Rachford and Rice5 integrated the partial differential equations for both linear and radial unsteady-state flow of ideal gas in porous media. Their published results were a substantial guide in this study of integration of the partial differential equation of linear flow with coefficients of the equation variable. The computing program was developed to treat effective permeability as being both distance-dependent and pressure-dependent. In this study all examples of effective permeability are assumptions designed primarily to aid in developing programs for giving effect to this and other variable coefficients. The accumulation of data for expressing the pressure dependency of the effective permeability is the objective of a concurrent investigation. FUNDAMENTAL DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATION OF FLOW Linear transient fluid flow is represented by the partial differential equation,

Jan 1, 1965

By Lincoln F. Elkins

The Spraberry Trend Field of West Texas was discovered in January, 1949. Drilling of 2.234 wells and production of some 45 million bbl of oil by January, 1953, indicated this to be an important field which will ultimately cover more than 400,000 acres. In addition to being the world's largest field in areal extent, the Spraberry has presented many problems in well completion and operation and has demonstrated unique reservoir performance characteristics. The pay section consists primarily of a few fine grained sandstone or siltstone members in a thousand-ft thick section of shale, limestone. and siltstone. Since porosity averages only 10 per cent and nearly all permeabilities are less than 1 md, conventional core analysis does not delineate the "pay" section. Mercury injection was used as a capillary pressure test adaptable to rapid routine use to select those intervals having low enough connate water saturation to contain commercially significant oil saturation. In the central area of the field this "pay" amounts to 16 ft of Upper Spraberry and 15 ft of Lower Spraberry sands. An interconnected system of vertical fractures, observed in cores, provides the flow channels for oil to drain into the wells hut most of the oil is stored in the matrix since the void volume of fractures is estimated to be less than 1 per cent of that in the sand. Initial potentials of wells range up to 1.000 B/D after fracture treatment which should be compared with estimated capacity of 5 to 10 B/D if oil had to flow into the wells through the sand itself. Without exception initial pressures of later drilled wells were significantly lower than initial pressures of earlier drilled nearby wells in a large area some 6 miles long. This means the earlier drilled wells had drained fluids from areas much greater than their 40-acre proration units. Since most of this performance occurred while the reservoir pressure was above the. saturation pressure it was analyzed by the compressible fluid flow theory. This analysis gave calculated initial presrilres which agreed within 30 psi of measured pressures of 60 per cent of wells in the area using 16-md permeabilitv corresponding to a fracture system substantially that indicated by cores and using combined compressibility of rock and its contained oil and water corresponding to the core analysis data. The most important feature of this analysis was the very close agreement between effective compressibility of the rock and its contained oil and water from the field performance and that from the core tests, because it meant there are no "islands" of low permeability reservoir rock left untapped in the inter-well area and thus no additional wells are necessary to insure that at least one well penetrates each "reservoir." Twenty-five of forty-four 40-acre spaced wells on three contiguous sections were wed in a four-month inteference test. Six shut-in well. were tested monthly for oil production. productivity index. gas-oil ratio and pressure buildup, and seven shut-in wells were tested for decline in reservoir pressure. Tests on 12 regularly producing wells gave comparative data for interpretation of shut-in test wells. Reduction in reservoir pressure, decline in productivity index. and increase in gas-oil ratio were found to be substantially the same in the shut-in test wells as those in the comparative regularly producing wells, meaning that the producing wells were depleting the

Jan 1, 1953

By C. S. Smith, W. M. Williams

THE quantitative study of grain shape in three dimensions has been a difficult one from the practical standpoint. Experiments on grain shape have usually been based on indirect observations of two-dimensional metallographic samples or on grains separated by intergranular corrosion or fracture. The latter method was used in particular by Deschl in 1919 in observing the shapes of ß-brass grains which had been separated with mercury. The similarity in detail of shape and packing that exists between bubbles in a soap froth and metal grains was noted by Desch who remarks that: "Surface tension has an. important share in determining the form of the crystal grains in a solidifying metal and such grains have a tendency to assume the shape of foam cells." More recently the relation between grain shape and grain growth has been studied by Harker and Parker.' The importance of surface tension forces in determining microstructure has been discussed by Smith," and an analysis of the topologi-cal rules that apply to such structures has been made by the same author.' Biologists have made valuable contributions to the general problem of shape and form in an assembly of space-filling cells subject to surface tension forces. Particular mention must be made of the extensive researches of Lewis5,8 on botanical and anatomical tissues, and Matzke and Nestler,7,8 who have made painstaking studies of bubble shapes in different kinds of foam. The present paper introduces a new method of studying individual grains in particular alloys—a method which may well find use in other metallurgical researches, and which seems also to be applicable in biological problems. Metal Grains in Three Dimensions The individual grains in a piece of metal or the bubbles in a soap froth must satisfy two conditions —they must be in juxtaposition so as to fill space and their interfaces must conform to the laws of surface tension. If a further restriction is imposed, namely, that only one type of polyhedron be used, then the sole possibility is an assembly of "minimal" tetrakaidecahedra," each polyhedron having eight doubly-curved hexagonal faces and six plane quadrilateral faces, the angles between adjacent faces being 120". Such cells could be stacked on a body-centered cubic lattice and would be quite stable. For many years it was believed that soap bubbles actually possessed the shape of this ideal body, and the hypothesis was the starting point for several investigations on the shapes of soap bubbles, biological cells, and metal grains. These investigations have shown, what might have been obvious in the beginning, that if the condition of symmetry is relaxed the tetrakaidecahedron is almost never found. If we retain only the rule that all contacts should be in local surface tension equilibrium and that the assembly should fill space, metal grains may then take a nearly infinite variety of shapes and no single one may be regarded as an archetype. These two conditions, which are fully discussed in one of the previously mentioned papers,' lead us to think of a typical grain in an annealed metal as having the following characteristics: 1—A grain may have any number of faces, but usually between nine and eighteen. The faces (each of which is, of course, shared by two grains in the assembly) may have any number of edges, but usually between four and six. Pentagons are by far the most frequent. 2—In the assembly of grains each edge is shared by three faces and by three grains. Each vertex is a point shared by four grains, six faces, and four

Jan 1, 1953