Search Documents

By A. C. Stevenson

THE hoisting equipment selected for use at the Hecla mine in 1907 was one of the first Ilgner type Ward-Leonard controlled hoists put into ser- vice. Development of the Hecla below the 2000-ft. level, which is the ultimate depth of the old equipment, and the additional service required by the "Star," necessitated additional hoisting capacity. It will be our endeavor to give j70u a brief summary of the engineering which determined the selection of the new equipment. The engineer in making his decision as to size, speed, etc., must .make a duty study approximating as closely &s possible the actual opePating conditions which will eventually result. Primarily, it is necessary to know what work must be done. The amount of material which must be hoisted under the maximum condition, as well as its distribution from the various levels of the mine was carefully studied by the mining engineering department, as the ore reserves, their relative-position and the probable method of working them had to be considered. The final figures reached, taking into consideration four different operating conditions, were, in tons per day, as follows:

Jan 1, 1930

By Granville S. Borden

TAXES in general are onerous and are not a pleas- ant subject for discourse. There are, however, some very cogent reasons why we should dedicate a part of our thoughts and services to the solution of our taxation problems: The happiness of every individual in the United States depends on the maintenance of a stable government. The collection of taxes is an indispensable function of a stable government, and the power to collect taxes depends in part on the question whether taxes are fair and equitable. As representatives of the mining industry it is our duty to .protect it from tax burdens that are discriminatory or- unjust, As citizens we find that taxation is incessantly endeavouring to unravel the knots in our purse-strings. Many members of the Institute may be consulted .by clients on taxation problems, especially the phases which re- late to depletion and valuation. It is the hope of the writer that the discussion of some of the current depletion and valuation problems will stimulate interest in the problems, and as a consequence of this new interest facts will be ventilated, opinions rendered and services performed, which will be a means of solving the problems to the satisfaction of all parties interested.

Jan 1, 1929

By AIME AIME

A COMMITTEE, representing a number of the larger oil companies in California, met in San Francisco on Nov. 12, 1919, for the purpose of discussing the depletion deductions to be allowed in computing the Federal income taxes under the Internal Revenue Act of 1918. It was the general sense of the committee that a system based on sound engineering principles should be formulated for estimating the value of oil properties, so that better and more uniform results might be obtained by the various companies in filing their returns with the' Bureau of Internal Revenue. A subcommittee, consisting of five petroleum engineers, was thereupon appointed to collect information and draft a report. The committee consisted of G. C. Gester, Carl H. Beal, J. H. G. Wolf, E. G. Gaylord, and E. D. Nolan; this report is the work of that subcommittee.

Jan 1, 1920

By Richard H. Olson, Lawrence A. Roe

Talc, when it can be isolated as a pure mineral, has a composition of 63.36% Si02, 31.89% MgO, and 4.75% H20. However, as an industrial commodity, talc rarely approaches theoretical purity. Nevertheless, the impure talc products of commerce find a multitude of uses and few substitutes in many industrial applications. Minerals commonly associated with and sold together in talcose mixtures are tremolite, chlorite, dolomite, mica, and magnesite. Steatite was originally a mineralogical name applied to pure talc. In today's commerce steatite generally refers to the massive variety of talc suitable for electrical insulator manufacture. Block steatite ore can be machined into various shapes. Impure varieties of massive or block talc are still commonly referred to as soapstone. Soft massive talc, suitable for crayon manufacture, has been referred to as French chalk. It has been common practice to discuss talc, soapstone, and pyrophyllite under the same general heading. In the case of soapstone and talc, this is a natural thing since many different types of platy, soft minerals exhibiting a high degree of lubricity (commonly referred to as "slip") have been grouped together and called soapstone or talc. Further, there is a mineralogical relationship of pyrophyllite to montmorillonite and talc to hectorite. Thus, talc and pyrophyllite are sometimes referred to as clay minerals. When finely divided talc or pyrophyllite are combined with water in proportions to make a slurry, the end product does have the appearance of a clay-water mixture. Soapstone for utensils and ornaments was mined by prehistoric Indians on Santa Catalina Island, CA. In the mid-1800s, soapstone from deposits along the western foothills of the Sierra Nevada mountains was used by white settlers for building and ornamental stone and in the linings and foundations of furnaces (Anon. 1956). Previous to 1916, the annual recorded production in California did not exceed 2000 tons. In the period 1912 to 1918, output rose sharply when the Talc City, Western, and Silver Lake mines were put into operation. From 1916 to 1935, the state's annual output of talc was in the range of 9,000 to 20,000 tons. In the mid-1930s, the use of talc in wall tile grew rapidly. In California, the production of talc grew to 63,000 stpy in 1943. The post-war building boom helped California production to grow to 120,000 tons in 1951. In 1968, production was 165,000 tons, 1969 was 145,000 tons, 1970 was 185,000 tons, and 1972 was 155,000 tons. Talc mining in New York state dates to about 1878 when a Colonel Palmer and associates opened the first commercial talc mine on the Nelson Freeman farm near Talcville, New York. In 1893, this operation was sold to the International Pulp Co. which changed its name to the International Talc Co. in 1944. The company was acquired by R.T. Vanderbilt Co., Inc. in 1974. A new talc operation, the Gouverneur Talc Co. owned by the R.T. Vanderbilt Co., began operations near Balmat, NY, in 1948. The initial capacity of the processing plant was 200 stpd (Gillingham, 1950). Subsequent expansions have increased capacity to over 600 stpd. The Vermont talc industry began just after the start of the 20th century. Talc was discovered in the area of the present Johnson, VT, talc mine in 1902. The American Minerals Co. initiated plant operations at Johnson in 1904. The Magnesia Talc Co. opened a plant at Waterbury, VT, in 1913 (Burmeister, 1963; Trauffer, 1964) and acquired American Minerals Co. in 1923. The Waterbury mine has a long history as a producer of talc crayons. The Eastern Magnesia Talc Co. was formed in 1924 by a merger of the Eastern Talc, Co.

