Search Documents

THE CHAIRMAN (G. C. STONE, New York, N. Y.).-It is a matter of history that the smelting of zinc ore in, a reverberatory furnace was the first process used in this country for making oxide of zinc, and it was abandoned for three reasons: it was very expensive, the recovery was, very poor, and the color of the oxide was worse. It is probable, however, that the Anaconda company is working in a totally different way. While Mr. Ingalls does not say so, my own suspicion is that they are using coal-dust firing, very much as they do in the copper furnaces, and by that method there is much more chance of getting good recovery.

Jan 1, 1918

By L. G. Jones

An approximate method of calculation is developed in this paper which allows duplication of radial unsteady-state gas flow computer results where Darcy's law applies, such as those reported by Aronofsky and Jenkins and Bruce, Peaceman, Rachford, and Rice. Moreover, the new calculation method can be used to obtain results for radial unsteady-state gas flow obeying the quadratic flow law proposed by Duwez and Green. Means are discussed for predicting well behavior at single or superimposed flow rates in finite or infinite reservoirs, determining reservoir rock properties from well-test data, reproducing and interpreting back-pressure test data, and determining the radial extent and reserves of gas reservoirs from well-test data. Example calculations are presented for gas flow following both Darcy's law and the quadratic flow law. INTRODUCTION Since the publishing of U. S. Bureau of Mines Monograph 7,1 most gas-well testing methods have been based on the equation q = y (pl2 - pw2)a, where q = production rate, pl and pw, are formation pressure and sandface pressure, respectively, and y and a are constants to be obtained from test data. These methods, used for predicting both deliverability and "open-flow" capacity of gas wells, have been useful and accurate in many cases but unsatisfactory in others. Even at best, however, they do not supply information about the formation or lead to an understanding of nonsteady-state gas flow in porous media. Many theoretically based studies2-4 of gas flow obeying Darcy's law have been made.' Since the partial differential equations which result from combining Darcy's flow law with the continuity equation are nonlinear, most of the published research consists of either numerical solutions or analytical solutions for linear approximating equations. Such solutions have been of limited value in field work due to their unhandy form and their failure to correlate most field data. There is evidence5,6 which indicates that Darcy's law is inadequate to describe gas flow at some flow rates of practical interest. A quadratic flow law: which reduces to Darcy's law at low rates, is more successful in accounting for experimentally observed behavior. This flow law has been applied successfully to a few hypothetical reservoir cases in work which has recently been published.10 However, these numerical solutions of the equations involved have been successful only on a special analog computer. Routine applications to field cases would be awkward and have not been attempted. The present paper describes an approximate method for computing nonsteady-state gas flow solutions which has been completely successful in predicting the results for both Darcy flow and quadratic flow obtained by elaborate numerical methods. The new calculation method allows determination of the observable variables in gas-well testing at constant rates. It is similar to the scheme of using a succession of steady states suggested by Muskat7 in that it makes use of steady-state and material-balance equations. It also is similar to the work of Aronofsky and Jenkins2 in that the new method includes Aronofsky and Jenkins' stabilized flow equation as a special case. It improves upon both of these calculation schemes in that it accurately describes all portions of reservoir history and suggests means of determining reservoir rock properties from well-test data. This paper deals only with production from a reservoir, in which case the rate is defined as being negative. The reservoir studied here is a homogeneous, disk-shaped porous body of uniform thickness, with all boundaries sealed except the inner radial boundary of the well.

By Y. Zimmels, I. J. Lin, I. Yaniv

The potential application of magnetic and electric separation techniques in beneficiation of fine-particle feeds are reviewed. Special features of these techniques are considered and compared with a view to increased amenability of such feeds to separation in electro- magnetic fields. The following techniques are discussed, with reference to those already developed and in use to-date: - Differential-susceptibility separation (DSF). - High-gradient differential-dielectric separation (HGDDS). - High-gradient magneto-gravimetric separation (HGMGS). In this context the new approach - whereby selectivity of solid- solid separation processes, in wet high-gradient magnetic and electric separation, can be improved - is described, including experimental results from DSF and HGDDS tests.

