Search Documents

By A. Kaufman

More than 3 trillion gal of water are used annually by the mineral industry. Of this, approximately 21/2 trillion gal are recirculated, the rest constituting intake water. The major users are natural gas processing plants and phosphate rock, sand and gravel, and iron ore producers. Water was used by the mineral industry for mining (6%), processing (64%), cooling and condensing (27%), and miscellaneous uses such as boiler feed and sanitary purposes (3%). Whereas total water use is dependent on the quantity of material processed and on the particular process requirements of an industry, recirculation is dependent on processing, as well as cooling and condensing requirements, quality of new water intake, and the necessity for treating new and discharged water. Consumed water, on the other hand, is dependent on the quantity of water recirculated, and temperature and humidity in the area. Based on this analysis, an increase in water use by the mineral industry of 21/2 times by 1985 is forecast. Wster intake, however, will only rise 62%, because of a substantial increase in recirculation. In one report of a special series concerning the water resources of the United States, the Senate Select Committee on National Water Resources, 87th Congress, estimated that water demand would double by 1980 and triple by 2000.' In view of the possible water deficiencies that might result from such expanded usage and the need for research guidance, the Bureau of Mines organized and carried out a statistical canvass of water use in the mineral industry for calendar year 1962. The data used in this paper, unless otherwise noted, are derived from that canvass.' The efforts in this paper are devoted toward summarization of the canvass and analysis of the data. SOME DEFINITIONS Intake: Water introduced from an external source for the first time into a given mine or plant regardless of quality. Intake water is also called new water, water withdrawn, or makeup water. Fresh Water: Water suitable for cooking and drinking. Saline Water: Water containing more than 1000 parts per million of dissolved solids. Contaminated Water: Water not suitable for domestic use, but excluding saline water. Recirculated Water: Water reused to conserve intake water. Solutions that are recycled primarily because of fixed metallurgical practices, such as copper leaching solutions containing sulfuric acid, are excluded. Gross Water Used: Recirculated water plus intake water. Also called total water used. Consumed Water: Water that is lost by evaporation, as well as water lost in product. Seepage and transferred water are not considered consumed. The use of either intake or recirculated water may result in consumption. However, because of difficulties in measuring consulmption, consumed water is defined as the residual between intake and water discharged from the mine or plant. Mineral Industry: For the purposes of this paper, mineral industry includes all metal and nonmetal surface and underground mines and their associated processing plants, as well as custom mills, coal washing plants and associated mines, petroleum and natural gas well drillers, natural gas processing plants, and secondary recovery operations. WATER USE IN THE MINERAL INDUSTRY Water used by the mineral industry constitutes a relatively minor fraction of the water withdrawn by all users. Data compiled by the U.S. Geological Survey in 1960 indicate that water withdrawals, exclusive of that used to generate hydropower, approximate 99 trillion gal annually.3 Our data indicate that water withdrawals by the mineral industry comprise only 1% of this total, or 2% of the water withdrawn by industrial users. Use by Industry: Water used by various mineral industries is shown in Table I. The table indicates that close to one-half of the gross water used by the mineral industry in 1962 was used by natural gas processing plants, followed by sand and gravel, phosphate rock, and iron ore producers. The largest aggregate users of water are also the largest users per dollar of product. For example, Fig. 1 shows that the natural gas processing and phosphate rock industries are very large users of water per dollar of product. Their recirculation per

Jan 1, 1968

By N. E. Brown, F. L. Versnyder

