Search Documents

By J. C. F. Randolph

This mill being now in successful work, it seems proper to record some professional memoranda concerning it. It is intended to handle an unusual ore, in a locality difficult of access, and its arrangement and method of procedure present some novel features of interest. It was prepared for the San Miguel property at Batopilas, State of Chihuahua, Mexico, owned by the Consolidated Batopilas Silver Mining Company of New York, at present under the management of ex-Governor Alexander R. Shepherd. This property contains a large series of more or less parallel veins, worked through a deep cross-cut tunnel 1300 feet in length. The vein rock of these veins is calcapar in diorite. While, near the stirface, ruby silver, black

Jan 1, 1882

By John C. F. Randolph

THIS mill being now in successful work, it seems proper to record some professional memoranda concerning it. It is intended to handle an unusual ore, in a locality difficult of access, and its arrangement and method of procedure present some novel features of interest. It was prepared for the San Miguel property at Batopilas, State of Chihuahua, Mexico, owned by the Consolidated Batopilas Silver Mining Company of New York, at present under the management of ex-Governor Alexander R. Shepherd. This property contains a large series of more or less parallel veins, worked through a deep cross-cut tunnel 1300 feet in length. The vein rock of these veins is calespar in diorite. While, near the surface, ruby silver, black

Jan 1, 1882

By Louis E. Toth, Alan W. Searcy

A modification of Le Claire's microscopic model for self-diffusion is developed in a form suitable for prediction of activation energies for diffusion in disordered substitutional solutions as well as in pure metals Bonding is considered as a localized interaction, and the energy of bonding between atoms of different types is taken as the arithmetic mean of the energies in the pure elements. The activation energies for vacancy formation and migration in substitutional alloys are shown to depend on the empirical constants developed for self-diffusion when the equations are adjusted for the mole fractions of the two elements. The differences in activation energies between alloys and pure metals are calculated and compared with the experimental differences. MaNY different correlations have been proposed for use in predicting activation energies for self-diffusion or diffusion in dilute alloys. Usually theories for self-diffusion cannot be applied to alloys nor can theories for dilute alloys be satisfactorily adapted to self-diffusion. We present here a correlation scheme which can approximate activation energies for diffusion, vacancy formation, and atom migration in substitutional-alloy systems of any concentration. The activation energy for self-diffusion QD has been correlated to the melting temperature Tm,1, 2 the sublimation energy L,,l the heat of fusion Lf,3 crystal structure, and valence.4 Working with mechanical parameters Le claire5 proposed that the activation energy for vacancy migration Qm could be related to appropriate shear moduli. He also proposed that the activation energy for vacancy formation 62, is proportional to Ls. THEORY We consider that Le Claire's expression for 4, is satisfactory. But an alternate expression for Qm in terms of the bulk modulus is just as satisfactory from a theoretical point of view for the prediction of the energy of vacancy migration in pure metals and can more easily be used for alloys. In order for Qm to be a minimum, not only the moving atom but also adjacent and nearby atoms must undergo distortions. These distortions are neither isotropic as is implied by a correlation with the bulk modulus nor unidirectional as implied by a correlation with a shear modulus, but are complex and multidirectional. A complete analysis in terms of the bulk modulus and all the shear moduli is not yet available. In a first approximation in which only one parameter is used, there is no obvious theoretical reason to prefer the shear modulus over the bulk modulus. The bulk modulus is a much more convenient parameter for use in predicting Qm. Shear moduli are known for less than half of the metallic elements while bulk moduli are usually available. Furthermore, because of the directional dependence of shear moduli, their estimation for alloys appears inherently more difficult than is the estimation of bulk moduli. We assume that the bulk modulus Bs reflects an average value of the shear moduli such that the strain energy at the point of maximum distortion of the lattice during vacancy migration can be expressed as

Jan 1, 1964

By Erwin Pink