Jan 1, 1980

By John B. Pearse

The works heretofore used in carrying out the Bessemer process have been constructed substantially as follows: The whole works or plant has been divided into three parts. 1. The division in which the metals for conversion are melted. 2. The division in which the metals are converted into steel and cast into ingots, or the converting department. 3. The engine and boiler department which supplies the hydraulic power required to move the machinery, and which furnishes the blast used for converting the crude metal into steel.

H. Van F. Furman, Denver, Colo. (communication to the Secretary): There are some statements in Dr. Carpenter's interesting paper which appear to require explanation, if not modification.

Jan 1, 1901

By John B. Pearse

THE works heretofore used in carrying out the Bessemer process have been constructed substantially as follows : The whole works or plant has been divided into three parts. 1. The division in which the metals for conversion are melted. 2. The division in which the metals are converted into steel and cast into ingots, or the converting department. 3. The engine and boiler department which supplies the hydraulic power required to move the machinery, and which furnishes the blast used for converting the crude metal into steel.

Jan 1, 1876

By Walter Crafts, C. M. Offenhauer

IN the course of study of the heat-treatment of low-alloy steels, the behavior of alloy carbides at subcritical temperatures was found to vary from that indicated by published investigations. In order to establish the ranges of formation of carbides in low-chromium steels, a series of steels containing up to 0.60 per cent carbon and 7.5 per cent chromium has been investigated with respect to the type of carbide formed after quenching and tempering at sub-critical temperatures. The types of carbide formed as products of austenite transformation to pearlitic and pseudomartensitic structures have also been determined. There have been numerous investigations of the carbides in chromium steels and diagrams of the iron-chromium-carbon system.' Westgren, Phragrnén and Ne-gresco2 identified the following carbides in chromium steels: Cementitc............ Fe3C (containing up to about 15 per cent Cr) Trigonal carbide...... Cr7C3 (containing up to about 55 per cent Fe) Cubic carbide........ Cr4C (or Cr23C6) (con- . taining up to about 25 per cent Fe) Orthorhombic carbide , cr3c2 (containing only a few per cent Fe) The most complete information on the iron-chromium-carbon system was supplied by Tofaute, Sponheuer and Bennek3 and by Tofaute, Küttner and Büttinghaus,4 By means of dilatometric and microscopic stzdies, sections of the diagram were established. This diagram showed the presence of the special carbide CI7C3 in steels containing about 2.5 per cent to 12 per cent chromium in low-carbon steels and from about 5 to 20 per cent chromium in I per cent carbon alloys from A1 to room temperature. The present investigation has shown that a carbide inversion from Fe3C to Cr7C3 occurs in chromium steels containing more than about I per cent chromium after hardening and tempering at subcritical temperatures above about 500°C. The type of carbide observed to be more stable in quenched and tempered specimens was also found to be produced by direct austenite transformation in corresponding ranges of temperature and composition. Experimental Procedure The steels used in these experiments were made in magnesia-lined crucibles in high-frequency furnaces and were poured into 2-in. ingots. The ingots were forged into I-in. round bars, and all heat-treatments were conducted on specimens of this section size. Sections of the bar stock 8 in. long were hardened in air, oil or water, depending on the composition, after heating to a temperature sufficiently high to dissolve the carbides. The low-carbon and low-chromium steels were not entirely martensitic after quenching, and contained some ferrite and pseudomartensite. No pearlite was observed in the microstructures after quenching. The quenched bars were cut into

Jan 1, 1942

By Walter Crafts, C. M. Offenhauer