THE apparent dependency of the electrical characteristics of hexagonal crystalline selenium on microstructure has aroused much interest in microscopical studies of selenium. Microscopic observations on the crystallization of selenium have been made by Escoffery and Halperin,' P. H. Keck,' and other investigators. It is the purpose of this paper to discuss the microstructural changes observed on polished cross-sections of single layers of selenium after various heat treatments. Observations were also made on crystallization of the free-surface layer of these deposits. In general, all of the transformations studied were either transformations of the vitreous selenium to hexagonal selenium or micro-structural transformation of the hexagonal selenium itself. Procedure The selenium used in this work was obtained from the American Smelting and Refining Co. and was approximately 99.96 pct pure. An intentional impurity of 1 part per 2,000 of bromine was added to the material prior to evaporation. A thickness of approximately 0.002 in. of this selenium was vapor deposited on an aluminum base plate. The maximum plate temperature during the vacuum vapor deposition was 140°C. Mounting of the cross-sectional specimens for metallographic study could not be done in plastic mounting media, as is customary, since temperatures in excess of 50°C would cause unwanted transformations. Consequently, a simple clamp-type device was used to mount the specimens for preparation. All grinding operations were then done carefully by hand in order that the specimen not become heated during this operation. Wet polishing was done on the conventional metallographic polishing laps, using successively finer grinding powders. An extremely careful polish is necessary, since observation and micrography of the specimens are done in the unetched condition under polarized light. The two observations of crystallization made on the free surface of vitreous selenium deposits (Figs. 4 and 5) were made on surfaces which were perpendicular to the cross-sections studied. These free-surface layers were examined directly, i.e., no pre- vious metallographic preparation, as obtained from the vacuum vapor deposition. Microscopic Observations A study was made of polished cross-sections of the vitreous selenium as-deposited. It was noted that in all cases there was columnar crystallization adjacent to the base plate, which appeared to occur during the vacuum deposition process. This observation has also been made by Keck? It also was observed that vagrant spherulitic crystallization occurred in the vitreous selenium. The term "vagrant" is used, since these spherulitic grains appear to crystallize at random throughout the vitreous selenium during the vacuum deposition process. Columnar crystallization at the A1-Se interface and a typical spherulite observed in a polished cross-section of "as-deposited" vitreous selenium may be seen in Fig. 1. Cross-sectional samples of vitreous selenium studied after heat treating individual samples for 20 min in 10" steps from 80" to 220°C revealed that crystallization—in this case, columnar crystal growth —proceeds from the aluminum base plate to the surface of the specimen (Fig. 2). Crystallization was microscopically observed to be complete after the 130°C heat treatment. Visual examination of the free surface of the specimen after the 130 °C heat treatment revealed the readily recognizable grey appearance of the completely crystallized selenium, in corroboration of the microstructural observations. No microstructural transformations then appeared to take place between 130" and 190°C. At 190°C the beginning of recrystallization appeared and proceeded until the columnar grain structure had been completely transformed to equiaxed grains between 210" and 220°C (Fig. 3). Naturally, the grain size of the recrystallized grains at the lower temperatures (190" to 210°C) was smaller than is illustrated in Fig. 3. In addition, polished cross-sections of deposits heat treated at 140°C for 10 min to cause complete crystallization and, subsequently, heat treated in 10" steps from 80" to 220°C for 20 min were studied. As expected, no microstructural transformations took place until the beginning of recrystallization was observed at 190°C. A comparison with the previously studied specimens revealed that recrystallization proceeded almost identically in the two experiments although in the first case the deposits were vitreous prior to the series of heat treatments and in the second case they had been crystallized by a previous heat treatment. By heat treating for longer times (180 min) at lower temperatures, the

Jan 1, 1956

By J. S. Aronofsky

The first step in a quantitative analysis ot the mechanism of oil displacement by water in a fractured reservoir is usually conceded to be the solution of the differential equation describing the saturation distribution of two immiscible fluids flowing in a porous medium, where the capillary pressure is taken into account. In such a system the production mechanism may consist of displacement of oil both by the flow of water due to natural or artificially imposed pressure gradients and by imbibition, which implies a flow of water not due to external pressure gradients. Owing to the presence of the two oil displacement mechanisms, the mathematical model given by the differential equation intended to describe the system may not properly represent the behavior of the physical system. In fact, in the reservoir the rate of water advance may be very slow, and in the case of a fractured reservoir with a great number of large fractures, the pressure difference determining the flow of water through the matrix may be much less than 1 lb/'psi over lengths of a few feet. In such a case, imbi-bition (the exchange between oil in the matrix and water in the fractures resulting from capillary forces) may become, with time. a significant element of the production mechanism. It occurred Lo the authors, however, that without going into a physical analysis of the process of production, it might be possible by means of simple abstract reasoning to throw some light on the variation of recovery with time under conditions occurring in a highly fractured oil reservoir with rising water table. The object of this paper is to present both the reasoning and its application to a reservoir of the highly fractured type. Specifically, the analysis given here was undertaken to try to explain the increase of recovery (as defined later) with time as observed in this reservoir, without having to assume unlikely variations in the reservoir parameters with depth. This attempt has been successful as will become clear upon comparison of the computed recoveries with the actual field data. ABSTRACT MODEL Let us consider a small volume. of porous matrix saturated with oil at time, t = 0. Let the process of oil displacement by water start at time, t = t,,. At some time, t, the process will have terminated. Then a volume of oil equal to or smaller than the original oil contained in the matrix will have been produced. The first basic assumption that describes the model and guides the forthcoming reasoning is that the oil production from the small volume. dv, is a continuous monotonic function of time and that it converges to a finite limit. Such an assumption is not inconsistent with the results of laboratory waterflood tests as well as with results of imbibition tests where this is, in fact, observed. The second basic assumption is that none of the properties which determine the rate of convergence change sufficiently during the process to affect this rate or the limit. Let it be assumed that the form of the function of time relative to production from the matrix volume. dv, is given by where V,(t) is the volume of oil produced up to that time t. R is the limit toward which the recovery converges, A is a constant giving the rate of convergence, and V,(t) is the volume of oil originally in place in the volume, dv. It should be noted that recovery at time I will be understood here to he It follows that r. tends to R as t tends to infinity. CONSTRUCTION OF RESERVOIR FROM ABSTRACT MODEL Let the reservoir consist of a series of identical blocks of porous matrix stacked vertically and separated by fractures. Let each of these blocks satisfy the conditions of our abstract model. These conditions are: (1) that recovery is a continuous mono-tonic function of time converging to finite limit, and (2) none of the properties that determine the rate of convergence change sufficiently to affect the rate or the limit. Let water be rising in fractures so that oil production from any part of the block starts when the water comes in conLact with it. TOTAL RECOVERY FROM THE RESERVOIR COMPARED TO RECOVERY FROM A SINGLE MODEL ELEMENT As stipulated in Eq. 1, in the case of the abstract model,

By J. Washburn, R. Drouard, E. R. Parker

Temperature dependence of the rate of recovery in zinc single crystals after a simple shear deformation at low temperature was investigated. Some tentative suggestions regarding the annealed and strain-hardened states of a crystal are discussed. RECOVERY may be defined as the gradual return of the mechanical and physical properties of strain-hardened metal to those characteristic of the annealed material; an increase in temperature increases the rate of recovery. The annealing process in strain-hardened polycrystalline metals is complicated by the inhomogeneity of strain which always exists in aggregates. Polygonization in bent regions of the crystals and growth of new almost strain-free grains starting at points of severe local distortion1-:' make it almost impossible to isolate and study the recovery process. Homogeneously strained single crystals, however, do not polygonize or re-crystallize and hence they can be used advantageously to study recovery. In such crystals strain hardening is completely removed by recovery alone. Since recovery is a process whereby certain lattice disturbances introduced by plastic flow are gradually reduced, a knowledge of the rate and temperature dependence of this process for various conditions of prestrain might be helpful in formulating a model of the strain-hardened state. For simplicity it seemed desirable to limit the type of prestrain to the simplest obtainable, i.e., simple shear strain. In the experiments to be described, recovery was studied by observing changes in the stress-strain curve of prestrained zinc single crystals held for various times at temperatures above that employed for straining. Single crystals were grown from the melt by a modified Bridgeman technique from Horse Head Special zinc 99.99 pct pure, and from spectrographically pure zinc 99.999 pct pure. They were grown as 1 in. diameter spheres and acid-machined' to the final specimen contour. The test section was a cylinder about 1/8 in. high and 3/4 in. in diameter. The conical sections adjacent to the test section were cemented into the grips so the load could be transmitted to the crystal as uniformly as possible. The specimens were oriented so that in testing the maximum shear stress was applied along one of the slip directions, [2110], in the (0001) plane. Details of the production and testing of such specimens have been presented.' Each test was carried out according to the following schedule: 1—The crystal was strained at — 50°C until it reached a maximum shear stress, ,,,. The strain rate was approximately 5 pct per min in all cases. 2—After straining, the crystal was unloaded before the temperature was changed. Unloading required about 3 min. 3—The temperature of the specimen was then increased from — 50°C to the temperature, T, of recovery. This change in temperature was completed in a time of less than 2 min. The specimen remained at temperature, T, for a time, t, which differed for the various specimens. 4—Thereafter the temperature was again reduced to — 50 °C in approximately 3 min. 5—While at —50°C, the stress-strain curve after recovery was obtained. 6—The specimen was then unloaded and annealed for 1 hr at 375 °C in a helium atmosphere to bring about complete recovery. Cooling to room temperature after anneal required 90 min. 7—The same crystal could be re-used for another test because the plastic properties after annealing closely duplicated those of the original crystal. The specimen was immersed during the test in a bath of methyl alcohol which, through a system of tubes, could be pumped through either of two heat exchangers to regulate the temperature; this was accomplished by circulating the liquid through coils immersed in a bath of acetone and dry ice for cooling or in a bath of warm water for heating. Test temperatures were thus maintained constant within ±1°C. The — 50°C temperature was low enough so that no measurable recovery occurred during unloading and reloading. The stress-strain curve continued after recovery along a path below, but approximately parallel to, the path of a curve obtained in an uninterrupted test. Fig. 1 shows some of the results from a specimen of 99.999 pct Zn. The amount of downward displacement of the curve due to recovery was a

Jan 1, 1954

By Kenneth Schellinger

THE terms surface tension and surface energy are well known when applied to liquids and are generally described by referring to the excess energy of the air: liquid interface as a result of unsaturated molecular forces surrounding the surface molecules of the liquid due to the presence of the air phase on one side. Such unbalanced forces produce the familiar water droplet of spherical form and are generally summed up as a surface tension measured in dynes per centimeter which can be shown mathematically to be equal numerically to a corresponding surface energy expressed in ergs per square centimeter. A specific surface energy, however, is best thought of as the energy necessary to produce one unit, of new surface on a substance. Hence, in producing a bubble in a flotation cell the impeller must supply surface energy corresponding to the air: liquid interfacial area on the interior of the bubble. Inasmuch as it is relatively easy to extend or contract the surface of a liquid, there are a number of successful methods for liquid surface tension, or energy, measurement based upon surface deformation. This happy state of affairs does not, however, extend to solids, which are considered to possess surface energies for the same reasons as do liquids, i.e. because of unsaturated ionic bonds at the solid: gas interface. As in the case of the flotation cell producing surface on liquids as new bubbles, it takes energy to produce new surface on solids as new particles. As every mill man knows, this surface is produced on mineral solids in a grinding mill by the action of a tumbling mass of iron balls. But here so much energy usually is wasted by the inefficient action of these balls that a large amount of heat is generated, and the surface energy production may be easily confused with the energy necessary to produce this ineffective heat. The tumbling balls and fracturing minerals ultimately take their energy from a rather large electric motor. It has been variously estimated that from 10 to 20 pct'" only of this energy from the motor does not appear as heat and may be presumed to appear as surface energy on the min- erals present. Such a production of new surface on the mineral phases is accompanied, of course, by a size reduction that is inevitable as more and more mineral interior molecules become surface molecules by the fracture exposure. This size reduction of mineral particles, although the most obvious feature and perhaps the sole object of the milling operation, is from this energy viewpoint the outward manifestation of the production of surface energy only. Measurement of the characteristic surface energies of pure minerals and their various mixtures in ores would be a step towards understanding of the energetics of the commercial grinding operation. In addition, the characteristic surface energy of a mineral is probably a physical property specific for that mineral, and therefore, from a scientific standpoint, should be measured. It is interesting to note that, in contrast to the large body of work on the surface tensions of liquid systems and biological systems, the field of solid surface energies has been neglected. Prior to 1920 it is difficult to find more than one or two references to work on solid surface energies in Chemical Abstracts. Since 1920 such references number somewhat less than 100, while those on liquid systems are numbered in the thousands. Much of this apparent neglect of the field of solid surface energies (the term is intended to be somewhat inclusive at this point and refers both to the solid: gas, and the solid:liquid interface) is because of the lack of a reliable method of measurement rather than any lack of scientific curiosity. It was, and still is, difficult to produce new surface on a solid without the simultaneous production of interior changes in the same solid which may consume part of the energy used. The extension in the surface area can be measured, but the interior crystal-

Jan 1, 1953

By Lawrence H. Van Vlack, Otto K. Riegger, Gerald I. Madden

The microstructural variations of periclase (MgO) in the presence of oxide liquids are examined under the steelmaking variables of: 1) temperature, 2) liquid composition, and 3) FeO additions under different oxidation levels. Attention is given to the distribution of the phases, both liquid and solid, and to the growth of individual crystalline grains. Silicate liquids penetrate more extensively between individual periclase grains than do liquids containing high percentages of Fe2O3 Higher MgO solubilities in the liquid and lower MgO contents of the solid favor more rapid grain growth. The presence of a second solid phase reduces the periclase grain growth rate and increases the amount of the solid-to-solid contact within the oxide microstructures at high temperatures. The service suitability of a refractory depends on many factors. Two are of major importance and include 1) the thermal resistance to melting, and 2) the mechanical resistance to loads at service temperatures. Neither is a simple consequence of the service temperatures because service conditions will alter compositions, produce partial melting, and induee phase changes. Consequently, the equilibrium phase relationships have been rather thoroughly studied and give a knowledge of the thermal resistance to melting, but do not give full information about the mechanical properties because two refractories with the same types and quantities of phase may have different microstructures. Although variations of microstructures with time, temperature, and composition have been subjects for extensive investigation in metals, only a limited amount of comparable microstructural work has been performed for refractory materials.' This study was an attempt to evaluate some of the consequences of service parameters upon the microstructures of refractories so that bases may be established for the analyzing of high temperature mechanical properties. Periclase (MgO) was chosen as the refractory oxide; variables included those which are encountered under steelmaking conditions such as 1) temperature, 2) liquid composition, and 3) FeO additions under various oxidation levels. Specific attention was given to the distribution of the phases, both liquid and solid, and to the growth of individual crystalline grains. The most closely related work on microstructures of polyphase materials is that of Van Vlack and Rieg-ger2 on the microstructure of magnesiowüstite [(Mg,Fe)O] in the presence of silica. In that work which pertained to solid solutions with less than 40 pct MgO, most of the quantitative work was performed on FeO microstructures. The chief conclusions concerning these relatively low-melting oxide solids were as follows: 1) the rate of crystalline grain growth is inversely proportional to the grain diameter, 2) grain growth proceeds more rapidly at higher temperatures but is slightly retarded by additional liquid content, and 3) a Silicate-containing liquid penetrates as a film between the individual magnesiowüstite grains independent of time, temperature, amount of liquid, or the MgO/ FeO ratio. The above observations are in contrast to prior work3 on the microstructure of silica in the presence of iron oxide-containing liquids where the liquid does not penetrate as a complete film between solid grains. The phase relationships for the compositions of the present work are shown in Fig. 1 which is a summary of the work of several investigators.4 Of importance is the fact that CaO forms a more stable structure with SiO2 and Al2O3 than do either MgO or FeO. The oxygen potential has little effect on periclase unless iron oxide is also present. The iron oxide is ferrous at moderately low oxygen levels, changing to ferric as the Oxygen potential is increased so the spinels, magnetite, and magesioferrite are formed.5 These two phases are relatively stable in air at steelmaking temperatures. I) EXPERIMENTAL PROCEDURE ractory were made with reagent grade Oxides. The magnesium oxide used was 99 pct MgO after ignition, and the iron oxide raw material had a minimum content of 99 pet Fe2O3. The CaO, SiO2, and A12O3 were also reagent grade raw materials. After mixing, the required compositions were pressed into pellets at a minimum pressure of 5000 psi to insure compaction of the raw materials and prevent excess void content. A silicon carbide element tube furnace was used with thermocouple control for sin-

