Search Documents

By R. E. Carter, C. Wagner, F. D. Richardson

By means of inert markers of radio-platinum, it has been shown that cobalt metal oxidizes by out-ward diffusion of cobalt atoms through the oxide. Oxidation rates have been measured at various temperatures and oxygen pressures and have been found to agree with the rates calculated from the Wagner equation and the authors' values for the diffusion coefficient of cobalt in the oxide. The distribution of radio-cobalt in growing oxide loyers has been accurately measured and found to be different from that predicted from the Wagner oxidation theory. Attempts have been made to measure the change of lattice parameter of the oxide with composition. SEVERAL investigators have shown that cobalt metal oxidizes according to the parabolic law,1-3 but no detailed study has been made of the oxidation mechanism. Cobaltous oxide has a metal-deficient lattice, but the defect concentration which is dependent on the oxygen pressure is only approximately known.' From the similarities in the structures of "CoO" and "FeO" and the rate laws governing the oxidation of cobalt and iron, Gulbransen- as suggested that cobalt, like iron," oxidizes by diffusion of cobaltous ions through the oxide layer. To substantiate this hypothesis, he has attempted to measure the oxidation rate of cobalt as a function of the oxygen pressure but has been able to report only that the rate depends on a "small" power of the pressure. On the other hand, Preece," Valensi,' and Arkharov' all have suggested from indirect evidence that oxygen diffuses through the COO to the metal-oxide interface. The present authors, in a previous paper,' have shown that the rate of diffusion of cobalt in cobaltous oxide is dependent on the oxygen pressure. In that investigation, they also found that cobaltous oxide disks made by complete oxidation of cobalt metal were porous at their centers (Fig. 5L, ref. 4). Both these facts strongly suggest that oxidation of the metal proceeds via outward diffusion of cations. In order to establish the mechanism unequivocally, oxidation experiments were made with inert radioactive platinum markers, and the dependence of oxidation rate on oxygen pressure was measured. In addition, an attempt was made to determine the form of the cation vacancy gradient, which the marker experiment showed must exist across a growing oxide layer. This was done by depositing radioactive Co on the surface of a cobalt disk and then measuring the tracer distribution across the oxide layer produced by partial oxidation of the metal. Finally, an X-ray examination was made in an attempt to find the expected change of lattice parameter of COO with oxygen pressure. Experimental The apparatus and materials used were substantially the same as those previously described by the authors.4 Marker Studies: The use of inert markers to follow oxidation reactions was first attempted by Pfeil" in 1929, and radioactive markers have recently been used by Davies et al.' In the investigation reported here, the markers were applied as an (NH4)2PtCI6 slurry in several thin streaks across both faces of four disks of cobalt metal 1 mm thick and 14 mm in diameter. The disks were then oxidized almost to completion at 1200°C in pure oxygen. After oxidation, they were mounted and sectioned and auto-radiographs were taken to locate the markers. Oxidation Rate: The oxidation rate of cobalt was measured at 1148°C at oxygen pressures ranging from 0.0055 to 1.0 atm and at 1000°, 1148", and 1350°C under 1 atm of oxygen. The rate was deter-

Jan 1, 1956

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 H. R. Brannon, R. J. Robinson, G. W. Kattawar

Accurate prediction of storm wave characteristics are needed for design of offshore structures. Statistical methods of random noise analysis provide techniques for predicting required wave properties. These techniques have been used to analyze and characterize storm wave profiles from Gulf of Mexico recording installations in 34, 65 and .98 ft of water. Correlations based on the results can be used to predict wave crest probabilities and wave shapes for a range of Gulf water depths and storm conditions. Example predictions of wave crest probability as a function of water depth for a particular set of storm conditions are given. INTRODUCTION Accurate predictions of wave crests and wave shapes are needed for design of offshore structures. With predictions of wave crests, platform deck elevations giving adequate protection from exposure to storm waves can be selected. Use of hydro-dynamic wave theories with the wave crest and wave shape predictions permits an estimate of the wave forces an offshore structure must withstand. Observation of ocean waves suggests use of statistical analysis for studying wave characteristics for water depths beyond the breaking wave zone; recent wave generation theories lend support to this approach. l-3 One of the most powerful methods of statistically analyzing wave profiles was introduced by Cartwright and Longuet-Higgins; most of the basic relations used were developed by Rice5.6 in studies of random noise theory. This method predicts wave crest elevation probabilities from two parameters characterizing a wave profile. The statistical approach also provides a means for predicting wave shapes. In the random noise model of a wave profile, surface elevation is represented as an infinite sum of sine waves with closely spaced frequencies and random phase angles. Power density is proportional to the sum of the squares of the amplitudes of sine waves having frequencies within a narrow band, and the distribution of power density as a function of frequency is called the power density spectrum, or power spectrum. The power spectrum characterizes an irregular sea, and hence finds use in motion studies of ships, barges and semisubmersible drilling platforms. From the power spectrum, a wave profile of any duration of time can, in principle, be calculated; this long wave profile depicts the many wave shapes to which a structure may be exposed. Thus, the statistical methods of wave analysis provide an approach to selecting wave shapes as well as wave crest elevations needed for design of offshore structures. For practical application of these techniques, power spectra and related quantities must be predicted from storm properties. Barber and Tucker 7 have reviewed correlations of wave properties and storm conditions. Their review summarizes work of Darbyshire,8 among others, in which a correlation of power spectra with wind intensity of the North Atlantic was developed. These results, plus more recent work by Pierson and Moskowitz and Kitaigorodskii,10 establish feasibility of correlation, but there is no theoretical basis for modifications extending the relations to other areas as remote as the Gulf of Mexico. Hence, to predict characteristics of Gulf of Mexico waves, a study of local observations is needed. This paper presents (1) a practical approach to computing wave profiles hat depict shapes of waves for use in force calculations, (2) a summary of relations for predicting probabilities of wave crest elevations, (3) correlations of parameters needed to apply these methods in the Gulf of Mexico and (4) examples of application of the techniques. THEORETICAL BACKGROUND The following sections summarize relations needed to calculate wave profiles and to estimate wave crest probabilities. These relations also identify quantities that must be related to storm

By J. Gurland