IN the course of study of the heat-treatment of low-alloy steels, the behavior of alloy carbides at subcritical temperatures was found to vary from that indicated by published investigations. In order to establish the ranges of formation of carbides in low-chromium steels, a series of steels containing up to 0.60 per cent carbon and 7.5 per cent chromium has been investigated with respect to the type of carbide formed after quenching and tempering at sub-critical temperatures. The types of carbide formed as products of austenite transformation to pearlitic and pseudomartensitic structures have also been determined. There have been numerous investigations of the carbides in chromium steels and diagrams of the iron-chromium-carbon system.' Westgren, Phragrnén and Ne-gresco2 identified the following carbides in chromium steels: Cementitc............ Fe3C (containing up to about 15 per cent Cr) Trigonal carbide...... Cr7C3 (containing up to about 55 per cent Fe) Cubic carbide........ Cr4C (or Cr23C6) (con- . taining up to about 25 per cent Fe) Orthorhombic carbide , cr3c2 (containing only a few per cent Fe) The most complete information on the iron-chromium-carbon system was supplied by Tofaute, Sponheuer and Bennek3 and by Tofaute, Küttner and Büttinghaus,4 By means of dilatometric and microscopic stzdies, sections of the diagram were established. This diagram showed the presence of the special carbide CI7C3 in steels containing about 2.5 per cent to 12 per cent chromium in low-carbon steels and from about 5 to 20 per cent chromium in I per cent carbon alloys from A1 to room temperature. The present investigation has shown that a carbide inversion from Fe3C to Cr7C3 occurs in chromium steels containing more than about I per cent chromium after hardening and tempering at subcritical temperatures above about 500°C. The type of carbide observed to be more stable in quenched and tempered specimens was also found to be produced by direct austenite transformation in corresponding ranges of temperature and composition. Experimental Procedure The steels used in these experiments were made in magnesia-lined crucibles in high-frequency furnaces and were poured into 2-in. ingots. The ingots were forged into I-in. round bars, and all heat-treatments were conducted on specimens of this section size. Sections of the bar stock 8 in. long were hardened in air, oil or water, depending on the composition, after heating to a temperature sufficiently high to dissolve the carbides. The low-carbon and low-chromium steels were not entirely martensitic after quenching, and contained some ferrite and pseudomartensite. No pearlite was observed in the microstructures after quenching. The quenched bars were cut into

Jan 1, 1942

By K. G. Davis, P. Fryzuk

WEINBERG,' in an investigation of solute distributions along unidirectionally solidified rods of dilute silver in tin alloys, concluded that, for rods of 2 mm diameter or larger, a very high degree of convective mixing takes place in the melt ahead of the advancing interface. Suggested causes for this mixing were (1) thermal convection, (2) convective currents set up by changes in fluid density associated with concentration gradients, and (3) currents set up by changes in volume during freezing. Evidence was presented that causes (2) and (3) were not primarily responsible for solute mixing. Further, convective mixing took place in horizontal rods with temperature gradients as low as 6°C per cm, putting doubt on the importance of thermal convection. In the present work, convective mixing is examined in melts with known initial solute distributions. The absence of a moving solid-liquid interface eliminates convection caused by volume changes on solidification. In the first experiment, three silica tubes 2 mm in diameter and approximately 20 cm long were filled with an alloy of 100 ppm Ag in tin (99.999 pct purity), placed adjacent to one another in a horizontal furnace, and solidified from one end at 1 cm per hr. A110, a Y emitter, was used as tracer. A few centimeters at the start of the rod were left un-melted, to determine the initial concentration Go. After solidification was complete, one sample was removed. The hot furnace was then brought back to

Jan 1, 1965

By Gordon R. Haworth, J. Terry Aimone

The capital-short mining industry, faced with weak metal markets, will need to devise new and novel methods to finance future expansions.

Jan 9, 1977

By R. J. Towner