Jan 1, 1983

By R. P. Alger, C. A. Doh, M. P. Tixier

The principle, the equipment and field operation of sonic logging are described. The tfio-receiver system produces logs independent of hole size and mud. Field experience is given and forms the basis for the interpretation of the log. The derivaiion of porosity values from measured velocities is discussed, according to the type of formations (hard formations, compacted sands, and unconsolidated formations). The time-average equation proposed by M. R. J. Wyllie, A. R. Gregory and L. W. Gardner is uscd as a basis for the computation of porosity in limestones, cemented sandstones and compacted sands. Variations in the lithologic character of limestones do not seem to change the porosity calibration markedly. The compaction of sands is related to the compaction of shale adjacent to them; and, thus, the shale velocity is used for the establishtnent of empirical relations for the computation of porosity in un-consolidated formations. Various forrnu.las are tentatively presented to account for shale and fluid content. Field experience demonstrates that considerable attenuation of sonic energy takes place in unconsolidated formations, particularly when gas bearing, and in fractrured formations. Unusually large attenuation produces skipped cycles, a feature easily recognized. The application of sonic logging to structura1 studies is featured, showing its possible integration with the dipmeter. Correlation and its application to the irrterpretation of seisrmic surveys are reviewed. The paper is illustrated with field examples. Sonic logging is the recording of the time required for a sound wave to traverse a definite length of formation. Sonic travel-times are inversely proportional to the speed of sound in the various formations. The speed of sound in subsurface formations depends upon the elastic properties of the rock matrix, the porosity of the formations and their fluid content and pressure. Below the "weathered" or low-velocity layer extending from 50 to 100 ft or so below the surface, sound velocities may range from about 6,000 ft/sec in shallow shales to as much as 34,000 ft/sec in dolomites. In hard formations (well cemented and/or compacted), the sonic log reflects the amount of fluid in the formations; hence, it correlates well with their porosity. In unconsolidated formations, which are usually of fairly high porosity, the sonic log gives an approach to porosity determination, when its readings are corrected for lack of compaction, shaliness, and fluid content. In such formations the sonic log may also indicate the presence of gas and may distinguish between oil- and water-bearing beds. Field experience with the sonic log in soft formations is more limited than that in hard formations. As a result, interpretation in soft formations is not as well developed at the present time. In all types of formations the sonic log is of considerable value for correlation and for the more detailed and accurate interpretation of seismic data. EQUIPMENT AND PRINCIPLE FIELD OPERATION The basic features of the sonic log" are a two-receiver system and two available short spacings. The log is recorded on film using a standard field unit. An SP curve is run along with the sonic curve to give a more interpretable log and to establish absclute depth control for comparison with other logs run in the same well. If desired, a gamma-ray curve may be recorded simultaneously with the sonic curve. Fig. 1 shows the sonic-gamma ray combination sonde. The sonic portion consists of a sound transmitter and three receivers mounted 3, 4 and 6 ft from the transmitter. When the first and second receivers are used, a sonic log of I -ft spacing (distance between receivers) is obtainea. Use of the first and third receivers provides a 3-ft spacing. The transmitter emits pulses at the rate of l0/sec. After the transmission of a pulse the first arrival of sound energy at each receiver triggers its response system. The differ-

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 J. Szekely, V. Stanek

