Search Documents

By John Wulff

In a previous paper Burwell and Wulff1 have shown by electron diffraction studies that allotropic transformations can be induced in 18-8 stainless steel by polishing to a depth of about 5 X 10-5 cm. The development of a successful etching technique, which removes accountable amounts of metal from the surface without pitting and without changing the root-mean-square profile as measured with an Abbot pro-filometer, led to the application of the whole technique and material to an investigation of various mechanical methods of surface preparation, such as grinding, lapping, superfinishing and cold-rolling. Ideally suited to such a study is 18-8 of low carbon and nitrogen content. The material can be heat-treated to produce a coarse-grained single-phase structure. On cold deformation there occurs not only a reduction in grain size but, at regions of maximum stress, a change of face-centered material to body-centered material, provided the metal is worked below about 200°C. When worked between 200"about and 415°C., diffraction studies evidence a reduction in grain size without the appearance of any new body-centered ferrite. Ferrite, which can be referred also to strains induced by carbide precipitation, is evident when the material is worked at temperatures between 4150 and 800°C. Stainless steel is advantageous because the usual surface-finishing operations do not readily oxidize it to any appreciable depth. In order to produce detectable oxide films on the surface of austenitic 18-8 in a time equivalent to that required to produce, for example, a fine ground finish on 1/2 sq. in., the sample must be held in a furnace or oxygen stream above 600°C. By using longer heating times, it is possible to achieve measurable oxide films on samples , whose surface temperatures exceed only about 250°C. Experimental Procedure The apparatus, material, and technique described in the previous paper1 were employed in this investigation. Samples of 18-8, I by I cm. by 4 mm., were carefully annealed in the austenitic range to produce a coarse-grained structure. They were then quenched and pickled to remove scale, On a lapping machine of the Graton2 type, both sides of each sample were polished to a a of 1/6 a wave length of sodium D light. Thereafter they were electrolytically etched to remove at least 1 X 10-3 cm. of the surface, a depth that is well below that of any defermations induced by surface preparation. samples of this material, carefully mounted in steel holders, were then sent to various industrial organizations known for special finishing techniques. In the instrue ment shops at the Massachusetts Institute of Technology check samples were finished by the same methods, as far as they were available. Abbot profilometer3 readings Were made on all the samples, which were then carefully degreased in double distilled benzine and alcohol, dried in a hot-filtered air blast, and kept in a desiccator. SUC-

Jan 1, 1941

By John C. McDonald

This report is a continuation of an earlier one with a similar titlell to which the reader is referred for such details of procedure as do not appear here. A brief summary will be given of the objects and methods of the research at this time. It was pointed out in the previous paper that properties of magnesium alloys, as reported by different workers, often differ. This was thought to be due to differences in the preparation of the specimens tested. The ideal way to compare alloys would be to prepare each of them in so many ways that reasonable assurance would exist that the best properties of each alloy had been developed. The alloys would then be compared on the basis of these properties. The treatment actually given the alloys of these studies is necessarily a compromise with completeness. Nevertheless, enough treatments have been given so that a considerable degree of assurance exists that the properties obtained are representative of the alloy as such, and not merely of the manner in which it was worked and heat-treated. Procedure The alloys were cast in the form of cylindrical ingots, which were extruded into bar before being rolled. Some of the ingots were made in 2-in. diameter molds whereas all of the alloys of the previous paper were cast in 4.in. molds. For the larger ingots, the extruded bar was 2 by 0.140 in.; lor the 2-in. ingots, it was I by 0.100 in. Preliminary work indicated that strip rolled from both sizes of the bars would have substantially the same properties. The bar was rolled at some temperature in the range 400' to 800°F., for a total reduction of 50 per cent. The temperature and reduction per pass were so selected that no cold cracking occurred, but were such as to leave the metal in a semihard condition. The hot-rolled metal was annealed at several temperatures in the range 400' to 800°F.; it was also cold-rolled, and given solution and aging treatments. Two standard test bars corresponding to each treatment were pulled in tension, and the average was obtained. Twelve treatments, in all, were given each alloy. The properties of the annealed specimen having the highest percentage of elongation were selected to represent the alloy in the annealed condition. The properties of the cold-rolled specimen having the highest yield strength were selected to represent the alloy in the cold-rolled condition. These properties were then plotted against percentage composition by analysis. Smooth curves were drawn among the points, thus graphically averaging the results over the entire range studied. The percentage of elongation in the cold-rolled state was not plotted, since it was always low, and varied in an irregular manner with composition. Alloys Studied Magnesium-calcium Alloys.—Fig. I presents the results on the alloys of the rnag-

Jan 1, 1941

By John Wulff

In a previous paper Burwell and Wulff1 have shown by electron diffraction studies that allotropic transformations can be induced in 18-8 stainless steel by polishing to a depth of about 5 X 10-5 cm. The development of a successful etching technique, which removes accountable amounts of metal from the surface without pitting and without changing the root-mean-square profile as measured with an Abbot pro-filometer, led to the application of the whole technique and material to an investigation of various mechanical methods of surface preparation, such as grinding, lapping, superfinishing and cold-rolling. Ideally suited to such a study is 18-8 of low carbon and nitrogen content. The material can be heat-treated to produce a coarse-grained single-phase structure. On cold deformation there occurs not only a reduction in grain size but, at regions of maximum stress, a change of face-centered material to body-centered material, provided the metal is worked below about 200°C. When worked between 200"about and 415°C., diffraction studies evidence a reduction in grain size without the appearance of any new body-centered ferrite. Ferrite, which can be referred also to strains induced by carbide precipitation, is evident when the material is worked at temperatures between 4150 and 800°C. Stainless steel is advantageous because the usual surface-finishing operations do not readily oxidize it to any appreciable depth. In order to produce detectable oxide films on the surface of austenitic 18-8 in a time equivalent to that required to produce, for example, a fine ground finish on 1/2 sq. in., the sample must be held in a furnace or oxygen stream above 600°C. By using longer heating times, it is possible to achieve measurable oxide films on samples , whose surface temperatures exceed only about 250°C. Experimental Procedure The apparatus, material, and technique described in the previous paper1 were employed in this investigation. Samples of 18-8, I by I cm. by 4 mm., were carefully annealed in the austenitic range to produce a coarse-grained structure. They were then quenched and pickled to remove scale, On a lapping machine of the Graton2 type, both sides of each sample were polished to a a of 1/6 a wave length of sodium D light. Thereafter they were electrolytically etched to remove at least 1 X 10-3 cm. of the surface, a depth that is well below that of any defermations induced by surface preparation. samples of this material, carefully mounted in steel holders, were then sent to various industrial organizations known for special finishing techniques. In the instrue ment shops at the Massachusetts Institute of Technology check samples were finished by the same methods, as far as they were available. Abbot profilometer3 readings Were made on all the samples, which were then carefully degreased in double distilled benzine and alcohol, dried in a hot-filtered air blast, and kept in a desiccator. SUC-

Jan 1, 1941

By D. T. Peters, S. Floreen

