Search Documents

By A. Satter, J. M. Campbell

Reported herein are the results of a careful and detailed study of the non-ideal behavior of pure gases and their mixtures. Included are: (1) new data on five ternary systems composed of methane, ethane and H2 S; (2) a simple compressibility factor correlation that is inherently superior to present correlations, particularly for gases containing H2S and CO2; and(3) a detailed study of combination rules and the effect of system composition on the choice thereof. This study makes use of the rather large mass of data already available in the literature. A complete re-examination of the data and ideas presented in the last 25 years was considered desirable as a prelude to our basic concern — the effect of diluents on gas behavior. A consideration of both the macroscopic and microscopic properties of gases provides a better insight which, in turn, gives a firmer basis for improved correlation techniques. Such a study has shown that expressing the compressibility factor Z as a function of acentric factor w, as well as reduced temperature and pressure, yields a correlation that is broader in scope. The study of various combination rules has shown that better results are obtained by "tailoring" the rule used to the system composition. To do so improves the basic reality of results by overcoming some of the anomalies often found when using Kay's rule alone. Tentative recommendations are made regarding the most reliable combination rule for use with a given class of gas. The data presented are useful for estimating the direction and magnitude of the expected deviation when using a given rule. Although more work is needed, particularly around the critical region and with CO2 mixtures, the advantage of the classification scheme proposed is apparent. INTRODUCTION When one attempts to write a PVT equation to fit the data for actual gases, greater precision is obtained by the use of a multiple number of empirical constants. This has lead to multiple-constant equations such as Benedict-Webb-Rubin, Beattie-Bridgman, Keyes, etc., which are capable of yielding very precise results for pure gases in a range for which data to get the constants are available. As a matter of practicality, though, the use of such equations for gas mixtures is limited. Because of the infinite number of gas analyses available, any attempt to compile the constants needed requires a prohibitive amount of experimental data. This could be overcome by the use of a combination rule, but there is no real advantage in doing so because the end result offers no practical impovement over the Z factor correlation. The most widely used method of predicting the volumetric properties of pure gases is based upon the "theorem of corresponding states". According to this theorem, "all pure substances have corresponding molal volume at corresponding temperature and pressure if the reference point of correspondence is the critical point". Generalized compressibility charts for gases were prepared first by Cope and associates1 in 1931 and later by Brown and co-workers2 in 1932. However, the most commonly used charts are those of Dodge,3 Nelson and Obert,4 Hougen and atsson: and Standing and Katz.6 The work of Katz and co-workers has provided us with basic data for the hydrocarbons most widely used today. Their original chart6 was compared with a relatively large amount of multi-component data for gases consisting almost entirely of normal paraffin hydrocarbons. A deviation of only + 1.2 per cent was obtained.39 In the 20 years following publication of this work it has been found that the behavior of most mixtures of paraffin hydrocarbons could be predicted by this correlation within at least 5 per cent. Where difficulty has been encountered it has largely involved one or more of the following circumstances: pressures above 4,000 psig, mixtures containing large amounts of heavy ends and/or aromatics, systems in the critical region and mixtures containing polar compounds and/or CO2. The abnormal error sometimes found with such gases, not too unexpected for this method, is

By H. Pew, J. H. Field

To alleviate pollution more restrictive legislation is being enacted, either limiting emission of pollutants or the type of fuel that can be utilized. The nature and magnitude of air pollution problems affecting the mining, preparation, coking and combustion of coal are described. Methods for combatting particulate emissions by use of mechanical separators and electrostatic pre-cipitors are discussed. Proposed methods to meet the problem of gaseous emissions currently receiving considerable attention are described, with special emphasis on methods to decrease pollution by sulfur oxides. Concern about air pollution goes back several centuries, but until very recently most effort has been aimed at coal smoke and other visible pollutants. The classic example of a 'successful' campaign for smoke abatement and control is the fruitful combined effort of the city of Pittsburgh and its surrounding Allegheny County, which eventually led to the reconstruction of downtown Pittsburgh at an estimated cost of one billion dollars. Historically, the city's downtown Golden Triangle district had been afflicted by pollutants evolving from steel mills, from a variety of other industries, and from railroad locomotives. Efforts to alleviate the situation prior to 1943 were virtually ineffective. In 1945, however, a comprehensive redevelopment plan was prepared and backed by state authority. Within a few years a clean, modern metropolis has evolved where once stood America's famous 'smoky city.' But the victory in Pittsburgh, as in various other American cities, has not solved the national problem. Current estimates indicate that 133 million tons annually of air pollutants from all sources still are emitted annually into the atmosphere above the United States. About 10% of this annual effluent is particulate matter so that most of the remaining pollution problems will be solved only when other effluents are reduced. Essentially, these are sulfur oxides, nitro- gen oxides, hydrocarbons, and carbon monoxide. Over the years, both states and local communities have tended to increase the restrictions on smoke and fly ash — problems mostly of concern in the combustion of coal. Prior to the middle 1950's, ordinances sometimes permitted emissions of smoke equivalent to as much as No. 3 on the Ringlemann scale. Since 1956, no ordinance has been passed which allows smoke of greater than No. 2. Under today's conditions of improved fuels, equipment and practice, a few communities have passed laws prohibiting emission of smoke of any density darker than Ringlemann No. 1. The majority of existing laws on fly ash emission in the U.S. limit emissions equivalent to 0.85 lb of fly ash per 1000 Ib of flue gas. In recent years, however, regulations which have been adopted give cognizance to the higher level of performance now obtainable with improved equipment. A comparison of the restrictions of five codes adopted since 1960 is given in Table I. The most stringent of these is the one for New York City which provides for a maximum emission of 0.6 lb fly ash per million Btu heat release (equivalent to roughly 0.51 lb/1000 lb of flue gas). The first comprehensive effort to restrict the emission of SO2 resulted from the passage of a 1937 law in St. Louis. This regulation stipulated that coal containing in excess of 23% volatile matter and 2% sulfur must be washed, thereby presumably producing some effective reduction in the input sulfur content. This was followed in 1949 by a Los Angeles County law which prohibited the emission of SO2 in concentrations greater than 0.2%. Most SO2-restrictive legislation passed since that date has been based on this limiting 0.2% SO2 by volume, although modifications are occasionally permitted under selected conditions, sometimes based on the fact that certain limiting ground concentrations are not exceeded — such as in the rules adopted by the San Francisco Bay area. To date, no legislation has been passed in the U.S. to limit the generation of nitrogen oxides from the combustion of fossil fuels. However, such oxides are considered to be of potential importance in air pollution control because of their possible detrimental effects on health and their reported role in the formation of photochemical smog. Interest in reducing oxides of nitrogen from powerplant and auto exhausts is increasing and regulations limiting their quantity can be expected in the future.

Jan 1, 1968

By R. T. Hukki

