Search Documents

By Lawrence A. Shepard, Richard H. Krock

A series of ductile, two phase W-Ni-Fe composites, sintered in the presence of a liquid phase, were tested in tension. Identical room temperature stress-strain curves were obtained for specimens containing from 80 to 92 wt pct W (58 to 75 vol pct W particles). The composites exhibited a maximum elongation of 29pct at room temperature, and 10.7 pct at 77 °K. The tungsten particles in the composite elongated by the same amount at these temperatures. Single-phase alloy specimens matching the composition of the composite matrix showed about one half the flow stress of the composites. The test results demonstrated that the mechanical properties of W-Ni-Fe composites are determined by the tungsten particles alone and are independent of matrix volume fraction or mean free path over the composition range studied. A general prediction of the plastic behavior of two-phase alloys or composites from a knowledge of the properties, size, shape, and dispersion of the individual components is not, at present, possible. The excellent theoretical and experimental exposition of the dispersion-strengthening problem applies only to cases where the stronger phase is present as fine particles and in volume fractions of a few percent. Materials in which the stronger phase constitutes a significant fraction of the volume, although of broad engineering interest, have received comparatively little analytical attention. This important group of materials is the concern of the present study. The degree of deformation of each phase during the plastic working of a ductile two-phase alloy has been established by Boas and Honeycombe,1 Clare-brough,2 and Clarebrough and berger.3 These workers empirically correlated the recrystallization temperature of the individual phases with the degree of cold work sustained by each. It was shown that the deformations of both phases are equal and therefore the same as that sustained by the composite whole when the volume fraction of the stronger phase exceeds 35 pct. At yielding, however, deformation begins in the weaker phase. A direct correlation of composite strength, as a function of deformation, with the strength of the individual components cannot be made. In the presence of stronger phase particles, the flow stress of the weaker matrix phase is enhanced by an additional hydrostatic component. The development and potential magnitude of matrix strength promotion upon straining has not even been experimentally established. Quantitative work in this field has been confined to ceramic-metal composites. These materials fracture at strains considerably below 1 pct by cracking of the ceramic particles.4-6 An indication of the potential constrained matrix strength capacity is given by the observation by Nishimatsu and Gur-land4 that the maximum tensile strength of a WC-Co composite is almost double the bulk strength of the cobalt matrix. The relationship between composite strength and geometry is ordinarily expressed in terms of the mean-free matrix path (referred to subsequently as MFP) between the stronger particles.* This factor takes into account both the volume fraction of matrix phase and the average size of the particles. unke17 concluded that composite strength varied directly with the decreasing log MFP, following the Gen-samer8 relationship. However, Gurland and ~radzil~ reported a strength maximum at an MFP of the order of 1 p in WC-Co composites. A ductile two phase metallic composite was chosen for the present investigation of composite mechanical properties with the view of avoiding the difficulties and ambiguities noted above. The W-Ni-Fe heavy alloy of Green, Jones, and pitkin9 contains 80 to 94 wt pct W, and nickel and iron in a weight ratio

Jan 1, 1963

By I. Javandel, P. A. Witherspoon

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

Jan 1, 1969

By G. A. Stamboltzis, N. Arbiter, C. C. Harris

Fracture under low-velocity free-fall and double impact and under slow compression have been investigated. The pattern of breakage and the size distribution of resulting fragments of sand-cement and glass spheres have been determined. Photoelasticity methods were used to simulate the stress distributions in free-fall impact in order to explain the observed patterns of breakage. Oblique fracture planes, occurring only in free-fall impact, develop along lines which coincide with the trajectories of maximum compression as determined by mathematical analysis. Breakage efficiencies for different modes of fracture were compared for both types of spheres. For the same specimen and loading system, static loading and low-velocity dynamic loading induce geometrically similar stress fields resulting in reasonably similar fracture patterns and shapes of fragments. This work had its origin in a desire to understand three fundamental processes involved in autogenous comminution: namely; free-fall impact, double impact, and slow compression. To simplify the investigation, it is preferable to reduce the entire multi-stage fracture process to its most elementary form; that is, to study products resulting from a single-stage operation, such as the breakage of a specimen under single fracture conditions. The goal of this research is to study the energy utilization in the single fracture of brittle specimens and to relate the pattern of breakage and the resultant fragment size distribution with the nature of the material, the specimen size, the manner of load application, and the rate of loading. Published works on single fracture have mainly been concerned with tests on large (> 1-in.) irregular mineral pieces, especially coal and coke. These were mainly friability tests employing a single blow on closely sized pieces. More recently, small mineral specimens in the sieve range and glass spheres have been investigated.5-7 Kick8 has reported free-fall and double impact tests with large large cast-iron, cement, and clay spheres. The establishment of a minimum breaking height independent of size in the free-fall impact tests, and the direct proportionality between the minimum work required for fracture and the volume of the specimen in double impact tests were used by Kick as an experimental proof of his classical law of "Proportional Resistances." In the present work, spherical shapes were chosen because of their simple body geometry and consequent impact and stress field symmetry. Because the study involves several physical principles in connection with brittle-fracture, it may be of interest in fields where the strength of materials is of importance. APPARATUS The equipment for the free-fall impact testing is illustrated in Figs. 1 and 2. This consisted of a massive (17 x 17 x 3-in.) hard steel plate onto which specimens were dropped and a specimen release mechanism which could be set at any desired height up to 10 ft above the plate so that the impact velocity could be varied up to about 25 ft per sec. Double impact breakage is characterized by two points of loading situated at opposite poles. For this mode of breakage the apparatus was modified (Figs. 3, 4 and 5) so that a falling mass provided an impact of predetermined magnitude on a specimen resting on the plate. A spring-loaded device operated immediately after fracture to arrest the falling mass and to record its residual kinetic energy. Precautions were taken to avoid secondary breakage. Slow compression tests were performed on a conventional hydraulic testing machine, the specimen being held between two hardened parallel bearing blocks. In the case of glass spheres, carbide bits were used. Photoelasticity studies were performed on two-dimensional models held in a loading frame equipped with a force gauge (Fig. 6). The models were viewed in a polariscope and isochromatic fringe patterns were recorded photographically. SPECIMENS Sand-Cement Spheres: A series of sand-cement spheres were prepared by molding in spherical glass flasks.

Jan 1, 1970

By R. W. Gurry