The age hardening behavior of an Fe-8Ni-13Mo alloy was studied after the matrix had been varied to produce either ferrite, cold u~orked ferrite, or nzassive nzartensite. The aging behavior of the cold worked ferrite and murtensite structures were very similar. The martensite aging kinetics were much different from those observed in earlier studies of aging of maraging steels, even though the martensite wzatri.r had the same dislocation structure as those found in maraging steels. The results suggest that the previously observed precipitation kinetics of maraging steels ?nay have been controlled by the nucleation be-haviov, which in turn were dictated by the alloy compositions and the resultant identities of the precipitating phases. IT is well known that the rate of precipitation from solid solution depends not only on the degree of super-saturation, but also on the density and distribution of dislocations in the matrix structure. These imperfections often act as nucleation sites, and may also enhance atomic mobility. 'Thus, the presence of dislocations is important since the type and distribution of precipitates may be determined by them. The precipitate density and morphology in turn affects the mechanical properties of the alloy. A number of studies have been devoted to the precipitation characteristics in various types of maraging steels.'-" These are iron-base alloys containing 10 to 25 pct Ni along with other substitutional elements such as Mo, Ti, Al, and so forth, that are used to produce age hardening. The carbon contents of these steels are quite low, and carbide precipitation is not believed to play any significant role in the aging reactions. After solution annealing and cooling these alloys generally transfclrm to a bcc lath or massive martensite structure characterized by elongated martensite platelets that are separated from each other by low angle boundaries, and that contain a very high dislocation den~it~.~~~~~~~~-~~ Age hardening is then conducted at temperatures on the order of 800" to 1000°F to produce substitutional element precipitation within the massive martensite matrix. Most of the aging studies to date have revealed several common traits in these alloys, regardless of the particular identity of the precipitation elements. Generally hardening has been found to be extremely rapid, with incubation times that approach zero. The agng kinetics, at least up to the time when reversion of the martensite matrix to austenite begins to predominate, frequently follow a AX/~~ = ktn type law, where x is hardness or electrical resistivity, t is the time, and k and n are constants. The values of n are frequently on the order of 0.2 to 0.5, which are well below the idealized values of n based on diffusion controlled precipitate growth models. Finally, the observed activation energy values are typically on the order of 30 kcal per mole, and thus are well below the nominal value of about 60 kcal per mole found for substitutional element diffusion in ferrite. The common explanation of these observations is that the precipitation kinetics are controlled by the massive martensite matrix structure. Thus, the absence of any noticeable incubation time has been attributed, after ~ahn," to the fact that the precipitate nucleation on dislocations may occur without a finite activation energy barrier. The low values of the activation energy are generally assumed to be due to enhanced diffusivity in the highly faulted structure. If this explanation that the precipitation kinetics are dominated by the matrix structure is correct then one should observe a distinct difference in lunetics between aging in a martensitic matrix and aging the same alloy when it has a ferritic matrix. Such a comparison cannot be made with conventional maraging compositions, but could be made with the alloy used in the present study. In addition, the ferritic structure of the present alloy could be cold worked to produce a high dislocation density so that one could determine whether ferrite in this condition would age similarly to martensite. EXPERIMENTAL PROCEDURE The composition of the alloy used in this study was 8.1 pct Ni, 13.0 pct Mo, 0.10 pct Al, 0.13 pct Ti, 0.012 pct C, bal Fe. The alloy was prepared as a 40 lb vacuum induction melt. The heat was homogenized and hot forged at 2100°F to 2 by 2 in. bar, and then hot rolled at 1900°F to $ in. bar stock. The aging lunetics were followed by Rockwell C hardness and electrical resistivity measurements. Samples for hardness testing were prepared as small strips approximately 2 by $ by 4 in. thick. Electrical resistivity was studied on cylindrical samples approximately 2 in. long by 0.1 in. diam. The method for making the alloy either martensitic or ferritic was based on the fact that the alloy showed a closed y loop type of phase diagram. At high temperatures, above approximately 24003F, the alloy was entirely ferritic. Small samples on the order of the dimensions described above remained entirely ferritic after iced-brine quenching from this temperature. In practice, a heat treatment of 1 hr in an inert atmosphere at 2500°F followed by water quenching was used to produce the ferritic microstructure. These samples were quite coarse grained and usually en-

Jan 1, 1970

By John C. Calhoun, R. C. Newman

A report is made of experimental work performed on the capillary retention of water within porous solid systems, the displacement being accomplished with air and various organic liquids. A portion of the experiments were designed to measure the lowering of vapor pressure of water within a porous solid with subsequent conversion of such vapor pressure data to capillary pressure values. The form of the capillary pressure curve at high capillary pressures has been elucidated from this data. Certain theoretical approaches are presented to indicate the correlation between work done by previous investigators. The present work is correlated also with theory. Surface area values are calculated for the various core systems studied. INTRODUCTION Leverettl in 1941 presented a paper which gave the essential concepts of capillary behavior in porous solids. His presentation included both theoretical and experimental aspects of the problem. He defined the term "capillary pressure" and applied it to an ideal porous system. His work included the definition of a dimensionless quantity which was a function of fluid saturation and which could be correlated with porosity and permeability for clean, unconsolidated sands. His experimental work was performed by the drainage of water from packed columns of unconsolidated sands. Subsequently, other authors have presented experimental techniques which permit the determination of capillary pressure data on small core samples. Data have been presented on porous systems other than unconsolidated sand and with the use of liquids other than water. Notable among these are the works of Hassler, Brunner and Deahl,2 Bruce and Welge,³ Amyx and Yuster, and Purcell. Others have indicated the applicability of such data to field problems without regard to the theoretical aspects of the capillary pressure saturation function. Within the recent published literature there have also been reported studies of the correlations between capillary pressure data and other fundamental properties of porous solids. It was the intent of the present reported work to investigate the behavior noted by water and various other liquids within porous systems in an attempt to amplify existing correlations between capillary properties, surface properties, and other fundamental characteristics of porous systems. The experimental work described here was performed as two different investigations but due to its nature it is reported as separate phases of the same general topic. CAPILLARY PRESSURE AND VAPOR PRESSURE LOWERING It has long been recognized that the vapor pressure above the curved surface of a liquid is a function of the curvature of the liquid surface." The capillary pressure is also a function of the curvature of the liquid surface. Both capillary pressure values and vapor pressure lowering are, therefore, functions of the liquid saturation of a porous solid. A quantitative relationship between the vapor pressure lowering and the capillary pressure can be simply developed by considering a porous solid system containing water in equilibrium with its vapor. The capillary pressure at any height within the system is defined as Pc= Pv -P1.......(1) where PI is the pressure in the vapor phase and PI is the pressure in the liquid phase. Consider that both the water and its vapor are continuous phases within the system and that at any point the two phases are in equilibrium. The gradients in the liquid and vapor columns within the porous system are: PI = plgdh dPv = pvgdh By integrating these gradients from the point where h = 0 to the height h, values of PI and PI at the point h are found. Let the point where h = 0 be that where Pc = 0. At this point also Pv is equal to Pve, the equilibrium vapor pressure of the liquid above a flat surface. From the first of these integrations since p is independent of height: P = p1gh + P, From the second, by substituting pv = PvM/RT, where M is the molecular RT weight: Pv M Pvo- RT gh By combining the two last equations to eliminate h, one obtains: P1 RT Pv Substituting, values for Pv and Pi in equation (1) gives PyRT RT Pv which is approximately Pc=RT/Mp1 In Pvo/Pv .......(2) This equation relates the vapor pressure above a curved surface with the capillary pressure across the curved surface in terms of measurable quantities. The equivalent of this equation is given by Freundlich."

Jan 1, 1949

By John C. Calhoun, R. C. Newman