It has been pointed out in a previous publication' that, for a given strain rate, serrated yielding can be observed only at certain test temperatures. Fig. 1 shows the results obtained for polycrystalline and re-crystallized tantalum, Va group in the periodic table of elements, where serrations in the course of deformation match the decrease of total elongation. However, there was no evidence of repeated yielding when testing molybdenum or tungsten, VIa group, at similar conditions, i.e., annealed and air-cooled samples (serrations have been found in quenched tungsten2). The reason for the observed serrations in the stress-strain curve are dissolved impurities on interstitial sites of the lattice. Of importance is the concentration of oxygen, carbon, and nitrogen, which are in random solution. Two methods to calculate the limiting concentrations for serrated yielding have been evaluated and applied with more or less success in estimating the role of carbon and nitrogen in mild steel.3 J. Friedel4 has given two equations, which determine the upper and lower strain rates i for the appearance of serrations at a certain test temperature: W denotes the binding energy of the interstitial atom, -1 ev for refractory metals;5 e' the portion of continuous strain in the discontinuous stress-strain curve, measured as -3x10-3; b the Burgers vector, c the concentration, see Table 11, in atom fractions; p the dislocation density, -2 X l09 per cm2 at an average stress ?, -15 kp per mm2 at medium strain. k stands for the Boltzmann constant, T for the absolute temperature, and D for the diffusion coefficient. Table I. Eq. [1], for a given strain rate, determines the minimum temperature, and [2] the maximum temperature, below or above which no "dynamic strain-aging" can occur. According to the calculation, which uses the given quantities, serrations should only be observed within the cross-hatched temperature regions in Fig. 1, and this is in good agreement with the experiment. The effect of varying impurity concentrations is indicated schematically in Fig. 2: only Eq. [I] is affected. Setting Eq. [I] equal to [2], one can calculate the limiting concentration cmin which first causes serrations. The second method uses the aging condition by Cottrell and Bilby6,7 \£i/ \ n. a / where n(t) is the number of impurity atoms, which segregate at a dislocation line of 1 cm length after the time t (t in our case is the time between "repeated yield drops" in the serrated curve, or, the time until the above mentioned strain E' is obtained, i.e., -1 sec for medium intensity or serrations). no is the number of impurity atoms in 1 cm3 and A = Qb4ep[(1 + v)/(l - v)], derived from equations given in Ref. 7. p is the shear modulus, 18,800 kp per mm2 for tantalum, 12,200 for molybdenum, AE the tetragonal distortion caused by interstitial impurities in a bcc lattice (approximately 0.4 for both tantalum and molybdenum) and v Poisson's ratio (0.35 for tantalum, 0.30 for molybdenum). From Eq. [3], combined with the condition for a pronounced yield drop as formulated by cottrel16 (at least one impurity atom segregated at the unit length of a dislocation), again a no min, and hence the concentration Cmin can be calculated.

Jan 1, 1970

By J. M. Campbell, T. S. Buxton

The most widely used methods of predicting the volumetric properties of gas are based on the principle of corresponding states, which asserts that the compressibility factor is a universal function of reduced temperature and pressure. Previous studies have shown that the acentric factor, as proposed by Pitzer, l is an important addition to reduced pressure and reduced temperature as factors affecting the compressibility factor. Results of this study indicate that, if the pseudocritical temperature and pressure used to determine the reduced conditions are adequately predicted, characterization of natural gas-carbon dioxide mixtures with the acentric factor will allow reliable determination of the compressibility factor. Comparisons of predicted and experimental compressibility factors have shown that the pseudo-critical constant rules of Stewart, Burkhardt and Voo 2 are satisfactory for hydrocarbon mixtures. However, these rules fail to predict the pseduo-critical constants for hydrocarbon-carbon dioxide mixtures. Based on graphically determined pseudocritical temperatures for binary hydrocarbon-carbon dioxide mixtures, a correlation which gives the required correction to the Stewart, Burkhardt and Voo rules was prepared and a compressibility factor prediction technique was proposed. To test the proposed technique, compressibility factors for live mixtures of methane, carbon dioxide and either ethne or propane were experimentally determined at 100, 130 and 160F and pressures up to 7,026 psia. The predicted and experimental compressibility factors for these live mixtures had an average absolute deviation of 0.55 percent. INTRODUCTION Two-parameter generalized correlations based on the principle of corresponding states (PCS) are presently available for reliably predicting the compressibility factor of lean natural gas, i.e., natural gas with low concentrations of hydrocarbons heavier than methane. When these same correlations are used for mixtures of natural gas and carbon dioxide, large deviations between actual and predicted compressibility factors are observed. This study was undertaken to provide a means for accurately predicting the compressibility factor for such mixtures. In the past, two avenues of approach have been employed to extend the applicability of the PCS. The first has been to introduce additional parameters. The second has been to develop combination rules for predicting pseudocritical constants which do not suffer from the limitations of the commonly used molal average constants. In this study, both of these approaches have been considered and employed to arrive at a method for predicting the compressibility factor of mixtures of hydrocarbons and carbon dioxide. SELECTION OF A THIRD PARAMETER Failure of the original PCS to predict the compressibility factor for all pure gases, regardless of their mass, shape or polar moment, has led to the introduction of additional parameters into the PCS. Parameters have been introduced to correct for quantum deviations,3-5 non-spherical or globular shapes6-9 and for polar moments.8-12 In addition to these parameters which are based on some microscopic property and intended to correct for some specific cause of deviation, more general parameters based on some bulk property have also been introduced.1913-16 To use additional parameters to correct for each of the possible causes of deviation among the individual constituents of a gas mixture would require a rather complex form of the PCS. To maintain the PCS in as simple a form as possible, it is more desirable to have a single third parameter which is based on some bulk property influenced by several of the factors causing deviations. The third parameter selected for use in this study is a factor originally proposed by Pitzer.9

By A. C. Emery

The physical and chemical properties for the metalliferrous feed materials to the Imperial Smelting Process blast furnace are defined and the means of achieving these properties by the sintering or briquetting routes discussed. Modern sinter plant performance data is included with a plant flow sheet showing draughting patterns. The necessity for ensuring an adequate SO2 content in sinter machine gases when roasting sulphide materials to enable conversion to sulphuric acid in a contact acid plant is described. The ability of the process to consume limited quantities of oxide materials, residues from other processes and other secondary sources is discussed. The argument for a briquetting route is developed out of this limitation and the merits of cold and hot briquetting examined. Development work on the hot briquetting route for this purpose is briefly described and several flow sheets outlining possible applications for hot briquetting to different raw material situations including the preparation of low lead charges for the furnace are presented.

Jan 1, 1977

By E. D. Glass, R. C. Hessing

The volume of gas in storage reservoirs may be compuled from estimates of hydrocarbon pore volume and gas density. However, both are dificult 10 estimate accurrately. Further, no adequate method has been presented for estimating reservoir performance during operation for Volumetric gas-storage reservoirs. Normally, gas-storage reservoirs exhibit small pressure gradients in the area conr(zilliitg the storage wells, even under conditions of maximum injection or withdrawal. A zone of low permeability usally surrounds, and is in communication with, tile permeable zone which causes the storage reservoir to exhibit a definite pressure lag, as evidenced by lack of pressure stabilization after extended periods of shut-in. An analysis is presented for determining pore volume and gas in place, and for predicting future reservoir behavior of such reservoirs. The mathematical development is presented and a solution is shown for an operating storage reservoir. INTRODUCTION Depleted gas reservoirs arc often used for gas storage. The reservoir usually contains a permeable area in which numerous storage wells arc located. The permeable zone permits high deliverability and injectivity with relatively small pressure loss. However, pressure equilibrium may not be attained even after long periods of shut-in because of the fact that the permeable zone is in communication with an adjacent tight section. Most of these gas-storage reservoirs are sand lenses, such as ancient beaches or offshore bars. The areal geometry of these reservoirs is normally such that flow in the tight zone will be approximately linear, and the solutions presented here should apply. DISCUSSION For purposes of this paper, a gas-storage reservoir is considered to consist of a volume of permeable rock in communication with tight rock, all of which is surrounded by an impermeable barrier. Pressure gradients across the permeable zone are small due to low flow resistance and well density. Hence, the permeable zone is considered to be at a uniform pressure which is only a function of time as gas is injected or withdrawn. The tight zone is considered to be linear with uniform permeability and porosity, and the pressure in the tight zone is a function of both position and time. l'wo solutions for flow in this system are used to define the reservoir and determine its properties. The first solution is for an infinite linear tight zone in communication with the permeable zone. The expression for constant gas-withdrawal rate is Eqs. I and 2 ore superposed to account for varying rates of production or injection and are used to fit gas-storage field reservoir history. The solution for the infinite case, Eq. 1, is first applied to field data for early times following a shut-in period. Constants BDa and c are ohtained by a least-squares fit of pressure and rate data. After the pressure wave from gas injection or withdrawal reaches the external boundary of the tight zone, the infinite solution will no longer apply and pressure behavior should be more nearly represented by Eq. 2 for a gas reservoir in communication with a finite linear tight zone. The values of the constants found from curve-fitting Eq. 1 to the early time data may then be used in Eq. 2 that only the coefficient- L containing the unknown

By J. B. Conway

RECENT studies of first-stage creep data have led to a special least-squares procedure for use in calculating the Andrade creep-equation constants. This procedure is easy to apply, uses only experimental data, and yields results which do not vary with the individual making the analysis. Taking logarithms of the Andrade equation I = lo(1 +ßl1/3)ekt leads to In l = In lo + In(1 + ßl1/3) + kt A residual can then be defined as the difference between the logarithm of an experimental length value and the logarithm of a calculated length value. Following the usual least-squares approach differentiation with respect to 10 and k of the expression for the sum of the squares of the residuals leads to two equations which cannot be solved directly for the three constants involved. However, values of ß can be assumed and corresponding values of lo and k can then be calculated. Then the sum of the squares of the residuals can be calculated. For a certain value of ß this sum will exhibit a minimum value to identify the least-squares values for ß, lo, and k. An illustration of the application of the above least-squares approach is obtained by employing the creep data for lead originally presented by Andrade' in the development of the one-third creep law. A plot of the sum of the squares of the residuals as a function of ß is shown in Fig. 1 along with the corresponding lo and k values. The minimum is seen to occur at ß = 0.0143 from which lo = 39.149 and k = 7.27 x 10-5. Based on these results the Andrade equation for these lead data is 1 = 39.149(1 + 0.0143 t1/3)e7.27 x 10-5t compared to the equation I = 39.13(1 + 0.01457/l/3)e6.95 X 10-5t developed by Andrade (calculation procedure not given but was probably a graphical solution). Calculating length values using the above equations and comparing with the experimental length values, residuals corresponding to the difference between experimental and calculated length values can be obtained. The sum of the squares of the residuals can then be calculated. For the case of the least-squares approach a value for this sum was found to be 2.37 x 10-3 compared to a value of 3.47 x 10" for the equation originally proposed by Andrade. In other words, while the constants proposed by Andrade provide an excellent fit for the lead data, a still better representation is afforded by applying the least-squares approach discussed above. A similar result was obtained for the lead data at 162oc.l The least-squares approach led to 1 = 9.866(1 + 0.0497 t1/3)el.011 x 10-2t compared to the equation I = 9.94(1 + 0.0433 l1/3)e1.06 X 10-2t given by Andrade. In this case, too, the sum of the squares of the residuals is slightly smaller for the equation based on the least-squares approach. It should be pointed out that a double trial and er-

Jan 1, 1965

By E. F. Sumner, C. K. Eilerts

A partial differential equation was developed and programmed to compute the cumulation and flow of the liquid phase of a gas-condensate fluid. It was found that such cumulation near the wellbore may attain saturations permitting liquid flow within periods of the order of 1 day. The radius at which this saturation obtains advances toward the external reservoir boundary with continued recovery, although it may move only 2 to 5 percent of that distance during the productive life of the well. This relatively high saturation can cause a serious dimunition in delivery capacity. The pressure dilferential required for flow of gas through a formation with such saturation may be as much as three times greater than the corresponding differential in a formation where saturation is less than 2 percent. Because of the radial geometry, total reservoir pressure differential required for flow with a zone of relatively high liquid saturation about the well was indicated to be two to three times greater than without the saturation. INTRODUCTION A computing program developed by the authors 1 to integrate the partial differential equation for transient radial flow of gas-condensate fluids required a prevalued function, W(Po) = kD(p0)/M(p0)z0(Po).......(1) related to cumulation of condensate in pores of the reservoir as pressure declines. Methods were available for relating mobility to saturation, saturation to liquid-vapor volume ratio, liquid-vapor volume ratio to pressure and, hence, mobility to pressure to obtain W(pD). However, certain assumptions had to be made in advance of program execution concerning movement of the liquid phase and change in radial distribution of liquid saturation as recovery progressed. Because reservoir fluid mobility and related pressure gradients are markedly dependent on liquid saturation, it was desirable to eliminate the necessity for these assumptions and in this respect improve the reliability of a result obtained by integration. Experimental research was begun to determine the relationship of variables that control the reservoir flow of low liquid-vapor ratio fluids. With the porous medium partly saturated, permeability to a two-phase flowing fluid could not be measured because the viscosity was unknown. Viscosity could not be measured without knowing permeability of the liquid-containing medium; but mobility, the ratio of permeability to viscosity, was measurable. This property, and liquid-vapor volume ratio of the flowing phases, have values determined by nature of the fluid and porous medium and, in addition, were found to be pressure, velocity and saturation dependent. A program developed by Kniazeff and Naville 2 for computing multicom-ponent, two-phase flow used polynomials to represent saturation dependence of permeabilities for gas and for liquid phases, and an effluent viscosity constant under reservoir conditions. The authors developed a program based on a partial differential equation of similar mass balance origin but designed to represent distribution of an unlimited number of components between gas and liquid phases utilizing the experimental information available as to variable dependence. Attention was given to liquid-vapor volume ratio of the flowing phases, a property which is essential as a liquid saturation boundary condition. A program for computing implicitly saturation increase at uniform intervals along the drainage radius was combined with an existing program for calculating pressures, and the two were used alternately in an iterative solution for both saturations and pressures. The function W(pD) needed for calculating pressures was replaced by the ratio of mobility to compressibility factor. BASIC EQUATIONS The partial differential equation for calculating

By H. H. Rachford, and J. D. Rice

The effectiveness of digital computing machines in making technical calculations depends on how well the work is arranged to utilize the capability of the machines. This note presents a particularly useful way of calculating hydrocarbon vapor-li.quid equilibrium in the flash vaporization (or condensation) system. The method is well suited to sequence-controlled computing equipment. It is not limited to equilibrium calculations and may be used for solution of most implicit equations in one variable. INTRODUCTION There is increasing interest in the use of electronic digital computers in research and engineering calculations. This is a fortunate and inevitable trend in view of both the increasingly extensive numerical work which is becoming a routine part of many daily production operations and growing demand for the overwhelming amounts of calculations required by newly developed numerical methods for solving heretofore unsolved problems. Machines are in many ways ideally suited to the task but of necessity present certain difficulties. for a particular poblem to be solved must often be formulated quite differently from the way it would be arranged for manual solution. This is done in order to take advantage of the inherent speed and precision of electronic computers and at the same time to limit the need for number storage to the capacity of the machine. Therefore, this note is submitted to present a general and quite powerful method of finding solutions of the frequently encountered implicit equation: F(x,y1,y2,...ya)=0.......(1) where the root, x, is to be found for a given set of y1, ... yn.. The procedure is well suited for use with computing machines, for it usually requires but little storage or programming beyond that necessary to evaluate F. HYDROCARBON EQUILIBRIUM The method is described in terms of the problem it was designed to solve, i.e., the calculation of hydrocarbon vapor-liquid equilibrium in flash vaporization. Given the composition of a hydrocarbon mixture and appropriate values of the equilibrium ratios K, to find the phase ratio and compositions in a closed system: let z, be the mol fraction of the i-th component in the mixture. If K, is the ratio yl/xl, where y, and x,

Jan 1, 1952

Apparatus and a method for determining the volumetric behavior of hydrocarbons at pressures up to 10,000 psia and at temperatures between 0° and 460°F are described. The equipment is suitable for measuring the total volume of the system and the volume of a liquid phase as functions of state within these ranges of temperature and pressure. Introduction During the past two decades there has been a marked increase in the ranges of temperature and pressure encountered in underground hydrocarbon reservoirs of commercial interest. In addition, greater interest in the intensive characteristics of the fluids found in these reservoirs has developed. Apparatus Several types of apparatus have been developed in the authors' laboratory for volumetric measurements of hydrocarbon systems. The earliest equipment described1 was characterized by constant volume, and the weight of material present was varied systematically in order to attain the different pressures desired in the study. Later another apparatus2 was built in which the effective volume of the cylindrical cell, which confined a hydrocarbon sample of constant weight was varied by injection or withdrawal of mercury. The effective total volume of the equilibrium cell was measured by ascertaining the position of the mercury surface but no measurement could be made which gave knowledge of the volumes of individual phases existing within the hydrocarbon system. A third apparatus3 was designed for use at higher pressures and temperatures than the earlier equipment would tolerate. This apparatus was of the variable volume type but was built in the general form of a U-tube, the mercury level being measured at the surface which was not exposed to the hydrocarbon system. This arrangement avoided the need of subjecting measuring equipment to the extremes of temperature at which studies were made, but it gave no provision for measurement of the volumes of individual phases. Although this last apparatus was satisfactory as to accuracy and convenience for many purposes, need was felt for measurements of the volume of the liquid phase of the hydrocarbon system. For this reason a fourth equilibrium apparatus was designed and constructed, and it is the purpose of the present paper to describe it. The equilibrium cell has a working range of temperature from 0° to 460°F at pressures up to 10,000 psi. Provision for measurement of both total volume and liquid-phase volume makes it particularly suitable for the equilibrium study of heterogeneous systems consisting either of simple mixtures or of complex field samples of hydrocarbons. The general arrangement of the ap-

Jan 1, 1948

Apparatus and a method for determining the volumetric behavior of hydrocarbons at pressures up to 10,000 psia and at temperatures between 0° and 460°F are described. The equipment is suitable for measuring the total volume of the system and the volume of a liquid phase as functions of state within these ranges of temperature and pressure. Introduction During the past two decades there has been a marked increase in the ranges of temperature and pressure encountered in underground hydrocarbon reservoirs of commercial interest. In addition, greater interest in the intensive characteristics of the fluids found in these reservoirs has developed. Apparatus Several types of apparatus have been developed in the authors' laboratory for volumetric measurements of hydrocarbon systems. The earliest equipment described1 was characterized by constant volume, and the weight of material present was varied systematically in order to attain the different pressures desired in the study. Later another apparatus2 was built in which the effective volume of the cylindrical cell, which confined a hydrocarbon sample of constant weight was varied by injection or withdrawal of mercury. The effective total volume of the equilibrium cell was measured by ascertaining the position of the mercury surface but no measurement could be made which gave knowledge of the volumes of individual phases existing within the hydrocarbon system. A third apparatus3 was designed for use at higher pressures and temperatures than the earlier equipment would tolerate. This apparatus was of the variable volume type but was built in the general form of a U-tube, the mercury level being measured at the surface which was not exposed to the hydrocarbon system. This arrangement avoided the need of subjecting measuring equipment to the extremes of temperature at which studies were made, but it gave no provision for measurement of the volumes of individual phases. Although this last apparatus was satisfactory as to accuracy and convenience for many purposes, need was felt for measurements of the volume of the liquid phase of the hydrocarbon system. For this reason a fourth equilibrium apparatus was designed and constructed, and it is the purpose of the present paper to describe it. The equilibrium cell has a working range of temperature from 0° to 460°F at pressures up to 10,000 psi. Provision for measurement of both total volume and liquid-phase volume makes it particularly suitable for the equilibrium study of heterogeneous systems consisting either of simple mixtures or of complex field samples of hydrocarbons. The general arrangement of the ap-

Jan 1, 1948

By Leigh J. Wardle

SUMMARY This paper describes a numerical method for calculating the three- dimensional stresses and displacements induced by tabular excavations in a rock mass consisting of parallel layers of distinct material properties. The method is based on the superposition of numerical solutions for elementary displacement discontinuities in a parallel layered rock mass. The method can be used to analyse complicated excavation layouts such as those encountered in room-and-pillar or longwall mining. The method can automatically take into account the effect of the earth's surface. The method offers an economic and efficient solution to stress analysis problems which cannot be adequately modelled using two- dimensional methods and for which the expense of three-dimensional finite element or general boundary element analysis cannot be justified. INTRODUCTION Underground mining often takes place in tabular or lenticular orebodies. In many cases, for example room-and-pillar and longwall mining, the excavation geometry cannot be adequately modelled by two- dimensional methods of analysis. In these cases resort must be made to three-dimensional methods of analysis. Although in principle three-dimensional finite element models can be used for these geometries, cost limitations have restricted their application to a small number of cases. For example, room-and-pillar systems of regular pattern have been analysed by numerically modelling

Jan 1, 1984

By R. E. Collins, Paul B. Crawford

Mathematical equations have been derived to show the effect of the slippage phenomomenon (Klinkenberg effect) on unsteady-state gas flow through porous media. The assumption is made that the slippage may be represented by an equation of the type It is pointed out that most of the existing data on unsteady gas flow may be corrected for slippage by a simple change in definition of the dimen-sionless groups. INTRODUCTION Several solutions to the various problems associated with the unsteady-state flow of gas through porous media have been recently presented in the literature.','" The cited publications have resented solutions for flow under ideal conditjons. The principal assumptions involved include: (a) the density of the gas obeys the perfect gas law; (b) the viscosity of the gas is constant; and (c) the gas permeability is constant. Few, if any, actual gases under reservoir conditions are believed to conform with the above, but the ideal solutions do provide the basis for making improved estimates on the performance of gas reservoirs. The purpose of this paper is to show the performance characteristics when assumption (c) is corrected for the slippage phenomenon or the Klinkenberg effect. This is believed to be a step nearer to reality. MATHEMATICAL DEVELOPMENT Case 1. Gas Flow — Ideal Conditions The equations describing the flow of gas under ideal conditions have been clearly presented in the cited literature; however, for completeness of this work they are repeated here. The basic differential equation describing the flow of gas is Equation (6). It may be derived as follows: The continuity equation, the perfect gas law and Darcy's law are combined to give: Manuscript received in the Petroleum Branch office Aug. 3, 1963 or If — is constant, Equation (5) is rewritten to form the basic Equation (6). or for the linear case A dimensionless equation is obtained by dividing through by (pl)z and making the substitutions shown below: The resulting equation is A similar procedure for the radial case leads to the following equations: where and Equations (11), (12), and (13) are those used by Bruce, et al., and Jenkins and Aronofsky for their computational work.

Jan 1, 1953

By Richard Moldenke

So far as the quality of the product is concerned, the history of the production of pig iron for foundry purposes is one of constant retrogression. The steps in this deterioration began with cold-blast charcoal pig iron, then anthracite iron, coke iron; then gradually warm to hot-blast coke iron, with the charcoal furnaces also heating their blast to get greater tonnages, and finally the present-day hot-blast coke irons with scrap additions to the ordinary burden that, according to one European report, have gone as high as 65 per cent. Parallel with this quality retrogression is the enormous increase in furnace tonnages from about 15 to over 1000 tons per day in excep tional instances of modern practice. The key to this situation is economic pressure. The effect is a growing differentiation between furnace production for gray iron and malleable foundries and for the production of commercial steel. As the foundry is often compelled to draw upon the pig irons made for steel purposes, for economic reasons, whereas the steel industry does not draw upon foundry pig irons proper, the situation is resolving itself into one of grave consequence for the life of the gray-iron castings industry, and hence the intensive present-day feeling in foundry circles for rigid specifications to govern the quality of pig irons sold for foundry use. In the discussion of a paper on the use of scrap in the blast furnace held recently in the Institute, a very pertinent remark was made by one of the participating furnacemen, who confessed that "the things found wrong by the foundryman with his iron were beyond his depth." It may be added that those things are also beyond the depth of most foundrymen who, however, know from their results that they are getting poorer castings today with many pig irons than they did with those same brands in former years. Under the former fracture buying, the foundryman at least could select the strongest metal, as judged by the sledge, from various brands of similar fracture and pig section, and his results proved

Jan 1, 1927

By G. A. Crowe, W. T. Bishop

WHILE the literature in the United States on general polarographic technique is extensive, very little appears to be directly applicable to the analysis of flotation tailings. For this reason, it seemed worth while to present a few polarographic methods of ore analysis that have been in constant use in the laboratories of the Hercules Experiment Station for the past few years. At the time of our initial investigation of the polarographic method, there was little or no mention of its specific application to flotation tailings. It was therefore necessary for us to establish suitable conditions for the analysis of our flotation tailings and to consider possible interfering constituents. The polarograph, or dropping mercury electrode, method was discovered in 1925 by J. Heyrovsky and M. Shikata,1 of the Charles University, Prague. This method has since received extensive development and the literature contains many articles on the use of the polarograph for the quantitative determination of metals as well as anions, oxygen in solution, and reducible organic compounds. Kolthoff and Lingane,2 by the publication of their book, "Polarography," in 1941 did much to advance the knowledge and use of the polarograph. In chapter 25 of that book, methods for the polarographic analysis of brass, steel and iron products, alloys, and zinc ore are given. Aluminum,3 cadmium,3.4 cobalt,4 copper,3-9 gold,10 iron,3,11 lead, 3,12,13 molybdenum," nickel,4.8,9 tin,15 and zinc3,4,5,7 in ores have been analyzed by the polarograph. However, most of the work on tailing samples has been published in Russian journals and translations are not yet readily available. In this paper we will describe the method that was developed and used in our laboratory for the determination of copper, zinc, and lead in flotation tailings. PRINCIPLE OF THE METHOD The polarograph is an instrument for applying a potential to a mercury pool and a dropping mercury electrode separated by a solution, and for measuring the resultant electrical current. The procedure is essentially a small-scale electrolysis, in which the voltage required to reduce or oxidize the material identifies it qualitatively and the current measures the concentration quantitatively. Different constituents in a sample may often be determined simultaneously provided they reduce at sufficiently different voltages. An example is the determination of copper and zinc in the same solution. The reduction potentials are quite far apart, as shown in Fig. I. The polarographic method is particularly useful in the range of concentration between those normally determined by spectographic and chemical methods. In

Jan 1, 1946

By M. J. Scoble, N. Rowlands, L. Piciacchia

Pronounced non-linear behaviour of mine backfill arises from the influence of confining pressure, cementing agent content, curing time, void ratio and water content. Past research into backfill behaviour under stress is reviewed and related to the needs for constitutive laws which can be integrated into numerical modelling for mine design where structures employ backfill support. A method for representing non-linear backfill behaviour is proposed, which is based upon numerical regression treatment of triaxial laboratory data. The non-linear behaviour observed from a program of laboratory triaxial testing of cemented hydraulic backfill is presented. The particular influence of confining stress and curing duration is high- lighted. The paper concludes by considering the alternative means of accounting for such behaviour in mine design.

Jan 1, 1984

By T. W. Peterson, J. W. White, E. E. Krist

In developing a minicomputer application for an industrial grinding circuit, the control program is usually coded and the compilation is debugged on a remote access computing system. Final debugging must wait for installation of the minicomputer. This work reports on an alternative. Here both compilation and execution bugs are eliminated before the hardware is installed by testing the control program against a dynamic simulation. Attention is directed toward development of a particular real-time control program for an industrial wet grinding circuit in a copper concentrator. Specific data were supplied by the Anamax Mining Co. for a grinding division consisting of a rod mill followed by two ball mills in parallel. This information was used to set up a plausible model structure.

Jan 1, 1977

By J. H. Swisher

HIGH-temperature anneals in wet hydrogen are widely used in research as a technique for decarburiz-ing Fe-C alloys. This technique is also used commercially in processing low-carbon steels for magnetic applications. The results of some calculations to be presented show that for temperatures up to 1100°C there is a critical temperature for maximizing the decarburization rate of alloys containing about 0.1 pct C. The low-temperature portion of the iron side of the Fe-C phase diagram1 is reproduced in Fig. 1. Preliminary calculations show that the critical temperature is between 727" and 905°C. In this range, the alloy contains two phases during heat treatment, ferrite and austenite. Decarburization therefore occurs by a process termed "diffusion accompanied by a phase change".' For intermediate stages of decarburization, the carbon concentration profile in a flat plate is shown in Fig. 2. The carbon content increases linearly from zero at the surface to a value, cB , at the two-phase boundary. cB is known from the phase diagram, see Fig. 1. At the interface where x = 5, the carbon content increases to ci, the initial concentration in the alloy. smith3 has derived an equation which describes the rate of movement of the boundary. With a rearrange-

Jan 1, 1969

By R. B. Gilmore

It has been the custom of the Railroad Commission of Texas to hold hearings each month with respect to existence and imminence of waste of oil and gas in Texas and its prevention, and to issue orders restricting the production of oil to the reasonable market demands and allocating production to the various fields in Texas. The practice of designating the actual days during the month on which wells are to be shut down was discontinued in the order for November 1942, and orders issued thereafter contain only the number of shut-in and producing days for the month. Orders issued by the Railroad Commission are divided into two main subdivisions, the first a general shutdown order of all oil wells in the various fields in Texas, and second, the exception of certain fields to this general order. Table I shows the producing days for all of the fields in Texas for 1944 under the first of these classifications : The fields that have producing schedules varying from the general order may be divided roughly into three subdivisions, as follows: Table 2 shows the producing days for certain fields that were prorated as a group throughout 1944. Other fields which were occasionally included with this group were also prorated under the general order and by special orders. The producing days for these irregular fields are shown in Table 4. In Table 3 are listed fields that were exempted from any shutdown days during December 1944, the majority retaining this exception throughout 1944. Certain fields included in the general order were added to this group on the first day of the month shown.

Jan 1, 1945