In a series of charts this paper presents the composition of all gas mixtures, composed of the elements carbon, hydrogen, oxygen and nitrogen, which at temperatures from 1000° to 1800°F are in equilibrium with (1) austenite of various carbon contents, (2) iron and iron carbide, (3) graphite, (4) iron and its oxide. THIS paper presents the results of calculations made to determine the proportions of CO, CO2, CH4, H2, H2O and N2 in gas mixtures which, at temperatures from 1000" to 1800°F, satisfy, individually, the following conditions: (1) equilibrium with carbon steel of specified carbon content throughout the stable 7 iron region; (2) equilibrium with iron and iron carbide, both above and below the eutectoid temperature; (3) equilibrium with graphite; (4) equilibrium with iron and its oxide. From these data a new type of chart has been prepared which permits one to determine, from a room-temperature analysis of the furnace gas in question, the effect of the gas upon the steel at furnace temperature. Thus a guide is provided as to the carbon content which will result and as to the possibility of oxidation or sooting. All of the atmospheres commonly used in the commercial heat treatment of steel are composed of some or all of the simple constituents CO, CO2, CH4, H2, H2O and N2, and also in some cases oxygen or higher hydrocarbons. As the gas mixture is heated certain alterations in the proportions of the constituents take place, and at the temperature of the furnace the composition is, in general, not at all that of the entering gas. If equilibrium in the gas phase is presumed a calculation involving the use of certain fundamental data predicts for a certain temperature the carburizing, decarburizing, or sooting characteristics of a given atmosphere, and, upon the additional assumption of equilibrium between gas and metal, also the oxidizing ability of the gas and the ultimate carbon content of the steel. The extent to which these predictions are fulfilled will depend upon how closely the assumptions of equilibrium are realized in any particular case. Although at lower temperatures equilibrium is often not reached, there is considerable evidence indicating that above about 1000°F the approach to equilibrium in the gas phase is generally substantially complete, particularly in the vicinity of the steel which acts as a powerful catalyst in promoting reactions in the gas phase. The extent of the gas-metal reactions is commonly affected essentially by the sojourn at temperature; equilibrium with respect to carbon, for example, is accomplished simply by continuing the treatment long enough. However close the approximation of equilibrium may be in any instance, the above mentioned calculation indicates the condition of the steel at equilibrium, and at the present time this is the only basis of control of a furnace atmosphere. The effect at elevated temperature of a gas of a certain composition at room temperature upon a given steel has customarily been judged by first calculating from gas equilibrium constants the proportions of the gas constituents at the temperature of the steel, then comparing certain of these with the proportions known to be in carbon or oxygen equilibrium. This process is tedious and time-consuming, in particular because the calculation generally involves a series of successive approximations; and since some knowledge and familiarity with this type of computation are required, not everyone needing the information has been, in the past, in a position to obtain it. By making a series of equilibrium calculations, combining with fundamental data for gas-metal reactions, and plotting the results in a special way, a chart has been prepared which enables anyone

Jan 1, 1951

By L. F. Mondolfo, B. E. Sundquist

The crystallographic orientation relationships resulting when lead is nucleated from the liquid by Ni, Cu, Ag, and Ge were determined. For each nucleating agent several definite orientatioz relationships were found. These relationships seemed to be controlled by good symmetry relations and low crystallographic disregistry between mating planes. For any given nucleating agent the under colling for nucleation was found fairly constant and independent of the orientation relationship and consequent disregistry. It was also found that, upon re melting and refreezing the Pb, the orientation relationship was changed. These findings prove that crystallographic disregistry is not the controlling factor in heterogeneous nucleation from the liquid. The results of this investigation tend to confirm the theory presented in a preceding paper that heterogeneous nucleation starts with the formation of an adsorbed layer of nucleated metal on the nucleat-ing impurity. Evidence is given that cavities in the nucleating agent act as centers of nucleation. IT has long been known' that solid extraneous particles are active in catalyzing phase transformations that occur in a system, particularly condensation and crystallization. It is well established that these heterogeneities act as catalysts by providing surfaces upon which nuclei of the precipitating phase can form with activation energies smaller than those required for homogeneous nucleation. Numerous investigations have shown that in this process of heterogeneous nucleation: a) the nucleus forms with one, or several, definite crystallographic orientation relationships with the nucleating phase2-4 and b) that there is a small range of undercoolings or super saturations characteristic of the nucleation of a given solid on a given Substrate.5-10 Turnbull and vonnegut11 have developed a theory based on theories developed by Volmer12 and Turn-bull and Fisher1= for heterogeneous nucleation from gases and liquids, that relates the super saturation or undercooling required for nucleation to the dis-registry between the lattices of the nucleus and the nucleating agent. This theory predicts that nucleation should occur with the orientation relationship between the nucleus and nucleating agent that minimizes the disregistry. Further, it predicts that the undercooling or super saturation necessary for nucleation should be a function of the disregistry. Numerous investigations have dealt with the orientation relationships resulting from the condensation of vapors onto crystalline solid substrates2,3 and a few with the nucleation of one phase by a second phase in solidification4,14. Others have dealt with the supersaturation8-10 and undercooling5-7 associated with nucleation in condensation and solidification respectively. However, there is virtually no report that gives both of these factors for the same system. In this investigation a study was made of the undercoolings and orientation relationships resulting when Pb is nucleated from the liquid by Ni, Cu, Ag, and Ge. It was the purpose of this investigation to check the Turnbull-Vonnegut theory, i.e., the importance of crystallographic disregistry between nucleating catalyst and nucleated metal. The results indicate that disregistry is not an important factor in nucleation and that the nucleation process is probably somewhat more complex than current theories suggest. EXPERMENTAL PROCEDURE Small single crystals of nickel, copper, silver, and germanium were prepared from materials of four to five nines purity, and the Pb used was also 99.999+ pet pure. Cu and Ag single crystals were prepared by sealing small chips of Cu or Ag in an evacuated quartz capsule and heating the capsule at 2000°F for 1 hr before cooling. Nickel crystals of 200 diam were also prepared in evacuated quartz capsules, but melting was done by heating the capsules in an oxy-acetylene flame for a few minutes. These spheres were invariably polycrystalline so

Jan 1, 1962

By Seymour G. Epstein

Using a modified capillary-reservoir technique, electromigration of dilute solutions of silver in liquid bismuth has been measured with relatively good precision. The effects of experiment duration, temperature. current density, and silver concentration have been investigated. Chemical diffusion coefficients of silver in liquid bismuth were also measured at temperatures ranging from 300°to 600°C. The effective valence of silver in bismuth was calculated as a function of temperature from a comparison of electric mobilities and diffusion coefficients. In an alloy containing 1 wt pct Ag, the electric mobility of silver increased from 1.65 x X sq cm per v-sec at 300°C to 2.38 x 10-4at 600°C; the diffusion coefficient in-creased from 3.0 to 17.5 x 10-5 sq cm per sec over the same temperature range. Calculated values of effective valence decreased from 0.28 at 300°C to 0.10 at 600°c. ALTHOUGH there has been considerable progress in the theoretical treatment of electron-transport properties in liquid metals, the mechanisms of mass-transport phenomena in liquid metals are still not clearly understood. This is especially true for the case of electro migration, alternatively termed elec-trotransport or electrodiffusion, which is mass transport induced by an applied electric field. A thorough knowledge of the mechanisms governing electromigration will considerably aid in understanding the electronic state of an impurity atom in a liquid metal. Electromigration has been long observed and studied in liquid metals and alloys. However, no satisfactory theory to explain the observed effects has been developed. The disparity in the present theoretical treatments, summarized by Verhoeven,1 has been attributed to the difficulty in describing the momentum transfer between electrons and ions and to the lack of precise measurements in liquid metals, especially at elevated temperatures. Prior to the last decade most of the electromigration studies were qualitative, in that only the direction of migration of alloy constituents was determined with any degree of reliability. These results have been summarized by schwarz,2 who was able to measure the electric mobility of several elements in mercury. More recently electromigration in liquid metals has been measured by a variety of techniques, with a varying degree of success. Mangelsdorf3,4 made precise measurements on dilute amalgams at room temperature, using a technique relating changes in composition to changes in alloy resistivity, but was unable to make measurements at elevated temperatures. Using a modification of this technique, Verhoeven and Hucke5 measured electromigration and resistivities in liquid Bi-Sn alloys. However, the resistivity technique is not sensitive to small composition changes and is impractical for very dilute alloys. Verhoeven and Hucke6 also measured electromigration of several solutes in liquid bismuth by adapting the capillary-reservoir method for determining diffusion coefficients in liquids. The reported precision of these measurements was poor, however, primarily due to the large uncertainties in the chemical analyses to determine concentration changes. Belashchenko7 has summarized the recent Russian investigations of electromigration, many of which, however, are of questionable accuracy. The present study was initiated to develop an experimental technique which would be both versatile and sufficiently sensitive to yield precise measurements of electromigration of solutes in dilute metal alloys at elevated temperatures. Using this technique, a comprehensive study of silver in liquid bismuth was made to determine the effects of experimental variables on the apparent electric mobilities. Chemical diffusion coefficients of silver in liquid bismuth were also measured over the same temperature range for comparison with the electromigration data. EXPERIMENTAL TECHNIQUE Electromigration Measurements. Electromigration of silver in liquid bismuth was measured by the capillary-reservoir technique, similar to that used by Verhoeven and Hucke,6 but modified to yield more