A report is made of experimental work performed on the capillary retention of water within porous solid systems, the displacement being accomplished with air and various organic liquids. A portion of the experiments were designed to measure the lowering of vapor pressure of water within a porous solid with subsequent conversion of such vapor pressure data to capillary pressure values. The form of the capillary pressure curve at high capillary pressures has been elucidated from this data. Certain theoretical approaches are presented to indicate the correlation between work done by previous investigators. The present work is correlated also with theory. Surface area values are calculated for the various core systems studied. INTRODUCTION Leverettl in 1941 presented a paper which gave the essential concepts of capillary behavior in porous solids. His presentation included both theoretical and experimental aspects of the problem. He defined the term "capillary pressure" and applied it to an ideal porous system. His work included the definition of a dimensionless quantity which was a function of fluid saturation and which could be correlated with porosity and permeability for clean, unconsolidated sands. His experimental work was performed by the drainage of water from packed columns of unconsolidated sands. Subsequently, other authors have presented experimental techniques which permit the determination of capillary pressure data on small core samples. Data have been presented on porous systems other than unconsolidated sand and with the use of liquids other than water. Notable among these are the works of Hassler, Brunner and Deahl,2 Bruce and Welge,³ Amyx and Yuster, and Purcell. Others have indicated the applicability of such data to field problems without regard to the theoretical aspects of the capillary pressure saturation function. Within the recent published literature there have also been reported studies of the correlations between capillary pressure data and other fundamental properties of porous solids. It was the intent of the present reported work to investigate the behavior noted by water and various other liquids within porous systems in an attempt to amplify existing correlations between capillary properties, surface properties, and other fundamental characteristics of porous systems. The experimental work described here was performed as two different investigations but due to its nature it is reported as separate phases of the same general topic. CAPILLARY PRESSURE AND VAPOR PRESSURE LOWERING It has long been recognized that the vapor pressure above the curved surface of a liquid is a function of the curvature of the liquid surface." The capillary pressure is also a function of the curvature of the liquid surface. Both capillary pressure values and vapor pressure lowering are, therefore, functions of the liquid saturation of a porous solid. A quantitative relationship between the vapor pressure lowering and the capillary pressure can be simply developed by considering a porous solid system containing water in equilibrium with its vapor. The capillary pressure at any height within the system is defined as Pc= Pv -P1.......(1) where PI is the pressure in the vapor phase and PI is the pressure in the liquid phase. Consider that both the water and its vapor are continuous phases within the system and that at any point the two phases are in equilibrium. The gradients in the liquid and vapor columns within the porous system are: PI = plgdh dPv = pvgdh By integrating these gradients from the point where h = 0 to the height h, values of PI and PI at the point h are found. Let the point where h = 0 be that where Pc = 0. At this point also Pv is equal to Pve, the equilibrium vapor pressure of the liquid above a flat surface. From the first of these integrations since p is independent of height: P = p1gh + P, From the second, by substituting pv = PvM/RT, where M is the molecular RT weight: Pv M Pvo- RT gh By combining the two last equations to eliminate h, one obtains: P1 RT Pv Substituting, values for Pv and Pi in equation (1) gives PyRT RT Pv which is approximately Pc=RT/Mp1 In Pvo/Pv .......(2) This equation relates the vapor pressure above a curved surface with the capillary pressure across the curved surface in terms of measurable quantities. The equivalent of this equation is given by Freundlich."

Jan 1, 1949

By J. E. Warren, J. H. Hartsock

Because of the extensive utilization of hydraulic fracturing for the stimulation of low-productivity wells, the two related problems of fracture design and evaluation have become economically significant and, as a consequence, helve motivated this investigation. The producing characteristics of horizontally fractured wells were studied to determine the fracture configuration that should be employed as the basis for the design of the treatment and to develop a method that can be used to establish the degree to which the design objectives have been achieved. The equations which describe the steady-state flow of a single-phase fluid into, and through. a finite-capacity fracture were solved numerically for an idealized reservoir-fracture model. The numerical results were used to obtain an apparent skin effect for each combination of the parameters considered. Based on the computed results, subject to the limitcitions implied by the assumptions that were made, the following general conclusions were drawn. 1. For a radius of drainage at least four times us large as the radius of the fracture, an apparent skin effect that is independent of the radius of drainage can be calculated. 2. The productivity of the hydraulically fractured system, relative to that of the unfructured well, can be determined from the apparent skin effect and can be used to establish design objectives. 3. In the evaluation of a fracture job, it is not Possible to determine both the radius of the fracture and its flow capacity uniquely from the apparent skin effect; an independent determination of one of the quantities is necessary. lNTRODUCTION Although hydraulic fracturing has been employed as a method for stimulating the productivity of literally hundreds of thousands of wells during the past 10 years, it is only in the last few years that improvements in fracture design1-6 and fracturing technique' have combined to increase the probability of obtaining a successful treatment to such an extent that the mechanics of the method may be considered to be standardized. From an economic point of view, however, two related questions must be satisfactorily answered before hydraulic fracturing can be used in the most profitable manner. The two questions are the following. I. For a particular well in a given formation, what are the optimum design specifications for the fracture treatment? 2. Have the design objectives been achieved by the fracture treatment? The significance of these questions has been recoguized, and some attempts to obtain answers have been made. Howard, et al,8 endeavored to determine the optimum treatment, based on maximizing profits, for any given formation; unfortunately, this work was based on a crude method for approximating the productivity of a well. Carter and Tracy9 utilized the same approximation to study the effect of fracturing on the behavior of a well producing by virtue of a solution-gas drive. Electrolytic models were used by van Poollen10 to investigate the variation in productivity due to fracturing; however, only a limited number of results were presented. Later, from the same model results, van Poollen, et al,11 attempted to justify an approximate expression for determining the productivity of a fractured well. It is quite apparent that there is a definite lack of the practical information necessary for specifying the optimum fracture configuration to be considered for design purposes. The only detailed attempt to develop a procedure for evaluating the result of a given fracture treatment appears to be that of van Dam and Horner.12 These authors described a technique for analyzing pressure build-down data, obtained immediately after fracturing, to determine the final fracture volume, the final fracture porosity, the fracture area, the fracture thickness and the in situ fluid loss of the fracturing fluid. While this approach should be useful whenever acceptable pressure measurements are available, it does not yield a value for the flow capacity of the fracture. Since the problems of fracture design and evaluation are inversely related, it should be sufficient to study the effect of the fracture configuration on the performance of a well. The primary objective of this investigation is to evolve a technique for computing the desired solutions. The secondary objective is to analyze these computed results in order to prescribe a method for evaluating fracture treatments. Because this study is exploratory in nature, its scope

By D. A. Koss, R. O. Gordon

Changes in the ultrasonic attenuation in steel specimens have been observed during tensile tests. Samples of AISI 1020, 1045, and 1095 steel quenched and tempered to a spheroidized condition have been used. Both attenuation and microstrain measurements fail to reveal preyield dislocation motion in the as-tempered specimens. Upon reapplication of the load to samples in the strain-aged condition there is a large preyield attenuation increase but no observable micvostrain up to 90 pct of the flow stress. These and other observed changes can be related to changes in the free disloca-tion density and the average dislocation segment length. The results show that the increase inflow stress in the strain-aged condition is caused by increased fric-tional stresses rather than by dislocation pinning. Strain aging is found to proceed through three distinct stages. The second of these follows the t2,3 law and is characterized by an activation energy of 20 kcal per mole. RECENT years have seen substantial progress in the understanding of the sharp yield-point phenomenon in metals; the influence of solute atom pinning of dislocations, of dislocation multiplication, and of the stress dependence of dislocation velocity on the stress-strain curve has been investigated. Evaluation of the relative importance of these factors in the various metals and alloys is not complete, however, because observation of the stress-strain curve alone is not sufficient to determine the yield mechanism. In the experimental study of the sharp yield point it would be desirable to be able to directly observe the dislocations in a bulk specimen. Ideally, the observations should be made during deformation because of the rapidity of recovery processes in lightly strained metals. While direct observation of dislocations under these conditions is not possible, certain physical properties sensitive to the dislocation structure can be observed while a test specimen is undergoing deformation. Of these properties, one of the most useful is the ultrasonic attenuation. The total observed attenuation, a, is due to a number of causes: acoustic diffraction, scattering by inhomogeneities in the structure, heat flow, magnetostriction, and dislocation damping. Satisfactory measurements of the absolute a are difficult to make because all of these factors cannot be evaluated in any given specimen, but changes in a due to changes in the dislocation contribution alone can be reliably observed by appropriate technique. Because dislocation damping is sensitive to small concentrations of long, free dislocation segments, observation of acoustic attenuation due to dislocations is particularly useful in the study of the early stages of plastic deformation. In the experiments to be described, acoustic-attenuation measurements have been made on steel samples as these were pulled through their yield in tension. While many details of dislocation-damping phenomena in the bee metals remain obscure, the main characteristics of this type of internal friction are established. At low strain amplitudes and at frequencies well below that of dislocation resonance (the conditions of these experiments) the dislocation contribution to a is given by1 where A is the density of freely oscillating dislocation segments, L is the average segment length, B is the damping constant for a moving dislocation, w is the frequency, O is an orientation factor, C is the dislocation line tension, and t is a numerical factor relating average segment length to the distribution of lengths. Because of its fourth-power dependence on L, a is highly sensitive to changes in dislocation free segment length such as might arise from unpinning, recovery, jog formation, or dislocation intersection. The above equation is not valid near the resonant frequency of dislocations, but the segment length in steel is believed to be much too short to cause over-damped resonance near 10 Mc; this will be discussed later. In the present experiments plastic-strain measurements with a sensitivity of 2 x 10-5 were also made during the early stages of loading in order to detect preyield strain. The general plan of the experiments is to observe