Steady state heat transfer and fluid flow rates are predicted for the counter-current contacting of gas and solid streams. For a .fixed pressure drop across the bed the convective heat transfer rate is affected by the gas flow rate, which in turn is affected by the temperature distribution in the system. Variations in bed voidage and in particle size were found to change the predicted temperature profiles quite markedly, in particular when the thermal capacities of the gas and solid streams were equal. IN Part I we examined the effect of step and pulse changes in the inlet gas temperature on transients in the flow of gases through fixed packed beds. Here we shall discuss heat flow and fluid flow interactions in the counter-current heat exchange between packed solids and a gas stream. Counter-current heat exchange in packed beds has been extensively studied. In his now classical paper, Furnas&apos; assumed constant property values and developed a simple analytical. expression for the temperature profiles in the gas and in the solid stream. A highly sophisticated analysis of counterflow systems, allowing for transients and interparticle conduction was developed by Amundsen. In preparing a very readable review of this topic Elliott&apos; also presented a number of elegant solutions, including the case involving heat generation. A great deal of work has been done on the application of these concepts to the modelling of the blast furnace. However, in all this previous work the mass flow of the gas was regarded as a quantity independent of the gas-solid heat exchange process. In the following we shall consider counter-current heat exchange between a gas stream and a packed bed of solids, under conditions where the pressure drop across the bed is specified, rather than the mass flow of the gas. FORMULATION Consider a packed bed, extending from y&apos; = 0 to y&apos; = L, in which moving solids and a moving gas are contacted in a counterflow arrangement. Let cold solids, with an initial temperature Tsi and mass flow GA, enter at y&apos; = L, and hot gas with an initial temperature of Tgi and mass flow Gg enter at y&apos; = 0. In contrast to the problems considered in Part I this system has a nontrivial steady state and we shall restrict our attention to the study of the pressure drop flow relationships and the resultant temperature distributions for steady state conditions. Under these conditions the formulation of the problem requires the statement of the following: 1) heat balance on the gas 2) heat balance on the solids 3) the mechanical energy balance equation and in addition the knowledge of the equation of state, and a functional relationship between viscosity and temperature and an expression for the heat transfer coefficient in terms of the other parameters. On neglecting temperature gradients within the solids and for a constant mass flow of the gas, the solution to this problem is well known.&apos;" The expressions for the temperature profiles of the gas and solid streams may be put in the following form: where Tg and T, denote the dimensionless gas and solid temperatures respectively, y is a dimensionless distance, and R is the ratio of the thermal masses. These quantities are defined as follows: For the particular case of equal thermal mass capacities by the gas and the solid, that is R = 1, the Eqs. [1]and[2] simplify to Eqs. (71 and[8], respectively: i.e., the resultant profiles are linear. As for most practical cases , solutions given by Eqs. [7] and [8] will contain only the inlet temperatures as parameters. For a fixed, predetermined value of C;b the above

Jan 1, 1970

By CHARLES E. W.

I. LOCATION. THE Mojave mining district is situated in a group of small hills centering around Soledad peak, in the Mojave desert, Kern county, Cal. These hills are about 4.5 miles SSW. of Mojave, a railroad town at the junction of the Santa Fé and Southern Pacific railroads. The town of Mojave is about 50 miles north of Los Angeles, to which it is connected both by telegraph and by rail. The distance on the Southern Pacific railroad is 102 miles, covered in 3 hr. and 45 min. by express train. A sketch-map of this region is shown in Fig. 1. The Mojave district is relatively isolated. The only other mining towns in this section of California are Randsburg and Johannesburg, in a district (containing the well-known Yellow Aster mine) at the northern end of a low ridge of hills, which begins about 20 miles NE. of Mojave, runs nearly N. for 35 miles, and varies in width between 3 and 6 miles. As will be seen later, the origin and formation of these hills are not the same as those of Soledad peak. II. CLIMATE AND VEGETATION. The climate and vegetation of the Mojave district are typical of the arid regions of the southwestern United States. During the winter and spring there is some rainfall, and the tempera-ture is generally from 60° to 70° F. in the day-time. The rainfall in the spring is often very considerable for short spaces of time, and cloud-bursts are not infrequent. In summer there is continual sunshine, with temperatures, as a rule, exceeding 100° F. in the day-time, while 115° F. in the shade of a house is not at all unusual; but at night the temperature falls considerably. No rain falls during the summer. The vegetation

Jan 1, 1906

By W. S. Hyler, W. L. Bruckart

The powder metallurgy and arc-cast types of wrought molybdenum stock were studied in rotating beam fatigue. Endurance ratios of unnotched specimens after 5x10 cycles were found to be 0.74 and 0.81, respectively. Fatigue notch factors (K,) of 1.36 and 2.05, respectively, were obtained on V-notched specimens having a theoretical stress-concentration factor Kt = 3.1. MOLYBDENUM and its alloys are rapidly gaining importance as engineering materials of construction for high temperature applications. AS a wrought material, molybdenum and its alloys are outstanding above 1800°F. It is necessary, however, with molybdenum to provide elevated temperature oxidation protection. This is currently being studied by a number of investigators. Thus far, there have been no published data on the fatigue behavior of molybdenum either at room temperature or at elevated temperature. The knowledge of fatigue behavior, particularly notched fatigue behavior, becomes increasingly important as the engineering importance of a material increases, especially from the design standpoint. The work described here was directed toward determination of the notched and unnotched rotating beam fatigue characteristics of molybdenum at room temperature. Also, it was desired to learn whether any difference exists between wrought molybdenum of the powder metallurgy and the are-cast types. Fabrication History of Bar Stock For this work, bar stock was procured from Climax Molybdenum Co. (are-cast type) and from Westing-house Electric Corp. (powder metallurgy type). The processing histories of these materials are listed in Table I. The first lot of are-cast molybdenum was deoxidized by aluminum additions. The second lot was deoxidized by carbon." Both lots of the powder Microstructure A microscopic examination of these materials, as received, revealed the powder metallurgy molybdenum to be free of inclusions and to be finely fibered. On the other hand, the are-cast C molybdenum possessed numerous carbide inclusions, and the arc-cast A appeared to have coarse grained fibers. These are shown in Fig. 1. Owing to the large apparent difference in fineness of fiber between the are-cast A and the powder product and are-cast C, a check was made of the recrystallized grain sizes to determine whether these differences persisted. A series of recrystallization anneals revealed mixed grain sizes in both the are-cast A and the powder product. Variations in grain diameter of as much as 50 times were found after full recrystallization at 1200°C for 30 min. There was no distinction between the two. Arc-cast C was not checked. Room Temperature Mechanical Properties Prior to subjecting the materials to fatigue tests, room temperature mechanical properties were determined. These data, together with the tensil: strength of the notched bars, data,togetherare listed in Table tensil: There were no significant differences in the tensile strengths of the two types of molybdenum. However, the powder metallurgy material was generally more ductile and yielded at a somewhat lower stress than the arc-cast material. Tensile strengths of the notched specimens indicated a slightly greater notch strengthening effect in the are-cast C than in the powder product, lot 2. Room Temperature Fatigue Properties Sample Preparation: Specimens were prepared for use in R. R. Moore rotating beam fatigue machines. These specimens had a continuous radius test section with a minimum diameter of 0.260 in. The specimens were prepared by grinding to shape and polishing to produce a smooth unworked surface. Polishing consisted of rotating a disk in the continuous radius section on which were attached, successively, emery belts of 240, 400, and 600 grit. The specimens were slowly rotated during the operation. The belt speed was high, and rotation was such that polishing scratches were longitudinal. Unnotched Fatigue Properties: Unnotched fatigue data were obtained from a number of specimens of each of the four lots of molybdenum. The results of the tests are presented in Fig. 2. From this figure,