Jan 1, 1967

By A. J. Ardell, O. D. Sherby

The elevated-temperature steady-state creep behavior of polycrystalline a Zr was studied in vacuo under constant tensile stress. The experiments were conducted from 660° to 845°C over the stress range 870 to 5400 psi. The steady-state creep rate, is, of a Zr was found to be related to the applied stress, a, and temperature, T, by the equation where C is a constant. Two interpretations of the data are considered, each requiring different values of the parameters ß, Q, and . It is shown that this formula can be derived using Gilman's phenomeno logical equation relating dislocation velocity to the shear stress for the case where dislocation mobility is limited by the nucleation of kink motion. From the fundamental dependence of strain rate on dislocation velocity and the density of mobile dislocations, the factor exp(ßo) can be related to the stress dependence of the dislocation density, and Q to the energy of kink formation. The two interpretations considered are as follows: Case I—the parameters are adjusted so that the equation above describes the data over the entire range of stress and temperature, implying that creep is controlled exclusively by dislocation glide; Case IIthe parameters are adjusted to allow for the possibility that creep is diffusion-controlled at high stresses but glide-controlled at low stresses. For both cases the values of the parameters are within the spread permitted by the data. Factors favorable to the observation of the unusual stress -dependence of i, for a Zr are discussed and it is shown that the anomalously rapid rates of self-diffusion and the low modulus on the principal glide system are primarily responsible for the anomalous creep behavior. Metal-lograflhic observations are presented and discussed with reference to the overall creep behavior. PREVIOUS investigations of the creep behavior of a Zr have been conducted from room temperature to about 500C The results of these investigations are primarily of engineering significance because of the variety of conditions under which they were carried out, and the low testing temperatures employed. Little information is available concerning the fundamental nature of the elevated-temperature creep behavior of a Zr. This investigation of the creep behavior of a Zr above half the melting temperature (the melting temperature, T,, of a Zr is estimated to be about 1860"K was initiated because there are two significant anomalies of the physical properties of a Zr: a) the rate of self-diffusion in polycrystalline a Zr is anomalously rapid: b) both Young's modulus for polycrystalline a Zr and the shear modylus as-sociated with the prismatic slip system { 1010)( 1120) (the principal slip system in a Zr 7) have anomalously large temperature variations The importance of these quantities to the understanding of creep behavior has been demonstrated, and it was hoped that this investigation would provide a critical test of the predictions based upon the analysis of creep data for other pure metals. However, as the investigation progressed it became increasingly clear that the unusual creep behavior of a Zr could not be analyzed in a conventional manner. 1) EXPERIMENTAL EQUIPMENT AND PROCEDURES Two grades of zirconium, designated A and B, were used in this investigation. The analyses, as supplied by the vendors, are presented in Table I. Material A, with a nominal purity of 99.84 pct, was supplied by the Wah Chang Corp. in the form of cold-rolled sheet from which specimens were machined without further fabrication. Material B, with a nominal purity of 99.95 pct, was purchased from the Nuclear Materials and Equipment Corp. in the form of iodide crystal bar.

Jan 1, 1968

By R. W. Fountain, J. F. Libsch

ONSIDERABLE time and effort have been ex- pended recently in research designed to provide a better understanding of the solid state transformations leading to the permanent magnet qualities of many commercial alloys. The knowledge derived from these investigations has resulted in the development of new magnetic materials as well as improvements in some of the older materials. The permanent magnet qualities of the newer alloys having high residual inductions and high coercive force values, have been attributed to the precipitation of a second phase or an order-disorder transformation in a solid solution. The ternary Co-Fe-V alloys provide an interesting group of magnetic materials having a wide range of magnetic properties. The low vanadium alloys' possess soft magnetic properties similar to the binary Co-Fe alloys, while at higher vanadium contents (about 10 pct) interesting permanent magnet properties are obtained.2,3 Previous work by Nesbitt' on the ternary Co-Fe-V permanent magnet materials, resulted in the conclusion that the precipitation of a second phase, ?, was responsible for the permanent magnet qualities of these alloys. However, unlike the conventional precipitation hardening alloys, the ? phase, which is stable at high temperatures, is precipitated in the low temperature a phase upon aging. Recent work by Geisler and Martin"." has established that the order-disorder reaction inherent in the binary Co-Fe alloys is also present in the ternary Co-Fe-V alloys. On the basis of their investigation, it was proposed that ordering of the a phase is also a possible explanation of the mechanical and magnetic hardening in these alloys. Thus, it appears that the Co-Fe-V alloys are rather unique in that two solid state transformations are taking place which may account for the properties of these alloys. In view of these facts, the present investigation was undertaken to study the mechanism of mechanical and magnetic hardening in the 10 pct V-Co-Fe alloy with particular reference to the role of the ? phase; and to correlate the results of this study with other magnetic properties. In addition, the influence of increasing vanadium contents on the transformation of Co-Fe alloys was also undertaken to supplement the main objective of the investigation. Experimental Procedure The Co-Fe-V alloys used for study consisted of: 1—a commercial lot of Vicalloy I, a 10 pct V alloy, and 2—alloys prepared by powder metallurgy and containing 52 pct Co and 0 to 14 pct V, balance iron. The commercial alloy was received in the cast and forged condition and was used primarily for the detailed study of hardness, electrical resistivity, and magnetic properties of the 10 pct V alloy. The powder metallurgy alloys were prepared from electrolytic iron, commercial cobalt, and ferro-vanadium powders. The processing cycle was as follows: 1—Elemental powders were weighed and mechanically mixed for 24 hr. 2—Mixed powders were pressed at 60,000 psi into the form of bars 5/8 x 5/8 x 5 in. 3—Bars were sintered in an atmosphere of purified dry hydrogen for 24 hr at temperatures ranging from 1380" to 1340°C, depending upon the vanadium content. The as-sintered bars were cut in half; and one half was machined to 7/16 in. diam for thermal analysis. The remainder of the bar was hot forged to approximately 0.050 in. for X-ray diffraction analysis. Additional details of sample preparation are presented subsequently according to the type of test performed. Thermal Analysis Measurements One of the simplest and most widely used tools for determining phase relations as they change with temperature and composition is to study the rate of change of temperature of a material as heat is supplied or extracted at a constant rate. The thermal analysis method used in this investigation is the inverse-rate method described by Smith.'

Jan 1, 1954

By D. S. Rowley, F. C. Appl