Jan 1, 1967

By Frederick C. Langenberg

IT has been shown1 that copper can be removed from iron-base alloys by solvent extraction with molten lead and sodium sulfide slags. In these processes copper removal is favored by a positive deviation of copper in iron from ideal behavior. Chou' determined the distribution of copper between liquid silver and iron and obtained a value for the activity coefficient of copper in iron at 1600°C, y°, of 9.12. Chipman" calculated the activity coefficient from the Fe-Cu phase diagram to be 12 at 1600°C. More recent work by Koros and Chipman" yields a value at 1600°C for y V equal to 8.0. In this latter work, the distribution of copper between silver and iron was measured at 1600°C, while the activity of copper in the silver phase was determined from calculations based on the Ag-Cu phase diagram. The present work presents another value of based on calculations from the Cu-Pb phase diagram, on the heat of mixing data for the Cu-Pb system, and on an experimental study of the distribution of copper between iron and lead at 1600°C. These experiments were part of a larger research on the removal of copper from iron-base alloys conducted at The Pennsylvania State University. Experimental Procedure—The distribution of copper between pure iron and lead at 1600°C was made in the furnace assembly shown in Fig. 1. The melts were made in high purity alumina crucibles under an atmosphere of prepurified argon. Temperatures were read with an optical pyrometer sighted into a %-in. hole drilled in the inductor. Extensive calibration of the optical readings with immersion readings had been carried out previously. The initial charges consisted of 100 g of electrolytic iron and varying weights of copper. The melts were brought rapidly to 1600°C and 75 g of lead added. The temperature was quickly readjusted to 1600°C, and the melts were held for 10 min at this temperature. A sample was then drawn from the iron layer in a fused silica tube. This sample was analyzed to determine the copper content of the iron layer at equilibrium, while the copper content of the lead layer was determined from a material balance. The validity of the material balance was confirmed by a chemical analysis of the entire lead layer for one heat in which the crucible was removed and quenched after the iron sample was taken. The calculated copper content of the lead-rich phase was 1.18 pct; while the entire lead layer analyzed 1.20 pct Cu. The distribution results are shown in Fig. 2. The distribution constant, N,.'r'"/N, is equal to 3.95. Discussion—From the Cu-Pb phase diagram' it can be seen that lead-rich liquid Pb-Cu solutions

Jan 1, 1957

By Frank G. Miller

This report presents developments of fundamental equations for describing the flow and thermodynamic behavior of two-phase single-component fluids moving under steady conditions through porous media. Many of the theoretical considerations upon which these equations are premised have received little or no attention in oil-reservoir fluid-flow research. The significance of the underlying flow theory in oil-producing operations is indicated. In particular, the theoretical analysis pertains to the steady, adiabatic, macroscopically linear, two-phase flow of a single-component fluid through a horizontal column of porous medium. It is considered that the test fluid enters the upstream end of the column while entirely in the liquid state, moves downstream an appreciable distance, begins to vaporize, and then moves through the remainder of the column as a gas-liquid mixture. The problem posed is to find the total weight rate of flow and the pressure distribution along the column for a given inlet pressure and temperature, a given exit pres5ure or temperature and given characteristics of the test fluid and porous medium. In developing the theory, gas-liquid interfacial phenomena are treated. phase equilibrium is assumed and previous theoretical work of other investigators of the problem is modified. Laboratory experiments performed with specially designed apparatus. in which propane is used as the test fluid, substantiate the theory. The apparatus. materials and experimental procedure are described. Comparative experimental and theoretical results are presented and discussed. It is believed that the research findings contributed in this * paper should not only lead to a better understanding of oil-reservoir behavior, but also should be suggective in regard to future research in this field of study. INTRODUCTION In recent years much time and effort has been consumed in both theoretical and experimental studies of the static and . dvnamic behavior of oil-reservoir fluids in porous rocks. Although lack of sufficient basic oil-field data, principally concerning the properties and characteristics of reservoir rocks and fluids, largely precludes quantitative application of research results to oil-field problems, qualitative application has become common practice. In effect. oil-reservoir engineering research is serving as a firm foundation for oil-field development and production practices leading to increased economic recoveries of petroleum. This province of research. however, still poses many perplexing problems. The thermodynamic behavior of two-phase fluids moving through porous media constitutes one facet of reservoir-fluid-flow research that has not received the attention it deserves. This report embodies a theoretical discussion of this subject and a description of a series of related laboratory experiments. The significance of the problem to oil field operations is indicated but in articular the report centers around a theory and method for analyzing the steady. macroscopically linear, two-phase flow of a fluid (a single molecular species) through a horizontal column of porous medium. For simplicity in showing how the thermodynamic behavior of two-phase fluids moving through porous media affects oil-reservoir performance problems, attention is focused temporarily on a particular well producing petroleum from an idealized water-free solution-gas drive reservoir, the reservoir rock being a horizontal, thin, fairly homogeneous sandstone of large areal extent confined between two impermeable strata. The flowing hydrocarbon fluid is considered to exist entirely as a Iiquid at points in the reservoir remote from the well; however. the decline in fluid pressure in the direction of the well causes vaporization of the hydrocarbon to begin at a radial distance r from the well. Upstream from r the fluid moves entirely as a liquid and downstream from r it moves either entirely as a gas or as a gas-liquid mixture depending on the properties of the hydrocarbon and on the thermodynamic process it follows during flow. The distance r would be variable under transient flow conditions. but for purposes of analysis the flow is considered to l~e steady at the particular instant of observation during the flowing life of the well of interest. If the flow were isothermal and the hydrocarbon a pure substance, the fluid would be entirely gaseous downstream from r. Thus, this isothermal flow process for a pure substance would require that the heat of vaporization be supplied at r. over zero length of porous medium, at the precise rate necessary to maintain the constant temperature. This means that the solid matrix of the porous medium (reservoir rock) and the surroundings (impermeable strata confining the reservoir rock) would have to serve as infinite heat sources. Heat-transfer requirements would be somewhat less severe for the isothermal flow of a multicorn-ponent hydrocarbon as bubble and dew points at the same temperature correspond to different pressures. In this instance isothermal conditions would be sustained without complete vaporization of the fluid over zero length of porous medium. Nevertheless. as the flow is in the direction of decreasing

Jan 1, 1951

By Frank G. Miller