Jan 1, 1956

By G. Derge, Ling Yang, John Henderson

Self-diffusion coefficients of aluminum have been measured by the capillary reservoir technique in liquids of the CaO-SiO2-Al2O3 system containirg equimolar portions of CaO and SiO2, in the temperature range 1400o to 1520oC. For the melt containing 6.0 mole pct Al2O3 the results can be expressed by the equation D = [(4.3 * 0.3) x 104] exp [(-85,000 * 17,500)/RT]cm~sec-&apos;; for the melt containing 12.5 mole pct Al2O3, by the equation D = (5.4 i 0.2) exp [(-60,000 * 10,00O)/RT] cm2sec-: A concept of the constitution of these melts has been developed which proposes that the ratios O/(Si +Al) and Al2 O3/CaO determine the dominant species present. In recent years considerable work has been done on the system CaO-SiO2-A12O3 with a view to elucidating the structure of liquids of this system. Conductivity,&apos; electromotive force, 3 density,4 and activity5 measurements have all contributed to an understanding of the problem, but as yet there is no entirely satisfactory theory that can account for all the observed experimental results. In principle, measurements of self-diffusion coefficients should give further information which is not obtainable in other ways. Self-diffusion coefficients for calcium, silicon,6 and oxygen7 have been measured. The present investigation extends the series of self-diffusion coefficient measurements to include aluminum and develops a concept of constitution for liquids of this system which is compatible with all the observed experimental results. EXPERIMENTAL The technique used in the measurement of the self-diffusion coefficient of aluminum was the capillary-reservoir method,8 A126 being incorporated in the melt contained in the capillary and inactive aluminum in the reservoir melt. This method was chosen because it minimizes convection currents&apos; in the diffusion zone. It has the further advantage that the boundary conditions are such that a slight modifica- tion of Anderson and saddington&apos;s8 solution of Fick&apos;s Law yields an equation which requires measurements only of the length and the initial and final average specific activities of the solidified melt within the capillary. where t is the diffusion time, Co is the average specific activity of the melt contained in the capillary at t = O, C is the average specific activity of the melt contained in the capillary at t = t, 1is the length of the capillary, and D is the self-diffusion coefficient. The equation is of somewhat different form than those previously applied to high temperature diffusion studies because the specific activity of the melt in the reservoir was negligible both before and after diffusion. The diffusion apparatus is shown schematically in Fig. 1. It consisted of a silicon carbide resistance furnace containing the inactive melt in a graphite crucible, internal dimensions 7/8 by 4 in., enclosed in a 1 1/4- by 1 1/2- by 30-in. McDanel tube. A 3/8-in. McDanel thermocouple sheath to which graphite radiation shields and a graphite protection tube were affixed, was used to hold the capillaries. By moving this sheath through a rubber seal at the top of the McDanel tube the capillaries could be lowered into or raised out of the melt. Temperature was measured with a Pt-Pt 13 pct Rh thermocouple which was calibrated periodically against a standard couple. The furnace temperature was controlled to ± 5°C through the thermocouple immersed in the melt. Temperature variation over the depth of the melt was less than ± 3°C. The inactive melts were prepared from pure native quartz, assaying better than 99.9 pct SiO2, A.R. grade A12O3 and A.R. grade CaCO,. The same quartz and CaCO3 were used in preparing the active melts, the

Jan 1, 1962

By William Hurst, A. F. van Everdingen

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst&apos; shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.&apos; The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy&apos;s equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

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

By A. F. van Everdingen, William Hurst

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst&apos; shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.&apos; The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy&apos;s equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By M. S. Hunter, D. L. Robinson