Jan 1, 1963

By Kenneth Schellinger

THE terms surface tension and surface energy are well known when applied to liquids and are generally described by referring to the excess energy of the air: liquid interface as a result of unsaturated molecular forces surrounding the surface molecules of the liquid due to the presence of the air phase on one side. Such unbalanced forces produce the familiar water droplet of spherical form and are generally summed up as a surface tension measured in dynes per centimeter which can be shown mathematically to be equal numerically to a corresponding surface energy expressed in ergs per square centimeter. A specific surface energy, however, is best thought of as the energy necessary to produce one unit, of new surface on a substance. Hence, in producing a bubble in a flotation cell the impeller must supply surface energy corresponding to the air: liquid interfacial area on the interior of the bubble. Inasmuch as it is relatively easy to extend or contract the surface of a liquid, there are a number of successful methods for liquid surface tension, or energy, measurement based upon surface deformation. This happy state of affairs does not, however, extend to solids, which are considered to possess surface energies for the same reasons as do liquids, i.e. because of unsaturated ionic bonds at the solid: gas interface. As in the case of the flotation cell producing surface on liquids as new bubbles, it takes energy to produce new surface on solids as new particles. As every mill man knows, this surface is produced on mineral solids in a grinding mill by the action of a tumbling mass of iron balls. But here so much energy usually is wasted by the inefficient action of these balls that a large amount of heat is generated, and the surface energy production may be easily confused with the energy necessary to produce this ineffective heat. The tumbling balls and fracturing minerals ultimately take their energy from a rather large electric motor. It has been variously estimated that from 10 to 20 pct'" only of this energy from the motor does not appear as heat and may be presumed to appear as surface energy on the min- erals present. Such a production of new surface on the mineral phases is accompanied, of course, by a size reduction that is inevitable as more and more mineral interior molecules become surface molecules by the fracture exposure. This size reduction of mineral particles, although the most obvious feature and perhaps the sole object of the milling operation, is from this energy viewpoint the outward manifestation of the production of surface energy only. Measurement of the characteristic surface energies of pure minerals and their various mixtures in ores would be a step towards understanding of the energetics of the commercial grinding operation. In addition, the characteristic surface energy of a mineral is probably a physical property specific for that mineral, and therefore, from a scientific standpoint, should be measured. It is interesting to note that, in contrast to the large body of work on the surface tensions of liquid systems and biological systems, the field of solid surface energies has been neglected. Prior to 1920 it is difficult to find more than one or two references to work on solid surface energies in Chemical Abstracts. Since 1920 such references number somewhat less than 100, while those on liquid systems are numbered in the thousands. Much of this apparent neglect of the field of solid surface energies (the term is intended to be somewhat inclusive at this point and refers both to the solid: gas, and the solid:liquid interface) is because of the lack of a reliable method of measurement rather than any lack of scientific curiosity. It was, and still is, difficult to produce new surface on a solid without the simultaneous production of interior changes in the same solid which may consume part of the energy used. The extension in the surface area can be measured, but the interior crystal-

Jan 1, 1953

By R. F. Stewart, A. W. Hall

A method has been developed for continuously measuring the moisture content of coal. The method is based on the thermalization of fast neutrons by hydrogen in the coal. Neutrons from a small radio-isotope source penetrate the coal, are scattered by hydrogen, and measured by a thermal neutron detector. The number of thermal neutrons counted can be directly correlated with the moisture content of coal. In a pilot-scale system, moisture was measured continuously within 0.2% in coal moving at rates up to 20 tph. The method is adaptable in industry for continuously measuring the moisture content of coal at high tonnage flow rates. Such an application would permit continuous recording of moisture in coal without sampling and facilitate quality control. An automatic and continuous method of measuring the moisture content of coal is needed by the coal industry. Automatic control of the coal quality would reduce the cost of coal preparation, improve the product, and thus indirectly increase the use of coal. Moisture in coal can be determined by several methods, but the time required to obtain samples and analyze them by existing methods makes it difficult, if not impossible, to control the quality of the product. Both producers and consumers need a method for continuous and instantaneous measurement of moisture content without sampling in order to regulate process equipment and keep the moisture content of coal within specifications. At the Morgantown, W. Va., Coal Research Center we are developing a nuclear method for continuous measurement of moisture in coal. This method is based on the thermalization of fast neutrons by hydrogen in the water and organic matter of coal. Neutrons from a small radioisotope source penetrate the coal, are scattered by hydrogen, and are measured by a thermal neutron detector. The number of thermal neutrons counted can be directly correlated with the moisture content of coal. Design of a moisture meter based on neutron thermalization depends on many variables, any or all of which can affect the sensitivity of the meter. These factors include those related to the nuclear aspect; type and size of neutron flux and source, type of detecting device, and background count; and those related to the coal being tested: rank, particle size, and ash content. A survey was initiated to eliminate the relatively insignificant factors and to ascertain the magnitude of the major effects. Such information was necessary to fully evaluate the technique and to establish design criteria. Coal contains a relatively large amount of hydrogen in the organic coal substance and the water of hy-dration of the shaly material as well as in the moisture. To apply this concept of moisture measurement to coal requires that the organic substance in coal from any one seam of a particular mine be uniform in hydrogen content. The difference in total hydrogen content of wet and dry coal is relatively small, so that a moisture measurement based on this concept requires a measurement between two large numbers to a high degree of precision. Thus, it was necessary to develop a highly precise instrumentation system for continuous measurement and to obtain a physical arrangement permitting measurement of moving coal with a minimum effect from density variation. EXPERIMENT WITH TRAYS OF COAL Tests were conducted with metal trays containing SO to 100 lbs of coal to develop an instrument system of high precision. A scaling system with a maximum instrument error of 0.2% was used to test different types of thermal neutron detectors. The most suitable type of detector was a boron-10-lined proportional counter tube. While this type of detector showed satisfactory stability, extensive testing disclosed a long-term count reduction probably due to some type of deterioration in the detector. However, development of an electronic system using dual detectors eliminated this deterioration as a serious problem. (The second detector would be used to measure a reference drum of dry coal — the difference in count rate between the wet coal and dry reference coal being a direct measure of moisture content.) Table I, column 1, shows typical results with a 1-curie plutonium-beryllium neutron source and a thermal neutron detector beneath a tray of coal and illustrates the precision of measurement. Consecutive measurements (indicated in Table 1, columns 2-5) of thermal neutrons at various times and positions be-