Assuming that rock behavior, during cutting with a single diamond, may be approximated by that of a rigid, Coulomb, plastic material, a theory of single diamond cutting action has been developed. Using this theory, the stresses on the diamond cutting surface and the components of the cutting force have been determined. Theoretical results agree reasonably well with available experimental data. The theory of cutting of single diamonds may be used as the basis for subsequent study of the perfonnance of surface set diamond drill bits and other surface set diamond cutting tools. INTRODUCTION For more than twenty years, surface set diamond bits have been used successfully in the mining and petroleum industries to drill and core rock formations that range from medium hard to the hardest known formations. Most extremely hard rock formations cannot be economically drilled or cored by any other means. In view of the growing use and importance of surface set diamond drill bits, it is reasonable to study the way in which these tools cut, and the factors that influence wear of the diamonds. Such knowledge will contribute to better bit design and will improve operating procedures. To analyze the cutting performance of these tools it is appropriate to consider the way in which an individual diamond cuts. Once this is established, the over-all performance of the bits can be determined by considering the total effect of the individual diamonds. In recent years there has been a growing interest in rock mechanics and there now exist several publications dealing with the deformation and failure of rocks. Notable experimental work has been reported by Bredthauer,1 Handin and Hager,2, 3 and Gnirk and Cheatham.4 An extensive collection of data pertaining to rock properties at low confining pressures has been given by Wuerker.5 Cheatham and Gnirk have recently presented an excellent review of the state of knowledge regarding analytical and experimental rock deformation and failure.6 Generally rocks are extremely weak in uniaxial tension, and weak or only moderately strong in uniaxial compression. Both tensile failure and compressive failure are generally brittle in nature. This means that there is little or no deformation or "flow" of the material before shear or tensile fracture occurs. Thus, until fracture occurs, rock behaves very nearly like an elastic material in uniaxial stress. Significant changes in behavior of rock occur under conditions of triaxial compressive stress. There is generally a marked increase in compressive strength and in ductility. It has been demonstrated that the failure becomes plastic (in a macroscopic sense) at sufficiently high values of confining pressure. Sharp wedge chisel indentation studies conducted by Gnirk and Cheatham indicate that many rocks behave in a plastic manner at relatively low triaxial compressive stress states.4 If "failure" under triaxial compressive conditions is considered to be either complete rupture or macroscopic plastic flow, then the general failure characteristics of a rock can be reasonably well represented by the commonly known Mohr envelope. This certainly is not a perfect failure criterion (since the effect of the intermediate principal stress is neglected) but it has been found to be adequate for most engineering purposes. Mohr envelopes usually are constructed graphically using experimental triaxial test data. For purposes of analysis it is convenient to represent the envelope in equation form. Among the several approximate forms appearing in the literature, a parabolic form is suggested by Cheatham,7 and an exponential form is suggested by Paone and Tandanand.8 The exponential form has been adopted for this study, and based on published data, appropriate parameters have been established for 19 rock types. A typical plot of the data compared with the approximate envelope, and a summary of the rock properties and

Jan 1, 1969

By John C. Sawyer

A brief review is given to show that catastrophic oxidation of stainless steel in the presence of lead oxide is similar in many respects to catastrophic oxidation as caused by MoO3 and V2O5, Experimental data are presented to show that lead oxide accelerates the normal oxidation process of stainless steel by chemically altering the normally protective oxide film formed on the base metal. Data are given to show that similar reactions occur during oxidation of stainless steel in the presence of MOO3,and V2O5. THE mechanism by which pure metals oxidize at high temperatures has been the subject of intensive research for many years and much has been accomplished towards understanding the processes involved. The mechanisms by which alloys oxidize have also been investigated, but due partly to the fewer number of investigations and partly to the complexity of the oxide systems present, available knowledge is not as extensive as with pure metals. When oxidation processes are further complicated by atmospheres contaminated with sulfide, halide, or oxide vapors, the knowledge of fundamentals becomes meager or altogether lacking. Thus, it is not surprising to find that the catastrophic oxidation of stainless steels in metal oxide vapors of MOO, and V2O5 is not a well understood phenomenon despite the fact that a number of investigations have been undertaken. The literature in dealing with the subject of catastrophic oxidation is primarily concerned with the effects of MOO3 and V2O5, but does mention that other metal oxide vapors, one of which is lead oxide, appear to promote a similar type of accelerated oxidation. While this latter effect is virtually unknown, the deleterious effect of lead oxide is of considerable commercial importance, as it is believed to be one of the fundamental reasons for exhaust-valve failure in gasoline engines. For this reason studies have been in progress for many years attempting to better understand the mechanism involved. Modern exhaust-valve alloys"2containing 20 to 25 pct Cr plus nickel and/or manganese are in effect high-grade austenitic stainless steels exhibiting ex- cellent oxidation resistance in air at exhaust valve operating temperatures of 1100° to 1700°F; yet these same alloys are prone to corrode during engine service because of the conditions attending their operation. These conditions are such that the valves become coated with deposits of complex lead compounds mainly derived from gasoline additives of tetraethyllead, ethylene dibromide, and ethylene dichloride. Literature searches3-7 and discussions with various investigators show that considerable difference of opinion exists as to the nature of the lead-bearing deposits found on the surface of exhaust valves and the atmospheres to which they are subjected; yet, all seem to generally agree that corrosion, a gradual transformation of alloy to nonmetallic products of corrosion, is a consequence of the presence of lead compounds. This majority opinion seems to be substantiated by laboratory engine tests which show that corrosion does not occur to any appreciable extent when lead-free fuels are used. The first observations of catastrophic oxidation reported in the literature were probably made by pfeil,8 but it was Leslie and Fontana9 who undertook the first detailed study of the phenomenon. Their work was chiefly centered on the acceleration of oxidation which took place when stainless steels containing more than 3 pct Mo were heated in air with restricted circulation. These investigators reported that the presence of other oxides near the surface would further accelerate the rapid oxidation which occurred. While they recognized that the effect was attributable to the presence of molybdenum, it was Rathenau and Meijering" who emphasized the fact that the presence of molybdenum in the alloy was incidental to the presence of MOO3 vapor in the atmosphere surrounding the specimens. Evans,11 as well as Cunningham and Brasunas 12 studied similar

Jan 1, 1962

By P. H. Turnock, R. D. Pehlke

The effects of mutual interactions between elements on the solubility of nitrogen in liquid iron alloys containing chrormium, columbium, molybdenurn, nickel, or silicon have been determined. The equilibrium nitrogen solubilities in the liquid iron alloys were measured by the Sieverts method. The second-order effects caused by the presence of chrormium, columbium , or silicon in the melt were found to be significant. The solubility of nitrogen in liquid iron alloys containing several alloying elements has also been measured as a function of melt temperature and nitrogen pressure. The heat of solution of nitrogen in microcomponent iron alloys has been found to be a function of the logarithnz of the activity coefficient of nitrogen, irrvespectilje of the composition of the melt. The heat of solution of nitrogen in pure liquid iron was determined to be 1265 ± 450 cal per g-atom of N over the temperatxue range 1600° to 1800°. Sieverts ' law was obeyed for all melt compositions studied in the pressure interval 0.4 to 1.0 atm. The solubility of nitrogen in liquid pure iron and in many dilute binary liquid-iron alloys has been the subject of many investigations. Pehlke and Elliott' made a comprehensive investigation on this subject in which they summarized and compared previous researches. However, the practical application of these solubility studies requires prediction of the nitrogen solubility in the multicomponent alloys encountered in operating practice. The interaction parameter as defined by wagner2 has been used to predict the behavior of a solute in a complex alloy, and a modification of this approach has met with reasonable success in several instances.'1"5 The necessity for extending the nitrogen solubility data on binary alloys is caused by the lack of direct measurements of the solubility of nitrogen in molten alloy steels. A direct comparison between measured solubilities and those predicted for multicomponent alloys based on data in binary systems has been made in only a few instances, and most of these were reviewed by Langenberg,6 who developed a graphical technique for predicting nitrogen solubilities in multicomponent iron alloys. The solubility of nitrogen has also been studied over the entire Fe-Cr-Ni ternary system by Humbert and Elliott.? Both studies compared the solubilities computed from binary iron alloy data with those measured in ternary alloys in the liquid state at 1600°C. Most of the nitrogen solubility measurements referred to above were for iron alloys at 1600°C; thus, the ability to predict nitrogen solubilities in molten iron alloys has been restricted primarily to alloys at 1600°C. Recently, Nelsson8 devised an approximate method for calculating the change in nitrogen solubility in molten iron alloys as a function of temperature, but an assumption of regular solution behavior limits the applicability of the method. Chipman and Corrigan9,10 have developed an empirical correlation. based on the activity coefficient of nitrogen in the iron alloy at 1600°C, which facilitates the estimation of nitrogen solubilities at higher temperatures. This method also has limitations, due to the necessity of calculating the activity coefficient for nitrogen in a multicomponent iron alloy at 1600°C from the nitrogen solubility data for binary iron alloys. The shortcoming of such a straightforward approach to the problem of predicting nitrogen solubilities is the disregard of mutual interactions between elements in a complex iron melt. This investigation was conducted in an attempt to provide a more reliable method of predicting nitrogen solubilities in multicomponent iron alloys. The effects of the mutual interactions between chromium, nickel, silicon, molybdenum, or columbium in iron were studied, using the Sieverts' technique1,7 to measure the effects of alloy composition, temperature, and pressure on the solubility of nitrogen in ternary liquid iron alloys.