An extremely fine subgrain structure found in aluminum and aluminum alloys is shown and a method for revealing this structure is described. The appearance and some of the characteristics of this structure are described and the possible significance of subgrain structure in terms of chemical, electrochemical, and metallurgical processes is considered. GRAINS and grain boundaries, which constitute the principal manifestations of the crystalline aggregates known as metals, play such an important part in determining the characteristics and properties of metals and alloys that their nature and behavior have long been of much concern to the met-allographer. Many investigators have speculated on the possible existence of a finer structure within the usual grain structure and, through X-ray diffraction and metallographic techniques, have revealed evidence that such finer structures exist. Among these structures are the mosaic blocks and domains described by Burgers1 and others in their contributions on crystalline structure and the dislocation theory, the small units associated with polygoniza-tion described by Cahn,2 and the intragranular blocks revealed more recently by Lacombe and Beaujard.3 his paper describes some novel methods for clearly revealing a very fine granular substructure within the usual grains of pure aluminum and aluminum alloys. It also compares the characteristics and behavior of this subgrain structure with those of the ordinary grain structure and gives some examples of the manner in which this substructure may influence the characteristics of aluminum alloys. Metallographic investigations of aluminum and its alloys have dealt extensively with the functions and behavior of grains and grain boundaries and have evolved a definite concept of the relation between grain structure, orientation and physical, mechanical, chemical, and electrochemical behavior. During the examination of etched samples, a roughening or mottling of the body of certain grains has occasionally been observed. Such roughening obviously represented differential attack by the etching solution, but its true significance was not at first realized. It is now known that this roughening probably rep- resented a subgrain structure, although the particular etching treatments used were not capable of revealing this subgrain structure clearly. Also, in practically all cases, the subgrains were apparently far too small to be resolved by the light microscope. Recently, a method has been developed for revealing clearly the subgrain structure of most aluminum alloys. This method is somewhat unique in that it eliminates much of the usual tedious mechanical polishing and the attendant danger of a surface layer of flowed metal and at the same time is simple, rapid, and effective. It involves the use of the Alcoa R5 Bright Dip, a patented chemical polishing treatment licensed by Aluminum Co. of America. This treatment produces a bright, highly polished surface at least the equal of any mechanical polish and also etches the subgrain boundaries as well as the ordinary grain boundaries to reveal both the grain and subgrain structure simultaneously. Metallographic specimens to be examined are taken through the customary grinding and polishing operations used for aluminum alloys.&apos; The usual care required in the final polishing operation to prevent a flowed surface layer and to obtain a high polish is not necessary, however, because any flowed layer is eliminated and a high polish is produced by the final immersion in the Alcoa R5 Bright Dip. Depending on alloy type, this dip takes from a few seconds to about 3 min, after which the sample is rinsed well, blown dry, and is ready for examination. The substructure can be revealed by some of the usual metallographic methods such as electropolish-ing or mechanical polishing and etching but no method has been found which is the equal of the new chemical polishing treatment. Mechanical polishing followed by the usual metallographic etching is not entirely satisfactory for revealing this substructure, probably because of the flowed surface layer resulting from the mechanical polish. If the usual etching methods are prolonged sufficiently to remove this layer, the specimen surface becomes quite rough and the substructure is not defined clearly.

Jan 1, 1954

By H. Matlock

Beam-columns with continuous or discontinuous transverse and angular loads and elastic restraints are represented mathematically in a manner corresponding to finite-element mechanical models. Solutions for linearly elastic cases are accomplished by direct recursion-equation elimination of unknowns, and for cases involving nonlinear loads and supports, by repeated trial-and-adjustment of complete elastic-type solutions. Five practical applications are presented by which various capabilities of the general method are illustrated. The examples include simulation of a pipeline launching and the bending and collapse of a conductor pipe in deep-water drilling operations. INTRODUCTION The diversified and complex operations being conducted in support of the production of offshore oil have produced a wide variety of new problems in structural mechanics. Many of these problems cannot be solved by conventional analytical or experimental approaches and have therefore stimulated inquiry into new methods of solution. It is the purpose of this paper to give a brief description of an analytical method for solving nonlinear beam-column problems and to illustrate the generality of the method by applying it to several problems related to offshore operations. Although the data for all of the examples are artificial, the cases discussed will be recognized as being based on actual situations. DESCRIPTION OF METHOD The analytical approach which is to be described stems from research on offshore foundations, in particular, the study of lateral-load behavior of foundation piles. The method represents a generalization of a technique which was suggested by Gleser&apos; and extended by Focht and McClelland2 and Reese and Matlock." To maintain generality and therefore versatility in the method, a considerable variety of input data is considered. Fig. 1 illustrates the various types of loads and restraints which might be applied to a short segment of a beam. Each of the terms illustrated in the figure is shown acting in a positive sense. Capital letters are used to designate concentrated quantities whereas the corresponding lower-case symbols indicate effects which are distributed in some way along the length of the beam. The transverse loads and the elastic restraints may be divided into two classes: (1) loads Q and q and springs S and s, which act normal to the axis of the beam, and (2) couples T and t and rotational springs R and r, which act in an angular sense. In addition, there is an axial force P which may be constant or may vary along the length of the beam. The bending stiffness of the beam (the product EI of the modulus of elasticity and the moment of inertia) is designated by the symbol F. The assumptions of conventional elementary beam theory provide the basis for the present method. Among these are the neglect of shearing and axial deformations and the limiting of consideration to straight beams of symmetrical cross-section. Lateral deflections are considered to be small compared to original dimensions. A deformed element of such a beam is shown in Fig. 2a. The curvature of the beam is approximated by the second derivative of the deflection y with respect to distance x along the beam. It is assumed that the material of the beam is linearly elastic and therefore the relation of Eq. 1 between curvature and bending moment M is the ordinary expression from elementary beam theory. A generalized beam element is shown in Fig. 2b. The element has been deflected in a positive direction by an amount y and tilted through a small positive angle dy/dx. The various loads and elastic restraint reactions are detailed in the figure. Temporarily, only distributed effects are considered; provision for concentrated loads and restraints will be made subsequently. Eq. 2 in Fig. 2b is obtained by summing the forces and moments acting on