The effects of variations of carbon content on the constitution and properties of sintered WC-Co alloys were studied. The grain growth of tungsten carbide during sintering was measured and it was shown that it takes place mainly by solution and reprecipitation of the carbide in and from the binder at the sintering temperature. It was found that the decomposition of W3C03C by carburization during sintering is responsible for the formation of very large grains. THE properties of sintered WC-Co alloys are critically dependent upon their final composition and structure. Slight changes of cobalt and carbon contents and small variations of grain size result in marked changes of mechanical properties and performance of products. Previous investigators have studied the microstructure as a function of composition' and have dealt with the influence of cobalt upon properties.' The present work is concerned with the effect of carbon content upon constitution, microstructure, and grain growth during sintering, as well as upon the physical properties of sintered compacts. Phase Boundaries at the Sintering Temperature The WC-Co section through the ternary diagram W-Co-C has been studied by Wyman and Kelley,' Takeda,V andford and Trent," and Rautala and Norton." Two phases, namely WC and a cobalt-rich solid solution (P), appear in the sintered structure if the carbon composition corresponds very nearly to the theoretical carbon content of WC (6.12 pct C). But small deviations from the ideal carbon content will bring about the occurrence of either graphite or the double carbide W,Co3C (7) . The 7 phase was long believed to be metastable, this conclusion being based mainly on the work of Takeda. More recent studies have shown that it is a stable phase of the ternary system,"5 but the actual phase boundaries have not been established and there has been some discussion as to whether or not the 7 phase occurs at some stage of sintering in the fully carburized alloys." In view of this uncertainty it was thought advisable to undertake an investigation of the phases present during the sintering of WC-Co alloys. Commercially available tungsten, cobalt, and carbon powders were used in the present study. A typical analysis of impurities in a sintered compact is reported in Table I. The desired carbon compositions were attained by varying the amount of carbon added to the tungsten during carburization and controlling the composition during sintering by packing the compacts into beds of mixed alundum-carbon powders of suitably selected ratios. The alloys were sintered in a hydrogen tube furnace at 1400°C, then brought to equilibrium at the temperature under investigation and quenched in water. The high temperature phases were identified by metallographic and X-ray diffraction techniques. The samples were etched with a 10 pct solution of alkaline potassium ferricyanide. The phase boundaries, established by metallo-graphic examination, are shown in Fig. 1, which is a vertical section through the ternary diagram at 16 pct Co. The temperatures of the invariant planes, at 1298" and 1357°C were established by Rautala and Norton. The vertical section shows that the three-phase field, WC + 7 + L, extends into the two-phase field, WC + L, but does not include the stoichiometric composition of WC. In practice this means that the 7

Jan 1, 1955

By E. H. Edwards, J. Washburn, E. R. Parker

The mechanism of strain hardening was discussed in connection with some recent observations on the stress-induced motion of dislocation boundaries and on the simple shear deformation of zinc, cadmium, and copper crystals. THE nature of work hardening accompanying plastic deformation of metals has long been the subject of speculation. As early as 1925, Becker' suggested a thermodynamic process to explain slip based on the statistical probability for the occurrence of local glide steps. The theory that plastic glide was due to the presence of lattice defects known as dislocations was first introduced in 1934 by Taylor: and Orowan, and Polanyi.4 Taylor presented a rather complete analysis of the behavior and interaction of dislocations. He assumed a certain type of stress field to be associated with a single dislocation, and postulated that work hardening would result from plastic flow in two ways. First, dislocations following one another across a slip plane might encounter a barrier which would impede their movement, thereby gradually diminishing and eventually stopping plastic flow on that plane. The second mechanism involved the concept of dislocations traveling in opposite directions on nearby planes attracting one another to form a meta-stable lattice. The interaction of dislocations in such an array was shown to be large enough to require an increase in stress before plastic flow would continue. Kochendorfer5,6 extended the concept of work hardening involving a back stress due to an accumulation of like dislocations on a slip plane. He suggested that the formation of a new dislocation is hindered by those dislocations already present. The hardening resulting from the interaction of newly forming dislocations with bound dislocations was designated as "formative hardening," and indicated to be the only type of hardening necessarily connected with the slip process. In another detailed picture of work hardening, Mott and Nabarro7 assumed that sufficient dislocations are primarily present in the crystal as a consequence of growth irregularities (mosaic structure) to initiate slip. Work hardening resulted from the interaction of migrating dislocations with localized internal stresses produced during deformation. More recently, however, Mott8 as modified his treatment to include the generation of dislocations at Frank-Read sources and he accounts for the hardening by the formation of sessile dislocations. It has been widely accepted that metals of hexagonal structure, usually possessing only one set of glide planes, exhibit mechanical behavior markedly different from those having cubic structures. However, the differences in behavior of the two classes may be less fundamental than suggested. Rohm and Kochendorfer9 have shown that a single crystal of aluminum subjected to approximately simple shear gives the type of stress-strain curve associated with the hexagonal metals. Similar behavior has been demonstrated for gold and silver crystals by Andrade and Henderson.'" Thus when glide in cubic metals is restricted to one plane, the strain hardening characteristics closely resemble those of the hexagonal metals. It would appear, therefore, that an investigation of some of the strain hardening properties of single crystals sheared in simple glide might provide a more complete understanding of the phenomenon of work hardening. This report relates a number of recent observations made on single crystals of copper, zinc, and cadmium tested in simple shear. In addition, pertinent observations of effects accompanying stress-induced movement of dislocation boundaries are reported. Experimental Procedure and Results Details of the production and advantages of the type of single-crystal shear specimen employed for the following tests have been presented in a previous publication.n The method of testing makes

Jan 1, 1954

By Robert B. Hall, Charles W. Bauer

Alunite, KA13(SO4)2(OH)6, was used from the 15th century until early in the present century as a source of potash alum and aluminum sulfate. Since early in this century its use for this purpose has declined to virtually nil with acid treatment of bauxite or clay replacing alunite as a source of alum. During the First World War alunite served as an emergency source of potassium sulfate fertilizer in the United States and Australia. More recently, alunite has been investigated as a nonbauxite ore of aluminum, with potassium sulfate and sulfuric acid as recoverable byproducts. The Soviet Union established a commercial-scale plant in the mid-1960s producing cell-grade alumina from alunitic ore and recovering byproduct potassium sulfate and sulfuric acid. Interest in nonbauxite ores, including alunite, is expected to continue in countries that presently are dependent on imports of bauxite or alumina to support their aluminum industries. The US has the largest aluminum producing and fabricating industry in the world but imports more than 90% of its aluminum-bearing raw materials. This chapter reviews the status of alunite as a potential source of aluminum with byproduct fertilizer, and the resources available to support this industry, whenever this becomes economically feasible. Historical Background Alunite has received little attention until recently, yet for centuries it held a higher place than many of the industrial minerals considered essential today. The deposits of alunite at Tolfa, northwest of Rome, were mined almost continuously from 1462 until about the middle of the present century (DeLaunay, 1907; Lombardi, 1977). The potash alum made from the alunite, called Roman alum, was a valued commodity in international trade, being used in the making of paper, in tanning leather, and as a mordant in textile dyeing. The industry at Tolfa is said to have been controlled by the Vatican and was a source of revenue for the Papal State until fairly recent times (Lombardi et al., 1977). Since early in the present century, the preferred method for making alum has been by acid treatment of bauxite, clay, or other aluminous material, so that alunite has fallen into near-total disuse as a source of alum chemicals. Appreciable quantities of mixed alunite-kaolin rock still are quarried in the Tolfa district for use by the Italian cement industry (Lombardi and Mattias, 1979), but the role of this material in the making of Italian cement is not clear. A large deposit of alunite-bearing rock in Chekiang province, mainland China, is said to have yielded a "great quantity" of alum over a period of several centuries (Yih, 1931). During the 19th century, alunite deposits in France, Spain, Australia, and other countries were exploited as raw material for making potash alum and aluminum sulfate. Interest in alunite was revived during the First World War, not for alum, but as an emergency source of potassium sulfate fertilizer. Large veins of nearly pure alunite had just been discovered on Alunite Ridge, 11 km southwest of Marysvale in Piute County, Utah (Butler and Gale, 1912), and these were quickly brought into production for making potash fertilizer to replace supplies cut off from the famous deposits at Stassfurt, Germany. An estimated 225 kt of Marysvale alunite was mined and processed at that time (Callaghan, 1973). Similarly, Australian alunite deposits were mined as a source of potash fertilizer during both world wars (Hall, 1978).