Jan 1, 1968

By M. R. Tek

A new method for correlating the data on multiphase flow through vertical pipe is presented. The correlation is based on a "two-phase f factor" concept which was developed and successfully applied to horizontal multiphase flow by Bertuzzi, Tek and Poettmann. Field data previously published on several flowing and gas-lift wells have been used as the basis for the developed correlation. The application of the method to actual two-phase flow problems indicate that this method is capable of predicting the pressure distribution in vertical multiphase flow strings well within the accuracy range usually desired for common engineering and design calculations. A new working chart developed for calculation of two-phase pressure gradients and a graphical step-by-step procedure for the computation of pressure distribution are presented along with an example problem. INTRODUCTION Multiphase flow through vertical pipe is encountered in many engineering installations. In petroleum, chemical process, nuclear engineering and many other industries, problems associated with simultaneous flow of two or more phases through vertical pipe have been of interest for a long time. This interest has increased considerably during recent years due to applications to new processes in petroleum production and refining and to problems of steam generation and heat removal from nuclear reactors. One prominent example of vertical two-phase flow is provided by the gas-lift process where oil, water and gas flow simultaneously. If the pressure profile in a gas-lift well can be predicted within reasonable accuracy, it would be possible to get good estimates of the power required to lift the oil, the optimum depth, the pressure and the rate at which to inject gas. Furthermore, the effect of production rate and tubing sizes on these quantities can be evaluated before any design decision is made on the installation and operation of the flow string. The majority of experimental work available in the literature deals with two-component systems where individual phase flow rates in and out of the pipe remain constant. The general problem of prediction of pressure drops in multiphase flow systems is very complicated. The co-existence of numerous flow patterns of widely different geometry and mechanism, conditions of surface instability and the nature of force fields acting upon the system are among the major difficulties commonly encountered. The classical approach of fluid dynamics which would be based upon the formulation and solution of Navier-Stokes equations has been found by many investigators completely devoid of any hope of success-—not only because of inherent nonlinearities but also because of insurmountable analytical difficulties standirng in the way of setting up the boundary conditions. The presence and effect of interfacial forces on mu1tiphase flow systems further complicate the theoretical approach. For these reasons, many investigators choose to adopt semi-empirical if not purely empirical approaches in order to obtain solutions of engineering utility. A prominent practical solution has been given by Poettmann and Carpenter' in the form of an empirical correlation. In their paper, total flowing densities of fluids and solubility effects of gas in oil have been taken into account for the correlation of field data covering a wide range of operating conditions for oil wells. They treated the gas. oil and water as a single phase of combined propcrties and correlated the multiphase friction factor as a function of the product (pvd), i.e., density X velocity X diameter of the flow string. However, because the product (pvd) is dimensional while the frictior factor is dimensionless, the generality of their result's is somewhat restricted. It seems that the omission 01 viscosity effect may be one of the reasons for the scatter of data as shown in their correlation chart. The concepts of the "two-phase f factor" and the "two-phase Reynolds' number function" were recently developed and successfully applied to correlate horizontal multphase flow by Bertuzzi, Tek and Poettmann.' Recently, two new methods of correlation by Tek and Chan' have been presented on simultaneous flow of liquid and gas through vertical pipe. These two methods of correlation and the working equations in field units necessary for their application are included in this paper. The extension of the conccpt of two-phasc Reynolds' number function successfully developed for horizontal flow into vertical multiphase flow systems, development and evaluation of working charts permitting calculation of two-phase pressure gradients and

By J. P. Denny

The constitution of the Ga-In system was determined by thermal methods. An experimentally determined metastable equilibrium line (an extension of the indium-rich liquidus) was obtained. The various alloys were studied metallographically using polished samples obtained by a casting method. These low melting alloys required a special dry-ice assembly to maintain a suitable temperature. RECENT interest in alloys that are liquid at room temperature has led to rather extensive investigation of gallium-base alloys. Widely distributed over the earth, gallium could be produced in substantially larger quantities than at present, if a significant demand existed.' One study' has established its presence in 12 out of 14 zinc blends, in all of 15 aluminum ores, in 4 out of 12 manganese ores, in 35 out of 91 iron ores, and in all of 7 magnetite ores. It occurs as a rule in minute amounts, however, leading to high extraction costs. Recent quotations run from $2.50 to $7.50 per g. During the course of the present investigation, portions of the system Ga-In have been redeter-mined, and the results of this study are presented herein. Thermal and metallographic methods have been employed. Lecoq de Boisbaudran, the discoverer of gallium, conducted the first investigation" on Ga-In alloys in 1885. The temperatures of incipient melting, and of completion of melting, were determined at four alloy compositions. In 1936, Hansen' constructed a eutec-tic-type phase diagram for the system Ga-In, based on . work. The existence of a Ga-In compound was regarded as improbable by Hansen, and subsequent investigations are in agreement. French, Saunders, and Ingle3 conducted a more complete study of the system in 1938, using thermal methods. Their phase diagram is a eutectic type, containing a unique concave-upward liquidus. The solid-solution range of gallium in indium was reported as 9.5 pct by weight, and that of indium as less than 1 pct, at the eutectic temperature. The eutectic composition, determined as being bracketed by the compositions showing a true horizontal at the eutectic temperature (16°C), was reported as 76 5-0.5 pct Ga and 24 i 0.5 pct In. Experimental Procedure The preparation of Ga-In alloys is simplified by the low melting points involved. Various compositions were prepared by melting in pyrex tubes, using a cover of distilled water or parafin to prevent the alloys from wetting the glass wall. In all cases, the melts were homogenized. by stirring. Where possible, both cooling and melting curves were determined. The extensive undercooling of gallium was found to prohibit a satisfactory cooling-curve analysis of gallium-rich alloys, however, and transition points on the gallium side of the eutectic could be determined only by melting curves. The inverse rate method of thermal analysis proved to be most satisfactory and was used to a great extent. Various heating and cooling rates were used, ranging from 0.2" to 5.0°C per min. Low temperature melting analyses were conducted within a constant temperature bath, maintained at about 70 °C. The alloys were solidified (under water or paraffin) within a pyrex tube, using dry ice; the tube was then sealed within a cold Dewar flask, the unit transferred into the constant temperature bath, and periodic temperature readings taken. The high temperature melting-curve determinations and all cooling-curve determinations were made in a vertical tube furnace. At near-eutectic compositions, the furnace was placed within a refrigerated room held at —20°C. Accordingly, the furnace on cooling approached —20°C asymptotically and permitted the determination of those phase transitions occurring below room temperatures. Temperatures were measured with a 30-gage iron-constantan thermocouple, immersed directly in the alloy. To prevent contamination of the melt, the leads and junction were coated with Lucite, applied by painting with a solution of Lucite in ethylene dichloride. Electromotive force measurements were made with a Leeds and Northrup precision potentiometer, type 8662. The couples were calibrated against the boiling point. of water and, at lower temperatures, against a calorimeter thermometer having a Bureau of Standards certificate. The melting points of gallium and indium used in the present investigation were determined as 29.-77° and 156.1°C, in good agreement with previously reported values of 29.78°C° and 156.4"C.' The spec-