By B. Trafton, H. L. Hartman

An earlier paper1 introduced a rapid solution for mine ventilation network problems, employing the digital computer. In this paper, refinements and additions to the computer program which expand its scope of application and a solution to a ventilation network having a coal-mine configuration are presented. The revised program-written in FORTRAN for the IBM 7074 computer in Penn State&apos;s Computation Center-is applicable to any mine ventilation network employing natural splitting and having its fan located externally. It accomplishes the following: 1) Distributes an assumed total quantity of air flow among all airways of the network, performing iterations until an acceptable level of error is attained and the flow in each airway is known. 2) Computes the pressure drop (head loss) in all airways and the overall network and mine head losses. 3) Determines the operating point for the mine and a given fan, adjusting all air flows and head losses to the actual values. The basis of analysis used is the Hardy Cross method of successive approximation for solution of any fluid network, explained in detail in the earlier article. In brief, corrections to assumed air flows are computed for all the loops of a network; following distribution to the airways of each loop, the process is repeated until acceptable accuracy is attained. When a hand or desk-calculator solution is used, there is advantage to close approximation of the actual flows in the original assumption of air-flow distribution in the network. Since publication of the previous paper, in which assumption of initial quantities was also necessary in the computer solution, the program has been modified to make the first distribution itself, assigning the total quantity to one specified path in a split and then proceeding to correct the flows. The number of iterations required for a high-accuracy solution is little affected, at least in networks of few loops, 10 to 15 being sufficient to reduce the error to 1 pct or less. The principal limitation of the original program was that it contained no provision for handling a fan. Further, it did not determine the overall head loss in the network and was not designed for multiple parallel openings, such as occur in coal mines. The revised program was written to overcome these disadvantages and also to print out the adjusted flows for all airways, corrected for the system operating point. At the present time, the program can handle only an external fan. The procedure employed in solving a ventilation network is now as follows: 1) Input data consisting of airway characteristics and loop designation along with assumed direction of flow, for each branch of the network and any trunk airways, are fed into the computer. 2) Employing the resistance relationship, where R is resistance, K is friction factor, P is perimeter, L is length, A is area, and n is number of airways (where multiple intakes and returns are employed); values are computed for all airways. 3) An arbitrary, continuous route of flow through the circuit—from beginning to end—is selected and the path identified by airway numbers and read in. 4) Based on a total assumed quantity of flow Q (usually 100,000 cfm.) the computer employs the procedure described in the previous article in distributing quantities of flow among all the airways until a satisfactory level of accuracy has been attained (in the current program, the machine is turned off when the cumulative corrective error in cfm, is reduced to five times the number of airways). 5) Airway head losses H are calculated by the relation H = R Q2 and cumulated along the original designated flow route to find the overall head loss for the network. Trunk airways are then included to obtain the mine head.

Jan 1, 1964

By T. M. Morris

THE rate of flotation of solid particles determines the percentage recovery of these particles which can be obtained during a given time interval. It is an established fact that the recovery is greatest in the fine sand size range and falls off on either side of this size range. It is also a fact that the history of the development of concentrating devices has been marked by a desire to recover the small particles: tables were developed to treat material too fine to recover in jigs; flotation recovers material too fine to be recovered on tables. As our ores become more complex and require finer grinding, a greater proportion is in the slow-floating size range. It is therefore of economic importance that the effect of size of particle upon rate of flotation be studied. Many studies have been made along this line. The present investigation attempts to quantify the rate of flotation by means of a rate equation, and then compares the effect of size upon the rate of flotation. It was found that the flotation rate constant increases with increase in diameter of particle in the size range smaller than the optimum size particle. The flotation rate constant decreases with an increase in diameter of particle in the size range coarser than the optimum size particle. This suggested two different actions taking place, and these are discussed. Suggestions for further work in order to improve the flotation rate of both fine and coarse particles are made. The usual method of expressing rate is to plot cumulative percent recovery vs time. If we knew the equation of this kind of curve we could replot it as a straight line, and the slope of the line would be an index of the magnitude of the flotation rate. It is easier to compare straight lines than curves. The basis used for obtaining the equation of the recovery vs time curve is the law of mass action. This law states that the "rate of a chemical reaction is proportional to the `active masses&apos; of the reacting substances present at any given time." This, of course, implies a constant temperature. We can modify this for our use by using the term flotation reaction instead of chemical reaction and the term floatable particle instead of active masses. [ ] where C denotes concentration of particles, t denotes time, k denotes a proportionality constant. If the exponent of C is 1, then we have a first order rate reaction. If the exponent of C is 2, then we have a second order rate reaction. For a first order rate reaction, we obtain, after integration, [ ] where k is again a proportionality constant and has the dimension of the reciprocal of time. C, is the original concentration or weight or percentage of floatable particles. C, is the concentration or weight or percentage of floatable particles remaining in the flotation cell after time t. If all particles of the mineral to be floated were capable of floating we could take C, as 100 and C, would be 100 - R, where R is the cumulative percentage recovery of mineral at time t. But it is not usual for all of the particles to be capable of floating, because of surface oxidation, etc. Therefore, we must introduce a factor x which is the percentage of unfloatable mineral. Hence the first order equation is [ ] It has been proposed by various workers1-3 that this equation describes the flotation rate. More recently N. Arbiter4 states that a second order rate equation fits the recovery vs time curve more closely than does a first order rate equation. Morris&apos; has discussed the evidence presented for this statement. It has been concluded that the first