Jan 1, 1967

By F. A. Schaufelberger

METAL can be recovered from a leach solution either indirectly by precipitation as a compound that is later reduced or directly by electrolysis, cementation, or chemical reduction, for example, with hydrogen. The widely used electrowinning processes suffer from the inefficiency of converting fuel to low voltage dc and from low current efficiency when complete recovery is attempted. With the batch hydrogen reduction methods outlined here, a unit of fuel converted to hydrogen in modern gas reform plants will make two to three times as much metal as can be obtained with electrolysis. With continuous hydrogen reduction, four to six times as much metal is obtained. Precipitation by reducing a metal salt solution with H, has been known for almost 100 years.' Commercial use of the process, however, awaited the research and development program initiated by Chemical Construction Corp. Within the last few years this program has brought about construction of commercial plants, listed on page 545. Ideal hydrogen reduction will precipitate a pure metal from a solution obtained by commercial leaching methods at a rapid rate without excessive temperatures and pressures. The metal precipitate will be of desired size and density, less than a percent left in the vessel being deposited on the wetted parts. The present discussion will outline Chemical Construction Corp.'s early development program and will discuss the chemistry and mechanics of reducing copper, nickel, cobalt, and cadmium from solution by H2. Work on selective reduction of nickel from cobalt has been described earlier.' The article is not concerned with precipitation of copper by disproportionation of cuprous solutions where yield is limited3 to 50 pet.* It does not discuss use of gaseous reducing agents other than hydrogen, such as SO2,4 which may lead to contamination with sulfur, or C0,5 which is considerably more expensive than hydrogen and produces a gaseous reaction product, CO2.† The article does not include use of other nongaseous reducing agents,7 which are far more expensive than CO. Also, notes on the history, engineering design, and performance of commercial plants have been given elsewhere.' It will be shown that a pure metal can be precipitated by reduction with hydrogen from solutions obtained by commercial leaching methods when the solution composition is controlled and proper acidity, proper metal ion concentration through controlled complex formation and hydrolysis, and adequate agitation for suspension of metal and for transfer of reducing gas are maintained. Reasonable temperatures and pressures are used which are less excessive than those used by previous workers, although they are maintained at above equilibrium values because of process kinetics.' It has also been found necessary to induce nuclea-tion and to control growth and agglomeration in order to get a powder product of suitable physical properties. Review of the Literature Muller, Schlecht, and Schubardt of I. G. Farben10 claimed a successive reduction of silver, copper, nickel, cobalt, and zinc from ammoniacal solution by applying progressively higher temperatures and hydrogen pressures. Metals produced by this technique were not pure, since no attempt was made to adjust the solution composition between the different reductions. It is doubtful that the zinc product contained any metallic zinc. The major contribution to the literature" was made by the Ipatiews,12 who prepared platinum, iridium, copper, nickel, cobalt, lead, tin, arsenic, antimony, and bismuth at often extreme conditions, up to the critical temperature of water (373ºC), at 1500 to 7500 psi, for as long as several days. Even cadmium and zinc were claimed to have been observed in trace quantities precipitated together with basic salts. In the early experiments solutions were not stirred. Extensive formation of oxides and basic salts occurred, probably long before sufficient hydrogen could be supplied in the unstirred liquid to decrease the metal concentration in solution by reduction. The importance of hydrogen pressure was

Jan 1, 1957

By P. L. Raffo, T. E. Mitchell

Single crystals of Ta-Re alloys have been deformed in tension and compression as a function of solute concentration and temperature. The temperature dependence of the yield stress is decreased by alloying. At low temperatures there is a minimum in the yield stress (at about 2 at. pct) as a function of concentration. Three-stage hardening was not observed in the alloys. Linear hardening rates, where observed, were characteristic of stage II in pure tantalum. At low temperatures and high concentrations, work softening was often prominent at high strains. Electron microscopy revealed that dislocations are mostly of screw orientation, although heavily jogged and cusped, and often occurring in dipole pairs. The alloys become more prone to cleavage with increasing rhenium content, due mostly to an increasing yield stress rather than a change in cleavage strength. Alloying apparently decreases the Peierls stress of bcc metals at low temperatures by facilitating kink nucleation at solule atoms on screw dislocations. Higher solute concentrations further increase the yield stress. Quantitative calculations and comparisons with other bcc alloy systems are made to show that the predominant hardening mechanism is due to the interaction of dislocations with misfitting solute atoms close to their slip planes. The interaction is so localized that the core dilatation of dislocations could be important. CONSIDERABLE interest has arisen in recent years concerning the mechanical properties of unalloyed bcc metal single crystals.'-7 These and similar studies on polycrystalline materials have resulted in a considerable understanding of the macroscopic yielding characteristics8, 9 and the process of work hardening.'' This latter advance is due in a large part to the increased use of transmission electron microscopy in P I RAFFO junior Membar AlMF is Metallyureist Iewis such studies and the general observance of three stages of work hardening in the shear stress-shear strain curves. The properties of substitutional alloy single crystals of bcc metals have not been determined in as much detail as in the pure metals. Insufficient systematic studies have been made on the effect of alloying on the magnitude and temperature dependence of the yield stress or the stress-strain curves. The alloy systems which have been studied to some extent have been Ta-W,11, 12 Ta-Mo," Ta-Nb,12'14 Mo-Re,15 W-Re,16 and alloys of iron with nickel,17 platinum,17 vanadium,18 and aluminum.18 A large body of data exists on Fe-Si alloys, but these have mainly been a result of the experimental convenience of using this material in place of unalloyed iron, rather than a comprehensive study of the mechanical properties of the alloy system per se. The present authors have previously reported additional data on the properties of Ta-Re, Ta-W, and Ta-Mo alloy single crystals.1g Some of the data on Ta-Re alloys will be described later in this paper. Briefly, single crystals of Ta-Re (1.2 to 10.6 pct), Ta-W (6.2 and 10.2 pct), and Ta-Mo (2.9 and 9.0 pct) were grown by the floating zone method in vacuum. Compression tests were performed in the temperature range 77" to 623°K. Slip was found to occur in Ta-Re alloys at 77' and 298°K on the (~ll)[lil] system and slip lines were coarser and straighter than in unalloyed tantalum. No three-stage hardening was observed in any of the alloy systems and the curves were nearly parabolic. At temperatures above 175°K the yield stress* increased monotonically with in- creasing composition while at 77" and 175°K the yield stress initially passed through a minimum after which it began to increase again for concentrations greater than about 2 pct Re. Similar minima were observed in Ta-Mo alloys and to a lesser extent in Ta-W alloys. The minimum was interpreted as being a result of a