Jan 1, 1953

By J. D. Testerman

A statistical technique to identify and describe naturally occurring zones in a reservoir and to correlate these zones from well to well is described. The technique is particularly useful in describing a reservoir where cross flow between adjacent strata is important in determining reservoir behavior. Although it has been used primarily for permeability zonation, the technique is general and can be used to correlate any reservoir property or related data, such as the information contained in well logs. INTRODUCTION One of the first problems encountered by the reservoir engineer in predicting or interpreting fluid displacement behavior during secondary recovery processes is that of organizing and using the large amount of data available from core analysis. Permeabilities pose particular problems in organization because they usually vary by more than an order of magnitude between different strata. Due to the sheer volume, it is almost always necessary to group data and to use an average value to represent a number of measurements. Perhaps the most common method now used to group permeability data is the capacity-fraction technique, which ranks permeabilities in order of magnitude, regardless of the physical location of the permeabilities within the reservoir. The cumulative per cent capacity is plotted against cumulative per cent thickness. This plot is divided into an arbitrary number of zones, generally of equal thickness. Five zones (or averaged groups of data) usually are obtained, each of which is treated as homogeneous in subsequent calculations. The division so obtained has no physical meaning; strata in the same zone, calculation-wise, are usually not adjacent in the reservoir. Reservoir engineering techniques being developed will handle crossflow that occurs between adjacent communicating reservoir strata because of imbibition and gravity segregation. Since crossflow occurs between physically adjacent layers within the reservoir, a new zonation technique recognizing the actual location of strata within the reservoir is necessary. Similarly, the recognition of natural zones is important for predictions of oil recovery by processes involving diffusion. One such process is miscible displacement, where predictions of lateral diffusion within the reservoir must recognize the actual location of the invaded zones in relation to the rest of the formation. Natural zones must also be adequately recognized to account for heat transfer within the reservoir during thermal exploitation. Because of the complexity of the problem, statistics appear to offer the only practical hope of dividing a reservoir into physically-meaningful natural zones. This paper presents a statistical technique for identifying these natural zones and for ascertaining which ones are likely to be continuous between adjacent wells. The zones defined have minimum variation of permeability internally and a maximum variation between zones. The technique is general and can thus be applied to reservoir properties other than permeability. The method will guide the reservoir engineer in estimating which zones are likely to be continuous between wells. However, a statistical correlation based on permeabilities in two different wells is no guarantee that the zones so defined are, in fact, continuous. Rather, the assumption of continuity must be consistent with geological data concerning the depositional environment, as well as justified on the basis of engineering judgment in combination with statistics, just as judgment is required with conventional zonation methods. CALCULATION PROCEDURE The reservoir zonation technique is a two-step operation. The steps are individually described, and a sample calculation is presented in the Appendix. ZONATION OF INDIVIDUAL WELLS First, the set of permeability data at a single well is zone, into Zones. These zones are selected so that variation is minimized within the zones and maximized between the zones. The equations4,6 used to zone the data are where B = the variance between zones, , = the number of zones, i = the summation index for the number of zones, j = the Sumation index for the number of data within the zone, mi = the number of data in the ith zone, k,. = the mean of the permeability data in the ith zone, k . = the over-all mean of the data in the well, W = the pooled variance within zones, N = the total number of

By M. Garfinkle

LITERATURE data1'3 have suggested that the resolved shear stress necessary to initiate slip in a particular slip system in tungsten single crystals is strongly dependent on the tensile orientation of the crystals. This orientation dependence is contrary to the critical resolved shear stress criterion. Because the Schmid factors for slip on particular {110}( I l l), {112}(Ill), or {123}(111) slip systems in a (100) oriented tungsten crystal are virtually identical to the respective Schmid factors for a (110) oriented crystal, slip should initiate on the same system in both orientations, independent of the values of the critical resolved shear stresses. Thus, it would be expected that, for a critical resolved shear stress model to hold, the proportional-limit stress for these two orientations must be similar. The expectation of similar proportional -limit stresses has not been borne out by reported experimental observations. Rose et nl.' reported that the proportional-limit stress for crystals near the (110) orientation was 92,000 psi, approximately three times that of (100) oriented crystals (31,000 psi). Garlick and probst2 reported greater than a fourfold difference. Most recently, Beardmore and ~ull~ examined tungsten crystals with orientations along the edge of the standard stereographic triangle and showed a gradual increase in proportional-limit stress from 28,000 psi for the (100) orientation to 93,000 psi for the (110) orientation. Again, the difference was more than threefold. As has been emphasized by these investigators, the apparent gross variation of proportional-limit stress with orientation precludes any meaningful analysis of slip mechanics using a simple critical resolved shear stress relation. The purpose of the present communication is to report that, by use of refined strain-measuring techniques, the proportional-limit stresses for the (100) and (110) orientations are quite similar. Obviously, to determine accurately the value of the proportional-limit stress, which is by definition the stress at which nonlinear stress-strain behavior first occurs, the method of measuring strain must be capable of high accuracy. The investigators cited used as the measured elongation the crosshead displacement of the tensile machine. This method, however, introduces extraneous strain contributions from the load train. Because these extraneous contributions are considerable in the initial portions of loading (as has been well-documented in the literature) the proportional-limit stress measurement, dependent on the deviation of a line from linearity, can be significantly affected. During the course of a recent study, room-temperature strain was measured directly on several (100)-and (110)-oriented tensile specimens axially oriented within 5 deg of the desired direction. Each specimen was fitted with two foil strain gages mounted on opposite sides of the 0.08-in.-diam reduced section. The tensile specimens were loaded at a crosshead speed of 0.01 in. per min on an Instron Universal Testing Machine. Strain was monitored through an external bridge circuit in which the specimen strain gages comprised two arms. The bridge output was recorded against the applied load. A strain of 2 x 10"" could be readily detected by this method. Some specimens were subjected to conventional loading while the load for others was cycled between zero and continually increasing values until the first deviation from linearity was observed. The results of these methods were similar, and the stress at which nonlinear behavior occurred was considered the proportional-limit stress. The proportional-limit stresses measured in this study are approximately 5000 and 6500 psi for the ( 100) - and ( 110) -oriented crystals, respectively, a range of values almost a magnitude less than those reported previously. Representative curves are presented in Fig. 1. Because of the nature of the measurements made, the point of departure of a line from linearity is often debatable. While the load-strain curve for the (100) -oriented crystals made a relatively abrupt deviation from linearity, the curve for the (110) -oriented crystal was more gradual, so that the actual point of tangency may be at a lower load than reported. Thus, the difference in values for the proportional-

Jan 1, 1967

By J. W. Snavely

A practical mathematical treatment is presented for the determination and control of conveyor belt acceleration, particularly for conditions of starting where vertical curves are involved. A typical sample problem is analyzed, with required calculations, to clarify the procedure. THE part that acceleration plays in starting a belt conveyor and its effect on belt conveyor design are well understood in a general way. Its practical importance is easily overlooked, however, and under some conditions, it is absolutely necessary to give the problem of acceleration detailed study. Most handbooks on conveyor belting design adequately present basic data for the determination of acceleration values. This paper will only attempt to present practical thinking and a convenient method of treatment of acceleration in belt conveyor design. Mathematical Analysis In working out the various problems of conveyor belt acceleration, the starting point, as presented by the handbooks, is the familiar formula of "force of acceleration is equal to the mass times acceleration." By expressing these fundamental quantities in terms of belt conveyor design, it is possible to arrive at the unsuspected conclusion that the acceleration time for horizontal belt conveyors is independent of the load, and instead, dependent upon the belt speed, the type of drive arrangement and drive pulley, and the idler coefficient of friction. The mathematics leading to this conclusion are shown in Table I, which has been prepared to show this derivation. While at first the conclusion just given may not seem to be reasonable, further reflection indicates that obviously the type of drive pulley and the type of drive do affect materially the tension in the conveyor belt, and thus, as clearly shown, the time of acceleration is dependent upon the factors mentioned. Inasmuch as all of the factors except time are predetermined by the belt conveyor design, it becomes relatively easy to establish the accelerating time and to reduce further this time determination to a simple graph from which the time in seconds can be read directly. Such a graph is given in Fig. 1. The table appearing on Fig. 1 should be explained further. For a given belt speed, the time of acceleration can be expressed as a percentage of the belt speed. The time of acceleration is also dependent on the drive arrangement, and changes in the drive arrangement consequently change the time of acceleration. It further follows that for a given belt speed, the time expressed as a percentage of that belt speed also changes with the type of drive. Obviously then, it becomes possible to graph the percentage of speed for .each type of drive against the belt speed and accelerating time, after which, for a given belt speed and type of drive, the time can be read directly in seconds. Two constants were established for Fig. 1, the first one being the limiting of the maximum acceleration tension to 35 pct of the full load operating tension in the belt. The purpose of this is to limit the total tension imposed upon the belt during the acceleration period to 135 pct of the full load operating tension, 'which is the amount required to start or breakaway the fully loaded belt conveyor from rest. The other constant is the friction factor used for the idler equipment, which has been established as 0.022. For installations where it is necessary to establish the values of acceleration, invariably high grade idler equipment is used, and it has been established from field experience that 0.022 for the idler friction factor is amply conservative. The use of this friction factor for idlers must be tempered with judgment, of course, for occasions will arise where more power than indicated is required to start, even with the very best of equipment, such as low temperature operations that tend to congeal the grease in the bearings and thus produce additional friction drag. An inspection of the table in Fig. 1 affords a convenient rule of thumb method for determining the acceleration time, which conveniently can be 5 pct of the belt speed in seconds. The 5 pct of belt speed figure is close to the average for most types of drives. In using Fig. 1 it must be emphasized that it applies accurately to horizontal belt conveyors only.

