Search Documents

By John Chipman, Raymond L. Orr

Recently reported high -temperatlnre thermal data have been incorporated in a reuieu' of the therinody-r~ar,zic propesties of iron. The recent data permit mot corlsisterlt and reliuble choices to be made for many of the pvoperties, resnlling in bettev eslablished tabrtlations of the thermodynamic functions. Values of Cp, h°t- - H°298. s°t - S2sa, and (G°T -H298/T are tabulated for -the solid and liquid phases of iron in both stable and melastable regions. Consistency of the tabulations has been maintained to the precision necessr~vg to yield the Gibbs energy change belween the bcc (a, 8) and fcc (y) phases of iron at temperatures pertinenl to alloy studies. THE importance of the metal iron has stimulated many investigations of its thermodynamic properties at elevated temperatures. Reviews and summaries of the data have appeared from time to time. Among these the 1951 paper of Darken and smith' deserves highest recognition for the quality of its criticism and for the precision and thermodynamic consistency of its tabulations. These and the more recent table of Kaufman, Clougherty, and weiss2 reproduce the Gibbs energy of the (a, 6)--Y transformation with sufficient sensitivity for alloy studies. Other recent compilations of equal validity3-' make no attempt to do this. Since this Gibbs energy change is not greater than 22 cal per g-atom over the temperature range 1100° to 1800°K, it is obviously desirable to tabulate values of (G°t - H298)/T to four places of decimals even though the absolute accuracy would not warrant such precision. The necessity for reconsideration of the problem is emphasized by the recent publication of several important researches on the enthalpy and heat capacity of iron. The older work is well-covered in the reviews cited and will be mentioned here only when actually employed in the tabulations. GAMMA IRON Transition temperatures are taken from the selec -tions of Elliott and Gleiser,5 based on the work of Boulanger7 and earlier workers. For the a-y (A3) transition at 1184°K we retain the value 215 cal per g-atom adopted by Darken and smith.' It has received confirmation in the recent work of Dench and Kubaschewski,' and of Braun8 using adiabatic calori-metry, and of Wallace, Sidles, and Danielson10 using a pulse heating technique. The enthalpy data of Olette and ferrier" and of Anderson and Hultgren12 at higher and lower temperatures may be extended to 1184°K to give the following values which are adjusted to the required value of ? 1184 Ha,1184 Ha 298 = 8030 cal per g-atom [11 hr,1184-ha298 = 8245 cal per g-atom [2]

Jan 1, 1968

By B. M. Larsen, L. S. Darken