Jan 1, 1969

By G. R. Speich, S. A. Kulin

The isothermal formation of martensite at subzero temperatures has been studied in an austenitic stainless steel. The amount of martensite formed isothermally in a given time was found to follow a C-curve behavior with decreasing temperature. Since isothermal rnartensitic transformation occurs at temperatures above the true M. temperature of the olloy, the observed M. temperature is dependent on the cooling rate. It is believed that a general martensite theory must include the basic concepts of the strain embryo theory, and also the important role of thermal fluctuations. UNTIL recently the isothermal component of the martensitic transformation had been considered as a trivial effect in theoretical considerations of this type of transformation. One of the postulates used in defining a martensitic transformation was that the reaction was temperature dependent but not time dependent. Cohen and his associates1'2 were the first to report the isothermal formation of mar-tensite. They obsel-ved the time dependency of the martensite reaction at room temperature in plain carbon steels and certain tool steels. They found that a small amount of additional transformation occurred after the specimen reached the quenching temperature and that the rate of this isothermal transformation decreased rapidly with time. These authors considered this phenomenon to be a "tailing-off" of the main athermal reaction. Kurdjumov and MaksimovaL 5 studied the isothermal formation of martensite at subatmospheric temperatures in steels containing 0.6 pct C-6 pct Mn, 1.6 pct C, and in an iron-base alloy containing 23 pct Ni and 3.4 pct Mn. In each case it was found that both the initial rate and the total amount of isothermal transformation exhibited sharp maxima with decreasing temperature. They stated that it was possible to completely suppress the martensitic transformation in the 0.6 pct C-6 pct Mn steel by drastic quenching in liquid nitrogen. They interpreted their results in terms of a thermally activated process. The maximum in the rate of nucleation vs. temperature curve was assumed to be a result of the temperature dependence of the rate of growth of the martensite plates. It has been demonstrated," however, that martensite plates can form with great rapidity even at temperatures approaching absolute zero. In addition," attempts to confirm the reported suppression of the martensitic transformation by rapid cooling of a 0.6 pct C-6 pct Mn steel have proved unsuccessful. On this basis, and on other theoretical grounds, the validity of some of Kurd-jumov's results and conclusions is open to serious question. Cohen, Machlin, and Paranjpe' have proposed a theory of martensite formation based on the nucleation of strain embryos and a growth mechanism involving cooperative atomic displacements. According to this theory thermal fluctuations can make a strain embryo supercritical and thereby produce martensite isothermally. Das Gupta and Lement' have investigated the isothermal formation of martensite in a 17 pct Cr-0.7 pct C steel at temperatures below the ambient. They found that partial suppression of the martensitic reaction occurred on quenching into liquid nitrogen. That part of the reaction which can be suppressed is taken to correspond to the isothermal rather than the athermal part of the rnartensitic transformation. The initial rate of isothermal transformation for this steel was found to go through a maximum with decreasing temperature. This is in conformity with the findings of Kurdjumov and Maksimova for the steels used in their investigation. Das Gupta and Lement conclude that isothermal martensite formation occurs by the nucleation of new plates and is always preceded by some athermal transformation. An iron-base alloy containing 14 pct Cr-9 pct Ni was used in this investigation to study further the kinetics of the austenite to martensite transforma-

Jan 1, 1953

By S. S. Marsden, S. H. Raza

An experimental study of the flow of line-textured. aqueous foams through Pyrex tubes is described. The foams range in quality F (ratio of gas volume to total volume) from 0.70 to 0.96 and behave like pseudoplastic fluids. At lower flow rates they exhibit laminar flow ad have apparent viscosities which increase with quality, and which cover a range of 15 cp to 255 poise for tubes of 0.25- to 1.50-mm radius rj. At higher flow rates a plug-like type of flow is developed, the extent of which increases with both T and ri. When the same foams flow through either open or packed Pyrex tubes, remarkably high streaming potentials ?E are often generated. These can easily reach 50v if nonionic foaming agents are used, but are at least an order of magnitude less for ionic foaming agents. A linear relationship between ?E and the pressure differential ?p is observed; this usually extrapolates to positive values of ?p at ?E of zero. The slope of the line increases with both F and ri. An equation was derived to describe the streaming potential of non- Newtonian fluids in circular tubes and was used to correlate experimental results. The calculated < potentials are of the right order of magnitude. INTRODUCTION Foams are both unusual and inniguing in their physical properties, and have been the subject of many scientific studies. However, present knowledge of foams is still fragmentary, specific and often contradictory. Apparent viscosity of foam is the physical property of greatest interest to both rheo1ogists and engineers. sibreel reported that the apparent viscosity decreased with increasing shear rate in a manner similar to some non-Newtonian fluids. Penny and Blackman2 reported that fire-fighting foams had both a limiting shear stress and a tensile yield stress. There is little doubt that some foams at least behave like non-Newtonian fluids, and have apparent viscosities considerably higher than those of either constituent phase. The high apparent viscosity of foam with its concomittant effect on mobility ratio and sweep efficiency no doubt prompted several attempts by research groups to use foam as a displacing agent in porous media. Based on recent experience, most of these groups probably succeeded in completely blocking fluid flow in the porous media and then abandoned their efforts. Two groups apparently found the successful combination of experimental parameters at about the same time.3,7 Others have recently added to our knowledge of foam flow in porous media and its use as a displacing agent.8-13 An experimental problem encountered by Fried3 was a transient blockage of foam flow in porous media when distilled water was used to prepare the foam-producing solution. Fried surmised that this was due to an electrokinetic effect and he eliminated it by using electrolytes in preparing foaming solutions. He also measured the streaming potential of a number of foams in capillary tubes which he found to be appreciably higher than those obtained when the constituent liquid flowed under comparable conditions. This paper presents results of a more comprehensive study of the streaming potential generated by aqueous foam flowing in both open and packed Pyrex tubes. It also adds to knowledge of the rheology of these foams as deduced from their flow behavior in open tubes. APPARATUS AND PROCEDURE A diagram of the apparatus used is shown in Fig. 1. Details of its construction, testing and use are described elsewhere.14 Careful selection of materials, extreme cleanliness and rather elaborate electrical insulation and shielding were necessary to obtain reproducible results (±5 percent). Both streaming potential and streaming current were measured with an electrometer. The design of the foam generator developed for this work is novel (Fig. 2). It will produce uniform fine-textured foams at pressures sufficiently

By Charles Stein