Jan 1, 1953

By J. W. Snavely

A practical mathematical treatment is presented for the determination and control of conveyor belt acceleration, particularly for conditions of starting where vertical curves are involved. A typical sample problem is analyzed, with required calculations, to clarify the procedure. THE part that acceleration plays in starting a belt conveyor and its effect on belt conveyor design are well understood in a general way. Its practical importance is easily overlooked, however, and under some conditions, it is absolutely necessary to give the problem of acceleration detailed study. Most handbooks on conveyor belting design adequately present basic data for the determination of acceleration values. This paper will only attempt to present practical thinking and a convenient method of treatment of acceleration in belt conveyor design. Mathematical Analysis In working out the various problems of conveyor belt acceleration, the starting point, as presented by the handbooks, is the familiar formula of "force of acceleration is equal to the mass times acceleration." By expressing these fundamental quantities in terms of belt conveyor design, it is possible to arrive at the unsuspected conclusion that the acceleration time for horizontal belt conveyors is independent of the load, and instead, dependent upon the belt speed, the type of drive arrangement and drive pulley, and the idler coefficient of friction. The mathematics leading to this conclusion are shown in Table I, which has been prepared to show this derivation. While at first the conclusion just given may not seem to be reasonable, further reflection indicates that obviously the type of drive pulley and the type of drive do affect materially the tension in the conveyor belt, and thus, as clearly shown, the time of acceleration is dependent upon the factors mentioned. Inasmuch as all of the factors except time are predetermined by the belt conveyor design, it becomes relatively easy to establish the accelerating time and to reduce further this time determination to a simple graph from which the time in seconds can be read directly. Such a graph is given in Fig. 1. The table appearing on Fig. 1 should be explained further. For a given belt speed, the time of acceleration can be expressed as a percentage of the belt speed. The time of acceleration is also dependent on the drive arrangement, and changes in the drive arrangement consequently change the time of acceleration. It further follows that for a given belt speed, the time expressed as a percentage of that belt speed also changes with the type of drive. Obviously then, it becomes possible to graph the percentage of speed for .each type of drive against the belt speed and accelerating time, after which, for a given belt speed and type of drive, the time can be read directly in seconds. Two constants were established for Fig. 1, the first one being the limiting of the maximum acceleration tension to 35 pct of the full load operating tension in the belt. The purpose of this is to limit the total tension imposed upon the belt during the acceleration period to 135 pct of the full load operating tension, 'which is the amount required to start or breakaway the fully loaded belt conveyor from rest. The other constant is the friction factor used for the idler equipment, which has been established as 0.022. For installations where it is necessary to establish the values of acceleration, invariably high grade idler equipment is used, and it has been established from field experience that 0.022 for the idler friction factor is amply conservative. The use of this friction factor for idlers must be tempered with judgment, of course, for occasions will arise where more power than indicated is required to start, even with the very best of equipment, such as low temperature operations that tend to congeal the grease in the bearings and thus produce additional friction drag. An inspection of the table in Fig. 1 affords a convenient rule of thumb method for determining the acceleration time, which conveniently can be 5 pct of the belt speed in seconds. The 5 pct of belt speed figure is close to the average for most types of drives. In using Fig. 1 it must be emphasized that it applies accurately to horizontal belt conveyors only.

Jan 1, 1953

By P. Majani, F. P. Kokesh, M. Grosmangin

Determination of the quality of cementation of casing in oil wells in the past has involved inflow and circulation tests to insure that the producing zones were adequately sealed off from the adjacent zones. Existing logging methods, such as temperature and radioactivity surveys, may detect the presence of cement behind the casing. However, the qualities of the cement (i.e., its hardness and particularly its bond to the casing) are not indicated. The new logging method described in this paper operates on the principle that the attenuation of a sonic pulse transmitted by a casing is greatly increased when that casing is bonded to an outer annulus of hard material (such as set cement) which has an appreciably smaller sonic-wave velocity than that of the casing. The down-hole tool contains a source of recurrent sound pulses which are detected by a receiver spaced a few feet from the source The amplitude of the detected casing-borne pulse is measured, and the resulting signal is transmitted to the surface where it is recorded vs depth. Since amplitude is a function of attenuation, the log is readily interpreted. Laboratory studies have shown straightforward relationship between attenuation and such variables as source-detector spacing and per cent of circumference bonded. It is shown that cement not set or not bonded to the casing has compara- tively little attenuating effect. Field examples show not only the cement top, but also the variation in cementation quality below the top. Further, the increase of bonding with time and after squeeze cementation is depicted. The detection of poor cement jobs is confirmed by production tests and by formation-test results. It is anticipated that the method will have wide application in evaluating cementation quality prior to formation testing in completions and recompletions. The analysis it affords may aid in further improving cementation techniques. INTRODUCTION The main purpose of oilwell cementing is to isolate a production zone from other undesirable zones. To investigate whether this purpose has been accomplished, several logging methods have been used such as temperature logs, radioactive tracer logs, etc. While these logs all respond to the presence of cement behind the casing, they do not indicate the degree of bonding of the cement to the casing. Early in the application of sonic logging, it was noticed that considerable attenuation of sound signals takes place in cemented pipe and is often made evident on the standard Sonic log by cycle skipping.' The development of a circuit capable of continuously recording the amplitude of the casing-borne sound signal has made possible an extensive series of laboratory and field tests, which gave the following results. The amplitude of a sound signal after it has traveled in a firmly cemented pipe is only a small fraction of that recorded by the same device in free pipe. This provides a wide spectrum of energy levels; for given local conditions, empirical values of the amplitude can be correlated with the quality of cementation. Interpretations made in this manner generally have been confirmed by production tests, circulation tests and squeeze cementation. The purposes of this paper are to give a general description of the new logging method and to present some laboratory and field results. CYCLE SKIPPING Early attempts to study the quality of the cement behind the casing were performed with a standard Sonic log, which measures the transit time At, and were based on the well known phenomenon of cycle skipping.' Cycle skipping normally is interpreted as being a manifestation of weak signals at the receivers. The log of Fig. 1 was run with the recording instrumentation adjusted to enhance cycle skipping. A mirror-image presentation was used for better visual interpretation. The transit time of sound in steel is about 58 microseconds/ft (corresponding to a velocity of 17,000 ft/ sec), and the portions of the log where this value is recorded are interpreted as zones with no cement bond. Where cycle skipping produces a higher value of At, weak signals (or a high rate of attenuation) are indicated, and a good cement bond may be present. The occurrence of cycle skipping, however, depends too much on instrument adjustment to give a uni-

By W. W. Stephen, J. P. Walsted, W. J. Kroll