SOME years ago we collated all laboratory data then available to us on the distribution at equilibrium of manganese and of sulphur between metal and simple slags, and used the results in setting up an equation, which proved to be in accord with determinations of the distribution ratio of manganese and of sulphur between liquid steel and open-hearth slag, over a wide range of operating conditions. We now present an outline of our calculations and inferences on these two closely related questions, which are of especial importance because of the advisability of conserving manganese wherever possible. In general, our results show that in the open-hearth furnace equilibrium between metal and slag, with respect to both manganese and sulphur, is substantially attained during the finishing period of the heat, or in about a half hour to an hour after any additions were made. This implies that the end result will be identical whether manganese is added as alloy or metal to the liquid steel or as oxide to the slag. The proportion of total manganese retained by the steel is favored: (I) by use of the minimum amount of lime, consistent with other operating limitations, required to give a lime-silica ratio of about 2.4, which usually corresponds to a minimum practical slag volume; (2) by a high concentration in the slag of MnO relative to FeO; (3) by high operating temperature. Absorption of sulphur by the slag from the liquid metal is favored: (I) by large slag volume; (2 by a high concentration of active (or "free") lime and of manganese oxide dis) solved in the slag; (3) by low concentration of silica, iron oxide, and phosphate dissolved in the slag. DISTRIBUTION OF MANGANESE BETWEEN SLAG AND METAL AT EQUILIBRIUM On the basis that the distribution of manganese between liquid steel and slag is controlled by the reaction [FeO + Mn = MnO + Fe [A]slag metal slag metal] the distribution at equilibrium, at a given temperature, is expressed by the equilibrium constant for this reaction, which is defined as [K" aMoO aFeI1 dFec) aMn] where a is the activity of the molecular species denoted by the subscript. The activity of each of the several components, although it is not known precisely, is approximately equal to the concentration of each component expressed as its mol fraction; that is, as the ratio of the number of mols of the given component to the total number of mols of all the molecular species present in the metal or slag, as the case may be. Moreover, since the activity, or concentration, of iron in the liquid metal is in fact substantially constant, the term [afe] may be combined with K" to give a new constant K'. Making these changes,

Jan 1, 1942

By R. F. Mehl, B. Lustman

THE formation of intermediate phase layers in cementation processes has been subjected to extensive qualitative investigation though to relatively little quantitative study; this work has recently been fully summarized by F. N. Rhines.1 It is necessary here only to enumerate the general conclusions drawn concerning the kinetic laws governing the growth of intermediate layers and to state what is known concerning the structural changes accompanying the generation of phase layers. I. When the intermediate phase is protective in type-i.e., when its density is less than that of the base metal, so that a continuous coating is formed and maintained the isothermal growth of the phase follows a simple parabolic equation: X2 = kt [I] where X is the thickness of the alloy layer, k the reaction rate constant, and t the time. The applicability of this equation to diffusion processes is tested by plotting the square of the thickness versus time, or by plotting the thickness versus the square root of time; the resultant achievement of a straight line in such. a plot has been con¬sidered sufficient proof of the validity of this equation. Exceptions have been found where liquid phases are involved in the reaction, where the diffusion layer is cracked,2 and where thin semimetallic oxide and sulphide coatings are formed on metals.3 In the first two cases the thickness increases linearly with time, and in the latter it is a logarithmic function of time. 2. The rate of thickening varies with temperature according to the Arrhenius equation k = Ae-Q/RT [2] where k is the parabolic rate constant of Eq. I, or some other rate-descriptive function, A is the action constant, e the base of natural logarithms, Q the heat of activation, R the gas constant, and T absolute temperature. This equation has been found to be generally valid except when transformations occur during reaction in either the base metal or in the reaction layer; in such instances Eq. 2 is still applicable with different values for the constants A and Q above and below the transformation temperature. 3. All phases that are stable at the temperature at which diffusion occurs will appear in the interdiffusion of two components. In some cases certain of the phases may be too thin to be observed and identified easily. Nonequilibrium phases have never been noted in thick reaction layers; i.e., layers of a thickness that may be discerned microscopically. The relative, thickness of each phase in a diffusion layer is proportional to the concentration gradient and to the rate of diffusion in that phase .4 The reaction layers tend to have a

Jan 1, 1942

By J. E. Chappelear

We are concerned with the problem of calculating the forces produced by passing waves on a submerged object (in particular a drilling barge) located on the ocean floor. We employ the linearized theory of waves (the Airy1 theory), since solutions are generally more easily obtained in such a theory. It will be indicated how an approximate solution to the problem can be obtained for waves of finite height by use of a linear combination of Airy waves. We shall assume the submerged object has an infinite length, thus neglecting diffraction and making the problem two-dimensional. We find that there is not a considerable difference between our solution and the usual engineering solution. The force, as a function of time, is plotted in Figs. 4 and 5, together with the usual solution for easy comparison. So far, no experiments have been found to compare with the results. As will be shown, the calculations can be extended to waves of finite height (Stokes2-Struik3 waves) in an approximate fashion. This can be done by calculating the force produced by each Fourier component of the wave, separately, and omitting the coupling between the various harmonics produced by Bernoulli's theorem. The problem and its boundary conditions will be briefly stated. We seek a function, +, the velocity potential, which is harmonic, ?2?=0...........(1) and which satisfies the boundary conditions ??/?n =0 .......(2) (the normal derivative of + is zero) on the ocean floor and the sides of the barge, and s??/?y = ?..........(3) on the ocean surface, where s = g/w2........(4) The physical reasons for these choices and an interpretation of ? are given by Lamb.' As previously stated, we assume also that the problem is two-dimensional. We set up our coordinates such that x = 0 is the center of the barge and y = 0 is the ocean floor. The coordinates are shown in Fig. 1, together with the location of the barge on the bottom. Our solution may be decomposed into a function which is symmetric about x = 0,?e, and one which is antisymmetric, +,, ? = ? + ?o...........(5) Thus directly above the barge we have as boundary conditions that for x = 0, ?o = 0, ...........(6) and ??e/?x = 0........(7) Now, far from the barge, we have that (for x positive and assuming the waves are incident from the left) ?e ? Ae cosh ky cos k (x + a) , . . . (8) and ?e ? Ao cosh ky cos k (x + ß) . . . . (9) where + means ". .. approaches ... for x large." Let us put Ae, Ao, and k equal to one, and remember that, in accordance with the derivation of Eq. 3, a sinusoidal time factor is to be supplied. Then we may rewrite Eq. 8 for ?e, and we find that

Jan 1, 1958

By B. M. Larsen, L. S. Darken

Some years ago we collated all laboratory data then available to us on the distribution at equilibrium of manganese and of sulphur between metal and simple slags, and used the results in setting up an equation, which proved to be in accord with determinations of the distribution ratio of manganese and of sulphur between liquid steel and open-hearth slag, over a wide range of operating conditions. We now present an outline of our calculations and inferences on these two closely related questions, which are of especial importance because of the advisability of conserving manganese wherever possible. In general, our results show that in the open-hearth furnace equilibrium between metal and slag, with respect to both manganese and sulphur, is substantially attained during the finishing period of the heat, or in about a half hour to an hour after any additions were made. This implies that the end result will be identical whether manganese is added as alloy or metal to the liquid steel or as oxide to the slag. The proportion of total manganese retained by the steel is favored: (I) by use of the minimum amount of lime, consistent with other operating limitations, required to give a lime-silica ratio of about 2.4, which usually corresponds to a minimum practical slag volume; (2) by a high concentration in the slag of MnO relative to FeO; (3) by high operating temperature. Absorption of sulphur by the slag from the liquid metal is favored: (I) by large slag volume; (2) by a high concentration of active (or "free") lime and of manganese oxide dissolved in the slag; (3) by low concentration of silica, iron oxide, and phosphate dissolved in the slag. Distribution 02 Manganese between Slag and Metal at Equilibrium On the basis that the distribution of manganese between liquid steel and slag is controlled by the reaction FeO + Mn -& MnO + Fe [A] slag metal slag metal the distribution at equilibrium, at a given temperature, is expressed by the equilibrium constant for this reaction, which is defined as K,, _ aMnO' afe aFeO ' aMn where a is the activity of the molecular species denoted by the subscript. The activity of each of the several components, although it is not known precisely, is approximately equal to the concentration of each component expressed as its mol fraction; that is, as the ratio of the number of mols of the given component to the total number of mols of all the molecular species present in the metal or slag, as the case may be. Moreover, since the activity, or concentration, of iron in the liquid metal is in fact substantially constant, the term afe may be combined with K" to give a new constant K'. Making these changes, 7

Jan 1, 1942

By B. M. Larsen, L. S. Darken

Some years ago we collated all laboratory data then available to us on the distribution at equilibrium of manganese and of sulphur between metal and simple slags, and used the results in setting up an equation, which proved to be in accord with determinations of the distribution ratio of manganese and of sulphur between liquid steel and open-hearth slag, over a wide range of operating conditions. We now present an outline of our calculations and inferences on these two closely related questions, which are of especial importance because of the advisability of conserving manganese wherever possible. In general, our results show that in the open-hearth furnace equilibrium between metal and slag, with respect to both manganese and sulphur, is substantially attained during the finishing period of the heat, or in about a half hour to an hour after any additions were made. This implies that the end result will be identical whether manganese is added as alloy or metal to the liquid steel or as oxide to the slag. The proportion of total manganese retained by the steel is favored: (I) by use of the minimum amount of lime, consistent with other operating limitations, required to give a lime-silica ratio of about 2.4, which usually corresponds to a minimum practical slag volume; (2) by a high concentration in the slag of MnO relative to FeO; (3) by high operating temperature. Absorption of sulphur by the slag from the liquid metal is favored: (I) by large slag volume; (2) by a high concentration of active (or "free") lime and of manganese oxide dissolved in the slag; (3) by low concentration of silica, iron oxide, and phosphate dissolved in the slag. Distribution 02 Manganese between Slag and Metal at Equilibrium On the basis that the distribution of manganese between liquid steel and slag is controlled by the reaction FeO + Mn -& MnO + Fe [A] slag metal slag metal the distribution at equilibrium, at a given temperature, is expressed by the equilibrium constant for this reaction, which is defined as K,, _ aMnO' afe aFeO ' aMn where a is the activity of the molecular species denoted by the subscript. The activity of each of the several components, although it is not known precisely, is approximately equal to the concentration of each component expressed as its mol fraction; that is, as the ratio of the number of mols of the given component to the total number of mols of all the molecular species present in the metal or slag, as the case may be. Moreover, since the activity, or concentration, of iron in the liquid metal is in fact substantially constant, the term afe may be combined with K" to give a new constant K'. Making these changes, 7

Jan 1, 1942

By W. J. Romanow

RECENT papers by Williams, Sherwood, and Boothby,1 Seybolt et al.2 Heikes; and Roberts4 attest to the increasing importance of the hexagonal ferromagnetic compound MnBi. Since the magnetic properties of MnBi crystals are highly anisotropic, it is useful and necessary to have appropriate aids for determining their orientations—a (00.1) standard stereographic projection and a corresponding table of crystallographic angles for this structure. The tables of crystallographic angles for the hexagonal system available in the literature5 are computed for c/a ratios greater than that of MnBi and cannot satisfactorily be used to find the crystallographic angles of MnBi. A table of angles between the crystallographic planes (HK.L) and (hk.1) was calculated for MnBi from the formula0 F=are cos Hh+Kk+ 1/2(Hk+hK) + 3/4 (a/c)2Ll/v(H2 + K2 + HK + 3/4 (a/c)2L2) (h2+ k2+hk + 3/4 (a/c)2l2) on an IBM Type 650 Magnetic Drum Data Processing machine. The lattice constants used are listed in Table I and are those determined by Willis and Rooksby7 at 20°C. More recently published values obtained by both X-ray and neutron diffraction"' show no difference from these significant to the calculation. From the crystallographic angles between the basal plane (00.1) and the general plane (hk.1) in Table I, an (00.1) standard stereographic projection was made for MnBi which is shown in Fig. 1. Acknowledgments The author wishes to thank R. G. Treuting and J. N. Hobstetter for their technical assistance, and Miss D. P. Booth for programing the IBM machine for the calculations. He appreciates the opportunity to read the paper of B. W. Roberts of General Electric Co. prior to its publication. References 1H. J. Williams, R. C. Sherwood, and 0. L. Boothby: Bulletin Amer. Physical Soc., 1956 vol. 1, p. 132. 2 A. U. Seybolt, H. Hangen, B. W. Roberts, and P. Yurcisin: AIME Trans., 1956, vol. 206, p. 605: Journal oF Metals, May 1956. 3R. R. Heikes: Physical Review, 1955, vol. 99, p. 446. .B. W. Roberts: A Neutron Diffraction Study of the Structures of Manganese B.W.Roberts:ANeutronDiffractionBismuthide, Physical. Review, 1956, vol. 104, p. 607. 5 A. Taylor and S. Leber: AIME Trans., 1954, vol. 200, p. 190; JouRNal of Metals February 1954. « C. S. Barrett: Structure of Metals, 2nd Ed., p. 635. McGraw-Hill Book Co. New York and London, 1952. 7B. T. M. Willis and H. P. Rooksby: Proceedings Physical Soc., London, 1954, vol. 67B, p. 290.

Jan 1, 1958

By W. Rosotoker

Using an expression for the free energy of a homogeneous phase as a function of composition, a relationship is derived which interrelates the phase boundaries extending from the allotropic transformation of the solvent metal. Comparisons between observed and calculated phase boundaries in the systems Fe-Ni, Fe-Mn, Fe-W, Ti-Mo, and Ti-Cb are presented. IT is the purpose of this paper to illustrate how the general free energy equations for two phases in equilibrium may be used to derive an expression which interrelates the phase boundaries extending from the allotropic transformation of the solvent metal. The free energy per gram atom of a disordered solid solution may, for many purposes, be represented by the equation:' F = ½NZ [C-V,, + (1-C) VbB + 2C (l-C)V] + K(T) +RT[C-lnC+ (l-C)ln(l-C) ] [1] where N is Avogadro's number; Z, the coordination number; C, the concentration of the solute element in atomic percent; the subscripts refer to concentrations at phase mixture boundaries; VAA, VBB, and V,, are the bond energies between like and unlike adjacent atoms, respectively, at the interatomic spacing characteristic of the phase. The bond energies are usually assumed to be independent of small changes in the lattice parameter occasioned by solute additions and variations in temperature; V = VAB ---------------—— ; K(T) is the sum of the heat con- tent and heat entropy; R, the molar gas constant; and T, the temperature in degrees Kelvin. The free energy of a homogeneous phase is there- fore the sum of an internal energy, a specific heat, and a configurational entropy term, respectively. The compositions of two phases in equilibrium at a given temperature are governed by the minimum total free energy principle. Graphically, this is represented by the points of tangency of the tangent common to the free energy vs. concentration curves of the two phases. Analytically, this may be expressed as: ( df) = (df) [2] dc dc Expanding both sides using eq 1, the following equation is obtained after convenient readjustment of the terms: RT In[(ca)(1-c)]-½N (Vaa-Vbb)(Za-Zb) +NV [Za(l-2Ca)-Z,(l-2C,)] [3] The following assumptions have been made at this point: 1—The specific heats of the two phases are independent of composition variations within the range under consideration. Accordingly, the K(T) terms disappear on differentiation. 2— (Vaa)o = (VAA); (VBB). = (VBB); Va = V, Theyjustification for these assumptions can lie only in the successful application of the resultant equation. It may be noted that Zener2 has used assumption 1 in the calculation of the boundaries limiting the (a+j3) phase mixture in the Cu-Zn system.

Jan 1, 1952

By O. L. Culberson, J. J. McKetta

INTRODUCTION The equilibrium constants for methane and for water, and for ethane and water have been calculated from experimental data for the two binary systems.2,3,11,12 These constants are for the two-phase systems for the temperature range of 100" to 340°F and the pressure range of 200 to 10,000 psia. The constants for ethane are greater than those for methane. Equilibrium constants for a natural gas from a natural gas-water system6 are about the same as those for methane. Equilibrium constants for water in all three systems are very nearly the same over the greater part of the range of temperatures and pressures studied. In ideal solutions, the vapor-liquid equilibrium constants for a component should be the same in all such solutions at a given temperature and pressure.5 Few solutions are ideal, and consequently the use of generalized constants often leads to error in situations where they do not actually apply. In decidedly non-ideal systems, it is desirable to have equilibrium constant data for the specific system. Vaporization equilibrium constants may be calculated either by thermodynamic methods or from experimental vapor-liquid equilibrium data. In the thermodynamic calculation. the equilibrium constant K may be shown to be" K = fP,/fl'..........(1) where fPi is the fugacity of the pure component i at the given temperature and at the vapor pressure of the pure component, and fp is the fugacity of pure i in the gaseous state at the temperature and total pressure of the system. At best the thermodynamic calculation results in a K applicable to ideal solutions. Equilibrium constants calculated from Equation (1) have been found to have only limited agreement with experimental values, even though the solutions for which they were calculated approached ideality."'" If the more volatile component is at a temperature above its critical, the accuracy of the equilibrium constant for the component is made more doubtful by the use of the vapor pressure term under such conditions. In order to make the thermodynamic calculation on a component above its critical temperature, one must resort to one of the empirical methods of extrapolation of the vapor pressure to conditions above its critical point.13 Application of the thermodynamic equilibrium constant then to the methane-water or the ethane-water systems would be risky since the systems are non-ideal, and the methane and ethane are above their critical in most petroleum operations. ETHANE-WATER SYSTEM In the temperature range of from 100° to 340°F and for pressures up to 10,000 psia, the vapor-liquid equilibrium data 'References given at end of paper. Manuscript received in the office of the Petroleum Branch March 16, 1951.

Jan 1, 1951

By R. B. Gilmore, Lewis W. MacNauhgton

The North Texas district, as' herein defined, includes the counties of Archer, Baylor, Clay, Cooke, Foard, Hardeman, Knox, Montague, Wichita, and Wilbarger. This area covers generally the crest and south flank of a system of buried mountains known as the Red River uplift. The oil and gas accumulations along this feature are in traps, which, although localized by structures incident to the regional uplift, usually are modified by stratigraphic changes in the sediments. Excepting those in southeastern Baylor, southern Archer, and southwestern Clay counties, all the fields within the district are on this Red River uplift. These exceptions, which have the same type of oil and gas accumulation as the other fields, are on the extreme north end of the Bend arch, which is a broad anticline plunging northward from the Llano uplift in central Texas to Archer County. The larger part of past oil and gas production has come from Pennsyl-vanian strata, with less important amounts from the Permian, and minor but increasingly important quantities from the Ordovician. Developments during 1942 Approximately 500 wells were drilled in the district during 1942, in contrast with some 1600 wells drilled during 1941, and only some goo oil wells were completed, a decrease from 1941 of more than 10 per cent in the number of successful drilling operations of the total wells drilled. Production increased slightly over the previous year in spite of a decrease in the number of producing wells. The average depth of discovery wells was some 600 ft. deeper than in 1941. Although more discoveries produced from horizons above the base of the Strawn series than those below the Strawn, the deeper horizons, including the Bend series of the Pennsylvanian, the Mississippian limestone, and the Simpson and Ellenburger of the Ordovician, furnished most of the addition to oil reserves. The largest single addition to the reserves was again in the Ellenburger formation in the K. M. A. field. During the year 47 wells were completed, bringing the total developed acreage in this horizon to some 5450 acres. The most active area was Montague County, in which eight fields, or potential fields, were discovered although only 39 oil wells were completed. K. M. A: Pressure Maintenance Satisfactory progress was made by the K. M. A. Pressure Maintenance Association in 1942, At the end of the year, 100 wells were being used for input and total monthly volume of gas returned to the reservoir reached a high of 430,971 M cu. ft. in December, with an input ratio of 587 cu. ft. of gas returned for each barrel of oil produced during the month. This represents a utilization of 50 per cent of the gas available for return during that month. Comparison of the present average bottom-hole pressure by means of a decline

Jan 1, 1943

By Lester Guttman

The connection between short-range order and thermodynamic activities in binary solid solutions has been pointed out by Birchenall1 who calculated approximate values of the energy of interaction of Zn and Cu and of the degree of short-range order in alpha-brass at one composition and two temperatures. The purpose of this note is to call attention to a somewhat more accurate theory than that used by Birchenall and to extend the calculations to a wider range of temperatures and compositions as a means of testing its applicability. This is an example of the well-known general method of comparing measured thermodynamic functions of a system with those calculated from a microscopic model by means of statistical mechanics. In this case, mathematical approximations must be used which make difficult any decision about the validity of the model. Birchenall's statistical method resembles closely that calletl by Fowler and Guggenheim2 the "quasi-chemical" method which has been shown to lead in certain cases to results identical with those obtained by the method of Bethe.3 The quasi-chemical method has been extended by Takagi4 to give general formulas for the free energy of binary solid solutions based on several lattices at all compositions and temperatures. Specifically, he treats the four systems typified by: I. ß—brass; 2. a— brass; 3. AuCus; 4. AuCu. When the temperature is so high that no long-range order exists, these formulas reduce to the same form for all lattice types, namely: F - FO = (E - EO) - T(S - S°) = ½N(c — x)(2vAB — vAa - Vbb) + M1/2KTzN[x log x + 2(c - x) log (c - X) + (I - 26 + X) log (1 - 2c + x)] -kT(z - I)N[c log c+(I-c) log (I -c)] Eq 1 Here the symbols have the following meanings: F—Helmholtz free energy E—nergy 5—ntropy k—the gas constant per molecule z—coordination number, (e.g. I2 in a-brass). N—total number of atoms c—atomic fraction of A Vaa, vbb, vab—nergy per bond of the corresponding bonds log means log, throughout The superscript ° denotes the thermodynamic fuliction in the standard state, which has been chosen to be pure A or B with the same crystal structure as the solid solution. The quantity x is the ratio of the number of A-A bonds to the total number and for minimum F satisfies the equilibrium condition Vaa + Vbb — 2Vab _ (c — x)2 kT l0g x(i - 2C + x The free energy may also be written in terms of /La and B, the chemical potentials or partial free energies per alom:

Jan 1, 1949

By Lester Guttman

The connection between short-range order and thermodynamic activities in binary solid solutions has been pointed out by Birchenall1 who calculated approximate values of the energy of interaction of Zn and Cu and of the degree of short-range order in alpha-brass at one composition and two temperatures. The purpose of this note is to call attention to a somewhat more accurate theory than that used by Birchenall and to extend the calculations to a wider range of temperatures and compositions as a means of testing its applicability. This is an example of the well-known general method of comparing measured thermodynamic functions of a system with those calculated from a microscopic model by means of statistical mechanics. In this case, mathematical approximations must be used which make difficult any decision about the validity of the model. Birchenall's statistical method resembles closely that calletl by Fowler and Guggenheim2 the "quasi-chemical" method which has been shown to lead in certain cases to results identical with those obtained by the method of Bethe.3 The quasi-chemical method has been extended by Takagi4 to give general formulas for the free energy of binary solid solutions based on several lattices at all compositions and temperatures. Specifically, he treats the four systems typified by: I. ß—brass; 2. a— brass; 3. AuCus; 4. AuCu. When the temperature is so high that no long-range order exists, these formulas reduce to the same form for all lattice types, namely: F - FO = (E - EO) - T(S - S°) = ½N(c — x)(2vAB — vAa - Vbb) + M1/2KTzN[x log x + 2(c - x) log (c - X) + (I - 26 + X) log (1 - 2c + x)] -kT(z - I)N[c log c+(I-c) log (I -c)] Eq 1 Here the symbols have the following meanings: F—Helmholtz free energy E—nergy 5—ntropy k—the gas constant per molecule z—coordination number, (e.g. I2 in a-brass). N—total number of atoms c—atomic fraction of A Vaa, vbb, vab—nergy per bond of the corresponding bonds log means log, throughout The superscript ° denotes the thermodynamic fuliction in the standard state, which has been chosen to be pure A or B with the same crystal structure as the solid solution. The quantity x is the ratio of the number of A-A bonds to the total number and for minimum F satisfies the equilibrium condition Vaa + Vbb — 2Vab _ (c — x)2 kT l0g x(i - 2C + x The free energy may also be written in terms of /La and B, the chemical potentials or partial free energies per alom:

Jan 1, 1949

By J. J. McKetta, O. L. Culberson

INTRODUCTION The equilibrium constants for methane and for water, and for ethane and water have been calculated from experimental data for the two binary systems.2,3,11,12 These constants are for the two-phase systems for the temperature range of 100" to 340°F and the pressure range of 200 to 10,000 psia. The constants for ethane are greater than those for methane. Equilibrium constants for a natural gas from a natural gas-water system6 are about the same as those for methane. Equilibrium constants for water in all three systems are very nearly the same over the greater part of the range of temperatures and pressures studied. In ideal solutions, the vapor-liquid equilibrium constants for a component should be the same in all such solutions at a given temperature and pressure.5 Few solutions are ideal, and consequently the use of generalized constants often leads to error in situations where they do not actually apply. In decidedly non-ideal systems, it is desirable to have equilibrium constant data for the specific system. Vaporization equilibrium constants may be calculated either by thermodynamic methods or from experimental vapor-liquid equilibrium data. In the thermodynamic calculation. the equilibrium constant K may be shown to be" K = fP,/fl'..........(1) where fPi is the fugacity of the pure component i at the given temperature and at the vapor pressure of the pure component, and fp is the fugacity of pure i in the gaseous state at the temperature and total pressure of the system. At best the thermodynamic calculation results in a K applicable to ideal solutions. Equilibrium constants calculated from Equation (1) have been found to have only limited agreement with experimental values, even though the solutions for which they were calculated approached ideality."'" If the more volatile component is at a temperature above its critical, the accuracy of the equilibrium constant for the component is made more doubtful by the use of the vapor pressure term under such conditions. In order to make the thermodynamic calculation on a component above its critical temperature, one must resort to one of the empirical methods of extrapolation of the vapor pressure to conditions above its critical point.13 Application of the thermodynamic equilibrium constant then to the methane-water or the ethane-water systems would be risky since the systems are non-ideal, and the methane and ethane are above their critical in most petroleum operations. ETHANE-WATER SYSTEM In the temperature range of from 100° to 340°F and for pressures up to 10,000 psia, the vapor-liquid equilibrium data 'References given at end of paper. Manuscript received in the office of the Petroleum Branch March 16, 1951.

Jan 1, 1951

By R. B. Gilmore, Lewis W. MacNauhgton

The North Texas district, as' herein defined, includes the counties of Archer, Baylor, Clay, Cooke, Foard, Hardeman, Knox, Montague, Wichita, and Wilbarger. This area covers generally the crest and south flank of a system of buried mountains known as the Red River uplift. The oil and gas accumulations along this feature are in traps, which, although localized by structures incident to the regional uplift, usually are modified by stratigraphic changes in the sediments. Excepting those in southeastern Baylor, southern Archer, and southwestern Clay counties, all the fields within the district are on this Red River uplift. These exceptions, which have the same type of oil and gas accumulation as the other fields, are on the extreme north end of the Bend arch, which is a broad anticline plunging northward from the Llano uplift in central Texas to Archer County. The larger part of past oil and gas production has come from Pennsyl-vanian strata, with less important amounts from the Permian, and minor but increasingly important quantities from the Ordovician. Developments during 1942 Approximately 500 wells were drilled in the district during 1942, in contrast with some 1600 wells drilled during 1941, and only some goo oil wells were completed, a decrease from 1941 of more than 10 per cent in the number of successful drilling operations of the total wells drilled. Production increased slightly over the previous year in spite of a decrease in the number of producing wells. The average depth of discovery wells was some 600 ft. deeper than in 1941. Although more discoveries produced from horizons above the base of the Strawn series than those below the Strawn, the deeper horizons, including the Bend series of the Pennsylvanian, the Mississippian limestone, and the Simpson and Ellenburger of the Ordovician, furnished most of the addition to oil reserves. The largest single addition to the reserves was again in the Ellenburger formation in the K. M. A. field. During the year 47 wells were completed, bringing the total developed acreage in this horizon to some 5450 acres. The most active area was Montague County, in which eight fields, or potential fields, were discovered although only 39 oil wells were completed. K. M. A: Pressure Maintenance Satisfactory progress was made by the K. M. A. Pressure Maintenance Association in 1942, At the end of the year, 100 wells were being used for input and total monthly volume of gas returned to the reservoir reached a high of 430,971 M cu. ft. in December, with an input ratio of 587 cu. ft. of gas returned for each barrel of oil produced during the month. This represents a utilization of 50 per cent of the gas available for return during that month. Comparison of the present average bottom-hole pressure by means of a decline

Jan 1, 1943

By B. M. Larsen, L. S. Darken

SOME years ago we collated all laboratory data then available to us on the distribution at equilibrium of manganese and of sulphur between metal and simple slags, and used the results in setting up an equation, which proved to be in accord with determinations of the distribution ratio of manganese and of sulphur between liquid steel and open-hearth slag, over a wide range of operating conditions. We now present an outline of our calculations and inferences on these two closely related questions, which are of especial importance because of the advisability of conserving manganese wherever possible. In general. our results show that in the open-hearth furnace equilibrium between metal and slag, with respect to both manganese and sulphur, is substantially attained during the finishing period of the heat, or in about a half hour to an hour after any additions were made. This implies that the end result will be identical whether manganese is added as alloy or metal to the liquid steel or as oxide to the slag. The proportion of total manganese retained by the steel is favored: (I) by use of the minimum amount of lime, consistent with other operating limitations, required to give a lime-silica ratio of about 2.4, which usually corresponds to a minimum practical slag volume; (2) by a high concentration in the slag of MnO relative to FeO; (3) by high operating temperature. Absorption of sulphur by the slag from the liquid metal is favored: (I) by large slag volume; (2-by a high concentration of active (or "free") lime and of manganese oxide dis) solved in the slag; (3) by low concentration of silica, iron oxide, and phosphate dissolved in the slag. DISTRIBUTION OF MANGANESE BETWEEN SLAG AND METAL AT EQUILIBRIUM On the basis that the distribution of manganese between liquid steel and slag is controlled by the reaction [FeO + Mn ;7± MnO + Fe [A] slag metal slag metal] the distribution at equilibrium, at a given temperature, is expressed by the equilibrium constant for this reaction, which is defined as [K" aFPI OFeO ? aMn] where a is the activity of the molecular species denoted by the subscript. The activity of each of the several components, although it is not known precisely, is approximately equal to the concentration of each component expressed as its mol fraction; that is, as the ratio of the number of mols of the given component to the total number of mols of all the molecular species present in the metal or slag, as the case may be. Moreover, since the activity, or concentration, of iron in the liquid metal is in fact substantially constant, the term aFe may he combined with K" to give a new constant K'. Making these changes.

Jan 1, 1942

By R. D. Vaughn

The limiting cases of non- Newtonian fluids flowing inside a concentric annular duct are developed without using a model of the fluid behavior. The solutions provide limits with which to test the various models of fluid behavior such as the power law and Bingham plastic models. The results of previous theoretical work are discussed in terms of limiting cases. This limiting case study also shows that experimental work on flow of non-Newtonian fluids in annular ducts should be confined to ducts for which the ratio of the radius of the inner wall to that of the outer wall is less than 0.3 and preferably less than 0.2. INTRODUCTION During the last 10 years the problem of laminar flow of non-Newtonian fluids in concentric annuli has received much attention large1y because of its application to the hydrodynamics of the wellbore. Recently solutions utilizing the power law1,2 and Bingham plasticl,3,4 models have been published. In this paper the method of limiting cases, which has been successfully applied to laminar-flow heat transfer516 will be applied to the problem of flow of non-time dependent, non-Newtonian fluids through annuli. This method permits solutions for the limiting cases to be made without using a model of unknown validity. The solutions, therefore, provide limits with which to test the various models which have been or will be proposed. A pertinent conclusion concerning the region of experimental work is also provided. DEVELOPMENT OF LIMITING CASES The limiting cases for the axial flow of fluids in concentric annuli may be defined with reference to Fig. 1. It is possible to define two limiting cases which pertain to the physical dimensions of the annulus. First, the annulus must degenerate to a circular pipe as the radius of the inner wall decreases or, as k = (k R/R) ? 0. Second, the annulus must approach the limit of parallel plates of infinite extent as the spacing between the inner and outer tubes becomes small in comparison with the radius R of the annulus, or as k?1. It is also possible to ascertain three limiting cases which pertain to fluid behavior. With reference to Fig. 2, as a fluid becomes progressively more pseudoplastic, the shear stress-shear rate relationship progressively approaches the indicated horizontal line more closely. At this limit the shear stress becomes independent of the shear rate. At the other extreme of increasingly dilatant behavior, the vertical asymptote is approached. Intermediate between these two limiting cases lies the case of the Newtonian fluid. Fluids which exhibit a yield shear stress also approach the limiting case of "infinite" pseudoplastic behavior. The remainder of the theoretical development is concerned with the description of and the solutions

By R. S. French, W. R. Hibbard

FOR tensile deformation, if the stress value is defined by the ratio of the load to the actual area, and the strain value by the natural logarithm of the ratio of the immediate length to the original gauge length of the sample, the resulting data when plotted with logarithmic coordinates result in a linear relation extending from the initial plastic yielding through the maximum load and in some cases, such as for copper and brass, to fracture'.'. The normal form of the curve, therefore, is a parabola, satisfying the following equation: S = Kem [1] where S and e represent true stress and strain, K is a constant, and m is a coefficient evaluating the slope of the curve and designated therefore as the strain hardening coefficient. Hollomon² has reported that this coefficient is a function of the yield strength as effected by grain size in alpha brass, and the carbon content and heat treatment of steel. In the case of steel, the value of the proportionality constant is a function of the carbon content. It was further shown in that paper2 and also by Low and Prater -hat, if the power relation between stress and strain remains valid, the strain hardening coefficient m, is equal to the strain at the maximum or necking load. Thus, an added significant feature is that this coefficient may serve as an index of the capacity for deformation prior to localized necking and subsequent fracture, e.g. drawability. The constant K may be defined as the strength of the material at unit strain. Together the K constant and the m coefficient describe an alloy's physical properties as far as strength and ductility are concerned for a particular yield strength and, in the case of steel, a particular carbon content. The work of Lacy and Gensamer4 gives data on alloy ferrites from which very close estimates of the strength of ternary and higher iron alloys have been made. It was the purpose of this investigation to develop similar data for commercial copper base alpha solid solution alloys. Preparation of Specimens Commercial binary alloys of copper with aluminum, beryllium, cadmium, nickel, silicon, tin and zinc were investigated. In addition, laboratory castings were made of a 1 pct silicon alloy and a 2 pct aluminum alloy to complete these particular series. Two types of copper were used, oxygen free OFHC and tough pitch FC. The analyses of these alloys are shown in table I. Standard ASTM sheet tensile specimens were used to determine the usual physical properties of each alloy and to obtain stress data at strains from 0.0005 to 0.01. Small (1/4 in. x 3 in. x tk.) samples

Jan 1, 1951

By Carl J. McHargue

Mechanical twins were produced in electron-beam melted columbium by high-speed impact at room temperature and by slow or fast compression at -196°C. The composition plane of the twins was { 112} and the shear direction was <111>. Notches in the twin bands often corresponded to traces of {110) of the matrix and appeared to be untwinned regions. Markings within the twin bands were interpreted as resulting from {110} slip in the twins. THERE has been much work in recent years concerning plastic deformation by glide, and the dislocation theory relating to glide has reached a relatively high degree of development. On the other hand, there have been fewer studies of deformation by mechanical twinning, and understanding of this process is far from satisfactory. This method of deformation is of interest for at least two reasons. First, it provides a mechanism in addition to glide for the relief of stresses, and, in the bcc and hexagonal close-packed metals may result in significant amounts of plastic flow. Secondly, there is the possibility that twins may act as barriers for dislocation movement, resulting in pile-ups which could nucleate cracks. As might be expected, the bulk of the literature on mechanical twinning in the bcc metals is concerned with iron. A good summary of the work done prior to 1954 is contained in the book by all.&apos; Recently the refractory bcc metals have become increasingly important. Limited studies have shown that tantalum,2,3 molybdenum,4,5 vanadium,6,7 tungsten,&apos; and columbium9-11 deform by mechanical twinning under some conditions. Alloys of molybdenum with rhenium and tungsten with rhenium show extensive deformation by twinning at room temperature.I2-l4 Most of these studies have dealt primarily with mechanical properties at low temperatures or have shown the existence of twins, and there is only a small amount of information concerning the conditions under which they form. The subject of the present paper is the formation of twins under stress in columbium with a consideration of their morphology. EXPERIMENTAL PROCEDURE The material used for these studies was taken from an ingot of columbium which had been melted twice by the electron-be am-method. The analysis of the ingot was (in ppm): B < 1, C = 10, Fe < 100, The cast ingot contained very large grains, and it was possible to obtain single-crystal prisms which measured from ¼ to 3/4 in. on a side. A few experiments were conducted on polycrystalline plate which was prepared by rolling material from the same ingot at room temperature and annealing at 1000 in a dynamic vacuum of 10-6 mm Hg. This gave a plate in which the grains had an average diameter of 3 mm. After the specimens were cut from the ingot, the six faces were metallographically polished and elec-tropolished to remove all traces of cold work. Most of the observations were made on the surfaces of the deformed specirllens without further treatment. Occasionally, etching after deformation was desirable. In these cases, an etchant of the composition 50 parts H2O, 5 parts HNO3 25 parts HF, and 10 parts H2SO4 was found to delineate the twins very well. Unless considerable care was taken to ensure the removal of all disturbed metal left by the mechanical polishing, etching failed to reveal many of the features discussed in this; paper. The specimen&apos;s were deformed either by impact or slow compression at 77°K (liquid-nitrogen coolant), 198°K (dry ice and acetone coolant), and 298°K. The impact load was delivered by a hammer except in one case where the load was delivered by a bullet. Slow compression was carried out on a hydraulic testing machine equipped with a chamber to hold the coolant. EXPERIMENTAL RESULTS It has been generally believed that the conditions favoring the formation of deformation twins are large grains, low temperature, and impact loading. In fact, Barrett and Bakish2 found twins in tantalum only after impact deformation at 77°K, and Adams, Roberts, and Smallman10 observed twins in columbium only at 20 For these reasons, the initial experiments of this study used impact loading. Hammer blows caused many bands resembling twins in single crystals a.t 77" but not at 198°K. Only a few slip lines were observed on any of the single-crystal specimens of this study—essentially all the deformation occurred by twinning. The appearance of the twins on the as-deformed surface is shown in Fig. 1. Although both Figs. 1(a) and l(b) are photomicrographs of twins taken at the same magnification and from the same crystal, they are startlingly different in appearance. Fig. 1(a) was taken from the crystal face approximately perpendicular to the shear direction, whereas Fig. 1(b) was taken from

Jan 1, 1962

By Kenric C. Owen

INTRODUCTION Situated 23 miles east of Pretoria the Premier Mine started diamond production in 1903. Two years later it produced the largest diamond yet discovered, the 3 106 carat Cullinan stone. For the period to 1932 when the mine was closed down due to economic circumstances, open cast mining was practised using inclined rope haulage traction and by this time had produced 48.5 million carats of predominantly industrial quality diamonds from 148 million tons of ore. In 1945 it was decided to re-open the mine using underground mining methods. A system using long hole benching was devised (Hodgson and Sewel 1960). This system was later modified but still forms the basis for approximately fifty percent of production. Following the successful introduction of block caving on associated diamond mines in the Kimberley area (Gallagher and Loftus 1960) this mining method was introduced to the Premier Mine in the late 60&apos;s and now accounts for nearly 50% of production. This paper will discuss block caving at Premier Mine in the light of our experience and discussion will be directed to the main geological and structural features of the orebody and host rocks and the related constraints imposed on the mining methods. Reference is made to the long hole benching mining method for comparative purposes. GENERAL DESCRIPTION OF OPERATIONS The ore from the block caving and bench mining areas gravitates via ore passes to a twin haulage, 500m below surface. Electric 13 tonne locomotives, each hauling a train of 10 Granby type cars, deliver the ore to two 42" x 48" (1.06m x 1.22m) jaw crushers. The ore is reduced in size to minus 0.15m before being hoisted to the treatment plant on surface. Hoisting is done with 12.5 tonne bottom discharge skips in a five compartment rectangular shaft using two 3 240 H.P. semi automatic Ward Leonard winders. A small single drum service winder operates in the other compartment. Men and materials travel in a separate shaft in a 5,4m x 2,8m cage operated by a Koepi winder. These two shafts are also the main intake airways to the mine. The two main extraction fans capable of 250 m3/s at 3,2 kPa are situated on surface and are connected to the underground workings by an incline and a network of air passes and return airways. GEOLOGY Although there are numerous occurrences of kimberlite bodies in the district, the Premier Mine kimberlite pipe is the only economic orebody. It is roughly oval in shape with surface dimensions of 900m on the long axis and 450monthe short axis. The surface area is 320 000 m2. The kimberlite has intruded a massive body of felsite and norite which is intersected by a number of faults. The contact between the kimberlite and host rock is clear cut and dips inward at an average angle of 80°. Cutting across the pipe and the country rock is a 75 metre thick gabbro sill, the top of which intersects the pipe below the 347m level and dips at 15&apos; to the north west. The kimberlite has been metamorphosed a distance of some 20m both above and below the sill contacts. Age measurements on biotite from the gabbro date the sill as 1 115 million years, thereby providing a minimum age of the pipe. (Pre Cambrian). All other pipes in South Africa are of Cretaceous age (60 million years). A simplified geological plan and section of the Premier Mine orebody is shown in Figure 1. It is thought that the kimberlite intruded in at least three distinct phases. The kimberlites of these different phases can be distinguished most easily by their characteristic colour. BROWN KIMBERLITE. This is the oldest kimberlite and forms a crescent shape in plan on the south eastern side of the pipe. In depth it increases in relative area. Although the brown kimberlite carries the highest diamond grade in terms of carats per ton all the various kimberlites in the Premier pipe are economic to mine.

Jan 1, 1981

By Appendix by A. S. Goldoni, R. E. Tate, E. M. Cramer

The diffision coefficient for self-diffision of plutonium in the temperature range 350" to 440°C has been measured by using puZ3 as the tracer isotope. Autoradiopaphic techniques were used to inzlestigate the possibility of grain boundary diffusion, hut only evidence for volume diffusion was found. The least-squares fit of the data gives the .following equation for the diffision coefficient: The computer least-squares technique for fitting the nonlinear equations is outlined. DETERMINATION of the self-diffusion coefficient of the fcc 6 phase of plutonium is in principle a straightforward experimental task. However, the chemical reactivity, the intense @ activity, and the toxicity of plutonium put limitations on the experimental techniques. The technique selected included roll bonding for preparation of the diffusion couples and pulse-height analysis of the @-particle activity to determine the distribution of the PU tracer in the diffusion couples. EXPERIMENTAL PROCEDURE Couple Preparation. Cylinders about 0.5 in. in diam were cast from two special stocks of plutonium, one of which had been enriched in puZ3&apos;, as shown in the isotopic analyses listed in Table I. For each roll-bonded composite sheet, a cylinder 0.437 in. in diam and 0.190 in. thick was turned from each kind of plutonium on a lathe in a 98 pct He atmosphere. The two freshly machined cylinders were positioned face to face in a tube of commercially pure aluminum which had been sealed at one end by welding and the assembly was evacuated on a vacuum manifold overnight. The elapsed time between machining the faces of the plutonium cylinders and evacuating the loaded tube was about 15 min. After overnight evacuation of the assembly the indicated vacuum was 1 X 10"5 torr or better. The aluminum tube was then warmed and pinched off with a cold-welding tool. The pinched-off weld was also fusion-welded as an additional precaution against leakage of air into the evacuated assembly. The assembly was immediately heated for 30 min in a 250°C furnace and reduced in thickness by being passed through a rolling mill with rolls heated to 200°C. The rolling schedule (four passes of 100 mils each with a 10-min reheat after two passes) reduced the thickness of the assembly to 100 mils. The rolled assembly was allowed to cool normally in air. The rolled assembly was sheared at the edges and the aluminum peeled from the composite plutonium sheet. The elliptical sheet was about 0.060 in. thick and usually three 0.383-in.-diam disks could be punched from its central portion. The sheet was heated on a hot plate to the ductile low 0 range (140°C as measured by temperature-indicating pellets) before each disk was quickly punched. The quality of the bond in the disks was evaluated by metallographic examination of the scrap sheet adjoining the hole left by the punch. Only well-bonded specimens without oxide in the interface (as indicated by metallography) were considered satisfactory for further use. More than half of the specimens so examined contained sufficient oxide in the interface to be rejected. Diffusion Anneal. Each disk to be diffusion-annealed was wrapped in 1-mil-thick tantalum foil and sealed within a Pyrex capsule evacuated to 1 x 10~5 torr or better. This capsule was then sealed within another Pyrex capsule at a similar pressure. The diffusion anneals were carried out at temperatures between 350" and 440°C in Marshall furnaces adjusted to have a temperature gradient of not more than *1/2"C over a 5-in. length. This gradient was then further smoothed by using a nickel tube as a liner in the furnace. The liner was divided into three longitudinal cavities by a septum of nickel sheet to which two calibrated Chromel-Alumel thermocouples were attached. One thermocouple was used for controlling the furnace temperature by means of a Brown Pyrovane controller equipped with a Capaciline anticipation circuit; the second thermocouple was monitored twice daily with a Leeds and Northrup K-2 precision potentiometer.

Jan 1, 1964