Tensile properties of Aluminum Powder Metallurgy (APM) products at 600°, 800°, and 1000°F were lowered by cold reductions of as little as 10 or 20 pct. The detrimental effect of cold work on yield strength at elevated temperatures increased with temperature. Heating cold-worked products at 850° to 1150°F removed the structural effects of cold work through the process of recovery and restored, at least partially, the strength at 600°F to the level of hot-worked material. Material containing 5 pct A1203 was recrystallized by cold reductions of 97 pct or more followed by annealing at 1100°F. Recrystallization lowered sharply the strength at room temperature, but raised the strength at 600°and 800°F relative to that of the cold-worked state The high strength of Aluminum Powder Metallurgy (APM) products at elevated temperatures'-4 is well-known. In this paper, some conditions of processing that can affect this important characteristic are presented. The following nomenclature will be used to identify different alloys: All of the above alloys are of the aluminum-aluminum oxide type made from milled powders of commercially pure aluminum with the exception of M486, which is fabricated from atomized A1-Fe alloy powder. The products made from milled powders are hardened by a dispersion of aluminum oxide par- ticles in the aluminum matrix, while the intermetal-lic compound FeAL is the dispersion hardening agent in M 486.' The process used for making all APM products involves the application of heat and pressure to form a compact from powder which is fabricated into an extrusion, sheet, strip, or wire by standard techniques. Probably the most common cold-worked APM product is sheet. Cremens, Bryan, and rant' determined the room-temperature tensile properties and the 600°F stress-rupture properties of XAPOOl and XAPOO3 sheet and reported that there was little difference between the as-extruded, cold-worked, and annealed at 1200°F conditions. On the other hand, Emmons6" found that the tensile and stress-rupture properties of XAPOOl and XAPOO3 sheet at 600" and 800°F were lower than those of extrusions. The tensile properties of XAPOO3 sheet were also determined at elevated temperatures by Doyle' and ~ayward.' According to Friske" and Kendal et a1.," cold working by swaging, rolling, and tube drawing reduced the tensile and creep-rupture strength of XAP 001 Recrystallized aluminum powder metallurgy products are rare because the dispersion in the aluminum matrix is usually too fine, uniform, and concentrated to permit recrystallization. However, Lenel12 was able to recrystallize wires made from aluminum powders containing about 1.75 and 3.0 pct oxide and found that the tensile properties at room temperature remained higher than those of annealed commercially pure aluminum because of the dispersed oxide particles. Westerman and Lenel13 found the kinetics of recrystallization of the powder products were very different from those of cold-drawn commercially pure aluminum wire. The alloy containing 3.0 pct oxide recrystallized to a very coarse grain size was more creep-resistant at 400" to 600°C (752" to 1112°F) than in the as-extruded, fine-grained condition, according to Ansell and weertman.14 Ansell and aenel" further studied creep at 500" to 620°C (932" to 1148°F) of recrystallized 1.75 pct oxide wire and confirmed the creep model of Ansell and Weertman. The present work is being reported for the fulfillment of the following objectives: 1) to show the effect of cold work on tensile properties at elevated temperatures;

Jan 1, 1964

By G. Ansell, J. Weertman