Jan 1, 1983

By Ihor A. Kunasz

Introduction Lithium minerals occur predominantly in pegmatites which contain mineral assemblages derived from the crystallization of postmagmatic fluids or from the metasomatic action by residual pegmatitic fluids. They have been the traditional sources of raw materials for ceramic and chemical industries. With the discovery of Searles Lake, brines became new sources of lithium. Indeed, the 1960s may well be referred to as the "brine decade" because several brine bodies of significant economic importance were discovered. In the United States, the brines of Clayton Valley, Nevada, are presently exploited by Foote Mineral Co. The brines of the Great Salt Lake of Utah, the Smackover formation, and the Imperial Valley geothermal field, have been defined as important resources of lithium. In South America, new occurrences have been identified in Bolivia and Argentina. At the Salar de Atacama, Chile, new production capacity is scheduled for 1984. Following the termination of the AEC purchase program in 1959, several companies such as Maywood Chemical Co., American Lithium Chemicals, and Quebec Lithium Corp. were forced to close their operations. Nevertheless, with the development of new applications in glass ceramics, air conditioning systems, synthetic rubber, and metallurgy (aluminum potlines), the lithium industry grew steadily, and with the potential application of lithium metal in batteries and nuclear reactors, a healthy growth can be foreseen. New products and new applications are likely to stimulate exploration for expanded, reserves as well as the search for new sources from which lithium could be extracted economically. Geochemistry of Lithium Lithium is the third element in the periodic table. It is the lightest of all the metals, having an atomic weight of 6.938, an ionic radius of 0.68Å and a charge of +1. The geochemistry of lithium has been extensively studied and has been summarized by Rankama and Sahama (1950), Goldschmidt (1937), and Hortsman (1957). The distribution of lithium in igneous rocks is controlled by its size and its charge, and by the (MgO+FeO) / Li2O ratio. In the early stages of crystallization of a magma, that ratio is very large. Consequently, both magnesium and iron are removed by ferromagnesian minerals in preference to lithium which is then concentrated in the residual magma. This results in an enrichment of lithium in silicic rocks and pegmatites (Strock, 1936). Pegmatites are coarse-grained igneous rocks formed by the crystallization of postmagmatic fluids. Minerals within pegmatites may also form by metasomatism (Jahns, 1955). Genetically, the pegmatites are associated with neighboring intrusives. Mineralogically, granitic pegmatites contain feldspar, quartz, and mica as the main constituents and a variety of exotic elements such as lithium, beryllium, tantalum, tin, and cesium, which may or may not occur in economically significant concentrations. Detailed studies by numerous investigators (Cameron et al., 1949, 1954; Hanley et al., 1950; Jahns, 1953, 1955; Page et al., 1953) indicate that pegmatites often exhibit an internal zonal arrangement, with each zone containing a specific suite of minerals. The lithium minerals are usually found in the intermediate zones, and although as many as 13 zones have been recognized by Cameron, et al. (1949), a complete zonal arrangement is rarely found. Zoning of pegmatite bodies has also been observed on a regional basis. The regionally zoned pegmatite sequences exhibit mineral assemblages and complexity according to their re-

Jan 1, 1983

By F. W. Warnaars

The Los Bronces copper deposit is located on the west side of the Andes Mountains in central Chile about 54 km (34 miles) northeast of Santiago. The deposit consists of a hydrothermal breccia complex that was emplaced in the western part of an earlier porphyry copper mineral system. This system exhibits propylitic, sericitic, silicic to potassic hydrothermal alteration with disseminated and stockwork copper-molybdenum sulfide mineralization over an area of about 12 km2 (5 sq. miles). The Rio Blanco Mine, about 800 m (0.5 miles) east of Los Bronces, is an underground block caving mine located in the central part of this system. The mineralization at Rio Blanco is characterized by one of the highest primary copper grades known for porphyry copper systems. Present mining at Los Bronces is from an open pit located in the northern part of a hydrothermal complex that is composed of at least seven different copper-bearing tourmaline breccia pipes. Together these pipe form one large contiguous body about 2 km (1.2 miles) long and 0.7 km (0.4 miles) wide at the surface. This body crops out at elevations of between 3.5 and 4.1 km (2 and 2.5 miles) above sea level. Each breccia pipe has a characteristic shape, matrix, type and degree of mineralization and alteration. Quartz, tourmaline, specularite, anhydrite, pyrite, chalcopyrite, sericite, chlorite, and rock flour are present in significant quantities in most breccia matrices. The clasts consist mainly of quartz monzonite, andesite, and minor amounts of quartz latite porphyry. The seven different breccia pipes are identified, from oldest to youngest, as Ghost, Central, Western, Infiernillo, Fine Gray, Anhydrite, and Donoso. The breccias coalesce, interfinger, or display gradational contacts. In places they form sharp contacts with each other and with the surrounding intrusive or volcanic host rocks. Early breccias are disrupted or cut by subsequent breccia events. Explosive emplacement was followed by collapse after the release of hydrothermal pressures. The Los Bronces breccia complex has sharp contacts with the surrounding intrusive stock and andesites. These contacts dip inward on the north, west, and south margins. The eastern contact is nearly vertical and/or dips slightly outward and suggests a westward tilting of the body of about 15° after emplacement. The breccia complex is not fully explored at depth, but is known to extend at least below the 3 km (2 mile) elevation in its central eastern part, about 1 km (0.6 miles) below the highest breccia outcrop. Primary copper sulfide minerals are most abundant in the Donoso and Infiernillo breccias where the matrix consists mostly of tourmaline, quartz, sulfides, specularite, and/or chlorite with minor amounts of rock flour, sericite, and clay. The volume of the matrix material is mostly smaller than in the other breccias. Secondary enrichment enhanced the primary copper grade in the southern two-thirds of the Los Bronces breccia complex and in much of the surrounding porphyry system. Primary grades of 0.3% to 0.5% copper have been enriched to approximately 0.9%. This enrichment may extend from surface outcrops to a depth of more than 500 m (1, 640 ft) in more permeable parts of the breccia complex. A leached and oxidized capping of variable thickness covers a large part of the enriched layer and exceeds 300 m (984 ft) at one location. The shape and depth of the enrichment blanket and the capping suggest that the enrichment process is still active.