The effects of several different prior treatments an the creep behavior of 99.9995 pct aluminum at 260°C and 1000 psi canstant stress are compared with annealed specimens. Quenching from 538oC, irradiation with 2 mev electrons, and tension-compression room temperature fatigue damage were used to change the substructure of the specimens prior to their creep testing. The quenched and the irradiated specimens showed a larger primary and transient stage contribution to the creep curve than those specimens which were only annealed prim to creep testing. The specimen receiving prior fatigue damage at room temperature showed no first stage or transient creep when tested under identical conditions as the above specimens and had an average creep rate seven orders of magnitude lower than that of either the annealed, the quenched or the irradiated creep specimens. In a previous electron transmission microscopy investigationl of the substructure developed during the creep of pure aluminum at 0.57 Tm, it was noted that a large concentration of vacancy loops and dislocation loops were present in and near the subboundaries while the interior of the subgrains contained few of these defects and had a low dislocation density, see Fig. 1. Yim and Grant2 and Hazlett3 have shown that the presence in nickel of a substructure developed by prior cold work effectively reduces the primary and transient contributions to the creep curve. Hultgren,4 McLean,5 Gervais, Norton, and Grant,6 Chang and Grant,7 and others have shown the same effect with the higher stacking fault energy material, aluminum. However, the specific defect responsible for the change in the creep behavior developed by the prior cold work was not established. Specifically, what effectively interfered with dislocation motion in these materials at elevated temperatures? Was it jogs on dislocations produced by their interaction with other dislocations or with vacancy or dipole loops, see Fig. 1, or did the subgrain walls determine the mean free path for glissile dislocations? This would be a realistic possibility only if subgrain boundaries are effective barriers to dislocation motion. The ability of a subgrain boundary to act effectively against glissile dislocations depends on the number and the arrangement of dislocations in the subboundary,8 which in turn is a function of the creep strain.&apos; This paper compares the creep rate of specimens possessing a large vacancy concentration to that of annealed specimens and with specimens having a stable subboundary wall containing Lomer locks produced by prior fatigue damage. EXPERIMENTAL PROCEDURE Vacancy Loops. Aluminum creep specimens, 99.9995 pct pure, were machined from zone refined rod to a 2-in. gage length and a cross-sectional area of 0.025 in. These specimens were annealed at 538oC, 1000°F, for 10 min and quenched into ice water. They were then reheated to the creep test temperature of 260°C, 500°F, and held at this temperature for 18 min to al-low for vacancy condensation and loop formation. The specimens were subsequently creep tested at 1000 psi constant stress at 260°C. The temperature along the gage length was monitored by three Pt-Pt, 13 pct Rh thermocouples, two embedded in the upper and lower shoulders of the specimen and one attached to the mid-point of the gage length- The temperature gradient was held to ± l°C or less throughout the Creep test. Elongations were measured by two LVDT&apos;s &apos;connected in an averaging &apos;On-figuration, with quartz extension arms placed against the specimen at the extremes of the 2-in. gage length. i K 3J?--&apos; * *H - >^> - &apos; ^L r 77 . *y J OJL, J Fig. 1—Transmission electron micrograph of pure aluminum strained 8.2 pct at 260°C and 1000 psi constant stress. Note the presence of numerous loops near the light dislocation tangle

Jan 1, 1970

By M. R. Tek, F. F. Craig, J. O. Wilkes, C. S. Goddin

The waterflood performance of a water-wet, stratified system with crossflow is computed by a finite difference procedure. The effects of five dimensionless parameters on tile oil displacement efficiency, water saturation con-tour.7 and crossflow rates are evaluated in the absence of gravity forces. Crossflow due to viscous and capillary forces is shown to exert a significant effect on oil recovery in a field-scale model of a two-layered. water-wet sandstone reservoir. The crossflow is at a maximum in the vicinity of the front advancing in the more permeable layer. Under favorable mobility ratio conditions, the comparted oil recovery with crossflow always is interrnerliate between that predicted for a uniform reservoir and that for a layered reservoir with no crossflow. INTRODUCTION The important erects of reservoir heterogeneity on waterflood performance are commanding increased attention in the technical literature. Much of this attention is centered on two categories of layered reservoirs: those in which layers are non-communicating and those in which crossflow of fluids occurs between the layers. In the first category, the reservoir is assumed to consist of discrete layers, each uniform within itself and differing from the others only in such properties as thickness, porosity and absolute permeability. The performance within each layer is calculated by one-dimensional flow theory, and the performance of the total reservoir is obtained by summing individual layer performances. Capillary and gravity effects usually are not considered. Representative publications dealing with thi5 type of reservoir are those of Stiles,&apos; Dykstra and Parsons,&apos; Hiatt,3 Warren and Cosgrove&apos; and Higgins and Leighton." Prediction of performance for reservoirs in the second category is considerably more difficult since viscous, capillary and gravitational forces all play important roles in causing crossflow between layers. A number of authors have investigated the simpler problem of two-dimensional displacement flow in a stratified system with a mobility ratio of unity and negligible capillary and gravity effects.&apos;; Others have considered two-dimensional, non-steady-state flow of a single, slightly compressible fluid in a stratified reservoir. A limited number of laboratory oil displacement tests in layered models with crossflow have been reported. Miscible floods (with resultant zero capillary forces) in layered five-spot models were conducted by Dyes and Braun," who studied the effect of mobility ratio with zero gravity forces, and by Craig et al. 12 who studied the effect of gravity forces at constant mobility ratio. Waterfloods in layered five-spot models (with cross-tlow due to capillary, viscous and gravity forces) were conducted by Gaucher and Lindley,"&apos; who showed the effect of gravity forces in causing underrunning of the injected water and by Carpenter, Bail and Bobek, 14 who demonstrated the reliability of Rapoport&apos;s" dimension-less parameters for scaling layered systems. Waterfloods in rectangular layered models were conducted by Richardson and Perkins."&apos; who investigated the effect of velocity at constant mobility ratio and with zero gravity forces, and by Hutchinson," who studied the effects of varying mobility, layer permeability and layer thickness ratios. The differential equations which rigorously describe waterflooding in a heterogenous porous medium are non-linear and do not facilitate analytical solution. By using finite difference approximations it is possible to obtain a solution to any desired degree of accuracy. Such a solution, using an alternating direction implicit procedure (ADJP), is described by Douglas, Peaceman and Rachford.18 In the present study, a computer program using ADJP explores systematically the effects of important parameters on waterflood performance of a two-dimensional, two-layered, field-scale model of a water-wet sandstone system. Particular attention is given to evaluation of the water saturation contours and crossflow rates at the interface between layers to gain improved understanding of the crossflow mechanism. PROCEDURE BASIC: FLOW EQUATIONS The basic flow equations for two-dimensional, two-phase, immiscible, incon~pressible flow in a porous medium are:

Jan 1, 1967

By H. C. Brassfield, R. E. Fryxell, E. C. Duderstadt, E. A. Aitken, P. K. Conn