Jan 1, 1952

By J. H. Jackson, P. D. Frost, C. H. Lorig, L. W. Eastwood, A. C. Loonam

R. S. BUSK*—I wish first to congratulate the authors of this paper both for the work done and the presentation of that work. We have also been working on this type of alloy development, but any technical contribution I might have is better reserved for a separate paper which is planned in the near future. However, I would like to insert a word of caution. The alloys described in this paper are still very much in the laboratory stage. There still remains the quite vexing problem of property stability at slightly elevated temperatures which must be solved before the properties quoted can be effectively utilized. These alloys are extremely interesting and I am confident that the problems can be solved. However, the present existence of those problems must not be overlooked. I have one question with regard to the work-hardening test. If the work harden-ability of magnesium is compared to aluminum using a more standard test such as the Meyer Analysis the magnesium is found to work harden the more rapidly. Is the test developed for this work truly a work-hardening test or a measure of something else ? J. H. JACKSON (authors&apos; reply)—I am sure we appreciate the comments of Mr. Busk and we are pleased that the Dow Chemical Co. has also undertaken work in this field. It is a big field. We, of course, are continuing our research for the Navy Bureau of Aeronautics, since, as clearly shown in the paper, the alloys are not completely developed. With regard to Mr. Busk&apos;s comments on the relative work-hardenability of aluminum and magnesium, the comparison we cited in the paper was intended merely to show that the results obtained by our method of testing were independent of the yield strength of the materials. We believe that this test method afforded a rough indication of the work-hardenability of the magnesium-lithium base alloys. It is not known whether this method affords greater or less accuracy than the Meyer Process, which is based upon a hardness-test impression. It is quite possible, since work-hardenability is a function of metal purity, that our results, obtained for commercial aluminum and magnesium, would be different if high-purity metals had been used. A. C. LOONAM—Lithium has some very interesting effects on magnesium. The paper showed that even small percentages increased ductility. These small amounts changed the axial ratio of the hexagonal magnesium from 1.624 in the pure metal to 1.606 in the saturated alpha solution, a value close to that of titanium, which has a ratio of 1.601. You will recall an article in the January, 1949, issue of the Journal of Metals where it was stated that titanium can be cold-rolled over 90 pct without any significant edge cracking. At 4.9 or 5 pct lithium, a change begins in crystal structure—from a hexagonal close packed to body centered cubic. There are two phases over a short range of lithium content, but at 10 pct the structure becomes completely cubic. The melting points of alloys in this range of lithium content are relatively high. Alloys even up to the six ratio have melting points above 500°C; this represents a relatively small reduction in melting point from pure magnesium itself. The magnesium-lithium system has a number of very interesting characteristics. This is only a very small part of the story. J. H. JACKSON—I would like to remind you that Mr. Loonam is one of the authors of this paper, and the fact is that this work is based on many of his ideas. I am sure we appreciate his additional comments in this discussion. R. LEITER*—I am interested in Fig 19 in which Mr. Jackson compares

Jan 1, 1950

By M. T. Hepworth, T. P. Papazoglou