FOR the production of carbides, various furnace types are available, especially those using arc, resistance, and high-frequency heating. Selection of a specific means of heating depends primarily on the material to be treated and the physical properties of the carbide produced. In the present case, zirconium carbide had to be prepared on an industrial scale as a raw material for the production of anhydrous zirconium chloride. Considering that a rather expensive pure oxide was to be used, the arc-furnace treatment recommended for zircon sands in a previous publication' was ruled out because of the considerable volatilization and dust losses caused by the blast of the arc. For this reason, either high-frequency or resistance heating seemed to offer more promise. Since there was not enough capacity of the former available, resistance heating was chosen. It was first thought that the Acheson silicon carbide furnace would be suitable for the present purpose, but the voltage in such a furnace, in which the current passes through the batch, varies from 220 to 75 v from the start to the end of a run. This variation is so great that a special tap transformer would have been required. Trouble was also expected by local melting of the carbide. Pure zirconium carbide melts at about 3527°C, but the melting point is brought down to 2427°C, according to Agte,2 when an excess of 6 pct C is present in the carbide. This we found confirmed by experiments in a high-frequency furnace. Excess carbon is needed in the batch to obtain a complete reduction. Fusion of the charge would cause great difficulties in an Acheson-type furnace because of the good electrical conductivity of the carbide as compared with that of the loose batch. Also, fused carbide is much more difficult to chlorinate than the spongy product that can be made in the radiation furnace described below. It was apparent that, to obtain a good-qual-ity zirconium carbide, the heat input would have to be well-controlled. The hairpin-resistor principle seemed to offer possibilities in this regard, and a furnace of this type was therefore developed. The advantages of the hairpin-resistor radiation principle have been discussed in previous publications, and a split-tube graphite-resistor furnace," now increasingly used in various laboratories, as well as a centrifugal quartz melting furnace4 of this type, has demonstrated the usefulness of this heating method. The hairpin-heater element has the following definite advantages over a straight resistor of the type used, for instance, by Georges:5 Its resistance is four times greater; it can expand freely; it is sturdier because of the larger diameter, and it has a larger radiation surface; there are no hot contacts that might wear out or overheat; only one clamp is used which permits assembling all electrical leads at one side of the furnace, making the other sides easily accessible to the operator. The shorter element and its larger diameter permit greater concentration of heat. The furnace developed is shown in fig. 1. The box (I), made of 2 1/2-in. graphite plates, has inside dimensions of 23x17x16 in. It contains the briquet-ted batch (2). The box is embedded in lampblack (3) up to the cover plate. The cover plate contains an opening for the gas escape (4) and for the observation hole (5), which permits measuring the temperature with an optical pyrometer. The cover plate is embedded in charcoal (6). The lampblack is contained in the insulating brick lining (7), held in the 1/4-in. sheet steel box (8). The graphite box is set on two rows of triangular graphite bars (9). The hairpin-heater element (10), the dimensions of which are given below the main drawing, extends horizontally in the graphite chamber and radiates freely on the batch. A graphite tube (11) keeps the lampblack from falling into the slot. The split electrode, which in reality is turned 90" against the drawing, is so arranged that the slot is vertical. The water-cooled packing gland (12) is insulated by an airgap from the heater element. A thin pipe (13)

Jan 1, 1951

By M. R. Tek