Jan 1, 1983

By Nicholas J. Grant, Rolf Nordheim, Bill C. Gissen

The system Cb-Re was examined in detail utilizing pure metals, careful melting techniques, and heat treatments. Metallographic and X-my methods were utilized for phase identification. In addition to soli-dus determinations, the composition limits and mode of formation of the intermetallic compounds 0 and X were determined. The binary constitution diagram Cb-Re has been the subject of a number of investigations.'-7 Two intermetallic phases have been described and tentative diagrams have been proposed. It was the aim of this investigation to obtain a complete and more accurate diagram through the use of purer alloys and improved techniques. The diagram worked out by Knapton7 was published after our experiments were concluded. EXPERIMENTAL METHODS The starting materials were Re powder of 99.95 pct purity, supplied by the Chase Brass and Copper Co., and Cb rondelles of 99.6 pct purity, supplied by the Electro Metallurgical Co. Typical analyses are given in Table I. The alloys were prepared from compacted Re powder and Cb rondelles; for critical compositions master alloys were used. The melts varied from 5 to 50 g and were melted nonconsumably in an Heraeus arc furnace, under ti-tanium-gettered laboratory grade argon. It was necessary to invert the buttons and remelt, followed by crushing and further remeltings, the number of which depended on the nature of the fracture of the alloy. The weight of each alloy was carefully checked after each melting operation. It was found that significant weight losses occurred only when loose Re powder was used. These losses amounted to a maximum of 3 wt pct in such instances, but did not occur after the charge was once molten, regardless of the composition. Typical loss of weight after 4 remelts, including crushing, amounted to 0.2 wt pct. Oxygen pickup during successive melting cycles did not appear to be significant, since the weight gain values never exceeded 0.01 wt pct. Control of the composition of the alloys was not too difficult; however, the problem of homogenization was a far greater one. Optimum alloying conditions were finally achieved by a combination of master alloy production and repeated crushing and remelting cycles. The compositions of the alloys utilized to determine the constitution diagram are listed in Table 11. The weight balance method was used ultimately as the sole analytical method, since repeated checks indicated that it was accurate to 0.2 pct, if the weight losses were small. The results were double-checked by a metallographic study of the homogeneity of the samples, since this was presumed to be the more serious factor. Heat treatments at 1200°C or lower were accomplished by sealing the specimens in Vycor tubes, heating in a Globar furnace, followed by air cooling. This was considered to be an adequate quench because of the very low temperatures involved. For temperatures of 1200" to 1900°, alloys were homogenized and heat treated in a tungsten-filament resistance vacuum furnace. From 1900" to 2750" a tantalum-tube resistance furnace was utilized, operating under a vacuum of 105 mm Hg. The temperature measurements up to 1200" were made by means of platinum-platinum 10 pct Rh thermocouples, with an accuracy of 5°C. For temperature measurements greater than 1200°C, optical pyrometers were used. They were calibrated against a standardized instru-

Jan 1, 1962

By E. H. Rennhack, J. F. Libsch

Definite transition behavior was found in unalloyed titanium at 0.13 pct 0 equivalent. The addition of 0.5 Sn, 1.0 Al, 0.5 Al, and 1.0 Sn lowers the tvansition temperature of titanium at oxygen equivalents exceeding 0.22, 0.29, 0.40, and 0.51, respectively. Metallographic evidence indicates that progressive amounts of oxygen restrict deformation at low temperatures resulting eventually in brittle fracture. The presence of tin or aluminum is thought to alter the atomic distribution of oxygen within the titanium lattice in such a manner as to facilitate some degree of slip and twinning. 1 ITANIUM-like iron and zinc experiences a transition from ductile to brittle fracture with decreasing temperature. The presence of sufficient quantities of the interstitial elements, carbon, oxygen, and nitrogen in solid solution, has been shown to progressively lower the room-temperature impact resistance of titanium with increasing interstitial content.1"1 Ogden,et al,2 reported finding no transition from ductile to brittle fracture with decreasing temperature in iodide titanium. However, a definite transition was found to occur at 0.36 pct 0 equivalent in an alloy containing 0.47 pct carbon in solid solution. The oxygen equivalent in this case was obtained by assuming the effects of carbon, oxygen, and nitrogen in combination on the impact properties of titanium to be additive on the same basis as their separate effects on tensile strength, namely: Oe = 2/3 (pct C) + pet 0 + 2 (pet N) They concluded that the tendency toward transition behavior is dependent on the total interstitial content in solid solution and not on any one particular interstitial element. Rennhack,4 however, recently found transition behavior in iodide titanium at an oxygen equivalent of 0.17 pct. Ogden, etal3, reported that the strengthening effects of carbon, oxygen, and nitrogen in titanium decrease with the addition of tin. They showed that tin additions increase the tolerance of titanium for the interstitials, that is, 10 pct Sn permitted individual additions of up to 0.4 wt pct of the interstitial elements without any significant loss in ductility. Aluminum tends to decrease ductility and impact strength of titanium for a given interstitial ~ontent.3 Van Thyne and Kessler5 found that the solid solubility of carbon in titanium doubled in the presence of 10 pct Al. The precise effects of tin and aluminum in the nrespncp of rorvfen nn the impact properties of titanium had previously not been determined. The high oxygen tolerance of titanium containing tin may provide a suitable means for improving the impact properties of titanium containing the interstitial elements, particularly oxygen. Unlike tin, however, aluminum has a high affinity for oxygen and will reduce rutile titanium dioxide (TiO2)by a thermit-type reaction.6 Thus, the atomic distribution of oxygen in titanium might be different in the presence of aluminum. The present work was undertaken to: 1) Determine the effect of oxygen on the transition behavior of unalloyed iodide titanium. 2) Determine the effect of small amounts of tin and aluminum on the transition behavior of iodide titanium in the presence of oxygen and other interstitial~. EXPERIMENTAL Three different lots of iodide titanium were employed to prepare the alloys studied in the present work. Their chemical composition is presented in Table I. The compositions of the tin and aluminum used for alloying are shown in Table 11. Oxygen additions were made using calcined rutile titanium dioxide (TiO2)of 99.7 pct purity. Alloy Preparation—Four titanium-base alloys containing separate additions of 0.5 and 1.0 pct each

Jan 1, 1960

By Howard Martens