This report presents developments of fundamental equations for describing the flow and thermodynamic behavior of two-phase single-component fluids moving under steady conditions through porous media. Many of the theoretical considerations upon which these equations are premised have received little or no attention in oil-reservoir fluid-flow research. The significance of the underlying flow theory in oil-producing operations is indicated. In particular, the theoretical analysis pertains to the steady, adiabatic, macroscopically linear, two-phase flow of a single-component fluid through a horizontal column of porous medium. It is considered that the test fluid enters the upstream end of the column while entirely in the liquid state, moves downstream an appreciable distance, begins to vaporize, and then moves through the remainder of the column as a gas-liquid mixture. The problem posed is to find the total weight rate of flow and the pressure distribution along the column for a given inlet pressure and temperature, a given exit pres5ure or temperature and given characteristics of the test fluid and porous medium. In developing the theory, gas-liquid interfacial phenomena are treated. phase equilibrium is assumed and previous theoretical work of other investigators of the problem is modified. Laboratory experiments performed with specially designed apparatus. in which propane is used as the test fluid, substantiate the theory. The apparatus. materials and experimental procedure are described. Comparative experimental and theoretical results are presented and discussed. It is believed that the research findings contributed in this * paper should not only lead to a better understanding of oil-reservoir behavior, but also should be suggective in regard to future research in this field of study. INTRODUCTION In recent years much time and effort has been consumed in both theoretical and experimental studies of the static and . dvnamic behavior of oil-reservoir fluids in porous rocks. Although lack of sufficient basic oil-field data, principally concerning the properties and characteristics of reservoir rocks and fluids, largely precludes quantitative application of research results to oil-field problems, qualitative application has become common practice. In effect. oil-reservoir engineering research is serving as a firm foundation for oil-field development and production practices leading to increased economic recoveries of petroleum. This province of research. however, still poses many perplexing problems. The thermodynamic behavior of two-phase fluids moving through porous media constitutes one facet of reservoir-fluid-flow research that has not received the attention it deserves. This report embodies a theoretical discussion of this subject and a description of a series of related laboratory experiments. The significance of the problem to oil field operations is indicated but in articular the report centers around a theory and method for analyzing the steady. macroscopically linear, two-phase flow of a fluid (a single molecular species) through a horizontal column of porous medium. For simplicity in showing how the thermodynamic behavior of two-phase fluids moving through porous media affects oil-reservoir performance problems, attention is focused temporarily on a particular well producing petroleum from an idealized water-free solution-gas drive reservoir, the reservoir rock being a horizontal, thin, fairly homogeneous sandstone of large areal extent confined between two impermeable strata. The flowing hydrocarbon fluid is considered to exist entirely as a Iiquid at points in the reservoir remote from the well; however. the decline in fluid pressure in the direction of the well causes vaporization of the hydrocarbon to begin at a radial distance r from the well. Upstream from r the fluid moves entirely as a liquid and downstream from r it moves either entirely as a gas or as a gas-liquid mixture depending on the properties of the hydrocarbon and on the thermodynamic process it follows during flow. The distance r would be variable under transient flow conditions. but for purposes of analysis the flow is considered to l~e steady at the particular instant of observation during the flowing life of the well of interest. If the flow were isothermal and the hydrocarbon a pure substance, the fluid would be entirely gaseous downstream from r. Thus, this isothermal flow process for a pure substance would require that the heat of vaporization be supplied at r. over zero length of porous medium, at the precise rate necessary to maintain the constant temperature. This means that the solid matrix of the porous medium (reservoir rock) and the surroundings (impermeable strata confining the reservoir rock) would have to serve as infinite heat sources. Heat-transfer requirements would be somewhat less severe for the isothermal flow of a multicorn-ponent hydrocarbon as bubble and dew points at the same temperature correspond to different pressures. In this instance isothermal conditions would be sustained without complete vaporization of the fluid over zero length of porous medium. Nevertheless. as the flow is in the direction of decreasing

Jan 1, 1951

By Elmer R. Kaiser

COMBUSTION of coal and transfer of heat from flames and gases to boiler surfaces continue to be of great interest to engineers here and abroad. Numerous investigations have been in progress to improve furnace and boiler performance and economy. The importance of better understanding of the processes and opportunities for improvement is apparent when it is remembered that heat from at least 500 million tons of coal a year the world over is being transferred to boiler water at efficiencies ranging mostly between 50 and 90 pct. Even slight gains in efficiency, economy, and labor saving become very significant when multiplied by the enormous quantity of fuel consumed. Also the competitive position of the large coal, oil, and gas industries in satisfying the fuel consumers is greatly affected by the achievements made through technical progress with each fuel. This paper is part of a continuing activity of Bituminous Coal Research, Inc., to extend the knowledge of coal utilization for steam generation and to seek promising directions for future research and development in cooperation with others. Particularly in the latter regard, numerous interviews were held during the last three years to seek the experience and advice of boiler and combustion-equipment manufacturers, electric-utility executives, and fuel engineers. A wealth of published information was also reviewed, which together with the interviews pointed to the advisability of further work on ash and sulphur control. For the present purpose a number of factors important to efficient heat liberation and recovery have been grouped as follows: 1—combustion, temperatures, and rates of heat liberation; 2—radiation, convection, and furnace and boiler configuration; 3—sponge ash, slag, and hard-bonded deposits; 4— low-temperature deposits and corrosion (cooling flue gas below dew point and air-pollution control); 5—the limitations of coal cleaning and boiler size and cost as related to fuel characteristics; 6—future possibilities and conclusions. The development of combustion apparatus for power boilers is progressing at a lively pace. There has been no letup in improvements in design of pulverized-coal-fired boilers, and there is a strong trend at present toward improving dry-bottom units. Spreader stokers with overfire jets and dust collectors as standard equipment are gaining favor. Less than 10 years in commercial use, cyclone burners are going into numerous installations here' and abroad.' Underfeed and traveling-grate stokers have long since been developed for heavy-duty operation, yet new developments in overfire jets and humidification of air blast have improved their performance. A water-cooled vibrating-grate stoker of German origin is being introduced into the United States and Canada." The primary objectives of an ideal coal combustion device are: capacity to burn the variety and sizes of coals likely to be economically available during the life of the unit; capacity to burn the coals automatically for a wide load range and rapid load fluctuations and to burn the coals completely to CO2, H2O, and SO2, which means without smoke and cinders, or carbon in the refuse; capacity to control and discharge all the ash in final granular form without ash adhesion to walls or tubes, and without flue dust; minimum furnace volume; minimum labor and maintenance; low initial and operating cost. Regardless of the method of burning, the gaseous products of coal combustion are N2, CO2, O2, H20, and SO?. By way of illustration, the coal analyses in Table I is assumed from an installation described by E. McCarthy.' When coal is burned with 20 pct excess air (theoretical air, 9.23 lb per lb of coal), the quantities of combustion gas shown in Table II are produced. In addition, the gases carry particles of fly ash, unconsumed cinders, soot particles, and small but significant amounts of vaporized oxides and sulphates of sodium, potassium, lithium, phosghorous, iron, and other metals. In recent years, germanium, one of the rare metals found in coal, has been shown to oxidize and vaporize at combustion temperatures and to be concentrated by reconden-sation at lower temperatures." Pulverized coal and cyclone flames" have peak temperatures of 3000' to 3500°F. Temperatures in fuel beds of large underfeed stokers reach maxima of 3000°F, sufficient to fuse almost any ash and to volatilize some of it. These peak temperatures are above the optimum necessary for rapid combustion, but they hasten heat transfer for ignition as well as boiler heat absorption. Furnace and gas temperatures increase with combustion air preheat. Low excess air has the same effect. Fine coal pulverization and highly turbulent combustion shorten the distance for fuel burnout, increase flame temperature, and speed up heat transfer. Rates of combustion of pulverized coal exceeding 200,000 Btu per cu ft per hr have been demonstrated in atmospheric gas-turbine combusters,

Jan 1, 1955

By B. K. Larkin, H. R. Bailey