General equations relating the pressure drop necessary to sustain the flow of a fluid through a porous matrix at a given rate have been developed. The results indicate that at high values of flow rate the pressure-flow behavior may not necessarily satisfy the usual Darcy equation. The mathematical analysis, carried through the micro-pore geometry and extended through the macro-reservoir scale, indicate that Darcy's law, of limited applicability to certain ranges of Reynolds numbers, can be generalized through the inclusion of some additional parameters. The "generalized Darcy equation" has also been formulated in dimen-sionless form permitting the evaluation of its predictive accuracy with regard to literature data. A comparison between predicted and experimental values indicates that the generalized Darcy equation predicts the pressure drops with good agreement over all possible ranges of Reynolds numbers. INTRODUCTION The limits and the nature of validity of Darcy's law' has been a subject of every-day interest to the industry for many years. It is well known that as the Reynolds number, characteristic of the fluid flow through porous media, becomes large, Darcy's law gradually loses its predictive accuracy and ultimately becomes completely void. For the last 20 years much has been said and written on this subject. Unfortunately little has been accomplished to bring about a satisfactory agreement, at least on the nature of the threshold of validity of Darcy.'s law. Fluid dynamists, geo-physicists, and engineers all had their individual views, explanations, interpretations and concepts on the subject. To some, a mechanistic analogy with pipe-flow proved a satisfactory explanation.' To others,' turbulence, in its random character, was incompatible with the geometric structure of consolidated porous systems. To some,4 turbulence merely represented a factor influencing the permeability measurements and again to others5,6,7 em-pirical or semi-empirical correlations proved satisfactory from an engineering viewpoint. Deviations from Darcy's law at high flow rates have been studied by systematic experiments by Fancher, Lewis, and Barnes.' In an article on the flow of gases through porous metals, Green and Duwezs conclude that the onset of turbulence within the pores appears unsatisfactory to explain deviations from Darcy's law. This view is held by many others. While the subject remained controversial for many years, the development of vast natural gas reserves throughout recent years further justified considerable interest on this problem from the standpoint of gas reservoir behavior. As large amounts of field data became available from the operation of many gas fields, it became evident that the steady-state behavior of gas wells was not, in general, in agreement or compatible with Darcy's law. This suggested a careful reconsideration of all mechanisms which may account for pressure drops in addition to viscous shear. In a series of articles9,10 . Hou-peurt indicated that deviations from Darcy's law may be explained on the basis of kinetic energy variations and jetting effects without resorting to assumptions on turbulent flow conditions. Another article by Schneebeli11 indicates that special experiments by Lindquist clearly demonstrated that the onset of turbulence does not necessarily coincide with conditions of deviation from Darcy's law. This view is also held by M. King Hubbert.12 Starting with the basic pressure-flow relations suggested by Houpeurt, the derivation, development and extension of analytical expressions to -supplement and generalize Darcy's law has been the objective of this work. MATHEMATICAL ANALYSIS Derivation of Dimensionless Pressure-drop, Flow-rate Relations In considering the flow of a fluid through a porous matrix geometrically represented by a succession of capillary passages in the shape of truncated cones,810 an approximate expression may be derived relating viscous and inertial, i.e., total pressure drop to the physical properties of the fluid, geometric properties of the rock matrix and the rate of flow: ?P/?r = µ/k V [ 1 + c(m4 - 1) p V/16n" mµ w] ..........(1) Let us formally set: c (m4 - 1) / 16n" m = a d ......(2) Such a representation is equivalent to assert that the term [c(m4 — 1)/ 16n"m], variable with various porous media and probably highly variable within a given porous medium, may be macroscopically defined as equal to a lithology factor times the aver-age grain diameter d. In view of the usual grain and pore size distribu-

Jan 1, 1958

By T&apos Ke, ing-sui

NUMEROUS investigators have observed that internal friction accompanies cold-working of metals and the effect of annealing is to reduce this internal friction.1,2 However, - most of the experiments were made at high stress amplitudes and the principal purpose was to study the increase of internal friction as a result of the applied cyclic stress during measurement. In order to study the internal friction introduced by cold-working applied prior to the measurement, the stress level applied during the measurement of internal friction must be sufficiently small. The results of measurement are significant and can be used for a base of comparison only when the applied cyclic stress is so small that the internal friction is independent of stress amplitude. Internal friction of cold-worked metals under small stress level has been studied by a number of workers8 -" he internal friction was measured around room temperature with a frequency of vibration of the order of kilocycles per second. The purpose of this paper is to report a study of the change of internal friction when severely cold-worked aluminum was annealed at successively higher temperatures until it was completely recrys-tallized. The measurements of internal friction were made over a range of temperature extending from room temperature up to the temperature of prior anneal. The frequency of vibrations used was about one cycle per second. The apparatus used for the internal friction measurements to be reported in this paper was a torsion pendulum with the specimen in wire form as the suspension fiber. The description of this apparatus and the method of measurement have been previously given.7,8 The applied stress was sufficiently small SO that the magnitude of internal friction is independent of stress level at each temperature range concerned. Corresponding to this stress the maximum shearing strain on the surface of the specimen is of the order of l0-5 and lower. The in- ternal friction (Q-1) is reported as 1/p times the logarithmic decrement. Internal Friction Versus Temperature of Anneal: Fig. 1 shows the internal friction measurements performed upon 99.991 pct aluminum subjected to 95 pct reduction in area. The final diameter of the wire is 0.033 in. This figure gives a general survey of the effect of temperature of anneal and of temperature of measurement. The internal friction of the cold-worked specimen was first measured at room temperature. It was then annealed at 50°C for one hour and the internal friction measured at 50°C and at room temperature. The same wire was successively annealed at higher temperatures for one hour and measurements were taken at the annealing temperatures and lower temperatures as before. Such a procedure was followed in order to stabilize the internal friction at the temperature of measurement so that during the measurement which generally takes about half a minute, there is no detectable change in internal friction. This series of measurements .was made up to 450°C. After each annealing a short test piece of the specimen, which had received the same past thermal and mechanical treatments, was taken out for metallographic examination. It is seen from fig. 1 that up to the annealing temperature of 250°C we have the following observations: for any given temperature of measurement, the internal friction is lower the higher the temperature of prior anneal. When the annealing temperature is 290°C or higher, the internal friction at the annealing temperature drops abruptly to a value which is much smaller than that for the previous curve. Metallographic examinations showed that the recrystallization of the specimen was completed after the annealing at 290°C. Fig. 1 shows that, as far as internal friction is concerned, there is no abrupt transition between the processes of recovery and recrystallization. Averbach has also reached the conclusion that recovery may be a process analogous to recrystallization on the basis of X ray extinction measurements in brass." The effect of annealing temperature upon the internal friction at room temperature is shown by curve I of fig. 2. In this figure the internal friction at room temperature was plotted as a function of annealing temperature. It is seen that the internal friction decreases rapidly at first with an increase

Jan 1, 1951

By A. V. Bradshaw, R. I. L. Guthrie

The behavior of large bubbles in the size range 4 to 25 cm3, rising through molten silver, has been studied. It was found that rising velocities were equivalent to those in aqueous systems of low viscosity. Mass transfer coefficients for oxygen bubbles dissolving in silver were found to be 0.036 ± 0.007 cm sec-1, being close to those predicted for transfer through the front surface of the spherical cap bubble only. It is suggested that the surface active nature of oxygen in silver could account for the relatively low coefficients obtained. MANY metallurgical processes involve interactions between gas bubbles and liquids. Examples include the removal of carbon monoxide in Open Hearth Steelmak-ing, the removal of sulfur by blowing air through copper matte during converting, and the removal of hydrogen from steel during vacuum degassing or inert gas flushing. The steps involved in such refining processes include; transport of the dissolved species to the bubble interface, adsorption and chemical reaction of the species at the interface, desorption of product molecules from the interface, and transport of product gas into the bulk gas phase of the bubble. It has been concluded1 that all the interfacial steps involved proceed so rapidly at steelmaking temperatures that transport of the solutes, present in the metal, become the important rate controlling factors provided nucleation phenomena are not restrictive. The O-Ag system was chosen for the investigation into gas bubble-molten metal interactions due to the relatively high solubility of oxygen that enables rates of oxygen transfer to be measured from changes in bubble volume. Other advantages of this system include the absence of a stable oxide phase at an oxygen pressure of 1 atm and the relatively low melting point of the metal which permits the use of a metallic container, providing that it is resistant to oxidation. In those metallurgical processes where bubbles have an important influence, bubble volumes are usually greater than 5 cm3. For this reason the present study relates specifically to single large bubbles of oxygen rising in silver. These bubbles adopt the characteristic spherical cap shape similar to that shown in Fig. 1 for a 30 cc bubble rising in water. After an initial investigation to determine the velocities of inert (nitrogen) bubbles rising in molten silver, experiments were carried out with oxygen and the rates of mass transfer between the oxygen bubbles and the silver were measured. EXPERIMENTAL Apparatus. The apparatus, Fig. 2, for containing molten silver, was constructed from "Nimonic 75" Alloy (75 pet Ni, 20 pet Cr, 5 pet Fe, Mn) and provided for the release of single bubbles from an hemispherical cup, situated at the bottom of the column. The cup was turned by translating the rotation of the drive shaft through 90 deg. This was accomplished by the use of a bevelled gear system, and a smooth drive was provided by the lubricating action of the silver on the gears. Since reliable high temperature seals at 1000°C were found to be impracticable, the filling and drive shaft tubes were extended outside the 3.5 kw resistance wire tube furnace, where connections were made using easily accessible O-ring seals. The apparatus remained gas tight to the atmosphere at pressure differentials far in excess of those used. The filling tube was connected via a small bore tube to the differential pressure transducer. Gas could be bubbled into the inverted cup from two i-in. tubes which passed down the inside of the column to preheat the gas. The temperature of the silver was maintained at 1020°C during all experiments. Measurement of Bubble Volume. In order to calculate mass transfer rates, it was necessary to obtain a continuous record of the bubble's volume during its passage through the column of molten silver. The method adopted for measuring the bubble volume involved closing off the top gas space to the atmosphere prior to each experiment, and recording the variation in gage pressure of this space during the formation and rise of the bubble. Since any change in bubble volume results in an equal change in top space volume, Boyles Gas Law may be applied (for isothermal con-

Jan 1, 1970

By L. Stone, H. Margolin

The Ti-C-N and Ti-C-O systems were investigated in the temperature range from 500° to 1400°C and in the composition range up to 2 pct C and 5 pct N or 0. Characteristic isothermal sections at 800°, 900°, 1000°, and 1300°C are presented. The Ti-N-0 system was studied in the temperature range from 900' to 1400°C with alloys containing up to 6 pct total alloying content. Characteristic isothermal sections at 1000° and 140O°C are presented. Melting-point data for all three systems are also included. THIS paper reports on one of a series of investi-gations which have been conducted on the phase diagrams resulting from interstitial alloying with iodide titanium. The other investigations involved delineation of the binary systems with carbon,' nitrogen and boron,h and oxygen. The Ti-0 binary system has also been investigated by Bumps et al.' In varying degrees, each of these interstitial elements has been shown to stabilize the low temperature a modification of titanium1-5 and each forms a face-centered cubic TiX compound (henceforth designated 6). In addition, the Ti-N and Ti-0 systems reveal a low temperature tetragonal phase (6) formed by a peritectoid reaction between a and TiX Experimental Procedure The development of experimental techniques for the study of titanium alloy systems has, to a large extent, become standardized. In this investigation, the equipment and procedures described in detail by Cadoff and Nielsenl have been used. Arc Melting: In general, binary alloys with carbon, nitrogen, and oxygen, prepared in the composition range of interest in this investigation, show negligible composition changes during arc melting. However, the possibility of the formation of some gaseous combination of alloying elements such as CO, CN, or NO during the preparation of these ternary alloys was considered. Calculations showed that the evolution of only 0.05 gram of such a gas would be detectable as a pressure change in the closed system used during preparation of these alloys. Such pressure changes were not observed. Consequently, nominal compositions have been used in plotting the data. The compositions of the materials used in the preparation of the alloys are shown in Table I. After melting for 3 to 5 min at 275 to 350 amp, the alloys were checked for homogeneity by microstruc-tural examination. Alloys containing up to 1 pct C were homogeneous in the presence of less than 3 pct N or 0. At higher alloying contents, some inhomo-geneities in the carbon distribution became evident. Alteration of the melting procedure toward longer times and higher currents did not improve the homogeneity of these alloys. Ti-N-0 alloys were homogeneous in the range to about 3 or 4 pct total alloying addition. Beyond this, almost all of the specimens showed as-cast microstructures consisting only of the phase. Consequently, inhomogene-ities could not be detected by examination of micro-structures. Ten alloys from each of the systems were analyzed for two of the elements present (oxygen being omitted in all cases and titanium being omitted in the Ti-C-N alloys). In all cases the analyses were found not to be sufficiently precise to serve as criteria for the total composition of these alloys. On the basis of phase distribution in heat-treated alloys, however, it appears that carbon is distributed throughout the alloys most uniformly, with oxygen and nitrogen following in that order. Heat Treatment: Specimens for heat treatment were wrapped in titanium sheet before sealing in the argon-filled quartz capsules. Heat-treatment times varied from 100 hr at 800°C to 0.5 hr at 1400 °C. After heat treatment the specimens were quenched by breaking the capsule in water. With the exception of alloys in the low composition region, heat treatment did not have an appreciable effect on the as-cast microstructures. Metallography: Following heat treatment, the specimens were prepared for metallographic examination by grinding on emery paper and electrolytic polishing. For the majority of the specimens a 10-sec etch with Remington "A" agent (25 pct HNO3, 25 pct Hf, and 50 pct glycerin) adequately

Jan 1, 1954