HARDNESS behavior of commercial titanium alloys following various heat treating processes has been studied for some time. However, the hardness of such alloys following a definite measured cooling rate from the single phase region has not been reported. Therefore, a series of experiments was conducted on the three alloys: 150A, 130A, and 130B. The manufacturers gave the following compositions for these alloys: Ti 150A: 0.045 pct C, 0.069 pctN, 1.32 pct Fe, and 2.68 pct Cr; RC 130A: 0.13 pct C, and 7.9 pct Mn: RC 130B: 0.11 pct C, 3.5 pct Mn, and 3.2 pct Al. A method similar to that described by Greninger was used, and the cooling rate varied from 6" to 9200°F per sec. The cooling curves were recorded on a Speedomax or a recording oscillograph, and the hardness was measured on a Tukon hardness tester using a Knoop indentor with a 500 g load. The curves of Fig. 1 show the variation of the Knoop hardness with the cooling rate for these three alloys. The alloys 150A and 130B behave in a similar manner, as shown by the curves. In neither of these alloys was it; possible to retain the high temperature ß phases, even by the most rapid cooling rate. The slow cooling rates produced a typical ß structure in these two alloys. As the cooling was increased, these phalses became more finely dispersed and the hardness increased. This increase in hardness with cooling rate continued for the 150A alloy until a rate of approximately 1000°F per sec was reached. At this cooling rate, the structure showed the result of the dif-fusionless transformation of the ß solid solution to the supersaturated a solid solution which is referred as a'. At higher cooling rates, the structure showed no change and the hardness showed no great change, as indicated by the dotted portion of the curve in Fig. 1. The increase in hardness with cooling rate for the 130B alloy continued until a rate of approximately 300°F per sec was reached. At this rate the structure was produced and the maximum hardness was reached. Higher cooling rates caused no marked changes in hardness or structure. The alloy 130A behaved in a different manner, as shown by the curve of Fig. 1. The structure produced by the slow cooling rates was a typical a+ß structure, which became finely dispersed as the cooling rate increased. The hardness increased with cooling rate and reached a maximum at approximately 120°F per sec. At this maximum hardness, the structure showed the result of the decomposition of the ß phase into the ß + phase structure. As the cooling rate was increased, smaller amounts of the at phase were formed and the alloy became softer. For cooling rates of 3000°F per sec and above, the high temperature ß phase was retained and the hardness of the alloy did not change markedly. Acknowledgment This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory, California Institute of Technology, under Contract No. DA-04-495-0rd 18, sponsored by the Dept. of the Army, Ordnance Corps.

Jan 1, 1957

By Arthur Taggart

Recent jig practice has shown such marked departures from the pronouncements of the textbooks, particularly as to particle size recovered and size range of feed, as to make it desirable to reexamine the mechanism of jig action. This paper presents the results of such a study. The general conclusions of the Paper are: I. A jig bed is a bi-functioned as well as bi-parted body. The upper layer comprises a plastic stream moving slowly in a substantially horizontal direction from feed lip to tail board. It serves to rough out tailing and light middling and to pick up heavy middling and lift it out 01 the jig box. The lower layer acts as the ultimate separating means by absorbing concentrate and rejecting heavy middling. 2. The separating function of the lower layer is performed in three definitely diflerent ways, depending upon the relative sizes of the grains presented to it and those composing it. Experimental The experimental work consisted of close and long-continued visual observations of operating beds in a plunger jig, eccentric driven, and in determinations of the settling rates of various metal spheres through a quartz bed under various conditions of operation. Falling velocities were measured by an ingenious method devised by A. 6. Dorenfeld, a student in the Mining Department at Colum- bia. Spheres of different metallic alloys, 3/8-in. diameter, were cast with a short piece of fine copper wire frozen in to Serve for suspension-A short piece of black cotton thread was looped to the top of the wire to serve as an indicator. The set-up for a test is shown in Fig. I. The thread should be boiled, if necessary, to make it freely water absorbent. The piece of clean plate glass D is supported rigidly at the proper height in a vertical position. The white scale E is similarly supported behind the glass. At the beginning of a test, D is filmed with water by holding it (best with clean rubber finger cots) under running tap water. It is then clamped in approximate position and C, already wet, is laid across it and drawn straight and vertical with the upper end at or near the top of D, and the ball A slightly above the top of the expanded bed. The water film holds the string to the glass. For a test the string is released, its upper end position noted just after the top of the ball becomes submerged, when the corresponding time is read, and thereafter simultaneous readings of position and time are taken as desired. Results of the visual studies of the bed are given at the pertinent places in the following discussion. Results of the settling tests are summarized in Figs- 2 and 3. Definitions The name bed as used herein refers to the entire mixture of solids and water filling the jig box. A layer is a stratum of a bed composed of particles of substantially the same specific gravities. A granular product is long-range when it approximates a crushed or ground product limited as to upper size only; a short-rafzge product is one closely sized. A bed is loaded when operating with continuous feed and discharge of products; it is unloaded when pulsating norma11y but neither being fed nor

Jan 1, 1943

By J. R. Kyte, C. C. Mattax

Previous workers have developed differential equations to describe oil displacement by water imbibition, but have not explicitly defined the relationship between recovery behavior for a single reservoir matrix block and its size. In the present work, imbibition theory is extended to show that the time required to recover a given fraction of the oil from a matrix block is proportional to the square of the distance between fractures. Using this relationship, recovery behavior for a large reservoir matrix block is predicted from an imbibition test on a small reservoir core sample. The prediction is then extended to analyze recovery behavior for fractured- matrix, water- drive reservoirs in which imbibition is the dominant oil-producing mechanism. Experimental data are presented to support the basic imbibition theory relating matrix block size, fluid viscosity level and permeability to recovery behavior. INTRODUCTION Imbibition has long been recognized as an important factor in recovering oil from water-wet, fractured-matrix reservoirs subjected to water flood or water drive. l-3 Recently, two approaches have been published which might be used to predict imbibition oil-recovery behavior for reservoir-sized matrix blocks. Graham and Richardson4 used a synthetic model to scale a single element of a fractured-matrix reservoir. lair,' on the other hand, used numerical techniques to solve the differential equations describing imbibition in linear and radial systems. This latter method requires auxiliary experimental data in the form of capillary pressure and relative permeability functions. These two approaches, i.e., synthetic models and numerical techniques, have been used to study a variety of reservoir fluid-flow problems. One pur- pose of this work is to present, with experimental verification, a third method for predicting imbibition oil recovery for large reservoir matrix blocks. This method uses scaled imbibition tests on small reservoir core samples to predict field performance. The imbibition tests are easier to perform than the capillary pressure and relative permeability tests required to apply the numerical method. Furthermore, when suitably preserved reservoir-rock samples are used, the properties of the laboratory system are the same as those of the field. This offers an important advantage over the use of synthetic models because there is usually some question as to how accurately reservoir-rock properties can be duplicated in such models. Based on the recovery behavior for a unit matrix block, an analysis is presented to predict oil recovery for a fractured, water-drive reservoir made up of many such unit blocks. In the analysis, it is assumed that the flow resistance and the volume of the reservoir fracture system are negligible compared with that of the porous matrix. These assumptions are generally consistent with observed characteristics of many fractured-matrix reservoirs, and have been employed in previous studies.2,3,6-8 further assumed that the effect of gravity on flow in the matrix blocks is negligible. On first thought, the latter assumption might appear to seriously limit application of the method. However, in a fractured reservoir, the effect of gravity on flow in a matrix block will be restricted by the height of the block. Furthermore, matrix permeabilities are often very low (10 md or less) in such reservoirs. This means that capillary or imbibition forces will be large, thus tending to minimize the relative importance of gravity. For these reasons, imbibition should be the dominant oil-recovery mechanism in many fractured-matrix, water-drive reservoirs. The predictive method presented in this report is applicable to such reservoirs. SCALING IMBIBITION THEORY The "scaling laws" applicable to the operation.