Permeation rates of hydrogen through are-cast tungsten were measured at temperatures from 1300" to 2600°C with hydrogen pressure differentials of 1 and 0.1 atm across isothermal membranes. Rates were nearly Proportional to square root of pressure, were inversely ProPortional to membrane thickness, and were temperature-dependent according to the equation: p = permeation coefficient in cu cm (stp)-mm per sq cm-min-atmlr2, R = gas constant in cal per mole-OK, T = absolute temperature in OK. Permeation rates of oxygen through tungsten were deternlined in the range 1965" to 2300°C and at calculated oxygen (02) pressures of 3.3 x 10-4 to 2.2 x at?n using sealed capsules containing Kmown quantities of tungsten oxide. Permeation coefficients obtained in this manner had a temperature dependency in the range of 40 to 45 kcal per mole of permeating species. Permeation coefficients for oxygen were higher than those for hydrogen in the temperature range studied. This is probably the result of oxygen solubility being much greater than hydrogen solubility at these temperatures at a gzven pressure. The purpose of this investigation was to determine experimentally the permeabilities of hydrogen and oxygen through polycrystalline tungsten at temperatures ranging from 1300° to 2600°C. Permeation rates of hydrogen through tungsten have been measured at temperatures up to 925°C,but no previous measurements have been reported for the temperature range covered in this investigation. Some diffusion measurements3&apos;&apos; and solubility measurements3 have been made for hydrogen at temperatures below 2200°C. No previous permeability studies of oxygen through tungsten are known to the authors; however, efforts have been made to determine diffusion rates5-&apos; and solubilities PERMEATION THEQRY Permeation has been defined by ortton" as the overall steady-state flow process from the gas phase on one side of a membrane or wall to the gas phase on the other side and includes the steps of adsorption, solution, diffusion, exsolution, and desorption; in a given temperature range, any of these processes may be rate-controlling. For diatomic gases, which dissociate prior to diffusion, the quantity of gas which permeates through a membrane or wall at a given temperature is directly proportional to the surface area of the membrane, the length of time of permeation, and the difference between the square roots of the pressures on each side of the membrane, and it is inversely proportional to the thickness of the membrane when diffusion-controlled, thus: q = total amount of gas permeating a membrane, p = temperature-de pendent permeation coefficient, A = area of membrane exposed, / = time, p i = gas pressure of high-pressure side, p2 = gas pressure of low-pressure side, d = thickness of membrane. When the surface reactions are not rate-controlling steps and the solution of the gas in the metal follows Henry&apos;s law, the permeation coefficient, P, is equal to the product of diffusivity and solubility. The permeation coefficient is also temperature-dependent according to the equation: Po = permeation constant, Q = activation energy, R = gas constant, T = absolute temperature. MATERIALS The tungsten materials used in this investigation were greater than 99.9 pct pure. The interstitial impurities ranged from < 1 to 50 pprn C, 7 to 35 ppm O, less than 15 pprn N, and undetectable quantities of hydrogen. Previous high-temperature studies of are-cast tungsten have shown that thermal treatment similar to that which the hydrogen permeation cells experienced reduces the carbon content from starting levels of 50 to 60 ppm to about 30 ppm. Metallic impurities were in the low (1 to 50) ppm range. The argon gas used in this investigation contained about 8 ppm N. None of the following gases were

Jan 1, 1968

By C. B. Manula, H. Gezik

In planning modem-day mining operations, management needs to pass from the area of subjective decision-making to an area of objective decision-makirlg. Planning procedures currently being practiced by many mine managers depend primarily on common sense and trial-and-error methods which oftentimes never reach optimum solutions. For this reason, an attempt is made to provide an understanding of linear programming to solve production planning problems. More specifically, a linear programming model is proposed to solve the production scheduling problem as it was found to exist in the crushed stone industry. Optimal production scheduling is of prime concern here, since the industry has to cope with seasonal demands for its products. The solution, as determined through the application of the model, provides management not only with an optimal production schedule but also with the following information regarding the improvement of existing facilities: (I) the location of bottlenecks which limit production and profit; (2) the directions in which the company i fixed equipment should be expanded or modified; and (3) the evaluation of unused resources. Future trends in stone production will quite likely parallel those in the construction industry. Continued record expenditures for structures such as dams, highways, and buildings are the reasons for this accelerated growth. Production forecasts for the next decade are predicted to be around one billion tons with an approximate 7% annual increase. The problem of increasing production from 570 million tons in 1966 to 1 billion tons by 1976 offers no small challenge. This situation is further compounded by other problems and economic factors. Population expansion, government regulation, labor agreement, and changes in technology wd cause managers and operators to be confronted by many more complexities than characterize present-day operations. Some of these. problems and economic factors are: (1) In some metropolitan and other urban areas, easily accessible sources of stone are almost depleted. Transportation costs will be a prime concern of the industry. (2) Population growth has forced residential expansion to a point where many quarry operators must conform to urban and suburban regulations. (3) Government regulation in regard to air and water pollution, noise, and reclamation as well as zoning developments necessitate that crushed stone producers keep abreast of changing conditions and requirements. (4) The crushed stone industry operates under pressure of a highly competitive situation. Every company&apos;s objective is to retain a favorable marketing position by exercising great attention to production and marketing costs. This is rather difficult, since the crushed stone industry is one segment of modern mining which has to cope with a seasonal demand for its products. Among the problems posed here, the one of exercising great attention to production and marketing costs is of current interest to management. Under seasonal demand conditions, this problem creates further difficulties for managers who need to plan and forecast. STATEMENT OF THE PROBLEM Various methods of Operations Research have been developed to handle a problem whose genesis is explained above. It is very seldom, however, that these quantitative means are used in practice by the mining industry. This paper describes the use of Operations Research methods to solve a production planning problem when demand is seasonal. The seasonal production planning problem is an important one because of the substantial costs associated with changing production, employment, and inventory levels to adjust to seasonal sales patterns. Most companies try to operate within reasonable levels of employment changes and hold inventories within storage capacities. In general, little or no attempts were made to develop minimum cosl schedules or to provide objective means to compare various production policies. When demand varies within a fixed time period, 11:-oduction requirements may be met by following one

Jan 1, 1970

By R. E. Collins

The existence of fluid migration across fixed boundaries in oil and gas reservoirs has been known for many years. Several techniques have been developed in the past for estimating The rate of migration across fixed boundaries as an aid in planninng field development and in the valuation of oilfied properties. The principal de/errc7nt to the use of these techniques lies in the rather extensive reservoir and field data required for ca1cula1ion.s of fluid migration. For this reason, a new, simplified procedure has been developed which makes possible the calculation of fluid migration with a minimum of field and resrvoir data. This new technique is based on certain. solutions of the differential equations describing flour in the reservoir which assume, for the portion of the reservoir of interest, that the formation can be approximated by a homogeneous rock of uniforrn thickness, that only a Single mobile fluid phase exists, and {hut fluid production at the well is solely a result of expun.rioti of rile resersoir fluds. The results of the, present work are compiled in a set of curves. These curves can he used to calculate both rule and cummulative fluid migration when the abovc assumptions ore ju.r/ifiril. The only data required for such calcu1ations are the production hisrories of all wells in the field, the permeability and porosity of the reservoir., the, compressibility and viscosity of the fluid. and the shape of the reservoir, INTRODUCTION The migration of fluids across fixed boundaries in oil and gas reservoirs has long been recognized as an engineering and economic problem, but no methods for estimating the extent of migration during the development stage of a reservoir have been published. The investigation reported here has been directed toward providing a simple technique for estimating the extent of such migration and has therefore been restricted to the simple case of a homogeneous reservoir of uniform thickncss containing compressible liquids. It is assumed that only a single mobile fluid phase exists and that fluid production at the wells is solely by expansion of the reservoir fluids. Obviously, few, if any. reservoirs conform to these assumptions during their entire productive life. On the other hand, most reservoirs approximate fluid expansion reservoirs during their initial stages of primary production. Consequently, the present work should be applicable during the field development of most reservoirs, regardless of whether or not they ultimately are gas cap, dissolved-gas, or water-drive reservoirs. The mathematical analysis of this problem has been outlined in the Appendix, while the results of the investigation are presented in the following discussion. DISCUSSION OF CALCULATION PROCEDURE SLNGLE Well, Ini-iimite Reservoir The procedure for calculating fluid migrations is based on solutions of the partial differentia1 equations describing the flow of a homogeneous compressible liquid in a homogeneous reservoir of uniform thickness. This differential equation has the form of the well-known diffusion equation. The equation, as applicable to our problem, is where p is fluid density, is porosity, EL is viscosity, c is compressibility,* and k is permeability; .r, y and t are the space coordinates and time. Derivation of this equation is given in the Appendix. The analytical solution of this equation for a single well of negligible radius producing at a constant rate. q, from a reservoir of infinite areal extent for a time, 1. can he readily obtained. From this solution the rate of flow and also the cumulative flow across a straight line segment of length, 1, at a perpendicular distance. 6. from the well can be computed. One end of the line must be on the perpendicular joining the well to the