G. Ansell and J. Weertman (Rensselaer Polytechnic Institute and Northwestern University, respectively) —The great increase in creep strength that you found in your indium-glass composites is quite striking. You ruled out a decrease in the dislocation source density as an explanation of the reduced creep rates. Could you tell in greater detail why you feel this is an inadequate explanation of the phenomenon. It appears to us that it is a reasonable one. Consider the spacing between the glass spheres in your composites. It is of the order of 10-3 cm. The spacing between dislocations in a (single-phase) annealed metal with a dislocation density of the order of 106 to 107 cm per cm3 also is of the order of 10-3 cm. Therefore it is not reasonable to expect that the arrangement of the dislocation lines is the same in the composite material as it is in pure indium. If the dislocation arrangements are different the density of the dislocation sources likewise will be different. Another way of saying the same thing is the following. The presence of the glass spheres reduces one of the dimensions of the metal phase to 10-3 cm. Thus the metal phase, in effect, is equivalent to a metal whisker or a thin metal film. These are strong because they have so few dislocation sources within them. T. D. Gulden and J. C. Shyne (authors' reply)—It is to be expected that the arrangement of dislocations in the indium-glass composites would be quite different than that observed in annealed pure metals. The dislocations probably exist primarily as segments connecting the dispersed particles. These segments and discontinuities at the very large area of particle-metal interface could act as dislocation sources. For this reason it was considered unlikely that the very high observed creep strengths were due to a low density of dislocation sources. A more likely explanation is that the distance a dislocation can glide is limited by the interparticle spacing. The particles averaged about twenty in diam and were less than 10 µ apart in the 40 pet composite. It seems reasonable to assume that a dislocation loop in a volume of metal surrounded by these relatively large particles would not bow out around them but would be restricted to the region between particles. A general expression for creep rate is given by the product of the dislocation source density, the average rate at which loops are generated from a source, and the average area swept out by a loop times the Burgers vector magnitude. Assuming that the area swept out is limited by the interparticle spacing and the rate of generation of loops is controlled by the rate of climb of loops over second-phase particles, as suggested by Ansel and Weert-men,9 a calculation of the dislocation source density may be made from the experimentally observed creep rates. The result comes out 10" to 10" per cm3 which seems reasonable in view of the assumptions involved. Thus it is concluded that the anomalously low creep rates that were observed were due not to a low dislocation source density but rather to a limitation on the area swept out by each dislocation loop.

Jan 1, 1964

By Rudolf Gahl

SOME time ago the Detroit Copper Mining Co. had to decide the question whether it would pay to re-treat slime-tailings, and several machines were tested in order to ascertain the type of construction which would give the greatest saving. In previous tests on the ore of this company, several tables proved far inferior to the Frue vanner, and for this reason they did not receive any serious consideration. In the late tests, however, one table was found to surpass the vanner in the percentage of saving made. The Frue vanner used in these tests had been adjusted approximately to the slope which the slime-vanuers in the concentrator of the Detroit Copper Co. had at that time, about 2 in. between the inside of the wooden posts, or 0.272 in. per ft. It had not been determined, however, if this slope was the most economical one for treating the slime-feed in the mill, or for re-treating the tailings from the slime-vanners, and I was instructed to ascertain by test-runs, under different conditions, the best adjustments for the work to. be done. In order to determine if a change in the adjustments of the machines would improve the results, it seemed sufficient to continue running the tests between the vanner and the table before mentioned, keeping the feed and all other conditions the same, but. varying the adjustments of the vanner .from run to run. I employed a vanner-man of long experience to operate these machines, who expressed the opinion that a slope between 2 3/8 and 2.5 in. between the posts (0.323 and 0.340 in. per ft.) would give the best results. Some runs made this way gave a considerable improvement. While for a slope of 2 in. (0.272 in. per ft.) the saving had been only 40.5 per cent. of that of the table, it .was 49.2 per cent. for 2.5 in. (0.340 in. per ft.) slope. A decrease of

Sep 1, 1909

By T. F. Archbold, P. G. McCormick

Silica fiber 11arkevs and hardness irzdentatiozs have been used to study the thermal diffusion of racancies in zinc placed in temperature gradients of 460o and 706oK per Cm All markers mored toward the higher temperatures, and the data show that the heat of trarlsport Q is greater than the eizthalpy of z,acancj3 forrjalion E by an amount equal to 3.5 kcal per mole. These results ore not in agreement with those obtained in either of the two previous studies with zinc. The problem is discussed in terms of experimental differences and racancy source and sink efficiencies. THE general interest in thermal-diffusion phenomena over the past several years has prompted a number of marker movement studies using pure metals. In order to understand better the relationships between thermal and mass flow, a general goal of the experiments has been to obtain values for the net enthalpies carried by the diffusing atoms. For self-diffusion by means of the vacancy mechanism, the determination of relative marker positions before and after an anneal in a temperature gradient permits a calculation of a value for (Q - Ef). Q is the transferred enthalpy, and Ef is the enthalpy of vacancy formation in the metal being investigated. The initial experimental study of thermal diffusion in pure metals was accomplished by Shewmon who found essentially no movement of surface-hardness indentation markers in single and bicrystals of zinc held for 20 days in a temperature gradient of 310°C per

Jan 1, 1967

By R. Stefanko

The use of strain gauges and compression pads for an underground bolt installation in a central Pennsylvania coal mine revealed that satisfactory anchorage could not be maintained in the shale stratum that constituted the immediate roof. Comprehensive laboratory tests resulted in a new shell designed with an increased area of contact and smaller surface serrations. Retightening of the bolt was also found to be beneficial for maintenance of high load values. Roof bolting in coal mines has progressed very rapidly in the last decade, and today this method of support is in general use. There are many theories concerning the functions of roof bolts, but the importance of beam building is indisputable. Proper tension in the bolts minimizes shearing along bedding planes and creates a thick, competent beam from individually weak layers. Bolting must be accomplished immediately upon exposure of the roof, and bolt tension must be maintained if roof failure is to be prevented. Seams are being mined today with roof bolt support where conventional timbering was not effective.' The widespread use of continuous mining equipment was made possible by roof bolts which permit long supported spans free of obstructions. Since adequate anchorage is the key to successful bolting, better methods for determining the efficiency of anchorage are continually being sought. Attempts have been made to measure bolt loads, thus assessing anchorage efficiency by employing sensing devices between the bolt head and roof. 3-s None of the techniques appears to be very accurate, and all have at least one serious deficiency. The greatest efforts are being directed toward the standardization of the pull test.6" Nevertheless, the author believes that the pull test is only one criterion for establishing the efficiency of anchorage and not necessarily the most important one. The conventional pull test involves the short-term application of forces, whereas the success of roof bolt installations is dependent upon the ability of anchorage to sustain relatively long-term forces. Therefore, it would appear that the time-dependent function of the anchor is the single most important parameter for the success of roof bolting in many types of rocks. UNDERGROUND TESTS Strain gauges and compression pads were employed to measure bolt loads in a central Pennsylvania coal mine.8 Bolting had been employed only sporadically in the shale top since dubious benefits were obtained with this type of support. Table I shows some typical results from this investigation. Anchorage losses were extremely variable among the bolts, but all showed a considerable degree of bleed-off. The tests indicated that anchorage efficiency was extremely low in this type of shale. The results confirmed that bolting under these conditions was practically worthless, and thus a waste of material and labor. Probably even more important, bolting constituted a hazard when reliance was placed upon it for support. Certainly the complete lack of confidence in roof bolting, which was apparent in the mine personnel, was justifiable since the bolts did not accomplish their purpose. PRELIMINARY LABORATORY BOLT TESTS The extreme losses in bolt tension that were detected in underground installations prompted a comprehensive study to assess the causes. However, preliminary studies were first made to determine the combined stresses in a roof bolt and to prove conclusively that the observed bleed-off phenomena were due to anchorage: slippage and not to some other physical manifestation. In all of the ensuing tests

Jan 1, 1962

By W. S. Ayres

By A. S. Emanuel

A correlation has been developed for calculating the phase behavior of, gas-condensate systems at reservoir conditions. The correlation is based on the principle of corresponding states and has been coded for an IBM 7094. Experimental K-values were determined for several gas-condensate systems at reservoir conditions to evaluate various semi-empirical parameters of the correlation. The approximate range of application of the correlation is 150 to 300F and 1,500 to 6,000 psi. INTRODUCTION The rapid development of digital computers during the past several years has made feasible the calculation of hydrocarbon phase behavior by methods based on rigorous thermodynamic principles. Good correlationsl have been developed for low to moderate pressures, but these techniques have not yet been extended successfully to reservoir fluids at high pressures. Consequently, the determination of phase behavior of oil and gas systems at reservoir conditions is still based almost entirely on generalized data correlations or on experimental studies of the fluid in question. While these methods have been used successfully many times, they do have inherent limitations that restrict their applicability. Generalized correlations, such as the NGShlA K-charts, are limited to the range of pressure, temperature and components for which the data were determined. The accuracy of these correlations is often questionable because the effect of total system composition is not well defined. Experimental studies offer a reliable method for determining phase behavior, but usually the studies are costly and time consuming. Recently, Leland and co-workers presented a new approach to calculating phase behavior from the principles of corresponding states.6-8 Corresponding states methods determine the thermodynamics properties of a given system by comparison with a reference substance whose properties are known. The accuracy of this approach depends on close chemical and structural similarity between the reference substance and the system in question and between components within the system itself. For high accuracy, it is usually necessary to correct for chemical and structural dissimilarities. In principle, however, the corresponding states method should be no less accurate at high pressures than at low pressures, provided reference substance properties are known. This paper describes an empirical modification of the basic correlation proposed by Leland, et al.6 for the specific purpose of calculating the phase behavior of gas-condensate fluids at reservoir conditions. The modified correlation, which has been programmed for an IBM 7094, may be used for either approximate or precise determination of fluid behavior depending on the amount of analytical and, experimental data available for the system. BASIC THEORY The basic theory of the corresponding states phase equilibria correlation was first published by Leland, Chappelear and Gamson.6 Subsequently, Leland, Chappelear and Leach published methods for improving the accuracy of the original theory.7-9 The aim of the correlation is to calculate the K-value of each component of a given system as a function of pressure, temperature, and over-all composition, where Once the K-values are known, the phase behavior may be determined directly by an appropriate flash calculation. The basic equation for calculating component K-values was taken from the work of Joffe. 4 For any component i of a mixture, the K-value is given by

Jan 1, 1969

By R. M. Asimow

Both smith1 and wells2 have shown that the dif-fusivity of carbon in austenite increases rapidly with increasing carbon concentration. Smith has shown, in addition, that part but not all of this increase can be attributed to the concentration dependence of the activity coefficient of carbon in austenite. From the work of Eshelby', it is known that, because of the image displacements, there is a uniform dilatation of the material in a body containing random centers of dilatation. It is thus tempting to explain the enhanced diffusivity as resulting from this apparent negative pressure caused by the interstitial carbon atoms.* That is, the activation free energy for diffusion is not directly a function of concentration if we consider that the alloy is subjected to a concentration-dependent pressure. To proceed with a quantitative solution, it is necessary to determine the activation volume AV for the diffusion process. An exact calculation of AV is quite a difficult problem and does not appear to have been solved as yet. A simpler procedure is to use the semiempirical, macroscopic theory of seyes' which relates the activation free energy AG to AV. Keyes has shown that for many types of diffusion processes the following relationship is valid Where p is the Gruneisen constant and x the compressibility. This relationship is based on the assumption that AG can be identified with strain-energy terms. In general there will also be contributions to AG due to a change in the energy of the valence electrons of the metal; thus the relation must be considered largely empirical. Surprisingly, for self-diffusion where it is well-known that the electronic contribution is quite important,' the Keyes relation works quite well for many systems. For interstitial diffusion where it might be expected that the relationship is on a firmer theoretical basis, the limited experimental data available is not in particularly good agreement as can be seen from Table I. However, it must be remembered that besides the uncertain basis of the Keyes formula there is considerable uncertainty in some of the experimental values of AV. The only real justification for using the Keyes formula in the present analysis is that it gives a result which satisfactorily solves the problem at hand and, further, there are no more reliable procedures for calculating AV available. The coefficient 3 in Eq. [I] is determined such that the room-temperature value for x is to be used. The activation entropy AS is estimated as 8.95 x 10"3 kcal per OK mole using the values in Table I1 and the Zener" relation AS = PAHIT, where p is a dimensionless thermal coefficient of Young's modulus which is to be evaluated at low values of T/T,, aH is the activation enthalpy, and T, is the melting temperature of the alloy. The resulting value of AG at 1000°C is 25.6 kcal per mole. The Gruneisen constant (p = 3cup/XC,), where p is the density, is evaluated from the room-temperature data in Table I1 after calculating C, from the relation. A value of p = 8.32 g per cu cm based on the value of a. in Table I is used. The result is C, = 5.0 x lo6 erg per g, H = 2.32, and AV = 3.41 cu cm per mole at 1000°C and infinite carbon dilution. It is assumed that the activation volume is not a function of carbon concentration. The theory could be refined to consider AV a function of concentration, but due to the uncertainties in the parameters there seems little justification for this. Carbon dissolved in austenite forms a nonideal solid solution, so that in writing an equation for. the diffusivity of carbon it is necessary to take into account the change in activity coefficient with concentration. The resulting equation is16

Jan 1, 1964

By Ralph H. Sweetser

OF SPECIAL importance in the design and construction of an iron blast-furnace plant are tile raw materials to be employed. Obviously the iron must come from some ore of that metal, but the many kinds available lead to important differences in plant arrangement, and in the design of furnaces and accessory equipment. For instance, the bulk of the gangue in the ore is subject to wide variation as is that of the ash in the fuel. A gradual increase in the bulk of raw material necessary to produce a ton of pig iron has been noted, but this has been more than offset, so far as costs of operation are concerned, by improvements in the machinery used for handling large quantities of ore. coal limestone, and slag. The mechanization of coal and iron-ore mining. the building of big ore boats and limestone carriers, of ore docks for fast loading and unloading, and the general use of heavy railroad equipment have all contributed to greatly reducing the number of man- hours required to mine, transport, and manufacture the raw materials into a ton of pig iron. This looks like a big saving and it is. If it were not for these big labor saving machines it would be impossible to produce the tonnage of iron and steel that is made each year in this country. On the other hand, the total capital investment necessary to bring down the man-hours per ton of pig iron to a 2.3 hours average at the blast furnace and to a total of only 12 man-hours from ore to metal (see Table 1) is so great that even the present high price of pig iron is probably not far above the true cost of manufacture when all the fixed charges are properly included.

Jan 1, 1939

By H. S. Price, J. C. Cavendish, R. S. Varga

This paper presents a new technique for solving some of the partial differential equations that are commonly used in simulating reservoir performance. The results of applying this technique to a simple problem show that one obtains accurate pressure values near wells, as well as accurate pressure pdients, which can be explicitly calculated. The method is completely rigorous in that convergence of the discrete numerical solution to the continuous solution for both pressure and pressure gradient is established. High-order, piecewise-polynomial approximations are used near the wells where pressure gradients are steep, while low-order, piecewise-polynomial approximations are used elsewhere to reduce greatly the calculation time. This combination is shown to give a uniformly good approximation to the solution. These approximations, obtained by using a Galerkin process with suitable Hermite subspaces, are shown to be theoretically and numerically superior to the usual approximations obtained from standard finite-difference techniques. Not only are much greater accuracies obtained, but computer times are also greatly reduced. The application of this technique to multiphase flow problems (e.g., single well coning problems) would have considerable practical interest, but such extensions of this technique with full mathematical rigor have not been made as yet. However, the numerical methods presented here are general, and in principle extend to multidimensional, multiphase flow. Moreover, the preliminary results given in this paper are sufficiently encouraging that we feel the effort in attempting these extensions is justified. INTRODUCTION The problem of obtaining accurate pressure distributions and pressure gradients around wells is of considerable importance in the numerical simulation of reservoir performance. The most common approach to solving this problem is to use finite difference techniques (see McCarty and Barfield9 or Peaceman and Rachfordl0). This approach, however, has many disadvantages, the major one being that many grid points are generally necessary for accurately describing the pressure distribution and the pressure gradient around wells. This need for a fine grid results in large computer times and often in prohibitively high costs. Besides investigating the method of finite differences, some authors, such as Welge and weber17 and Roper, Merchant and Duvall,13 have considered a combination of analytical and numerical techniques with some success. These approaches, however, are all nonrigorous and quite often cannot even be applied. In this paper, we present a numerical formulation of high-order accuracy, based on the Galerkin method, for solving this problem. We treat here only the partial differential equation that describes steady-state, single-phase flow. However, the methods presented are general and in principle extend to multidimensional, multiphase flow. Specifically we treat special cases of the problem in two dimensions described by: where G is a rectangle (with sides parallel to the coordinate axis) with boundary ?G, ?/?n denotes the outward normal, and y and B are non-negative constants such that y + B > 0.

Jan 1, 1970