By Noel Jarrett, Allen S. Russell, Bernard M. Starner, Stanley C. Jacobs

A process has been developed to refine high-grade commercial aluminum to 99.99 pct purity. This enzploys precipitating titanium, vanadium, and zirconiu~ as borides. The upgraded liquid is partially crystallized under conditions that prevent massive freezing. The resulting crystals are compacted and washed free of adhernig impure liquid aluminum by melting the top crystals while draining off the liquid. Finally, the bottorn portion of the crystal bed is melted and removed as the pure product. Leakage from the crystallizat~on vessel is prevented by a layer of alumina granules that are not wetted by the aluminum. ALUMINUM of 99.99 pct purity is desired in a variety of applications. These include electrolytic capacitors, catalyst carriers, decorative trim for automobiles, aluminum powder, and high-strength alloys. This market was for many years filled by refining commercial-purity aluminum in a three-liquid-layer electrolytic cell of the type originally proposed by Hoopes. To meet the expanding need for lower-cost high-purity aluminum, a refining process based on chemical precipitation and fractional crystallization was developed. The availability of this super-purity metal was announced in November 1961, and three patents1-3 relating to the process were issued in 1965. An alternate process, comprising continuous fractional crystallization of a portion of a liquid-aluminum feed stream under conditions of strong agitation at the surface of the growing crystals, was published by Dewey4 in 1965. Earlier proposals by Larsen 5 and by Regner6 also emphasized the advantage of rapid agitation of the melt in the region of the crystallization zone. johnson 7 more recently proposed extracting purified aluminum from an alloy by mixing crystals of pure aluminum into the molten alloy, allowing them to grow, and removing the pure aluminum as by casting. Aamot's 8 recent patent describes a freeze refining process that involves controlled partial freezing against cooled moving retort walls. The Alcoa process to produce high-purity aluminum is schematically described in Fig. 1. Relatively pure aluminum is treated with boron in a furnace (a) to precipitate high melting point impurities, the supernatant liquid is transferred to a crystallization vessel (b), the liquid is partially crystallized with mechanized "tamping" (c), the remaining impure liquid is tapped out and the crystal bed partially and directionally re-melted to wash the underlying crystals (d), and the balance of the crystals is melted and withdrawn as high-purity aluminum (e). PRINCIPLES OF FRACTIONAL CRYSTALLIZATION Commonly encountered impurity elements form two types of systems with aluminum. Titanium, vanadium, and some other transition elements in periodic groups IV, V, and VI are the high melting point impurities. Initial solidification of aluminum containing small percentages of these elements will produce solid markedly enriched in them at temperatures slightly above the freezing point of high-purity aluminum. Most other elements, such as silicon and iron, produce eutectic type reactions with aluminum. In this case the first solid formed on initial solidification at temperatures slightly below the freezing point for high-purity aluminum is markedly depleted while the remaining liquid is enriched in these elements. The equilibrium change in composition during solidification of aluminum contaminated with silicon, iron, or titanium is illustrated in Fig. 2. Here, for clarity, the temperature scale is expanded beyond the

Jan 1, 1968

By Arthur Taggart

Recent jig practice has shown such marked departures from the pronouncements of the textbooks, particularly as to particle size recovered and size range of feed, as to make it desirable to reexamine the mechanism of jig action. This paper presents the results of such a study. The general conclusions of the Paper are: I. A jig bed is a bi-functioned as well as bi-parted body. The upper layer comprises a plastic stream moving slowly in a substantially horizontal direction from feed lip to tail board. It serves to rough out tailing and light middling and to pick up heavy middling and lift it out 01 the jig box. The lower layer acts as the ultimate separating means by absorbing concentrate and rejecting heavy middling. 2. The separating function of the lower layer is performed in three definitely diflerent ways, depending upon the relative sizes of the grains presented to it and those composing it. Experimental The experimental work consisted of close and long-continued visual observations of operating beds in a plunger jig, eccentric driven, and in determinations of the settling rates of various metal spheres through a quartz bed under various conditions of operation. Falling velocities were measured by an ingenious method devised by A. 6. Dorenfeld, a student in the Mining Department at Colum- bia. Spheres of different metallic alloys, 3/8-in. diameter, were cast with a short piece of fine copper wire frozen in to Serve for suspension-A short piece of black cotton thread was looped to the top of the wire to serve as an indicator. The set-up for a test is shown in Fig. I. The thread should be boiled, if necessary, to make it freely water absorbent. The piece of clean plate glass D is supported rigidly at the proper height in a vertical position. The white scale E is similarly supported behind the glass. At the beginning of a test, D is filmed with water by holding it (best with clean rubber finger cots) under running tap water. It is then clamped in approximate position and C, already wet, is laid across it and drawn straight and vertical with the upper end at or near the top of D, and the ball A slightly above the top of the expanded bed. The water film holds the string to the glass. For a test the string is released, its upper end position noted just after the top of the ball becomes submerged, when the corresponding time is read, and thereafter simultaneous readings of position and time are taken as desired. Results of the visual studies of the bed are given at the pertinent places in the following discussion. Results of the settling tests are summarized in Figs- 2 and 3. Definitions The name bed as used herein refers to the entire mixture of solids and water filling the jig box. A layer is a stratum of a bed composed of particles of substantially the same specific gravities. A granular product is long-range when it approximates a crushed or ground product limited as to upper size only; a short-rafzge product is one closely sized. A bed is loaded when operating with continuous feed and discharge of products; it is unloaded when pulsating norma11y but neither being fed nor

Jan 1, 1943

By H. S. Price, D. A. T. Donohue