Hydrogen diffusion in hcp and bcc titanium in the temperature range 610° to 900°C was measured by reacting cylindrical specimens of high-purity titanium at constant temperature with small constant pressures of hydrogen of the order of a few Torr. In the range 0 to 6.7 at. pct H in hcp titanium, the diffusion coefficient D in sq cm per sec is given by: where T is in °K and R in cal per g-atom OK. Only three measurements of diffusivity in bcc titanium were made: at 900&apos; : the diffusivity in ß titanium was found to be 1.8 x 10-4 sq crn per sec. McQUILLAN1 has studied the constitution of the Ti-H system and reported the concentration of hydrogen as a function of pressure of hydrogen and temperature. He found hydrogen absorption was reversible and rapid at temperatures above 500°C in both crystallographic modifications, the a (hcp) phase and the ß (bcc) phase, with the latter phase stabilized by hydrogen additions. Wasilewski and Kehl 2 measured the diffusivity of hydrogen in hcp and bcc Ti-H alloys. For bcc alloys, they reacted cylindrical specimens of titanium with hydrogen gas at 1 atm pressure. They measured the rate of absorption of hydrogen isothermally on specimens ranging from zero hydrogen content when above the transformation temperature to hydrogen contents sufficient to stabilize the bcc phase when below the transformation temperature. By integrating Fick&apos;s second law for a cylindrical geometry, assuming a diffusion coefficient independent of composition, they calculated diffusivities as a function of temperature. In the range 650" to 1000°C their results in bcc titanium yielded the expression: where D is in sq cm per sec, T in OK, and R the gas constant in cal per g-atom-OK. Their studies on the hcp a phase were conducted by exposing long cylindrical titanium specimens to hydrogen at a sufficiently high pressure to establish a surface layer of fl phase 80 to 120 u thick for the duration of an isothermal experiment. Concentration profiles were then determined by hydrogen analysis from chips obtained upon machining quenched specimens. The diffusivity of hydrogen was then calculated yielding the expression: The present investigation was undertaken principally to redetermine the diffusivity data on a titanium by a direct absorption from the gas phase at sufficiently low hydrogen pressures to avoid the presence of a surface bcc phase and also to check the diffusivity data at low hydrogen concentrations in the bcc stmcture. EXPERIMENTAL PROCEDURE High-purity electrolytic sponge titanium (99.9 pct) of BHN 91 was arc-melted, forged, rolled, and machined into cylinders of final hardness, upon annealing, of BHN 92. Cylindrical specimens of two sets of dimensions were used: 1.524 by 2.931 cm and 1.905 by 3.663 cm diam and length, respectively. The experimental apparatus used for annealing specimens and for admitting known amounts of hydrogen at constant pressure is shown in Fig. 1. In a typical experiment, a specimen was lightly abraded with fine-grain silicon carbide, washed with water and detergent, and then thoroughly degreased with reagent-grade tri-chlorethylene and acetone. The specimen was suspended by a Chrome1 wire in the center of a quartz furnace tube, without contact with the tube, and heated to 925°C (as measured by a thermocouple adjacent to the specimen, outside the tube) for a period of at least 4 hr under a vacuum of from 10-6 to -7 Torr as measured by an ionization gage. The furnace temperature was then adjusted to the temperature chosen for the particular diffusion study and the specimen allowed to equilibrate thermally. Hydrogen was then ad-

Jan 1, 1969

By T. M Morris

THE rate of flotation of solid particles determines the percentage recovery of these particles which can be obtained during a given time interval. It is an established fact that the recovery is greatest in the fine sand size range and falls off on either side of this size range. It is also a fact that the history of the development of concentrating devices has been marked by a desire to recover the small particles: tables were developed to treat material too fine to recover in jigs; flotation recovers material too fine to be recovered on tables. As our ores become more complex and require finer grinding, a greater proportion is in the slow-floating size range. It is therefore of economic importance that the effect of size of particle upon rate of flotation be studied. Many studies have been made along this line. The present investigation attempts to quantify the rate of flotation by means of a rate equation, and then compares the effect of size upon the rate of flotation. It was found that the flotation rate constant increases with increase in diameter of particle in the size range smaller than the optimum size particle. The flotation rate constant decreases with an increase in diameter of particle in the size range coarser than the optimum size particle. This suggested two different actions taking place, and these are discussed. Suggestions for further work in order to improve the flotation rate of both fine and coarse particles are made. The usual method of expressing rate is to plot cumulative percent recovery vs time. If we knew the equation of this kind of curve we could replot it as a straight line, and the slope of the line would be an index of the magnitude of the flotation rate. It is easier to compare straight lines than curves. The basis used for obtaining the equation of the recovery vs time curve is the law of mass action. This law states that the "rate of a chemical reaction is proportional to the &apos;active masses&apos; of the reacting substances present at any given time." This, of course, implies a constant temperature. We can modify this for our use by using the term flotation reaction instead of chemical reaction and the term floatable particle instead of active masses. dC Thus----— kc [1] where C denotes concentration of particles, t denotes time, k denotes a proportionality constant. If the exponent of C is 1, then we have a first order rate reaction. If the exponent of C is 2, then we have a second order rate reaction. For a first order rate reaction, we obtain, after integration, where k is again a proportionality constant and has the dimension of the reciprocal of time. C, is the original concentration or weight or percentage of floatable particles. C, is the concentration or weight or percentage of floatable particles remaining in the flotation cell after time t. If all particles of the mineral to be floated were capable of floating we could take C, as 100 and C, would be 100 - R, where R is the cumulative percentage recovery of mineral at time t. But it is not usual for all of the particles to be capable of floating, because of surface oxidation, etc. Therefore, we must introduce a factor x which is the percentage of un-floatable mineral. Hence the first order equation is It has been proposed by various workers1-&apos; that this equation describes the flotation rate. More recently N. Arbiter4 states that a second order rate equation fits the recovery vs time curve more closely than does a first order rate equation. Morris- as discussed the evidence presented for this statement. It has been concluded that the first

Jan 1, 1953