THE accepted basis of comparisons between mills of different diameter is the percentage critical speed. If n = actual mill speed, rpm, nc = calculated critical speed, rpm, np = calculated percentage critical speed, and D == inside diameter of the mill in feet, then n, In the following analysis capacity, T, is expressed in short tons per hour, tph, and power consumption, P, in kilowatts, kw. Accordingly power consumption per unit of capacity, P will be expressed in kilowatt hours per short ton, or kw-hr per ton. In all equations D refers to the inside diameter of the mill in feet and v to the peripheral speed of the mill in feet per minute inside the liners. ' Comparison between separate mills must be based on equivalent grinding conditions, i.e., same feed, same size distribution of feed, same size distribution of product, and same percentage of solids. In addition, comparisons between separate rod mills must be based on the same rods, same type of liners, and same percentage rod load. Comparisons between separate ball mills presuppose the same balls, similar liners, and same relative ball load. The practical np-range through which the equations apply varies, being narrower for fine grinding in ball mills and wider for coarse crushing in rod mills. The Relationship between Capacity and Speed It is the general belief that the capacity, T, of a tumbling mill is directly proportional to the speed of the mill, other things remaining constant.' Mathematically this is represented by the equation T - c¹ n tph [4] where c, is a factor related with the grinding characteristics of the ore, method of reduction, and the units chosen. It is proposed here that the general equation relating mill capacity and speed should be of the form T = c¹ nm tph [5] In other words, the capacity should be proportional to the mill speed raised to power m, the numerical value of the exponent being 1 5 m 5 1.5, depending on the circumstances. Eq. 5 can also be written in the following forms: T = c, (np)m tph, and [6] T=Ca vm tph, [71 where v = peripheral speed of the mill in feet per minute. If the observed capacity of a mill at speed n¹ is = T¹ tph, the capacity T² of the same mill at speed n² should be T² = T¹ (n²/n¹)tph [8] The Relationship between Power Consumption, Mill Diameter, and Speed The only well known theoretical deduction relating power consumption, P, and mill diameter appears to be the formula of duPont introduced by Gow, Guggenheim, Campbell, and Coghill.' According to duPont, the power required to operate a mill is a function of the mass of the balls, of the lever arm of the ball mass, and of the speed of the mill. The ball mass per unit of mill length is proportional to the square of the diameter, the lever arm is directly proportional to the diameter, and the critical mill speed or any percentage thereof is inversely proportional to the square root of the mill diameter. Following this reasoning, the original duPont formula is of the form P = c4D² c D • c6D-0.5 = c7D2.5 [9] If the mill speed in the above equation is expressed in terms of Eq. 3, the duPont formula may be written as follows: P=f1(D2) f2(D) f³(—np) or [10] vD P = c np D2.5 kw [11] Eq. 11 may also be derived from the mechanical principle of force, which is equal to mass x acceleration. Power necessary to operate a mill may be considered to be an homogeneous linear function of the force developed. Ball or rod mass per unit of mill length is a function of D2. The acceleration factor of the ball or rod mass is a function of the peripheral speed of the mill. Thus P = f4(F) = /x(D2) f5(v) Indicating that v = nDn, and n = c9 np /vD, the above equation becomes P = f2 (D2) -fa (D ca np/vD

Jan 1, 1955

By Charles Y. Clayton

THE general public does not know that there is such a thing as metallurgy and it is very seldom that you see the word metallurgy in print except in technical magazines. Perhaps it is more to the front now than in former years; the Dee. 30, 1929, issue of Time mentioned under the heading of "Science7' certain new achievements in metallurgy, which may be cited as a rare instance of the use of the word in anything but the technical press. Too often the part the metallurgist plays in the advancement of the different branches of engineering is overlooked. To say that the public does not know there is such a profession as metallurgical engineering seems to be a bold statement, but had there beep time a questionnaire could have proved this beyond doubt. Suppose a carefully-planned questionnaire had been sent to all high school students asking them, among other questions, to list the different branches of engineering, I believe that the several professions would have been listed in this manner: (1) Electrical engineering; (2) mechanical; (3) civil, and (4)' mining. One in each hundred might have listed metallurgy as an engineering profession. Even the Society for the Promotion of Engineering Education in its Bulletin No. 10-"Investigation of Engineering Education" lists Civil, Mining, Mechanical, Electrical, and Chemical, leaving out Metallurgical Engineering.

Jan 1, 1930

By W. Bonfield

The temperature dependence of the microscopic yield stress (the stress to produce a plastic strain of 2 x 10-6 in. per in.) and the stress-plastic strain curve of polycrystalline Cu 1.9 wt pct Be have been measured for the solution treated condition, an intermediate condition containing G.P. zones and ?' precipitate and the overaged ? precipitate condition, in the range from -58° to 200° C. A transition in micro -yield behavior and a large temperature dependence were noted for the intermediate condition, which are interpreted in terms of the interaction of glide dislocations with two differently sized zones. In comparison the microscopic yield stresses of the solution treated and overaged conditions were less sensitive to temperature variations and are satisfied by the Mott-Nabarro and dislocation bowing theories, respectively. A determination of the temperature dependence of the yield stress of a precipitation hardening alloy has provided a powerful tool for evaluation of the operative deformation mechanism. There is a marked contrast between the effect of temperature on the yield behavior of a metal containing coherent zones or intermediate precipitates, which can be "cut through" by mobile dislocations, and a metal containing a dispersion of noncoherent particles, through which dislocation "bowing out" is the dominant role of deformation.' These studies have in general been confined to single crystals, as it was considered that similar experiments on polycrystalline material did not produce good data because of the lack of sensitivity with which the yield stress could be determined. However, this objection has been removed by the introduction of mi-crostrain techniques, with which the yield stress in polycrystalline materials can be measured to a strain sensitivity of 10-6. Such measurements have not only shown that the deformation of polycrystalline precipitation hardening alloys can be examined with the same detail as single crystals, but also that some unexpected results are obtained.' In this paper the results obtained from a study of the temperature dependence of the microscopic yield stress (the stress to produce a plastic strain of 2 x 10-6 in. per in.) and the stress-plastic strain curve of a polycrystalline Cu 1.9 wt pct Be precipitation hardening alloy (Berylco 25) are discussed. The temperature dependence of the alloy was measured for three different conditions: 1) The solution treated condition (a supersaturated solid solution of a containing ~12 at. pct Be3) which is obtained by water quenching the alloy from 800° C. 2) The condition of y' intermediate precipitate, to- gether with some G.P. zones,' which is produced after an aging treatment of 2 hr at 315°C from the solution treated condition. (The alloy was cold rolled to 40 pct reduction prior to aging to minimize grain boundary precipitation effects.)4 3) The condition with equilibrium ? precipitate structure2 which is developed after an aging treatment of 24 hr at 425° C. EXPERIMENTAL PROCEDURE Tensile specimens of gage length 1 in. and with rectangular cross section of 0.18 by 0.06 in. were prepared from the solution treated, cold rolled alloy and were either resolution treated for 1 hr at 800°C, followed by water quenching, or aged for 2 hr at 315°C and 24 hr at 425° C to produce the desired precipitate structures. The microstrain characteristics of the aged specimens were determined at temperatures from —58" to 200° C and those of the solution treated specimens from -58° to 30° C. Each temperature was controlled to ± 0.2°C, which was a level of stability sufficient to eliminate thermal expansion effects from the measurements (~1.2°C temperature increase produced an extension of 2 x 10-6 in.). The microplastic behavior of the specimens in the temperature range below 82" C was measured with a standard Tuckerman strain gage,5 while at temperatures above 82°C a modified Tuckerman gage with a reduced strain sensitivity (4 x10-6 in. per- in.) was used. A load-unload technique was used to establish values of the microscopic yield stress. The specimen was strained at a constant cross head speed of 2 x 10-2 in. per min to a given stress level, at which the total strain was measured. Then the specimen was immediately unloaded at the same rate and any residual plastic strain determined. This procedure was repeated for an increasing series of stress levels until the microscopic yield stress was established by a direct measure of the stress to produce a residual plastic strain of 2 x 10-6 in. per in. (It should be noted that, as reversible dislocation motion occurs at stresses less than the microscopic yield stress,2 the plastic strain rate at this level was not constant.) In an ideal test, the microscopic yield stress would be determined from a continuous stress-strain measurement, rather than from a load-unload sequence, in order to eliminate mechanical recovery effects.6 However, it was found experimentally that mechanical recovery was negligible in Cu 1.9 wt pct Be at small plastic strains for all the temperatures investigated, as the microscopic yield stress was independent of the number of load-unload cycles employed (i.e., the values measured for specimens subjected to different numbers of cycles was within the experimental scatter determined for specimens tested in an identical manner). Therefore, it is reasonable to consider the microscopic yield stress determined in the load-unload

Jan 1, 1969

By Kenneth M. Bravinder, Jonathan E. Koogle, Dean H. Sheldon

THE solution of subsurface geological problems requires an analysis of vertical and horizontal dimensional relationships. For many, the ability to visualize structures in three dimensions is not easily acquired, and for areas where the structural conditions are complicated by faulting the need arises for illustrative models. This has led to the many forms of peg models, in which the usual course of construction includes a peg for the well bore, and the suspension of strings from peg to peg to illustrate the formation correlation between wells, thus developing the usual folds and flexures. Such models proved unsatisfactory for the analysis of fault problems. Consideration was then given to the use of transparent materials. A model of the Dominguez oil field was constructed in 1935, which consisted of a series of glass panels from well to well on which the formations were shown in color. This model, though serving its purpose at the time, proved to be impractical from the standpoints of cost, breakage hazard, and extreme difficulty in maintenance. The cost of materials was estimated to be approximately $300, and the construction required the services of one man for a period of from 10 to 12 months. Maintenance after actual construction was an additional expense. Consideration was then given to the use of celluloid, and after some experimental work it was found that the advantages in its use were: (I) a short construction time, (2) ease of handling the material, (3) sections are easily reworked, (4) low cost. CONSTRUCTION It would be of little value to recount here the difficulties encountered in the building of the first celluloid model. The conclusions drawn, however, suggest the following procedure for those who contemplate the construction of such a model: I. Determination of the Scale.-Accurate geological interpretation requires the use of the same scale horizontally and vertically. It should be emphasized that the limited size of the celluloid sheets (20 by 48 in. maximum dimensions), the area to be studied, the maximum depths to be plotted, and area of possible future development are the controlling factors in adopting a scale to be used. 2. Construction of the Base.-A satisfactory base may be constructed by using two pieces of three-ply wood separated by strips of i by 2-in. material. This type of construction prevents warping of the base. The size of the base is determined by the dimensions of the area to be studied plus 5 in. on all sides. A map of the same scale as that adopted for the model, indicating the lines of section to be used, should be fastened to the base. 3. Construction of the Frame (Fig. I).-In the construction of the frame the use of 1/2-in. o.d. pipe welded in a box shape is suggested, with two horizontal rods on each of the four sides, top and bottom. The top rod should be at least I in. higher than the

Jan 1, 1941

By Howard A. Meyerhoff

IN western Argentina, in the Province of Mendoza and the Territory of Neuquén, there is a series of solid bitumen deposits which are claimed to be the most extensive in the world. In a linear belt 500 km (311 miles) long, roughly paralleling the axis of the Andes but occupying a variety of positions with respect to the orogenic structures, more than 100 outcrops have been reported. Approximately 70 of these have been prospected, and some mining has been planned or attempted for about half this number. Limited tonnage, high sulphur content, inaccessibility, labor shortages, and disastrous fires have taken a heavy toll among these operations, and in 1947, only 4 deposits were being actively mined, but 8 more were being explored and developed with a view to placing them in production. Inasmuch as these deposits are currently yielding a substantial fraction of Argentina's limited domestic supply of solid fuels, a brief account of the deposits and of the small industry based upon them will provide a commentary on Argentina's fuel problem. CHARACTER OF SOLID BITUMENS The solid bitumens of Mendoza and Neuquén are petroleum derivatives, as is indicated by the following ultimate analysis: [CONSTITUENTS PER CENT Hydrogen6.8 Carbon82.4 Nit.rogenI.6 Oxygen2.7 Sulphur2.2 Ash4.3 100.0] In physical appearance they bear a close resemblance to coal and are commonly called carbón. Their specific gravity has an extreme range of 1.10 to 1.777 in 36 samples weighed, but the average is between 1.14 and 1.24. The material is black, lustrous to dull, with blocky to sub-conchoidal fracture. Though brittle, the purest material is softer than anthracite and exhibits slightly more cohesion. This bituminous material is generally known as asfaltita or asphaltite, but definitive chemical studies indicate that the term asphaltite or grahamite can be applied correctly only to the exceptionally pure material of the Auca Mahuida district in Neuquén. Some other name should be used for the chemically variable bitumens found in all other districts, and the designation "asphaltitic pyrobitumen" has been proposed by Carlos A. S. Piscione.1 At the present time these bitumens, are finding favor as solid fuel and as a gas-producing medium. They yield heat values of 8650 to 9750 cal,2 or 14,500 to 16,750 Btu per pound, and in consequence of their excellent heating qualities and the scarcity of imported coal, they are marketed as far away as Buenos Aires at a price of go to 95 pesos, or approximately $22 to 23, per metric ton. In burning, the material swells and decomposes into solid, fluid, and gaseous fractions which burn violently though nonexplosively. Its combustibility has been responsible for disasters and virtually complete destruction at Mina La Esperanza and Mina Santa Marta, and for minor mishaps elsewhere, as at Mina Escondida in August 1947.

Jan 1, 1947

By Howard A. Meyerhoff

In western Argentina, in the Province of Mendoza and the Territory of Neuquen, there is a series of solid bitumen deposits which are claimed to be the most extensive in the world. In a linear belt 500 km (311 miles) long, roughly paralleling the axis of the Andes but occupying a variety of positions with respect to the orogenic structures, more than 100 outcrops have been reported. Approximately 70 of these have been prospected, and some mining has been planned or attempted for about half this number. Limited tonnage, high sulphur content, inaccessibility, labor shortages, and disastrous fires have taken a heavy toll among these operations, and in 1947, only 4 deposits were being actively mined, but 8 more were being explored and developed with a view to placing them in production. Inasmuch as these deposits are currently yielding a substantial fraction of Argentina's limited domestic supply of solid fuels, a brief account of the deposits and of the small industry based upon them will provide a commentary on Argentina's fuel problem. Character of Solid Bitumens The solid bitumens of Mendoza and Neuquen are petroleum derivatives, as is indicated by the following ultimate analysis : Constituents Per Cent Hydrogen.................... 6.8 Carbon...................... 82.4 Nitrogen..................... 1.6 Oxygen...................... 2.7 Sulphur...................... 2.2 Ash.......................... 4.3 100. 0 In physical appearance they bear a close resemblance to coal and are commonly called carbon. Their specific gravity has an extreme range of 1.10 to 1.777 in 36 samples weighed, but the average is between 1.14 and 1.24. The material is black, lustrous to dull, with blocky to sub-con-choidal fracture. Though brittle, the purest material is softer than anthracite and exhibits slightly more cohesion. This bituminous material is generally known as asfaltita or asphaltite, but definitive chemical studies indicate that the term asphaltite or grahamite can be applied correctly only to the exceptionally pure material of the Auca Mahuida district in Neuquen. Some other name should be used for the chemically variable bitumens found in all other districts, and the designation "asphaltitic pyrobitumen" has been proposed by Carlos A. S. Piscione.1 At the present time these bitumens are finding favor as solid fuel and as a gas-producing medium. They yield heat values of 8650 to 9750 cal,2 or 14,500 to 16,750 Btu per pound, and in consequence of their excellent heating qualities and the scarcity of imported coal, they are marketed as far away as Buenos Aires at a price of go to 95 pesos, or approximately $22 to 23, per metric ton. In burning, the material swells and decomposes into solid, fluid, and gaseous fractions which burn violently though nonexplosively. Its combustibility has been responsible for disasters and virtually complete destruction at Mina La Es-peranza and Mina Santa Marta, and for minor mishaps elsewhere, as at Mina Escondida in August 1947.

Jan 1, 1949

By H. R. Peiffer

The increase in electrical resistivity, ?pT,at 78°K was measured as a function of elongation, E, at 78°K for a 2 pct (approximately) by weight finely divided alumina in silver material. The amount of increase in electrical resistivity (?pT) is a function of the amount of alumina present and can be represented by: ?pT = cen where c and n are constants. The constant c increases and the exponent n decreases with increasing alumina content. Approximately 33 pct of the resistivity increase is recovered upon annealing at room temperature. RECENTLY, dispersed phase alloys' have shown promise for use at elevated temperatures. These materials have a high creep resistance, a high-temperature strength considerably above that of the pure matrix, and a high electrical conductivity. These materials also show other properties such as retardation of re crystallization and resistance to flow even above the melting point of the matrix.' It is believed that high-temperature creep is retarded in these materials because the dispersed phase blocks the motion of dislocations. According to schoeck3 and also Weertman,4 the rate controlling process is the climb of dislocations around the dispersed phase particles. Any mechanical property of these materials, or any other material for that matter, can be explained on the basis of defect structures and their interrelations. It would certainly be interesting to learn something about the effect of a dispersed phase on dislocation densities and also their effect on point defect concentrations following cold work. This is the purpose of the research discussed herein. The method used was the measurement of the electrical resistivity as a function of elongation at low temperatures. This scheme has already been applied to the studies of pure metals.5 Electrical resistivity measurements conducted on pure metals such as copper, silver, gold, and so forth, have indicated that defects such as vacancies and dislocations are formed during deformation.5 It has been shown by van Bueren6 that the resistivity change, ?p, should increase with elongation, E, according to the following relation: Ap=Al/2 + B£3/2 where A and B are constants of the system; ?p is in µ? cm. The first term refers to the amount of resistivity change associated with line defects and the second term with the change associated with point defects. A later calculation yielded exponents of 3/4 and 5/4, respectively.6 Manintveld7 has shown that the resistivity of copper increases with strain according to a 3/2 power law, whereas van Bueren has shown that the resistivity of silver increases according to a 5/4 power law. Assuming then that van Bueren's calculations are correct, it appears that line defects contribute very little to the electrical resistivity. There is some doubt, however, as to the validity of van Bueren's expressions for all cubic metals.' In fact, the data on copper and silver are not conclusive. In this in-

Jan 1, 1962

By M. King Hubbert

In 1856 Henry Darcy described in an appendix to his book, Les Fontaines Publiques de la Ville de Dijon, a series of experiments on the downward flow of water through filter sands, whereby it was established that the rate of flow is given by the equation: q= - K(h2- h2)/I1 in which q is the volume of water crossing unit area in unit time, I is the thickness of the sand, h1 and h2 the heights above a reference level of the water in manometers terminated above and below the sand, respectively, and K a factor of proportionality. This relationship, appropriately, soon became known as Darcy's law. Subsequently many separate attempts have been made to give Darcy's empirical expression a more general physical fornzulation, with the result that .so many mutually inconsistent expressions of what is purported to he Darcy's law have appeared in published literature that sight has often been lost of Darcy's own work and of its significance. In the present paper, therefore, it shall be our purpose to reinform ourselves upon what Darcy himself did, and then to determine the meaning of his results when expressed explicitly in terms of the pertinent physical variables involved. This will be done first by the empirical method used by Darcy himself, and then by direct derivation from the Navier-Stokes equation of motion of viscous fluids. We find in this manner that: q = (N4(p/p)[g - (l/p) gradPI = WE, is a physical expression for Darcy's law, which is valid for liquids generally, and for gases at pressures higher than about 20 atmospheres. Here N is a shape factor and d a characteristic length of the pore structure of the solid, p and IL are the density and viscosity of the fluid, a = (Nd2)(pip) is the volume conductivity of the system, and E = [g - (l/p) grad p] is the impelling force per unit mass acting upon the fluid. It is found also that Darcy's law is valid only for flow velocities such that the inertial forces are negligible as compared with those arising from viscosity. In general, three-dimensional space there exist two superposed physical fields: a field of force of charucteristic vector E, and a field of flow of vector q. The force field is more general than the flow field since it has values in all space capable of being occupied by the fluid. So long as the fluid density is constant or is a function of the pressure only, curl E = 0, E = - grad F where F = gz + dp/p The field of flow, independently of the force field, must satisfy the conservation of mom, leading to the equation of continuity div pq = — f ap/at , where f is the porosity, und t is the time. For steady rnotion ap/at = 0. and div pq = 0 .

Jan 1, 1957

By Potts N. L., Eilerts C. K., Summer E. F.

The second-order, nonlinear, partial- dillerential equation representing the transient radial flow of gas - condensate fluids in reservoirs has been integmted by using finite -dillerence equations and electronic computers. Effect was given to pressure -dependent permeability, viscosity, and compressibility and to distance-dependent permeability. The influence of a second-degree velocity term in the Darcy equation was investigated. Implicit methods were used and practical, convergent solutions were obtained with material balance to less than 6 x 10-4 for recovery of one- hall the reserve at constant flow rate. Integration results provide the productive period of a reservoir for a given constant rate and the fraction of the fluid initially in place that can be recovered in that period. The properties of a lean and a rich fluid are represented in a set of integrations designed to demonstrate the effect of different constant-recovery rates and significant variables emphasized one at a time. INTRODUCTION A method is needed for computing the availability or reserve of a gas in a formation, making use of all technical and engineering information that is pertinent. As a step in this direction, a program for computing transient linear flow1 was developed that utilizes principles of earlier finite-difference computing for a similar purpose2-4 and gives effect to significant pressure- and distance-dependent properties of reservoirs and their contents. The present paper pertains to transient, radial flow of gas-condensate fluids. In the partial-differential equation of flow, effect is given to the variables viscosity and compressibility,5 and also to permeability. Included in the study is the quadratic form of the Darcy equation of flow that has been the subject of field tests6 and that has been applied to a gas with constant properties in transient flow computing. 7 Inclusion in the differential equation of variable coefficients to represent properties greatly complicates the higher derivatives of the equation. Because it was impractical to make a required improvement in certain of the finite-difference derivatives by Taylor-series expansion, five-term derivatives were used in the implicit computing. Related methods were developed that can improve general facility in manipulating finite-difference equations. BASIC EQUATIONS PARTIAL DIFFERENTIAL EQUATION The basic partial differential equation for transient radial flow in the direction of decreasing radius is l_ a(vpr) _? .......(1) r gr at where r = radius, L v = apparent velocity aIong radius, L/t p = density, M/L3 t = time, t Ø = porosity, dimensionless. In this equation, only the porosity is regarded constant. The density is where p = pressure, m/Lt2 M = molecular weight of fluid, m T = constant temperature of fluid, T z(p) = pressure-dependent compressibility factor, dimensionless R = gas Constant, mL2/t2t. The physical boundary of the reservoir is at re, and the producing boundary is defined to be at tq with rq = re /200. By definition,

Jan 1, 1966

By H. I. Aaronson

Isothermally formed plates of proeutectoid a, lengthen at approximately the rates predicted by the Zener-Hillert equation, indicating that their lengthening kinetics are controlled by the volume diffusion of chromium in the ß matrix. Solution annealing treatments in the vicinity of the transus introduce a second mode of transformation, in which the lengthening rate is at least 1.8 x 104 faster than predicted by this equation. Little or no volume diffusion appears to accompany the latter transformation. The Presence of oxygen is evidently a Prerequisite for the occurrence of this reaction. X HE effects of thermal history prior to the initiation of transformations which proceed by nucleation and diffusional growth are currently evoking considerable interest. The rate of nucleation (but not the rates of growth) of pearlitel and of bainite,2 and the over-all rates of aging reactions in Al-cu3 and Al-si4 alloys have been shown to be appreciably dependent upon thermal history in the solid-solution region"' or in the two-phase region prior to the onset of pre~ipitation.3,4 In none of these reactions, however, was the transformation mechanism changed by variations in the temperature-time path used to reach the reaction temperature. Preliminary experiments reported by M. K. McQuillan5,8 and a systematic investigation by Farrar and Margolin7 have disclosed that the transformation behavior of hypoeutectoid Ti-Cr alloys is considerably more sensitive to thermal history. Qualitative evidence strongly suggests that not only the rate of nucleation but also the rate of growth and even the mechanism of the ß — a reaction can undergo remarkable changes when a solution annealing treatment in the vicinity of the ß/a +ß, or transus temperature, is interposed between a high-temperature solution anneal (conventionallq at 1000°C) and the isothermal reaction treatment. The present investigation was undertaken in an attempt to ascertain the nature and origin of these effects. Study of the changes in the kinetics and morphology of the ß - a, reaction produced by varialtions in thermal history was the primary means employed for this purpose. EXPERIMENTAL PROCEDURE The compositions of the alloys used in this in~esti~ation are given in Table I. These alloys were homogenized for 5 days at 1050° C in evacuated (2 x 10 -5 mm of Hg), zirconium-gettered fused quartz capsules. The capsules themselves were enclosed in a porcelain chamber evacuated with a mechanical pump during the homogenization anneals. Upon completion of homogenization, the alloys were cooled to room temperature and then reheated to 550°C for 3 1/2 days in order to transform the alloys to a finely divided mixture of a and TiCr2, and thus provide an approximately constant initial mi-crostructure. Individual specimens were sealed in evacuated and internally gettered Vycor capsules and heat treated as previously described.8,9 EXPERIMENTAL RESULTS Ti-10.64 Pet Cr. This investigation was largely performed on the 0.64 pet Cr alloy. The "base"

Jan 1, 1962

By H. Brandhorst, C. G. Dunn

THE main purpose of the present note is to provide further information on the effect of orientation on strain-induced grain boundary migration in sheet material. A secondary purpose is to draw attention to the presence of a possible (hkl) dependent surface energy driving force, when the two grains of a bicrystal specimen have different (hk2) planes at the gas-metal interface, and to suggest means of eliminating this effect from the strain effect. If two grains A and B of a bicrystal have specific surface energies ? A and ? B respectively, there exists a driving force for grain boundary migration equal to 2(? B - ? A)/t, where t is the thickness of the specimen.' When t is quite small and the residual strain energy densities €A and €B of grains A and B, respectively, adjacent to the boundary are nearly equal, it is possible that the surface effect may be larger than the strain effect. he difference € B - € A gives a driving force per unit area of grain boundary due to an applied strain.)2 However, with increasing t the surface effect decreases while the strain effect presumably remains constant. Therefore, at some sufficiently large value of t, it is possible to ignore surface energy driving forces.* On the other hand, the surface effect is eliminated for all value of t when both grains have the same (hkl) orientation or have identical values of surface energy. Although the previous study provided information on the effect of orientation on strain-induced erain boundary migration,' no orientations with the (100) plane in the plane of the sheet were included. The present paper provides information of this kind. A seed crystal reorientation technique4 was used to prepare a large bicrystal of silicon iron (3 1/4 pct Si) consisting of two (100) grains in the (100) [001] and (100) [Oil] orientations. After a cold rolling strain (7 pct reduction to 3/4 mm final thickness), small specimens were prepared as before2 and then annealed in hydrogen. The boundary migration was such that the grain in the (100) [011] orientation always grew. After overnight anneals at 1100° or 1000°C, none of the initial (100) [001] grain remained. An anneal of 30 min at 950°c, however, provided information on the start of boundary migration. Fig. 1 shows both the initial position (lower line) and the final position (upper curved line) of part of the boundary, as revealed by thermal grooves. In a different area, part of the final boundary was not revealed by thermal grooves; the rate of boundary migration presumably was too high for the growth front to be defined by the time dependent thermal grooving process. In another area there was a vacated intermediate position connected to the initial position of the grain boundary. This indicates that only part of the boundary

Jan 1, 1962

By O. Jantti, Erkki Laurila, R. T. Hukki

INVESTIGATION of the methods of analyses for magnetite and ilmenite in the Otanmaki iron-titanium ore and respective mill products has resulted in certain improvements in the methods conventionally employed. The outcome of the investigations is summarized as follows: 1—Development of a diamond-drill core analyser for the determination of magnetic susceptibility of drill cores containing magnetic minerals. 2—Development of a magnetite-ilmenite analyser for pulverized products. 3—Development of a potentiometric method of chemical analysis for the determination of iron, vanadium, and titanium in succeeding steps from one test sample. 4—Proposal for a continuous and automatic magnetite assay method. Diamond-Drill Core Analyser The conventional method of assaying diamond-drill cores consists of splitting the cores, pulverizing them and running chemical analyses on representative samples. The physical properties of magnetic minerals, however, offer an entirely different approach to obtain similar final result by the measurement of respective magnetic susceptibilities in a properly designed instrument. This instrument is shown in Fig. 1. Essentially it consists of a magnetic unit, A, and of a meter unit, B. On top of A is a slot where a piece of diamond-drill core C is visible. The magnitude of disturbance in the original magnetic field caused by the magnetic minerals in the drill core can be observed from the microammeter. This reading is a function of the combined susceptibilities of all magnetic grains in the magnetic field. The magnetic unit is built of E-sheets of a shell transformer, as shown in Fig. 2. The primary coil and two identical secondary coils are wound around the central leg and end legs, respectively. The drill core under examination is placed into one of the air gaps as shown in Fig. 2. The number of ampere-turns of primary coil excited by 50 or 60-cycle ac give a maximum induction of about 600 gauss in the air gaps. The secondary coils are connected in series in the opposite direction. Thus their voltages compensate each other if no drill core exists in the air gap or if the susceptibility of the drill core is equal to zero. As soon as the susceptibility becomes greater than zero, the balance of the magnetic circuit is disturbed, resulting in increased voltage in the terminals of coil 1 and in a respective decrease of voltage in the terminals of coil 2. The voltage difference thus produced is a function of the effective susceptibility of the drill core. If the magnetic grains in the drill core are unevenly distributed, the voltage difference varies as the core is rotated around its own axis. The voltage difference is measured with a Bren-tan0 tube voltmeter. If extreme sensitiveness is desired, the established symmetry of the magnetic circuit does not allow reduction of the 0-voltage difference far enough. Therefore a compensatory voltage of a correct phase taken from the terminals

Jan 1, 1952

By L. F. Coffin

Monotonic tension and controlled cyclic-strain tests have been conducted on an annealed steel at various strain rates and temperatures in the weep range. Evidence of a strain-induced precipitate reaction abozle 500°C was determined from pronounced cyclic-hardening effects, increase inflow stress, inter granular fracture, drastic decrease in low-cycle fatigue resistance, and loss in tensile fracture ductility, as well as through electron-microscopy examination. It was found that strain rate as well as temperature had a pronounced effect on these observations. A parameter of the form P = T(7.95 - log ip) was determined, which was directly related to the ductility, defined either as the true strain at fracture in the case of simple tension, or as 2N 1/2 ? E p in the case of cyclic plastic strain. A highly satisfactory comparison of the cyclic life as measured and as computed from the ductility predicted by this parameter was obtained for a range of temperatures from 400° to 650°C and for crosshead speeds of 0.002 to 0.2 in. per min. MANY metals subjected to sustained loading at elevated temperatures derive resistance to creep from solid-state reactions of various alloying constituents. These reactions can be many and varied, and occur as a consequence of time and temperature according to reaction-rate principles. Strong evidence exists that many of these reactions are influenced by plastic strain, the argument being put forth that dislocations act as sites for nucleation and growth for the new phase.' Many examples can be cited where plastic strain during transformation influences the structure, and as a consequence affects the properties.2 The low-alloy steels serve as particularly good examples for investigating the effect which aging reactions have on resulting properties. Alloying elements such as carbon and nitrogen, and carbide formers such as manganese, chromium, molybdenum, and so forth, can, by the particular aging reaction which occurs, produce strengthening effects at various temperatures, and so endow strength to the alloy over a broad range of temperatures. These effects have been systematically studied by Glen both during creep3 and for constant strain-rate tension testing.4,5 The work of Glen is of considerable interest in connection with the present study, particularly in regard to the effect of various alloying additions on the true stress-strain curve and on fracture ductility at various elevated temperatures. He has shown that strain-age hardening occurred during plastic straining, which was quite specific both to the temperature of the test and to the alloying element involved. When strain-age hardening was noted at a particular temperature, a .minimum in reduction of area was also obtained. Provided the temperature for the minimum reduction of area exceeded 300°C, all minimum ductility fractures were found to be inter-crystalline in nature. For carbon steels containing molybdenum, two minima in reduction of area were noted, one occurring at 200° to 250°C, the precise temperature depending on the percentage of molyb-

Jan 1, 1964

By Walter R. Hibbard, Cyril Stanley Smith

IN 1929 one of the authors1 determined the constitution of copper-silicon-manganese alloys contailling over go per cent copper. Through a combination of circumstances the presence of the copper-silicon kappa phase was overlooked in both binary and ternary systems. In the binary system the omission has already been rectified.2,3 To establish the limits of this phase in the ternary system, the samples that had been heat-treated at the time of the original research have now been re-examined, using more careful technique. Full details of the compositions used and the temperatures and times of heat-treatment are to be found in the earlier paper.2 The etching reagent now used to distinguish the kappa phase had the following composition: hydrogen peroxide (30 vol.), 2o c.c.; water, 25 c.c.; 20 per cent potassium hydroxide solution, 5 c.c.; ammonium hydroxide (0.90 sp. gr.), 50 c.c. 'I'able I lists the newly observed structures. In virtually every case the constituents are the same as those previously recorded, except where kappa is now shown to be present as a result of the improved etching technique. X few other changes will be noted (No. 138 at 450°C., No. 147 at 550°C., Nos. 153 and 149 at 75o°C.). The change from liquid to beta in equilibrium with kappa in No. 137 at 800°C. is based upon newly heat-treated samples. The alloy is very near the solidus at 800°C. and probably there was a slight error in the old heat-treatment temperature. The isothermal phase diagrams shown in Figs. I to 6 were reconstructed on the basis of the new studies. The presence of manganese does not greatly change the composition limits of the kappa phase or cause the K-a + Y change to occur at a markedly different temperature. There is some uncertainty regarding the range of existence of beta, for alloys were not spaced very closely in composition. The proposed diagrams, however, are in accordance with the available experimental data and with the dictates of the phase rule. Theoretically necessary fields with no experimental points are delineated by dotted lines. The alpha boundaries at all temperatures are shown together in Fig. 7. This includes also some additional points on the alpha boundary at 300' and 350°C., which were obtained as a result of electrical-conductivity measurements. Two series of seven alloys with increasing manganese content but with silicon constant at 1.90 + 0.03 per cent or 2.42 ±. 0.05 per cent were cast in ingots of 3-in: diameter. These were hot-rolled, homogenized and drawn to wire, finished at 0.08-in. diameter with a 61 per cent reduction of area following annealing at 700°C. Conductivity determinations were made before and after annealing for 30 days at 300°C. or 21 days at 350°C. Figs. 8 and 9 show the change of conductivity produced by this annealing, plotted against manganese

Jan 1, 1942

By Walter R. Hibbard, Cyril Stanley Smith

IN 1929 one of the authors1 determined the constitution of copper-silicon-manganese alloys contailling over go per cent copper. Through a combination of circumstances the presence of the copper-silicon kappa phase was overlooked in both binary and ternary systems. In the binary system the omission has already been rectified.2,3 To establish the limits of this phase in the ternary system, the samples that had been heat-treated at the time of the original research have now been re-examined, using more careful technique. Full details of the compositions used and the temperatures and times of heat-treatment are to be found in the earlier paper.2 The etching reagent now used to distinguish the kappa phase had the following composition: hydrogen peroxide (30 vol.), 2o c.c.; water, 25 c.c.; 20 per cent potassium hydroxide solution, 5 c.c.; ammonium hydroxide (0.90 sp. gr.), 50 c.c. 'I'able I lists the newly observed structures. In virtually every case the constituents are the same as those previously recorded, except where kappa is now shown to be present as a result of the improved etching technique. X few other changes will be noted (No. 138 at 450°C., No. 147 at 550°C., Nos. 153 and 149 at 75o°C.). The change from liquid to beta in equilibrium with kappa in No. 137 at 800°C. is based upon newly heat-treated samples. The alloy is very near the solidus at 800°C. and probably there was a slight error in the old heat-treatment temperature. The isothermal phase diagrams shown in Figs. I to 6 were reconstructed on the basis of the new studies. The presence of manganese does not greatly change the composition limits of the kappa phase or cause the K-a + Y change to occur at a markedly different temperature. There is some uncertainty regarding the range of existence of beta, for alloys were not spaced very closely in composition. The proposed diagrams, however, are in accordance with the available experimental data and with the dictates of the phase rule. Theoretically necessary fields with no experimental points are delineated by dotted lines. The alpha boundaries at all temperatures are shown together in Fig. 7. This includes also some additional points on the alpha boundary at 300' and 350°C., which were obtained as a result of electrical-conductivity measurements. Two series of seven alloys with increasing manganese content but with silicon constant at 1.90 + 0.03 per cent or 2.42 ±. 0.05 per cent were cast in ingots of 3-in: diameter. These were hot-rolled, homogenized and drawn to wire, finished at 0.08-in. diameter with a 61 per cent reduction of area following annealing at 700°C. Conductivity determinations were made before and after annealing for 30 days at 300°C. or 21 days at 350°C. Figs. 8 and 9 show the change of conductivity produced by this annealing, plotted against manganese

Jan 1, 1942

By H. L. Talbot, A. J. Hall, R. S. Young

A method for the chemical differentiation of oxidized cobalt from sulphide cobalt in complex mine and concentrator products is of considerable interest to a number of mining companies throughout the world. At Rhokana Corporation, in Northern Rhodesia, for instance, cobalt is found chiefly in the sulphide form, as carrollite, but a cobalt oxide mineral also occurs irregularly underground, and relatively small quantities in the mill feed from time to time are difficult to detect visually at the concentrator. The oxide cobalt is not amenable to the sulphide flotation process used at Rhokana and nearly all of it passes into the tailing. The occasional period when high tailing losses occur with no appreciable change in the cobalt content of the mill feed has naturally been disturbing to the concentrator staff, who had no means of determining whether an influx of cobalt oxide ore or laxity on the part of the mill flotation operators as responsible for these results. A procedure for determining oxidized copper minerals, based on their solubility in dilute sulphuric acid saturated with SO2, has been used for many years in the copper industry. In the presence of cuprite the procedure fails, but otherwise it gives results quite satisfactory for control pur- poses. Empirical methods for determining oxidized forms of zinc and lead have been proposed, but they are generally unsatisfactory except in isolated cases of a very simple ore. No published procedure for determining oxidized forms of cobalt or siinilar elements has come to the authors'notice. Experimental Procedure There are obviously two alternatives in attacking a problem of this nature. One is to employ a solvent for oxide cobalt that will not dissolve cobalt sulphide. The other is to find a substance capable of oxidizing the cobalt sulphide to sulphate, which of course is water-soluble and can then be leached from the water-insoluble oxide cobalt. A representative sample of Rhokana cobalt oxide ore was secured, ground to minus 200-mesh, and carefully rolled to furnish a standard sample for this work. The analysis of this material gave the following percentages: Co, 0.55; Cu, 1.64; Fe, 2.98; S, 0.08; SiO2, 41.42; Al2o3, 11.20; CaO, 13.82; MgO, 3.80. As a source of cobalt sulphide a similar large sample of carrollite was prepared, analyzing Co, 23.36 per cent; Cu, 27.32; and S, 31.13. This sample of course was not pure carrollite, but was the best obtainable in quantities sufficient for a standard intended for considerable work. It was a pure sulphide sample, microscopic examination showing only carrollite, a little bornite, and gangue, with no visible sign of oxidized ore present. A wide variety of reagents as used in an effort to selectively dissolve oxide cobalt or to oxidize sulphide cobalt to

Jan 1, 1949

By Y. L. Yao

THE NaT1-type compound may be considered as the penetration of two diamond lattices in such a way that a superstructure of the bcc lattice is formed. Examples of the NaT1-type compound of lithium are LiZn, LiCd, LiA1, LiIn, and LiGa. Except for LiZn and LiCd, each molecule of these equiatomic compounds contains a total of four valence electrons. In LiZn and LiCd it was suspected (for example, see Ref. 1) that the fourth electron was withdrawn from the 18-shell of zinc or cadmium. An incomplete 18-shell would have to show strong paramagnetism. In 1955 Klemm and Fricke2 determined the magnetic susceptibility of LiZn, LiCd, LiA1, and LiIn (also NaTl and NaIn). The molar susceptibility of the outer electrons of LiZn and LiCd was found to be not strongly but weakly positive; that of LiAl and LiIn not negative but weakly positive. Thus they questioned the hypothesis about LiZn and LiCd in particular and the description of the NaT1-type structure in general. However, they did not determine the magnetic susceptibility of LiGa. This note is to report results for all five compounds of lithium. The starting materials are 99.9 pct pure lithium sticks containing less than 0.005 pct Na. 99.9 pct Zn, 99.9 pct Cd, 99.99 pct Al, 99.999 pct In, and 99.999 pct Ga. Weighed amounts of two components were wrapped in tantalum foil and placed in a graphite crucible. It was found that the tantalum sheath minimized the vaporization losses of lithium during melting. The sheath was easily stripped from the alloy rod and the tantalum in the alloy was barely detectable (less than 0.01 wt pct). The melting was performed in a resistance furnace under an argon atmosphere. Alloy specimens were annealed for 1 hr at 400°C in vacuum. The magnetic susceptibility at room temperature was determined by the Faraday Method, using the procedures previously described.' The results are shown in Table I. The interpolated values at 50 at. pct of Li and the average values at 48 to 52 at. pct of Li at 20°C determined by Klemm and Fricke2 are also shown. It is seen that the present value for LiZn and LiCd is slightly more positive but the present value for LiAl and LiIn is slightly less positive than their corresponding values. Also, it is noted there is an abrupt decrease of susceptibility of LiAl from 51.2 to 52.5 at. pct of Li. Klemm and Fricke2 found that the total susceptibility of LiZn, LiCd, LiA1, and LiIn is almost independent of temperature. slater4 gave a working expression for the diamagnetism of the atomic core as follows:

Jan 1, 1964

By D. T. Peterson, T. Carnahan

The diffusion coejTicients of carbon, nitrogen, and oxyget were determined in $ thorium over the tempernilcre range 1440" io 1715°C. The diffusion coyfiicir?zls are given by: D = 0.022 exp (-27,000/RT) jor carbo)~, D = 0,0032 exp(-l7,00Q/RTj for nitrogen, and D = u.0013 expt(-11,UOU/RT) for oxygen. Cavl~orz was found to increase the hardness of thoriunz nearly linearly with concentration over the range 100 to 1000Ppm carbon. ThORIUM has a fcc structure up to 1365°C and a bcc structure from this temperature to its melting point at 1740°C. Diffusion of carbon, oxygen, and nitrogen in bcc thorium was of interest in connection with the purification of thorium by electrotransport.' In addition, it was possible to measure the diffusion of all three of these interstitial solutes in the same bcc metal. Only in niobium, tantalum, vanadium, and a iron have all three interstitial diffusion coefficients been measured in a given bcc metal. Diffusion coefficients have been measured for carbon and oxygen in a thorium by Peterson2, 3 and for nitrogen by Gerds and Mallett.4 Activation energies for diffusion are reported by the above authors to be 38 kcal per mole for carbon, 22.5 kcal per mole for nitrogen, and 49 kcal per mole for oxygen. Values of the diffusion coefficients of carbon and nitrogen in 3 thorium have been reported by Peterson et al.' However, these were secondary results of their investigation of electrotransport phenomena in thorium and it was hoped that the present study could provide more precise data. EXPERIMENTAL PROCEDURE The specimens used in this study were the well-known pair of semi-infinite bar type. The couple was formed by resistance butt welding two 0.54-cm-diam by 3.0-cm-long bars of thorium together under pure helium, the concentration of the solute being greater in one cylinder than that in the other. The finished couple then contained a concentration step at the weld interface and diffusion proceeded only along the axis of the rod. The thorium used in this study was prepared by the magnesium intermediate alloy method.5 The total impurity content was less than 400 ppm. The major impurities were: carbon, 100 ppm: nitrogen, 50 ppm; and oxygen. 85 ppm. The total metallic impurity content was less than 150 ppm. The high solute concentration portions of the diffusion couples were prepared by adding the solute to the high-purity thorium in a non-consumable electrode arc melting procedure. Carbon and nitrogen were added in the form of spectroscopic graphite and nitrogen gas while a Tho2 layer was dissolved by arc melting to add oxygen. High-purity thorium formed the low concentration portions in the carbon and nitrogen couples. The low oxygen portions were obtained by deoxidizing high-purity thorium with calcium for 3 weeks at 1000°C according to a method reported by Peterson.3 The high C-Th contained 400 ppm C, the high N-Th contained 400 ppm N, the high 0-Th contained 220 ppm 0, and the low 0-Th contained 25 ppm O. The high O-Th was brine-quenched from 1500°C to retain most of the oxygen in solution at room temperature. These concentration levels were all below the solubility limits in 0 thorium at 1400°C. A resistance-heated high-vacuum furnace was used to heat the couples. The samples were mounted horizontally on a tantalum support which had small grooves near each end. Spacer rods of thorium, 0.4 cm in diam, were placed in these grooves to prevent contact between the sample and the tantalum support. This arrangement should have prevented contamination of the sample by contact with the support. In further effort to reduce contamination, the oxygen diffusion couples were sealed inside evacuated outgassed tantalum cylinders lined with thorium foil. Thorium rings around each end of the samples acted as spacers in this case. Pressure during diffusion runs was about 10-6 torr after an initial outgassing stage. Temperature measurements were made by sighting on black body holes in the sample support adjacent to the samples with a Leeds and Northrup disappearing-filament optical pyrometer. Temperatures were constant during a diffusion anneal to ±5C. The observed temperatures were corrected for sight glass absorption after each diffusion run. The pyrometer was checked against a calibrated electronic optical pyrometer and a calibrated tungsten strip lamp with the electronic pyrometer being taken as the standard. All temperature readings agreed to within ±3C over the temperature range 1450" to 1690°C. Time corrections due to diffusion during heating and cooling were necessary because of the short diffusion times. The diffusion times ranged from 6 min for the oxygen sample run at 1690°C to 90 min for the carbon sample run at 1500°C. A series of temperature vs time plots were made for heating and cooling of the samples to the various diffusion temperatures. This data was then used in a method according to shewmon6 to determine the time corrections. The corrections amounted to

Jan 1, 1970