The system of equations describing displacement of a hydrocarbon liquid by a hydrocarbon vapor in a porous medium where mass transfer takes place between the phases is solved numerically for a variety of gas injection processes. Even though the method of solution is quite general, only systems with three hydrocarbon components are considered. Computer simulations of displacement processes wherein mass transfer between phases is both considered and neglected are compared, and it is shown that neglecting mass transfer can give pessimistic displacement efficiencies. INTRODUCTION The role of the gas displacement process in the recovery of petroleum has been subjected to a series of detailed analyses; as a result, a number of predictive models have been published in the literature. However, because of major simplifying assumptions, most of these models do not completely represent the physical system. As a result, the effect of making the simplifying assumptions is unknown. Therefore, a complete representation of this process — one without major simplifying assumptions —should lead to a fuller understanding of the process, and perhaps to methods of improving it. The general method of developing a model for two-phase fluid flow in a porous medium is to solve simultaneously the continuity equation, the energy equations and the equation-of-state for each phase under the prescribed initial and boundary conditions. For an isothermal system, the energy equations reduce to the momentum equation. Darcy's law. However, since natural gas is the vapor state of the reservoir liquid, interphase mass transfer may take place with concomitant changes in both the intensive and extensive thermodynamic properties of each phase. It is this phenomenon that has often been omitted in previous mathematical models. An additional relation, then, which accounts for mass transfer between the phases, must be included with the other equations to specify a complete model. Completely formulating the equations to be solved is not a difficult task1 but obtaining their solution has been intractable up to now. Availability of large-memory, high-speed digital computers now makes an attack on this formidable problem possible. This paper presents a preliminary study of the problem. Since this investigation is intended to be exploratory, it is restricted to the linear, horizontal, isothermal, two-phase viscous flow of oil and gas in an oil reservoir. In the early development of predictive models of this process, the reservoir system was considered as a unit and various forms of the material balance equation were proposed.2-4 Pressure and saturation gradients were then added in the Buckley-Leverett model.5 The Buckley-Leverett formulation considered the fluids to be incompressible; thus, the mathematical model reduces to a steady-state system. In the 1950's, studies incorporating numerical techniques were being published.6-9 These mathematical models differed in the efficiency of finite difference techniques, the inclusion or exclusion of capillarity or the number of space dimensions considered. To solve these nonlinear, partial differential equations, each phase was considered to be homogeneous with time; therefore, mass transfer between phases was neglected. The effect of mass transfer on the gas displacement process was first reported by Attra. 10 He simulated the one-dimension flow system by a series of cells in each of which the fluids were equilibrated during a time step. In addition, the pressure throughout the system during each time step was predetermined and constant-phase velocities were calculated according to the Buckley-Leverett incompressible fluid flow model. Welge et al.11 developed a model for the displacement of oil by an enriched gas where composition is considered to be a dependent variable. Fluid

By J. M. Smith, D. Kunii

Effective thermal conductivities of sandstones filled with stagnant fluids were measured using a steady-state technique. Data were obtained for seven sandstone samples, taken from four different locations and ranging in permeability from 18 to 590 md. The measurements with gases (helium, nitrogen, air and carbon dioxide) covered a pressure range from 0.039 psia to 400 psig. Data were taken for four liquids — n-heptane, methyl alcohol, 79.8 weight per cent glycerol-water solution and pure water at atmospheric pressure. The experimental results were used to evaluate the theoretical equations for predicting stagnant conductivities developed earlier 4 The low-pressure measurements permitted evaluation of the consolidation parameter bpDp/ks (necessary to utilize the theory) for the various types of sandstones. Using these characteristic values, the tbeoretical equations correlated well with the experimental conductivity data for the several fluids and rock samples. INTRODUCTION An aspect of heat transfer in solid-fluid systems of considerable current interest is the effective thermal conductivity of porous media. The stimulus for study of the subject arises from the need for sound procedures for designing thermal methods of petroleum production. The general system occurs when there exists a flow of fluid through the pores of the solid material. However, a logical starting point in developing a theory for predicting the effective thermal conductivity in the general system is to attack the special case when the porous solid is filled with stagnant fluid. Since the flow rates anticipated in thermal production processes are very low, such stagnant conductivities k,O are also of practical significance. Kunii and Smith4 recently reviewed the existing data for stagnant conductivities and proposed a theory for heat transfer in porous rocks. This leads to equations for predicting the stagnant conductivity as a function of the properties of the fluid and solid phases. The existing experimental information1-3,5,7,8 covered a range of porosities and agreed well with the theory. However, the available data were insufficient to examine critically the effects of two important properties on koe, — the thermal conductivity of the fluid and the pressure. The objectives of the present study are twofold: (1) to measure the effect on koe of (a) the pressure for gases, and (L) the thermal conductivity of the fluid kg for bloth gases and liquids; and (2) to determine for rocks of different porosities and permeabilities the apparent solid conductivity ks and the consolidation parameter bpDp/ks required in applying the theoletical equations. Measurements extrapolated to zero pressure are helpful in evaluating ks and the consolidation parameter. However, the extrapolations must be based upon data obtained at low pressures because in small pores the molecular conductivity of the fluid may be unusually depressed. This occurs when the mean free path of the molecules is of the same order of magnitude as the pore diameter, that is, in small-diameter pores at low pressures. The phenomenon has been studied recently by Schotte6in beds of unconsolidated fine particles. With particles in the size range of 45 to 200 microns, a significant reduction of kg with pressure was found even at absolute pressures greater than 1 atm. In the sandstones employed in the present investigation, pore diameters lower than those experienced by Schotte are likely, so that the effect would presumably be greater. SCOPE OF MEASUREMENTS Stagnant conductivities were measured for four types of sandstones using a total of seven samples. The characteristics of the samples, including porosity and permeability values, are given in Table

By C. A. Hartnagel

During- the past 10 years the annual production of petroleum in New York has averaged close to 5,000,000 bbl., the total for the period being 49,881,000 bbl. In 6 of the 10 years, the production was slightly in excess of the 5,000,000 mark. The 1944 production of 4,697,000 bbl. represents a 7 per cent decline from the previous year. Of the total 1944 output, the Allegany County district, which includes two small areas of Steuben County, contributcd 3,614,065 bbl., a decline of 5 per cent from 1943, and the contribution of Cattaraugus County, which represents the northern extension of the Bradford pool of Pennsylvania into New York was 1,082,935 bbl., a decline of more than 14 per cent. Since March 1942, the posted price of New York crude oil has remained at $3.00 per bbl. On lug. I, 1944, a Federal subsidy of 75¢ per bbl. was granted to oil producers, thus making the price actually $3.75 per bbl. Since all of the oil produced in New York comes under the designation "stripper well operations," it is entitled to the Federal subsidy. During the last few years there has been little change in the number of wells drilled annually. In the Allegany field, which accounts for three fourths of the oil produced in the state, 1240 wells, including water-intake wells, were drilled during 1944 as compared with 1223 in 1943. There was no marked increase in the number of wells drilled after granting the subsidy of 75¢ per bbl. of oil on Aug. I, the increase being but 14 wells during the last half of the year over the preceding six months. In flooding operations it is not to be expected that even a substantial increase in the price of crude oil will be reflected quickly in an increase of the number of wells drilled or in a marked increase in production of crude oil. In flooded areas, oil is obtained through a program extending over a number of years. This program involves the drilling of the water or input pressure wells, drilling of the oil wells, building up the oil flood to maximum production, and finally the watering out of the well when the amount of oil produced in proportion to the water becomes so small that pumping is no longer profitable. The various operations may extend over a period of 10 years or more, during which time the price of crude may have fluctuated greatly. Although the rate of production can be controlled to a limited extent, it is evident that a flood when once begun must be continued regardless of the price of crude. During the past seven years New York crude has had a price range from $1.68 per bbl. in 1938 to the present price of $3.75 per barrel The field pattern now in general use for wells in a flood project is the "five-spot." Under this arrangement, the lease is divided into squares with a water well at each corner and an oil well in the center. In the all-over pattern of the lease each oil well is acted upon by four water wells and each water well supplies part of the pressure to four oil wells. Spacing between the input water wells and the oil wells usually varies from zoo to 225 ft. or, on an acre basis, there is one oil well and one water well for each 2 1/2 acres. In recent years there has been a tendency to increase