One of the most complicated and potentially one of the most promising secondary recovery methods is that of underground combustion. A number of field tests1,2,3 have been performed recently, apparently indicating that the underground combustion process is technically feasible. Relatively little has been published concerning laboratory or theoretical investigations of the underground combustion process. A few results of laboratory tube runs have been published4,5 and these data are valuable in understanding the process. The tube run data also emphasize the complexity of the thermal recovery process. Under various conditions, problems in heat transfer, multiphase fluid flow and combustion rates arise. For example, at low temperatures the combustion process proceeds so slowly that the front cannot be propagated. It is a problem in heat transfer to predict what injection rate of air is necessary to sustain combustion in a particular reservoir. Certain fluid flow aspects of underground combustion have been studied by Wilson, et a13. In the development of some of the flow equations, assumptions of adiabatic combustion were made. Vogel and Krueger1 have devised an analog model of the heat conduction process for a constant temperature cylindrical source moving radially outward in a homogeneous conducting medium. Jenkins and Ramey7, in a discussion of the Vogel and Krueger paper, present an analytical solution of a heat conduction problem related to problems considered by Vogel and Krueger. In the present paper an idealized model of the heat conduction problem is considered, which is more general than the model of Vogel and Krueger or Jenkins and Ramey. In particular, initial well heating, vertical heat losses and arbitrary frontal velocities are included. Equations for the temperature are obtained in terms of integrals. These integrals have been evaluated numerically for the case of frontal movement corresponding to a constant air injection rate and the case of constant frontal velocity. These results are presented graphically. From these graphs one may estimate conditions under which nearly adiabatic temperatures are obtained. The application of these results requires some knowledge of reservoir thermal properties, combustion rates and residual fuel content. Assuming no vertical losses it is possible to obtain explicit evaluations of these integrals in the linear con- stant velocity case and in the radial case corresponding to a constant air injection rate. PRESENTATION OF RESULTS The geometry of the radial system is illustrated in Fig. 1 where the reservoir is represented as a slab of thickness h. By making an energy balance for an infinitesimal volume, one obtains the following equation. T= x T+F(z.r).........(1) where F (z,r) =HWv1/pc (r-r1), if z h/2 = 0, if z >h/2 . In Eq. 1, V- represents the Laplacian operator and 8 is the Dirac delta function. Eq. 1 is the well-known diff-sion equation with a source term and may be applied to a linear system. For example, if the source is moving in the direction of the x axis, then the argument of the delta function is x - 1,. In deriving Eq. 1 the following approximations were made concerning the actual combustion process: (1) all heat transfer is by conduction in a homogeneous medium, (2) the formations bounding the reservoir have the same thermal properties as the reservoir, (3) the combustion fuel is consumed instantaneously as the front arrives, (4) there is a constant fuel concentration

By Leonard A. Wood, A. G. Winslow

Subsidence has occurred in several areas of the upper Gulf Coast region of Texas, although in most cases this is not evident without precise instrumental leveling. As referred to in this report, the region com- prises about 7200 sq miles of the Gulf Coastal Plain and includes all or parts of the following counties: Harris, Fort Bend, Waller, Brazoria, Galveston, Liberty, Chambers, Orange, and Jefferson (Fig. 1 ). Topographic relief is low, ranging from sea level at the Gulf of Mexico to 300 ft or more in northwestern Harris County.

Jan 10, 1959

By L. S. Darken

It has long been recognized that the driving force in an isothermal diffusion process may be regarded as the negative gradient of the chemical potential (partial molal free energy) of the diffusing substance. This equivalence is exactly analogous to the fact that a mechanical force is the negative gradient of the potential energy. Nernstl in 1888 expressed these driving forces in terms of osmotic pressure, which is directly related to the chemical potential and free energy. The relationship between the diffusivity, mobility and activity coefficient of electrolytes in aqueous solution has been developed in great detail by Onsager and FuossZ with results which have recently been confirmed experimentally by Harned and Nuttall. This interpretation of the driving force in isothermal diffusion has been discussed earlier by the author,4,5 who pointed out4 that "for a system of more than two components it.is no longer necessarily true that a given element tends to diffuse toward a region of lower concentration even within a single phase region"; departure from the behavior of an ideal solution may be so great that the concentration gradient and the chemical potential gradient, or activity gradient, may be of different sign, thus giving rise to "up hill" diffusion. Since, in the discussion of that paper,4 some skepti- cism as to the validity of this conclusion was expressed, it seemed desirable to adduce experimental evidence to support it. Evidence of such diffusion in nonmetallic systems has already been reported. For example, it was observed by Hartley" in a series of experiments in which sodium chloride diffused from a saturated solution in about 5 pct acetone in water into a solution virtually free of salt but with essentially the same acetone concentration. Specifically, a wine glass containing crystalline salt and saturated solution was totally immersed in a large glass cylinder containing the acetone-water solution without salt; the difference in specific gravity was such as to prevent convection. Diffusion of salt was accompanied by a marked depletion of acetone in the wine glass; after 3 days the acetone concentration decreased from 5.37 pct at the top of the glass, where the salt concentration was low (0.93 pct), to 1.84 pct at the bottom where the salt concentration was high (21.64 pct). This "uphill" diffusion was found to be such that the partial pressure (activity) of acetone was essentially equalized throughout the system when a steady state gradient of sodium chloride was maintained. In order to demonstrate the existence of such diffusion in metals, a series of four weld-diffusion experiments was made. In these measurements pairs of steel of virtually the same carbon content, but differing markedly in alloy content, were welded at the end and held at 10500 C for about two weeks. Subsequent analysis showed that carbon had diffused so as to produce an inequality

Jan 1, 1949

By J. G. Faller, L. P. Skolnick