Jan 1, 1945

By R. A. Ritter, J. P. Batycky

A numerical technique has been developed to permit estimating the pressure gradient associated with laminar flow of thixotropic liquids through long Pipelines. For this purpose the pipeline is divided into a number of radial and longitudinal increments within which rheological properties of the fluid may be considered as constant at any time. Then, provided only that the fluid flow curve is defined at every duration of shear, it is possible to predict the instantaneous pressure gradient at any cross-section along the pipeline for each desired flow rate and pipe size. The technique consists of an iterative integration of shear rate to establish the appropriate value of the wall shear stress at each location. Consistency of fluid in the increment is determined by the flow history of that increment, while the radial flow associated with variations in velocity profile is accounted for by adjusting the width and radial position of the increment. A number of pressure profiles, computed at each of several flow rates, provide a convenient basis for pipeline design and pump selection. INTRODUCTION In recent years considerable attention has been given to predicting pressure drop associated with the isothermal laminar flow of time-independent non-Newtonian fluids in pipes and annuli. 6,11,12 The approach generally has been to develop analytical relationships between flow rate and pressure drop based on simple constitutive models which hopefully provide an approximate description of the rheological properties of the fluid. Analytical solutions are highly desirable since the influence of all pertinent parameters can be readily determined. Unfortunately, however, this approach is restricted to simple flow geometries and- frequently leads to erroneous results due to inadequacies in the model. In certain cases a solution may be obtained through applying appropriate numerical techniques.l,14 For example, a digital computer program is available for predicting the velocity profile and pressure drop encountered by any Newtonian or time-independent non-Newtonian fluid flowing under laminar conditions in a cylindrical pipe or annulus.l In this paper the consistency behavior of the fluid need only be described in terms of basic rheological data. Analyzing flow systems involving fluids with time-dependent rheological characteristics is considerably more complicated since substantial changes in consistency may occur because of sustained shear action. This sensitivity to shear frequently persists for several hours. Consequently, variations in pressure drop and/or flow rate resulting from the aging process and addition of unsheared or partially sheared fluid to the system must be considered for purposes of pipeline design. This paper outlines a numerical method for predicting the transient and steady-state laminar flow behavior of a thixotropic liquid in a pipeline of arbitrary length (i.e., at a specified constant flow rate, the instantaneous pressure gradient may be determined at any time after start up and at any location along the pipeline). Several such pressure gradient profiles, computed at several flow rates, may be combined to produce a complete portrait of the system response. This flow portrait provides a reasonable basis for pipeline design and for selecting a suitable pump characteristic. TIME-DEPENDENT RHEOLOGICAL BEHAVIOR The most familiar time-dependent rheological properties are those exhibited by thixotropic liquids. Many of these materials, particularly thixotropic crude oils, generally display an apparent yield stress in that a finite pressure gradient is required to initiate flow. Then, under the influence of sustained shear at a constant shear rate, the consistency systematically decreases to some final limiting value. The shear stress obtained at

By Donald Ofte, L. J. Wittenberg

The absolute viscosity coefficients and activation energies for viscous flow of bismuth, lead, and zinc are reported. The viscosities were measured in an oscillating cup viscosimeter from nearly 1000°C to within 1 deg of the melting point of zinc; 2 deg of the melting point of lead; and in the supercooled condition of bismuth, 9 deg below the melting point. Two methods of relating the change in the viscosity coeflicient with temperature were compared. From the plot of the logarithm of viscosity (77) as a function of reciprocal absolute temperature (1/T), the equations of the lines which fit the data are: AS a part of the Reactor Fuels and Materials Research program at Mound Laboratory, precise measurements of the physical properties of liquid metals are of interest because the desirability for a low vapor pressure liquid coolant in high-temperature nuclear power plants has focused attention on the properties of the liquid metals.' Also, a clearer understanding of the liquid metal state can be achieved when better measurements are available. For these reasons, the viscosity values of certain liquid metals have been determined. The viscosities of the liquid metals, bismuth, lead, and zinc are measured in an oscillating cup viscosi-meter. In this apparatus the liquid is sealed in a right-circular cylindrical cup which is attached to a torsion fiber so that the system forms a torsion pendulum. The damping effect of the molten metal upon the normal oscillations of the torsion pendulum is a function of the viscosity of the liquid. This method is advantageous because of the small amount of sample required (8 to 10 cc), the ease of remote operation of the apparatus in a vacuum system, and the containment of the liquid metal in a hermetically sealed, nonreactive, tantalum crucible. With the liquid metal sealed in the crucible, loss of liquid by vaporization and chemical reactions with atmospheric gases are prevented. Additionally, an absolute viscosity determination is obtained which avoids the sources of possible error in a relative measurement. The difficulties in a relative method include the repetitious cup calibrations, and the uncertainty inherent in the practice of calibration of the apparatus at room temperature and extrapolation to the high temperatures necessary for the investigation of most liquid metals. The oscillating cup viscosimeter was first described by Meyer in 1891.2 The first complete mathematical treatment offering an absolute measurement of viscosity with this type of apparatus was proposed by Andrade and chiong3 in 1936 for an oscillating sphere filled with a liquid. This solution was modified to include a right circular cylinder by Hopkins and Toye4 in 1950, and R. Roscoe5 in 1958. Roscoe's treatment for the calculation of the viscosity of a liquid at a particular temperature from the logarithmic decrement and period of oscillation of the torsion pendulum was used in this investigation. A detailed description of the method of calculating the viscosity from these measurements has been reported earlier.' This calculation has been programmed for the IBM 704 and IBM 1620 computers for expedition of the otherwise laborious calculations. APPARATUS The cutaway view of the apparatus in Fig. 1 shows the suspension system within the high vacuum enclosure. The low residual gas pressure. maintained between 1 x 10-4 and 5 X 10-6 mm Hg during a determination, avoided both air drag on the pendulum and

Jan 1, 1963