The distribution of cobalt and vanadium over non-equivalent crystallographic sites in C14-type Zr(VCo), alloys has been investigated. An anomalous X-ray scattering technique developed by Skolnick, Kondo. and lavine7 by which the separation in the scattering factors of two similar atoms can be enhanced was employed. Six alloys spanning the pseudobinary section ZrV1.6Co0.4-ZrVO.6CO1.4 at 10pct steps showed a nonrandom compositionally dependent distribution. Specifically, at high vanadium content cobalt preferentially occupied sites of type (6h) and vanadium, sites of type (2a; at low vanadium content the reverse was observed. In addition to the distribution fraction the structural parameters x and z were obtained. There was no significant deviation of these parameters from those obtained in the ideal C14 structure. Certain suggestions are made to account for the observed nonrandomness in the distribution of atoms on the two types of sites. INTERMETALLIC compounds of formula AB2 iso-morphous with MgCu2, MgZn2, and MgNi2 are known as Laves phases. Because Laves phases exhibit high symmetry and coordination numbers, the highest possible for an AB2-type compound,1 they are among the most frequently observed compounds in nature. In recent years interest has centered about the purely transition metal Laves phases2-' in efforts to understand the function of atomic size and electronic structure in the formation of these compounds. It has been observed that pseudobinary Laves phase systems often show a variation of structure across the phase diagram. Such a system is the ZrV2-ZrCo2 in which the structure varies from cubic MgCu2 to hexagonal MgZn2 to cubic MgCu2.4 Some understanding about the conditions under which the second modification is stable can perhaps be gained by studying the distribution of cobalt and vanadium atoms on lattice sites in the MgZn2 modification of the system ZrV2-ZrCo2. In both the MgZn2 and MgNi2-types there exist nonequivalent positions open to occupancy by the B element, whereas in the MgCu2 prototype all sites are equivalent. Skolnick, Kondo, and La-vine7 have developed an anomalous scattering technique suitable for this type of investigation. Whereas the influence of size on the formation of a Laves phase is well recognized, no substantial evidence has been put forth in support of the size ratio dependence of a particular prototype. Berry and Raynor8 suggested that RA /RB ratio does indeed affect the type of structure that is chosen, MgZn2 compounds tending to cluster about 1.225 while MgCu2 compounds were found at larger deviations from this ratio. Dwight,3 however, from a study of 164 Laves phases does not believe this generalization to be justified. Electronic effects are certain to play a part in the stability of Laves phases in general and in the choice of a structure type in particular. For example, size along would favor the formation of Laves compounds of Ti, Zr, Hf, Ta, or Nb as the A element with nickel or copper as the B element. The absence of such is attributed to an unfavorable electron : atom ratio by Elliott and rostoker.4 Early experiments of Laves and witte9 with pseudobinary and pseudoternary systems of the three prototypes established the dependence of crystal structure upon electron: atom ratios. They observed that the MgCu2 structure dissolved elements of higher valency until the electron: atom ratio of =1.8 was reached; the MgZn2 likewise dissolved elements of lower valency replacing zinc. witte,6 from calculations of the electron : atom volumes of Brillouin Zones, obtained limits of stability for the prototypes MgCu2 and MgZn2. Elliott and Rostoker4 used these limits with considerable success in the all-transition element Laves phases they investigated. According to witte,6 compounds between the electron :atom ratios of 1.80 and 2.32 were of the MgZn2 type; those above and below exhibited the MgCu2-type structure. On the basis of these limits and an assumed valency of zirconium based upon the near tetra-valence of titanium, Elliott and Rostoker obtained valencies for the first transition series elements. For the Laves phases with which this investigation is concerned, ZrV2 and ZrCo2, the authors calculated electron :atom ratios of 2.54 and 1.56, respectively. These ratios are for the MgCu2-type structure and straddle the stability band of the MgZn2 modification. One could, therefore, predict that a pseudobinary system ZrV2-ZrCo2 should pass through the MgZn2 modification in traversing the composition diagram from one end to the other. Implicit in this assumption is a smooth change of the electron: atom ratio from 2'54 to 1.56. MOSS10 states that his finding the low temperature structure of ZrCr2 to be C15 instead of C14 alters greatly Elliott's valency of zirconium and hence the assumed valencies of the other metals. Such a quantitative correlation of structure with electron : atom

Jan 1, 1963

By Morris Muskat

A theory is presented for calculating the performance history of complete water-drive systems producing from idealized stratified formations. The general equations are applied to systems where the permeability stratification is either of the exponential or linear type. Calculations were carried through for different degrees of permeability stratification, but with special emphasis on the effect of the mobility ratio between the produced oil and the invading water on the resultant performance. These results are also expressed graphically as curves for the initial water breakthrough recovery, for the different degrees of stratification, as a function of the mobility ratio, and of the composition of the produced fluid stream as a function of the cumulative oil recovery. For several typical cases the latter has also been plotted as a function of the cumulative oil and water throughflow. The general result is that when the mobility of the oil is lower than that of the invading water the channelling tendency resulting from the permeability stratification becomes aggravated as the higher permeability zones become flooded out. Situations of this type would obtain when producing low gravity or highly viscous oils. Conversely, if the mobility of the oil is high compared to that of the invading water, the flooding of the high permeability zones will lead to a retarding and choking effect, and the gross bypassing phenomena will be partially suppressed. These conditions would correspond to those of flooding high gravity or low viscosity oils. A discussion is given of the various basic assumptions made in the analysis, including that of ignoring the stripping phase of the production history as implied by relative permeability concepts. INTRODUCTION The physical ultimate recoveries from oil reservoirs are basically determined and limited by the physical oil displacement processes associated with the reservoir producing mechanism. In practice, however, the economic ultimate recoveries are further limited by the mobility of the reservoir fluids and the uniformity and continuity of the producing formation. In fact, it is the differential depletion between the component parts of the composite reservoir which ultimately determines the total recovery at the time of field abandonment. While this observation applies to both the solution gas drive and gravity drainage mechanisms, in which use is made only of the energy contained within the original oil-bearing reservoir, it is of even more paramount importance under operations wherein the energy associated with extraneous fluids provides the ultimate oil expulsion mechanism. Whether the invading fluid is the water from an edgewater drive, water injected for pressure maintenance, gas injected for pressure maintenance, or gas returned to the formation in a cycling program, it is often the continuity and uniformity of the producing section which will control the economic efficiency of the operations. The importance of the problem of reservoir non-uniformity does not, of course, lessen its complexity or the difficulties of its solution. In fact, these are inherently such that the concept of a "general" solution is virtually meaningless. About all that can be reasonably hoped for is the analysis of specific and well-defined types of non-uniformity which may give some degree of approximation to actual reservoir conditions. Since variations in the nature of the reservoir which depend only on the position along the streamlines will not lead directly to major differential depletion development within the reservoir, the types of non-uniformity considered thus far have involved stratification assumptions. That is, the producing section has been replaced by a multi-layer "sandwich," each uniform areally, and differing from the others only in its basic physical constants as to thickness, porosity and permeability. The fluid motion in the composite system is thus approximated by a parallel superposition of the independent fluid movements in the individual strata. For the specific application to cycling operations the theory of the effect of permeability stratification has been developed for both discontinuous' and continuous types of permeability stratification. Among the latter, treatments have been given of systems in which permeability distributions are governed by exponential, linear' or probability3 functions. In all these studies complete dynamical equivalence was assumed between the injected dry gas and the displaced wet gas. The overall effective permeability of each stratum was therefore considered as constant and independent of the degree of invasion of the injected fluid. In the case of the displacement of oil by water, the assumption of dynamical equivalence between the water and oil will be strictly valid only by accident. Even if the oil viscosity should be the same as that of the water, the effective permeability to the oil in the presence of the connate water will in general be quite different from that of the water behind the water-oil interface flowing past the trapped residual oil. As a result the effective permeability for the stratum as a whole and rate of water invasion will change with time as the intrusion continues. The differential fluid motion in the individual strata will thus also vary with time. Qualitatively, it is easy to predict the resultant effects. If the permeability to viscosity ratio of the invading fluid exceeds that of the fluid displaced, the stratification and bypassing effect of the perme-

Jan 1, 1950

By Brahm Prakash, R. Schuhmann

Chemical conditions for flotation and nonflotation of quartz with oleic acid as collector and barium, calcium, aluminum, iron, and tin as activators were studied using a simple vacuum-flotation technique in glass-stoppered graduates. The detailed study of barium activation led to an interpretation based on ideal Langmuirian chemi-sorption. FLOTATION of quartz is of practical importance as something to be avoided in soap-floating many types of ores. Clean, unactivated quartz is not floated with fatty acids and soaps, such as oleic acid and sodium oleate, in the quantities normally used for flotation. However, data in the literature indicate that almost any multivalent cation will activate quartz if given an opportunity. Thus, a common problem is to prevent activation of quartz by the various inorganic cations inevitably present in flotation pulps. Wark and his coworkers1 have demonstrated the reversibility of the chemical reactions and adsorptions involved in the activation, depression, and collection of the common sulphide minerals. The procedure in much of their work was to bring a mineral surface to equilibrium with solutions of known pH, collector concentration, and activator concentration, and then to test the floatability of the mineral by contact-angle measurement. From the data, graphs were constructed with pH and reagent concentrations as coordinates. These graphs show fields of flotation and fields of nonflotation, separated by narrow transition regions whose locations are shown by so-called contact curves. From the shapes and locations of the contact curves, which roughly separate fields of flotation from fields of nonflotation, a quantitative understanding of the interaction of the reagents with each other and with the minerals often can be deduced. The study of quartz flotation to be described in this paper follows in broad lines the approach of Wark and coworkers. That is, pH, activator concentration, and collector concentration were varied to find equilibrium conditions of flotation and non- flotation, and the results are presented graphically by means of contact curves. However, instead of testing for floatability by measuring the contact angle on a polished surface, a simple vacuum flotation technique was developed and used. Purified oleic acid was the collector and terpineol the frother. Barium activation was studied in some detail, and exploratory studies were made of activation with calcium, aluminum, ferric iron, and stannic tin. Preparation of Materials Quartz: Large lumps of high-grade vein quartz were crushed dry in a cone crusher and rolls. The —20, +28-mesh portion was screened out and used in the subsequent steps. This material was passed through a high-intensity magnetic separator to discard iron, then leached twice with hot concentrated HCl and washed repeatedly with distilled water. The cleaned sand was then wet ground with porcelain balls in a porcelain pebble mill, deslimed repeatedly by settling and decantation to discard —800-mesh material, and again washed with hot HCl followed by distilled water. The resulting stock of quartz was stored under water. Chemical analysis gave 99.8 pct SiO2. Table I gives the size analysis of the quartz used for flotation tests. Calculations from these data, using shape factors given by Gaudin and Hukki9 indicate a specific surface of about 500 cm2 per g. Blank flotation tests in distilled water, and in water with added frother, showed the prepared quartz to be completely nonfloatable and thus indicated the absence of organic contamination. Oleic Acid: The preparation of oleic acid was based on fractional vacuum distillation of methyl oleate2,3 followed by regeneration of oleic acid, and finally fractional crystallization of oleic acid from acetone solutions at low temperatures." The pure oleic acid was stored in a refrigerator. The iodine number of the oleic acid was found to be 90.0 (theoretical 89.93). Oleic acid was used in the form of a dilute water solution of sodium oleate, after preliminary flotation tests showed no effects of form of addition and order of addition of reagents when an adequate conditioning time (that is, 30 min) was provided. Other Reagents: Sodium hydroxide solutions low in carbonate were prepared by first making 1:1

Jan 1, 1951

By Brahm Prakash, R. Schuhmann

Chemical conditions for flotation and nonflotation of quartz with oleic acid as collector and barium, calcium, aluminum, iron, and tin as activators were studied using a simple vacuum-flotation technique in glass-stoppered graduates. The detailed study of barium activation led to an interpretation based on ideal Langmuirian chemi-sorption. FLOTATION of quartz is of practical importance as something to be avoided in soap-floating many types of ores. Clean, unactivated quartz is not floated with fatty acids and soaps, such as oleic acid and sodium oleate, in the quantities normally used for flotation. However, data in the literature indicate that almost any multivalent cation will activate quartz if given an opportunity. Thus, a common problem is to prevent activation of quartz by the various inorganic cations inevitably present in flotation pulps. Wark and his coworkers1 have demonstrated the reversibility of the chemical reactions and adsorptions involved in the activation, depression, and collection of the common sulphide minerals. The procedure in much of their work was to bring a mineral surface to equilibrium with solutions of known pH, collector concentration, and activator concentration, and then to test the floatability of the mineral by contact-angle measurement. From the data, graphs were constructed with pH and reagent concentrations as coordinates. These graphs show fields of flotation and fields of nonflotation, separated by narrow transition regions whose locations are shown by so-called contact curves. From the shapes and locations of the contact curves, which roughly separate fields of flotation from fields of nonflotation, a quantitative understanding of the interaction of the reagents with each other and with the minerals often can be deduced. The study of quartz flotation to be described in this paper follows in broad lines the approach of Wark and coworkers. That is, pH, activator concentration, and collector concentration were varied to find equilibrium conditions of flotation and non- flotation, and the results are presented graphically by means of contact curves. However, instead of testing for floatability by measuring the contact angle on a polished surface, a simple vacuum flotation technique was developed and used. Purified oleic acid was the collector and terpineol the frother. Barium activation was studied in some detail, and exploratory studies were made of activation with calcium, aluminum, ferric iron, and stannic tin. Preparation of Materials Quartz: Large lumps of high-grade vein quartz were crushed dry in a cone crusher and rolls. The —20, +28-mesh portion was screened out and used in the subsequent steps. This material was passed through a high-intensity magnetic separator to discard iron, then leached twice with hot concentrated HCl and washed repeatedly with distilled water. The cleaned sand was then wet ground with porcelain balls in a porcelain pebble mill, deslimed repeatedly by settling and decantation to discard —800-mesh material, and again washed with hot HCl followed by distilled water. The resulting stock of quartz was stored under water. Chemical analysis gave 99.8 pct SiO2. Table I gives the size analysis of the quartz used for flotation tests. Calculations from these data, using shape factors given by Gaudin and Hukki9 indicate a specific surface of about 500 cm2 per g. Blank flotation tests in distilled water, and in water with added frother, showed the prepared quartz to be completely nonfloatable and thus indicated the absence of organic contamination. Oleic Acid: The preparation of oleic acid was based on fractional vacuum distillation of methyl oleate2,3 followed by regeneration of oleic acid, and finally fractional crystallization of oleic acid from acetone solutions at low temperatures." The pure oleic acid was stored in a refrigerator. The iodine number of the oleic acid was found to be 90.0 (theoretical 89.93). Oleic acid was used in the form of a dilute water solution of sodium oleate, after preliminary flotation tests showed no effects of form of addition and order of addition of reagents when an adequate conditioning time (that is, 30 min) was provided. Other Reagents: Sodium hydroxide solutions low in carbonate were prepared by first making 1:1

Jan 1, 1951

By S. W. Churchill

A digital computer was used to obtain an exact numerical solution for the transient behavior of the insulation and earth adjacent to an isothermal, submerged flat surface for a single set of parametric values. Comparison of the computed results with analytical solutions for limiting conditions revealed that a complete and general solution for all parametric values could be constructed from these limiting solutions. Complete and general solutions for insulated spheres and cylinders can also be constructed from limiting solutions. The procedure is illustrated in part for an insulated spherical tank. The heat flow was found to fall rapidly to a pseudo-steady state; after some time, it then decreased slowly to zero for a flat plate and cylinder, and to a low steady-state rate for a sphere. The accumulative heat flow during the initial falling-rate period may be a significant fraction of the heat flow during the entire first year. INTRODUCTION Underground storage of liquefied natural gas and liquefied petroleum gas has received considerable recent attention. The rate of heat flow from the earth to the storage cavity and the resulting temperature field in the earth are important factors in the technical and economic evaluation of such storage facilities. The objective of this paper is to indicate how complete solutions can be developed for the transient flux and temperature field for various geometrical configurations. The representative properties and dimensions, and the resulting parameters utilized in the illustrative results, are indicated in Table 1. The results presented herein are for dry earth. The effect of the latent heat of solidification of water in the soiI will be described in a subsequent paper. The temperature field in the insulation and earth is determined by energy balances, boundary conditions and initial conditions. The physical problem can be described mathematically with virtually no idealizations insofar as physical properties are known. It appears possible to solve the equations by numerical integration on a high-speed digital computer for any geometrical configuration and conditions. For complex situations, however, the computations are expensive and the results are highly specific. Analytical solutions have been developed for a few simple but important geometrical configurations and conditions, including one-dimensional heat transfer from (or to) earth at an initially uniform temperature to (or from) flat plates, spheres and circular cylinders. The solution for an insulated flat plate has also been derived but is in the form of a slowly converging series involving taklated functions. It was planned to use a computing machine to evaluate this series for a number of representative conditions. However, upon examination of the results of preliminary computations it was discovered that a complete, general and accurate solution could be developed by interpolation between several much simpler solutions for limiting conditions. This technique then was used to develop a a complete solution for an insulated sphere. The equations presented herein are derived in Carslaw and Jaegerl and other books on applied mathematics, or they are simple extensions of these previous results. Hence, all derivations are omitted. Emphasis is instead placed on the results and their