Search Documents

By G. P. Marino

The "conductivity ", Kr, of a spherical void in a linear temperature gradient due to radiation is calculated by integration of the general expression for diffuse radiation between two gray surfaces of arbitrary orientation. The spherical geometry allows an exact evaluation of the phenomenon without resorting to up-proximate expressions for the radiant heat flux. The resulting expression for Kr is three quarters of the approximate value obtained by Loeb. IN the determination of the thermal conductivity of porous materials such as sintered UO2 and ThO2 fuel pellets, the magnitude of the thermal conductivity in the pores due to radiative heat transfer across the pore, Kr, may be important in assessing the total conductivity of the pellets. Loeb' has calculated Kr for various pore geometries by approximating the net flux of heat between two surfaces of temperature difference AT as 4eT~AT and obtained: K,=4yoet3odp [1] where a is the Stefan-Boltzmann constant, E is the emissivity of the radiating surface, To is the average absolute temperature of the region across which radiation occurs, dp is the largest dimension of the gap in the direction of heat flow, and y is a geometric factor. That is: y = 1 for laminar pores and cylindrical pores with axes parallel to the heat flow direction, y = 2/3 for spherical pores, y = /4 for cylindrical pores with axes perpendicular to heat flow direction. It was the purpose of this work to obtain an exact expression for K, in spherical pores by integration of the general expression for the net heat transfer between two incremental areas on the surface of a spherical pore. This procedure yields a value of Kr 3/4 of that obtained by Loeb. I) CALCULATION OF Kr Fig. 1 illustrates the general geometry which defines the problem. The total emissive power for a gray body is given by: E = EoT4 [2] where E is the radiated energy per unit area of the radiation surface per unit time, T is the absolute temperature, o is the Stefan-Boltzmann constant (o = 5.67 X lo- '' watt per sq cm-"K1), and E is the emissivity of the surface. Although Eq. [2] gives the exact expression for the total emissive power of any gray body, all practical problems involve the net exchange of thermal radiation between a body and surrounding bodies. The two major considerations in the solution of such problems are a) shape and relative position of the surfaces and b) the emissivities of the surfaces. In this work it shall be assumed that all incremental areas on a pore have the same emissivity and, further, that little or no reflection occurs (i.e., the radiation hitting a surface is either absorbed or transmitted). Thus, the solution of the stated problem requires only that consideration a be satisfied. Consider two incremental surface areas dA 1 and dA2 of the same material in a general orientation to one another. Let dAl be at temperature TI and dAz be at T2. Fig. 2 illustrates the geometrical parameters involved. That is, 71, is the normal of dA, n2 the normal of dA2, r the line connecting the centers of each area, 81 the angle between nl and r, and Q2 the angle between 7,2 and r. Using Eq. [2] and the Lambert cosine law of diffuse radiation it is easily shown that:293 WdAf- dA, = ecr T*------^fpi------" dAidA2 [ 3a 1 and

Jan 1, 1970

By H. K. Najarian

This paper describes recent development work on the Wheaton-Najarian Vacuum Condenser for Zinc at Josephtown Zinc Smelter, wherein, to dissipate the heat of condensation, the "airlift" action of gases bubbling through the molten zinc in the condenser is utilized to circulate the zinc through an outside "cooling well" and cooling it continuously. AN earlier paper1 gave an account of initial development of the Wheaton-Najarian Vacuum Zinc Condenser at the Josephtown Zinc Smelter of the author's company. This paper is a progress report of further development work over the years directed toward improving the efficiency and increasing the capacity of the condenser. In this type of zinc condenser, zinc vapor with accompanying noncondensable gases from electro-thermic zinc smelting furnaces is passed upwardly through a mass of molten zinc confined in a refractory-lined U-tube. The zinc vapor condenses in the molten zinc and the noncondensable gases are drawn out of the condenser by suction applied over the surface of the molten zinc. The heat of condensation of zinc, amounting roughly to 1000 Btu per lb of zinc condensed, is absorbed by the molten zinc and, in turn, is dissipated continuously to keep the temperature of molten zinc in a range to permit effective condensation. At the time of publication of our first paper, we had at Josephtown Smelter five units of 69-in.-bore metal furnaces and one larger metal furnace having 96-in. bore. Each 69-in. metal furnace had a condenser with 22-ft overall length. The 96-in.-bore furnace had a condenser with overall length of 31 ft. Condensers were of circular cross section, with corrugated steel outer shell and carborundum lining. Diameter of these condensers was 32 in. inside the refractory lining. Molten metal in the con- denser was kept cool by spraying cold water on the outside of the condenser shell. Under good operating conditions, 22-ft condensers produced some 12 tons of metal per day; and the 31-ft condenser, about 22 tons. When the furnaces produced zinc vapor at higher rates, invariably the condenser temperature would rise to levels where excessive blue powder would form and accumulation of accretions in the condenser would be accelerated. We found that, with the system of cooling in use, we had to provide heat-exchange surface of some 150 sq ft to condense 1000 lb of metal per hr with reasonable efficiency. The wall of the condensers was made up of 5/8-in.-thick corrugated-steel shell having next to it several layers of mica and inner refractory lining of special carborundum shapes to fit the corrugation in the steel shell, and having a minimum thickness of about 2 in. Relatively low rate of heat transfer through the thick condenser wall and slow internal circulation of the molten zinc limited the rate of heat exchange between hot metal within the condenser and the cold water flowing over the outer steel shell. The thickness of the condenser wall was about minimum that we considered suitable for continuous operation. As we learned to operate the electrothermic furnaces with higher power input, producing zinc vapor at higher rates, it became apparent that we needed condensers having larger heat-exchange surfaces. Our 96-in.-bore furnace had 6 single-phase power circuits with transformers rated at 500 kva each and normally good for power input of 3000 kw.

Jan 1, 1959

By Harvey T. Kennedy, Charles R. Olson

In many investigations of the performance of petroleum reservoirs the assumption is made that the liquid, if below its bubble-point pressure, is at all times in equilibrium with gas. On the other hand, observations by numerous investigators have indicated that gas-liquid systems including hydrocarbon systems, may exhibit supersaturation to the extent of many hundred psi in the laboratory. Up to the present, there has been no reliable data on which ti judge the actual extent of supersaturation under conditions approaching those existing in petroleum reservoirs. The work reported here deals with observations and measurements on mixtures of methane and kerosene in the presence of silica and calcite crystals. Bubbles were observed to form on crystal-hydrocarbon surfaces in preference to the glass-hydrocarbon interface or to the body of the liquid. Statistically, it was found that the number of bubbles formed per second per square centimeter of crystal surface was a function of the supersaturation only. and the function was evaluated graphically. Supersaturations were observed up to 770 psi, under which condition bubbles formed quickly and with considerable violence. With decreasing degrees of supersaturation, the frequency of bubble formation became less, until at 30 psi supersaturation and lower, no bubbles were observed to form, even though the observation at 30 psi was continued for 138 hours. It was found that silica and calcite crystals had identical effects, within experimental error, in accelerating the formation of bubbles, and that small amounts of water and crude oil had no effect on the results. It is shown that the maximum supersaturation that can exist in a reservoir may be calculated from the data presented and from the area of the rock surface. It is also shown that the number of bubbles formed in the reservoir, in order of magni- tude, may be calculated for any rate of pressure decline imposed on the reservoir by production. The bearing of the number and distribution of bubbles on reservoir performance is discussed. INTRODUCTION A liquid system is supersaturated with gas when the amount of gas dissolved exceeds that corresponding to equilibrium at the existing pressure and temperature. The degree of super-saturation may be conveniently expressed as the difference between the bubble-point of the mixture and the prevailing pressure. Thus, if a mixture having a bubble-point of 1,000 psi at a given temperature exists in single liquid phase at 700 psi at the same temperature. it is supersaturated to the extent of 300 psi. There are many examples of high supersaturations, mostly in aqueous solutions, reported in the literature. Thus, Kenrick, Wismer and Wyatt1 showed that water may be saturated with oxygen, nitrogen or carbon dioxide at 100 atmospheres, and the pressure reduced to one atmosphere without producing bubbles immediately. When liquids are in a state of tension, they may be considered as supersaturated at least to the extent of the tension. The tensile strength of water has been reported as 30 atmospheres by Meyer,1 60 atmospheres bv Budgett. 30 to 50 atmospheres by Temperley and Chambers,4,5 200 atmospheres by Dixon6. and 223 atmospheres by Briggs.7 Vincent8,9 determined the tensile strength of a mineral oil as 45 psi. Gardescu1 maintained pressures for short times in a model reservoir at 115 psi below the bubble-point. It should be noted that the high supersaturations observed were obtained on systems carefully purified to remove particles or surfaces which might promote the formation of bubble; These "nuclei" were considered as contaminants which interfered with the determination of a property of the liquid. In

Jan 1, 1952

By M. W. Legatski, D. L. Katz

The best currently available description of the longitudinal mixing properties of a porous medium is an equation of the form which relates the effective longitudinal dispersion coefficient Dp to the molecular diffusion coefficient Do, the electrical resistivity factor F, the porosity <f> and a Peclet number. If the parameters dpa and m are determined for a porous medium of known porosity and electrical resistivity factor, then a dispersion coefficient may be estimated for a given flow rate and a given gas pair. A new method, featuring on-line gas analysis by thermal conductivity and on-line data reduction by analog computation, was developed and used to determine these mixing parameters for eight naturally occurring sandstones and two dolomite samples. The exponent m of the above equation was found to vary between 1.0 and 1.5. The characteristic length dpa in the above equation was found to vary between 0.25 and 1.9 cm, with an average value of 0.4 cm for sandstones. Measurements were made on two cores in which paraffin wax had been deposited by evaporation from a pentane solution. They indicated that the presence of an immobile phase such as connate water could increase the dispersion coefficients significantly. INTRODUCTION While the petroleum and chemical industries have studied the mixing of miscible liquids flowing in consolidated porous media and of miscible gases flowing in unconsolidated porous media, relatively little data have been presented to describe the mixing of gases flowing through consolidated porous media. Such data are of particular interest to the gas storage industry. For instance, the U. S. Bureau of Mines is storing large quantities of a rich helium-nitrogen gas in contact with a natural gas in a dolomite reservoir. Since the rich gas occupies only 15 percent of the total reservoir volume, it is essential that the extent of rich gas-natural gas mixing be predicted and understood as a function of rock properties, pressure and rate of movement. This investigation was concerned only with the determination of longitudinal dispersion coefficients. It is understood that a transverse dispersion coefficient, which characterizes mixing perpendicular to the direction of flow, may be an order of magnitude less than the coefficient characterizing mixing in the direction of bulk flow.5719 It should also be recognized that the use of any dispersion coefficient is in itself a simplification. It is necessary to assume that mixing in a porous medium may be characterized by the equation for flow in a single direction. A number of authors1 l8,19 have pursued the mixing problem, not in terms of the so-called "dispersion model" described by Eq. 1, but in terms of a "mixing cell model". This model supposes that a porous medium is constructed of a large number of small mixing chambers and that the concentration of the diffusing component within each mixing chamber is uniform. Fick&apos;s law (Eq. 1) assumes that there is no gross by-passing of one fluid by another, and that there are not stagnant pockets of gas in the system under consideration as discussed by Coats and Smith.8 These assumptions are not always valid for flow through porous media and it is important to recognize the limitations upon Eq. 1. The most serious limitation upon the assumption of Fick&apos;s law mixing in a porous medium is

By W. B. Gogarty, H. R. Bailey

Methods are presented for determining molecular diffusion coefficients by using data from, capillary flow experiments. These methods are based on a numerical solution (presented in a previous paper) of the partial differential equation describing the combined mechanisms of flow and diffusion. Results from this numerical solution are given and compared with the approximate analytical solution of G. I. Taylor. The numerical solution is valid over a much larger time range. These methods are applied to experimental results for the fluid pairs water-potassium permanganate solution and amyl acetate-orthoxylene. Both of these fluid pairs have approximately equal densities and viscosities. Graphical and numerical techniques are presented for deter mining diffusion coefficients from the flow data. Values obtained by these techniques are compared with values obtained by other methods. INTRODUCTION The molecular diffusion coefficient is known to be a variable in determining the amount of mixing in a miscible displacement process. The effect of molecular diffusion on dispersion in longitudinal flow through porous media -has been examined by different investigators. These investigators concluded that at low velocities of flow, the amount of dispersion is approximate1y proportional to the molecular diffusion coefficient. The influence of diffusion on fingering, channeling, and overriding has been mentioned by other investigators.4 Recent studies5,6 have been made on the effects of molecular diffusion in connection with the problem of gravity segregation. Many different methods have been developed for the experimental determination of molecular diffusion coefficients? These methods differ mainly according to boundary conditions selected and analytical procedures used. Nevertheless, all of these methods have the condition in common that the bulk fluids in which diffusion is occurring are stationary with respect to each other. In connection with a series of papers on mixing in capillary flow,8-10 Taylor suggested the use of a flow method for determining molecular diffusion coefficients. Additional studies 11-13 have been conducted on miscible displacements in capillary tubes, but the data from these studies were not used for the specific purpose of determining diffusion coefficients. The flow method proposed by Taylor results in a single value of the diffusion coefficient for the fluid pair used in the displacement experiments. This single value represents the true value for the fluid pair when the diffusion coefficient is independent of concentration. If the diffusion coefficient is a function of concentration, the single value obtained by the flow method gives an average value for the coefficient of the fluid pair. These average values are based on diffusion taking place over the entire range of concentration, i.e., from 0 per cent of one fluid to 100 per cent of that same fluid. In field applications of the miscible displacement process, gradients occur over the same range of concentration as are found in the displacements in capillary tubes. Molecular diffusion coefficients obtained from the capillary flow method should, therefore, be especially relevant to field operations. This investigation was undertaken to evaluate the feasibility of obtaining molecular diffusion coefficients from capillary -flow experiments. In making this evaluation, diffusion coefficients were first determined for two systems from data obtained in capillary flow experiments. These values of the diffusion coefficient were then compared to values obtained by other methods. MIXING IN CAPILLARY FLOW—THEORETICAL The theoretical basis for determining molecular diffusion coefficients from capillary flow experiments is the partial differential equation relating the mechanisms of flow and diffusion. For the case of laminar flow and with a fluid pair having equal densities and viscosities, the equation is given by

By M. Bevis, A. G. Crocker

A. G. Crocker and M. Bevis (Battersea College of Technology)— Richman and Richman and conard17 have recently published the very interesting results of their studies on the plastic deformation of Fe-Ni-C martensites and of an Fe-5 wt pct Be alloy. Other work on the range of Fe-Be alloys containing between 16 and 27 at. pct Be has also been reported by Richman." One of the fascinating features of this work is that the materials investigated deform primarily by means of twinning and that several twinning modes are operative in addition to the well-established (112 ) mode of bcc materials The composition planes of these unusual twins were determined using a single-surface analysis and reported as (0131, "{045)" or "{089)", and "{2,3,10)" or "(127)". (It is now thought1&apos; that the "{0,1,13)" habit reported for the Fe-Ni-C alloys&apos;&apos; is more likely to be near "{127)". Although the full twinning elements of these unusual twinning modes could not be determined directly from the experimental results, the authors attempt to deduce the elements by considering twin intersectionslQ and by using the projection technique of Jaswon and Dove.&apos; The basic assumption of their approach is that the (013) habits are truly rational planes but that the other habits are irrational and merely approximated to by the indices quoted above. They thus deduce that the twinning mode describing the (013) habit is either Type I or compound and that the modes describing the other habits are Type 11. This deduction must however be false as all twinning modes in cubic structures obeying the classical orientation relationships of deformation twinning must be compound. Type-I and Type-I1 modes can occur in all of the other crystal systems (indeed in the triclinic system compound modes are not allowed), but in the cubic system only compound modes can arise. This fact, which does not appear to be well-known, fol- lows immediately from the relations between the four twinning elements defining a twinning mode. Thus for a Type-I mode which is uniquely defined by the rational twinning elements K1 = hi and tj2 = vi (i = 1, 2, 3), the Kz qd 771 elements are given by: where the summation convention of the tensor calculus has been adopted. For cubic materials the metrics cij Iand cij reduce to the unit tensors 6ij and 6ij so that we have It is clear from these equations that both ki and ui must take on integral values; i.e., Kz, and ql must be rational, if hi and vi are integers, i.e., if K1 and 77~ are rational. Thus for cubic materials a Type-I mode must degenerate into a compound mode and the same result also arises for Type-II modes. The apparently irrational habits observed by the authors can thus not be explained in terms of Type-I1 twinning if the classical laws of deformation twinning are to be obeyed. An alternative explanation of the irrational habits is that they may arise from twinning modes with four irrational twinning elements. These modes, which do not obey the classical orientation relationships associated with deformation twinning, have been shown to arise as the result of double twinning processes in magnesium,22 but may also describe a simple shear on an irrational plane in an irrational direction. Indeed double twinning processes in bcc

Jan 1, 1963

By W. B. Gogarty, H. R. Bailey

Methods are presented for determining molecular diffusion coefficients by using data from, capillary flow experiments. These methods are based on a numerical solution (presented in a previous paper) of the partial differential equation describing the combined mechanisms of flow and diffusion. Results from this numerical solution are given and compared with the approximate analytical solution of G. I. Taylor. The numerical solution is valid over a much larger time range. These methods are applied to experimental results for the fluid pairs water-potassium permanganate solution and amyl acetate-orthoxylene. Both of these fluid pairs have approximately equal densities and viscosities. Graphical and numerical techniques are presented for deter mining diffusion coefficients from the flow data. Values obtained by these techniques are compared with values obtained by other methods. INTRODUCTION The molecular diffusion coefficient is known to be a variable in determining the amount of mixing in a miscible displacement process. The effect of molecular diffusion on dispersion in longitudinal flow through porous media -has been examined by different investigators. These investigators concluded that at low velocities of flow, the amount of dispersion is approximate1y proportional to the molecular diffusion coefficient. The influence of diffusion on fingering, channeling, and overriding has been mentioned by other investigators.4 Recent studies5,6 have been made on the effects of molecular diffusion in connection with the problem of gravity segregation. Many different methods have been developed for the experimental determination of molecular diffusion coefficients? These methods differ mainly according to boundary conditions selected and analytical procedures used. Nevertheless, all of these methods have the condition in common that the bulk fluids in which diffusion is occurring are stationary with respect to each other. In connection with a series of papers on mixing in capillary flow,8-10 Taylor suggested the use of a flow method for determining molecular diffusion coefficients. Additional studies 11-13 have been conducted on miscible displacements in capillary tubes, but the data from these studies were not used for the specific purpose of determining diffusion coefficients. The flow method proposed by Taylor results in a single value of the diffusion coefficient for the fluid pair used in the displacement experiments. This single value represents the true value for the fluid pair when the diffusion coefficient is independent of concentration. If the diffusion coefficient is a function of concentration, the single value obtained by the flow method gives an average value for the coefficient of the fluid pair. These average values are based on diffusion taking place over the entire range of concentration, i.e., from 0 per cent of one fluid to 100 per cent of that same fluid. In field applications of the miscible displacement process, gradients occur over the same range of concentration as are found in the displacements in capillary tubes. Molecular diffusion coefficients obtained from the capillary flow method should, therefore, be especially relevant to field operations. This investigation was undertaken to evaluate the feasibility of obtaining molecular diffusion coefficients from capillary -flow experiments. In making this evaluation, diffusion coefficients were first determined for two systems from data obtained in capillary flow experiments. These values of the diffusion coefficient were then compared to values obtained by other methods. MIXING IN CAPILLARY FLOW—THEORETICAL The theoretical basis for determining molecular diffusion coefficients from capillary flow experiments is the partial differential equation relating the mechanisms of flow and diffusion. For the case of laminar flow and with a fluid pair having equal densities and viscosities, the equation is given by

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

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

By M. C. Fuerstenau, A. M. Gaudin

A homogeneous suspension is viewed by X-rays. The radiation density seen, affords a measure of the extra absorption due to the solids contained. This radiation density, at a predetermined depth, varies with time. The chart of this relationship gives the desired determination. s it is frequently necessary to know the distribu-rt tion of particle sizes in various products used in industrial applications, a number of means have been devised for these determinations. Some of them include direct microscopic measurement, screening of product, extraction of samples from a settling slurry, air and liquid elutriation, adsorption of gases and liquids, and light absorption. The absorption of light as a means for particle size analysis has been known and employed for some time. X-ray absorption has not received the attention that light has enjoyed primarily because of the complex and costly instrumentation involved with X-rays. In 1954, Brown and Skrebowskil suggested the use of X-rays for this type of determination but did not supplement their suggestion with experiment. In 1958, the authors conducted size analysis experiments on suspensions of quartz with the Trans-viewer: an X-ray unit specifically adapted for sedimentation research, and Rossa conducted experiments in a similar manner in 1959 with a stable gamma emitter, AmZ4l, on suspensions of uranium oxide. To illustrate the applicability of X-ray absorption for analyzing particle size distributions, an essentially mono-energetic and well collimated X-ray beam was used as a source. The beam was made mono-energetic by placing a zirconium filter before the counting window together with feeding the pulses from the proportional counter into a reverter, a device capable of discriminating against pulses of various heights. The sketch in Fig. 1 shows the collimation arrangement in more detail. As shown, the X-ray beam was collimated by inserting a 3/8-in. ID pipe into one of the ports of the X-ray tube and further by two small steel plates with slits 3/32-in. wide and 1/2-in. long placed in front of the window of the Geiger-Mueller (G-M) tube. One of the plates was located 1/2 in. and the other 2 1/2 in. from the window of the counter; both plates were contained in an aluminum holder mounted on the box containing the G-M tube. This type and degree of collimation should eliminate any intensity due to scattered radiation; the transmitted intensity would then be solely a functior. of the absorption by the suspension at that level. Experimentally, homogeneous suspensions of finely ground quartz were prepared, Daxad (a sulfo-nated compound made from wood) was added as dis-persant, and the transmitted X-ray intensity measured as a function of time at a fixed depth in the suspension. In more detail, the experimental procedure involved the following: 1) measuring the transmitted X-ray intensity through the Daxad solution at a predetermined position or settling depth in the absence of solids, 2) placing a known amount of dry solids into the vessel and agitating thoroughly, 3) measuring the transmitted X-ray intensity at time zero (i.e., immediately upon cessation of stirring), and 4) measuring the transmitted X-ray intensity as a function of time (a 10-sec counting interval being chosen and counting being started 5 sec before the desired time). Since a linear relationship exists between percentage solid in a slurry and the logarithm of transmitted intensity for a homogeneous material, the cumu-

Jan 1, 1961

By T. Vasilos, J. T. Smith

THE mechanism of pressure sintering, or hot pressing, for ceramic materials, has been investigated by several researchers.1-8 Plastic flow has been suggested as the rate-determining mech-anism1,2 while grain boundary sliding3,4 and fragmentation processess have been found controlling in other studies. Vasilos and Spriggs5,6 have reported that pressure sintering occurs by a diffusion-controlled mechanism. For the pressure sintering of metals, investigations with lead7 and with tin6 have not identified the operative mechanism during consolidation of these powders. This note reports the results of a pressure sintering study initiated to determine the mechanism controlling the densification of spherical copper powder. Temperatures of 300°, 400°, 500°, and 600°C at 2500, 5000, 7500, and 10,000 psi pressures were investigated. The copper powder for this study was obtained from Metals Disintegrating Co., Elizabeth, N.J. The minimum purity was reported to be 99.9 pct. The powder size was -325 mesh and a lineal analysis showed the average size to be less than 1 p. An alumina die body and punches were employed. Graphite spacers were placed between the copper powder and the punches before pressing. It was expected that the graphite would form carbon monoxide with any entrapped air and form a reducing atmosphere within the die cavity. This condition was generally realized during pressing although a minor amount of oxidation occurred at 500" and 600°C. The oxide particles were located randomly in the pressed samples and were not considered to influence the densification process. Calibration tests showed the temperature of the copper pellet to be about 20°C less than the die wall temperature, probably resulting from heat losses through the punches. This difference was compensated during the experimental pressing operations. The hot pressing equipment employed for this work has been described by Spriggs et a1.9 and was calibrated with a Baldwin-Lima-Hamilton load cell. Travel of the punches was monitored with a dial indicator showing increments of travel to 0.0005 in. Relative density was calculated at selected time intervals for a known copper weight, punch travel, and die diameter. The accuracy of this calculation was *1 pct. Average grain size after hot pressing was obtained by lineal analysis. Typical densification curves are presented in Fig. 1 and Fig. 2 where relative density is plotted vs log time as a function of applied pressure for the 300" and 600°C pressure sintering temperatures, respectively. The curves for 400" and 500°C were similar to those at 300°C and consequently are not shown. The curves at 300°C as well as 400" and 500°C do not show an end-point or limiting density for a particular applied pressure. At 600°C, a density in excess of 95 pct of theoretical has been reached such that further densification is extremely slow for the 5000, 7500, and 10,000 psi pressures. It is suggested that further densification at these pressure levels may be retarded or prevented by discontinuous grain growth, by entrapped gases, and by other microstructural limitations as described by Cob1e.10 The curve for 2500 psi hot pressing pressure is a straight line over the time interval although a slight change in slope may be occurring at the last time intervals. The plastic flow models,1,2 derived for hot pressing, predict the existence of an end-point (or limiting) density dependent on the yield point of the material at the temperature of interest. The yield stress of copper at a 0.005 strain level has been

Jan 1, 1965

By A. M. Skov, L. F. Elkins

Inclusion of anisotropic permeability in mathematical analysis of pressure transients observed during development of the huge Spraberry field indicates a major fracture trend which is in good agreement with that observed by fluid-injection tests spread over a 12- by 17-mile area. Delineation of this trend is important in selecting a pattern of injection for the pending large-scale water flooding in this field. Determination of reservoir parameters yielding best agreement between calculated pressures and observed reservoir pressures in newly completed wells was made using an IBM 650 computer. INTRODUCTION The Spraberry field covering 400,000 acres is a tight sand of less than 1-md permeability cut by an extensive system of vertical fractures. Primary recovery dominated by capillary retention of oil in the fractured sand matrix blocks is less than 10 per cent of oil in place. Strong forces of capillary imbibition of water into the sand, coupled with water flow under dynamic pressure gradient, indicate considerable increase in oil recovery can be achieved through water flooding. Best results will occur if the pattern of water injection is selected to force the water flow across the grain of the major fracture system. Existence of an oriented vertical fracture system in the Spraberry, observed first in cores, was highlighted more recently by the 144-fold contrast in permeability along and at right angles to the major fracture trend required to match relative water breakthrough times in Humble Oil & Refining Co.&apos;s waterflood test there. Spraberry operators since have conducted two gas-injection tracer tests for further areal confirmation of the fracture trend. Re-analysis of early reservoir pressure transients for evidence of anisotropic permeability has permitted many more local determinations of major fracture trend without resort to further field tests. This paper is limited to updating analysis of reservoir pressure transients to include anisotropic permeability as a test for orientation of the major fracture trend in the Spraberry. The reader is referred to Refs. 1 and 2 for information about general Spraberry reservoir performance and to Refs. 3 and 4 for information about significance of fracture orientation in selection of the injection-well pattern for water flooding the Spraberry. RESERVOIR PRESSURE DATA—DRIVER AREA During early development of the Spraberry Driver area, Sohio Petroleum Co. made the extra effort to measure the initial pressure in each of the 71 wells in a 5-mile-long area immediately after completion. Progressively greater reductions in pressure ranging up to 400 psi were observed throughout the six-month development period. Detailed data are presented in Ref. 1. Since the reservoir oil was undersaturated some 300 psi initially, early reservoir performance involving 55 new well pressures is subject to analysis as flow of a single compressible fluid in a porous media. Assumption of uniform permeability in all directions yielded good agreement between calculated pressures and observed pressures of these wells in the earlier study,&apos; but subsequent, additional, mathematical development to include anisotropic permeability in the transient pressure considerations and present availability of electronic computers to perform the much more extensive arithmetical calculations now yield even better agreement. The previous analysis, assuming uniform permeability, consisted essentially of calculating pressure reduction expanding circularly around each producing well and summing these effects at the time and location of each newly completed well for comparison with the measured pressure reduction. Permeability, effective fluid and rock compressibility, and permeability X thickness were varied until the best match with measured pressures was obtained. The present analysis, assuming anisotropic

By N. E. Garner

Laboratory experiments were conducted to investigate the cutting of a single diamond on limestone and shale under simulated down-hole conditions. A high-pressure chamber was modified by adding a rock-rotating device so that planing tests could be run. Weight on the diamond, diamond geometry and differential pressure across the face of the rock were varied. It was found that (I) high differential pressure reduced the volume of rock cut by the diamond at a given bit weight by strengthening the rock and changing the mode of failure; (2) a finely powdered rock left in the bottom of the grooves reduced the volume cut by as much as 50 percent; and (3) orientation of a single diamond about the axis of a drill point had considerable influence on diamond cutting efficiency. INTRODUCTION Diamond bits consist of many diamonds held in a matrix provided with water courses for fluid passage to clean, cool and lubricate the diamonds as they drill formation. Each diamond contributes its share to the overall effectiveness of a bit, but a lack of understanding of the performance of a single diamond has hampered efforts to engineer improvements in diamond bits. The need for those improvements prompted this study. This paper describes the findings of a laboratory investigation of the outting action of a single diamond where the differential pressure between the wdlbore and formation, diamond geometry and vertical force on the diamond were independently varied. The tests were conducted in a pressure chamber that provided relative horizontal motion (planing) between a rock and diamond under simulated down-hole pressure conditions. EQUIPMENT HIGH-PRESSURE CHAMBER WITH ROTATING DEVICE The high-pressure chamber used in previous rock mechanics studies&apos; was modified for this study by addition of a rock-rotating device. The modified chamber (Fig. 1) differs from the old chamber in that the rock can be rotated to provide relative horizontal movement between the rock and the penetrator which may be a small element of a diamond or a drag bit. In this study the penetrator was a single diamond. The rock sample holder inside the chamber is mounted on a thrust bearing and is connected to the external drive mechanism by a stem that passes through two seal elements. The seal elements are needed to allow formation pressure to enter the bottom of the rock sample through a small hole in the stem. The Kapseal seal elements (Tefion boots backed by O-rings) provide low friction (hand free) even at differentials of 5,000 psi. The drive mechanism is a variable-speed transmission with a speed range from 0 to 30 rpm and an output torque of 2,800 in.-lb. A right-angle-drive gear box couples the drive to the stem. The new design has two independent pressure systems (formation and borehole) to simulate down-hole rock stress environment. Formation pressure enters from the bottom of the rock through a 6-in. diameter area, and borehole pressure acts over the remaining surface. A 14-in. layer of zero fluid loss oil-base mud poured over the top of the rock and an O-ring seal on bottom allow a differential between the borehole and formation pressures. Bayol 50 was used as the hydraulic fluid in the chamber. The differentials reported here were obtained by elevating the borehole pressure while keeping the formation pressure at zero.

By B. Wilhelmi, W. H. Somerton

Simultaneous measurement of pore and elastic properties of socks under a wide range of triaxial stress conditions may be made with equipment and methods developed in this paper. Tests were run on three outcrop sandstones at several confining pressures and axial stresses up to approximately 80 percent of failure stress. The greater reduction in porosities and permeabilities occurred in the hydrostatic loading portion of the stress cycle, but substantial reductions in both properties did occur upon application of deviator (triaxial) stress. Triaxial elastic moduli were calculated from strain gauge measurements. Both Young&apos;s moduli and shear moduli were found to increase with increasing confining pressure. Poisson&apos;s ratios increased with axial stress but showed no consistent change with confining pressure. Only very general relations were observed between pore and elastic properties. The lesser decrease in porosity and permeability at higher stress values is related to increase in rigidity of rocks at higher stresses. INTRODUCTION Many investigators have studied the effects of hydrostatic stress on the properties and behavior of rocks. A few have studied the effects of triaxial stress fields on individual physical properties of rocks. To our knowledge no one has measured simultaneously pore and elastic properties of rocks under a wide range of triaxial stress conditions. This is the purpose of this work. Rocks exist in the subsurface under what Jaegerl refers to as "polyaxial" stress conditions (all principal stresses unequal). When two of the three principal stresses are equal, this is referred to as triaxial loading. The triaxial method of loading is most commonly used in modern testing work because experimental difficulties in measuring physical properties of rocks are greatly reduced over polyaxial testing. Although triaxial loading may not reproduce subsurface conditions precisely, subsurface stress conditions are never known precisely anyway. Triaxial testing would seem to be more realistic than testing under hydrostatic loading (all principal stresses equal) for application of results to subsurface problems. The objectives of this work were threefold: (1) to develop methods of measuring simultaneously pore and elastic properties of rocks under triaxial stress conditions, (2) to determine whether changes in these properties are of sufficient magnitude to be considered in subsurface calculations and (3) to attempt to relate changes in pore properties and elastic properties since they will have been measured at the same time under identical loading conditions. REVIEW OF EARLIER WORK PORE PROPERTIES Several investigators have studied the effects of stress on the physical properties of rocks. Fatt and Davis2 showed the reduction of permeability of rock with hydrostatic loading. McLatchie et al.3 and Knutson and Bohor4 related permeability reduction to compessibility of reservoir rock, again under hydrostatic loading. Wble,5 and later Redmond,6 evaluated the effects of applied radial pressure on conductivity, porosity and permeability of sandstones. In this latter work, no axial stress was applied. Dobrynin7 reviewed much of the previous work done . on the effects of overburden pressure on physical properties of rocks. He concluded that changes in physical properties are controlled to a large extent by pore compressibility of the rock. Several sets of curves developed from empirical relations are presented. These curves permit estimation of change of properties with net overburden pressure, knowing the pore compressibility and porosity of the rock. Net overburden pressure is defined as (Pe - 0.85 Pi ) where Pe is the

By R. D. Nelson, F. E. Bowman

Evidence has been found for a (010)a?(111)B[102]a /,11[341]/B crystallographic relationship between the monoclinic a and B phases of plutonium during the a?B transformation. The results of an X-ray diffraction study, where [010] fiber textured a samples were used in the absence of suitable single crystals indicate that the (OIO)a planes become (111)B planes upon transformation under near-equilibrium conditions. The transformation strain along the fiber axis, which was the [010] direction, was measured dila-tomekically; the measured 8.5 pct expansion compares favorably with the 8.44 pct expansion in the in-terplanar spacing represented by the (010)a?(111)B relationship. From plots of the atom positions in each of the related planes, a correspondence was deduced between the [102]a, direction and the [341]B direction. Further analysis of the two lattice structures in the light of this relationship revealed a rotation about the [102] direction of 1 deg 43 min and a slight shearing in this same direction of about 30 min. THIS paper discusses the (B crystal structure of plutonium and the orientation relationship existing between the monoclinic a and B phases during the a?ß transformation. The kinetics of the a?ß and (3 — a transformations have received extensive consideration, and because of the crystallographic complexity of these phases, the question of the mechanisms by which these transformations occur is especially intriguing. Nelson and Shyne1 evaluated the existing experimental evidence and concluded, from the indirect evidence available, that both the ß ? a and a ? ß transformations occur by a martensitic process. Their data reveal the anisotropic nature of these transformations which suggests the existence of specific crystallographic relationships between the parent and product phases. In this paper, we are concerned with the confirmation of the existence of such a relationship resulting from the a ? ß transformation and the derivation of the crystallographic details necessary to specify it. EXPERIMENTAL RESULTS X-ray diffraction and dilatometric data were obtained from samples of a plutonium having a [010] fiber texture. These cylindrical specimens, approximately 1 cm in diam, were prepared by transforming the ß phase to a under the influence of a uniaxially applied compressive load.&apos; The [0l0] direction is parallel to the direction of applied stress; it is also perpendicular to the (010) planes. This correspondence between the plane normal and the direction having the same indices is, in the monoclinic a plutonium structure, restricted solely to this (010)[010] combination. Such a specimen was mounted in an evacuated, elevated temperature specimen holder on an X-ray dif-fractometer so that the diffracting surface was perpendicular to the fiber axis. The a diffraction pattern, Fig. 1, taken at about 115ºC, revealed the predominance of the (020) reflection. This monoclinic, P 21/rn space group yields no (0k0) reflections when k is odd. The temperature was raised slowly to 160°C to effect the complete transformation to the phase; then it was lowered to 120°C so that a diffraction pattern of the ß could be obtained near the a = ß equilibrium temperature for direct comparison with the a at the same temperature. The (3 diffraction pattern, Fig. 2, obtained under these conditions showed a (222) reflection relatively much more intense than in the random pattern. The other reflections appearing with sig-

Jan 1, 1970

By D. A. Flanagan, B. L. Landrum, P. B. Crawford, B. C. Norwood

Herern is described a heat conduction model which has been developed to predict the pressure distribution within a reservoir being subjected to rrnsteady-state depletion. Although potentiometric models have been widely used to study reservoir pressure distribution, their application is limited to the study of steady-state cases. Electrical analyzers may be used to study unsteady-state phenomena, but their application and simulation have certain limitations. Much of the now accepted theory of unsteady-state fluid flow in porous media has evolved from previous heat conduction knowledge. An analogy between the flow of fluid in an unsymmetrical reservoir and con-duct ion of heat in a similarly shaped unsymmetrical metal plate is made to show adaptability of the hear conduction principle to transient reservoir cases. A description of the model and necessary instrumenrarion is included. Data show influence of well location on pressure distribution during unsteady-state deple[ion. Application of results is made to a hypothetical reservoir and to a reservoir located in an irregular aquifer. INTRODUCTION Development of methods for predicting behavior of petroleum reservoirs has resulted partially from adaptation of electrical and thermal flow concepts to the flow of fluids through porous media. Extensive investigation of electrical and thermal flow problems has resulted in mathematical solutions applicable to many reservoir problems. The relative ease of measuring electrical flow has brought wide acceptance of the use of the potentiometric model to study certain field problems. Potentiometric models are restricted to the study of steady-state conditions because of the instantaneous response of electrical flow to changes in the applied potential. The electrical reservoir analyzer has been one of our most useful instruments for studying certain phases of unsteady-state behavior1. Other than approximation techniques, no method exists for determining pressure distribution in an irregular reservoir during primary depletion. Heat flow, like fluid flow, is inherently transient. The model described here was developed to use this transient nature of heat flow, measurable in the laboratory, to simulate transient flow of fluids in the porous medium of the reservoir above the bubble point. A metal plate shaped like the reservoir to be studied was heated by a "well" source. Temperature increase in the plate is analogous to pressure decline in the field if the field is being depleted, or pressure increase if fluid is being injected. This study is intended to determine the feasibility of using the transient nature of the heat flow to study the transient behavior of fluid flow. It is hoped that this new technique will aid in the understanding of one of the more complex phases of reservoir engineering. THE FLUID FLOW-HEAT CONDUCTION ANALOGY The similarity of the equations which describe the flow of fluids through a porous medium and the conduction of heat through a solid is well known2. Rate of flow of fluids is given by Darcy&apos;s law3. while the rate of flow of heat is given by Fourier&apos;, (I,, - KA 4T........(2) Combining Eq. 1 with the equation of continuity for a compressible liquid", i.e.,results in Expressing y in terms of the compressibility, i.e., pr = peey ,..........(5) and by expanding the exponential and neglecting higher orders of cp one can obtain

By W. T. Turrall, M. J. Cook

A discussion of modern developments in beneficiation of fine sizes of anthracite, this paper includes a description of the plant flowsheet, an analysis of operating results, and a summary of fundamentals of mineral separation relative to plant operation. DURING the past several years Lehigh Naviga-tion Coal Co. has installed equipment at Tamaqua and Lansford collieries for cleaning and sizing No. 4 buckwheat —3/32 + 3/64 in. and for recovering and cleaning flotation coal —20 + 200 mesh at Tamaqua colliery. Papers describing some of these installations have been presented to AIME. At Coaldale colliery, however, methods used to clean fine coal were inadequate and inefficient. Management decided upon a thorough investigation of all existing processes before building a plant to treat the combined three sizes mentioned. Investigation, carried out on company property and at other plants, entailed a study of hydrotators, hydroclassifiers, tables, spirals, flotation machines, methods of screening, and methods of dewatering. On the basis of this investigation a flowsheet was designed which necessarily excluded many processes in use, and the writers wish to emphasize that no criticism is implied by omission of processes or machines. A determining factor in selection of many machines was the policy of standardizing equipment wherever possible. Flowsheet and Operational Data The first principle to be considered in designing of any plant is efficiency of cleaning as related to laws of classification. Beneficiation of coal may be described as separation of two materials having different specific gravities. Although there is an inherent ash content even in the purest coal which slightly affects its specific gravity, it is recognized that many particles are true middlings. The percentage of these particles allowable in the final product depends upon the ultimate ash, or specifications of the consumer. Richards in his study of free settling and hindered settling classification&apos; determined that particles of equal size and different specific gravities had unequal settling velocities. He also determined a definite size relationship for equal settling velocities with free settling and hindered settling conditions. Considering the design of machines used to clean No. 4 and No. 5 buckwheat, it may be assumed that both free settling and hindered settling conditions exist and that ratio of size for practical operations is between the established values. The essential factor is that the more efficiently and closely sized the feed is maintained, the more nearly perfect the resulting separation. Flowsheet design was governed primarily by factors related to the 1200 to 1500 tons per hr feed to the main plant. It should be stressed that because of the many sources of supply from mine and stripping, the physical characteristics of this feed varied considerably. The most significant influence was varying percentage of fines, necessitating that all machines be capable of operating under maximum peak load during periods of surge. Since cleaning practice is a wet process, about 9000 gal per min containing 14 to 15 pct solids, or 350 tons per hr, must be treated. This slurry contained all the fines, or —3/32 in., but also percentages of oversize, up to 1 % in., dependent on spills or breaks in screen jackets, a not uncommon occurrence in preparation plant practice. The flowsheet developed is shown in Fig. 1. Circled numbers represent points of sampling. Table I gives operational data including size, ash by size, composite ash, rate of flow in gallons per minute, percent solids dry weight, and average tons of solids per hour. It should be noted that these results are subject to human error in sampling practice. Detailed information reported is a representative evaluation of plant operation.

Jan 1, 1954

By J. A. Hren, G. Thomas

A high-temperature stage has been used to study precipitation and dissolution in thin foils of Al-20pct Ag alloy directly in the electron microscope. The results show that the aging characteristics are similar to those observed for bulk materials, and in addition, direct evidence is presented for clustering at the solid solution temperature. The limit of solubility of G. P. zones is-350°C. The heterogeneous precipitation, the growth and dissolution of ?&apos;have been studied in detail. Growth of ?&apos; parallel to the basal planes, occurs by continuous motion of inter -facial partial dislocations. Growth stops at intersecting particles due to formation of sessile a/6 <110> interfacial dislocations. Growth normal to the habit plane occurs by the spontaneous nucleation of partial dislocations on alternate (111) planes and the thickening rate is --10-3= times the lengthening rate. The precipitates dissolve in layers and the final product is a whole a/2 <110> dislocation . The nucleation -growth -dissolution processes are thus completely reversible. Kinetic data obtained from the dissolution rates of ?&apos;particles enable the diffusion coefficient to be determined and the activation energy is found to be 37.5(±4) kcal mol-1, which is that expected for dzffusion of silver in aluminum. The velocity of the interfacial dislocations is thus dependent solely on the rate of diffusion of silver atoms. RECENTLY it was demonstrated that a high-temperature stage could be effectively used in the electron microscope to follow directly the processes of precipitation and dissolution in thin foils of A1-4 pct Cu alloy.&apos; These experiments also provided kinetic data which enabled the diffusion coefficient to be estimated. This new development of electron metallography opens up many fruitful areas of research and it was decided to investigate the aging behavior in thin foils of A1-20 pct Ag alloy, particularly the mechanism of heterogeneous precipitation (i.e., precipitation on lattice defects). This system was chosen because the aging characteristics of bulk specimens are now well known,2"7 and because we wished to obtain more details of the role played by partial dislocations and stacking faults in the nucleation growth and dissolution of the hcp ?&apos; phase. This is of interest for the general case of fcc-- hcp transformations. Part of the A1-Ag phase diagram is shown in Fig. 1. The use of a high-temperature stage also enables a study to be made of the clustering process (formation of G. P. zones) and to determine the structure of the solid solution prior to the usual quenching and aging procedures. It is also possible to obtain kinetic data from measurements of the rates of dissolution of precipitate particles. Unlike the results of previous experiments on A1-4 pct CU 1 where much of the initial precipitation occurred on the foil surfaces, the results on the A1-20 pct A1 alloy were identical to those reported for bulk specimens6,7 and no surface precipitation was ever observed. Thus, it was possible to follow the aging sequence: G.P. zones — ?&apos; — ? directly in the thin foils. This paper is concerned with the formation and dissolution of zones and the ?&apos; phase, and no attempt has been made so far to follow the mode of discontinuous precipitation of ?. I) EXPERIMENTAL DETAILS The high-purity A1-20 pct Ag alloy was prepared by the Aluminum Co. of America in the form of 0.001-in.-thick sheets. These were cold rolled with intermediate anneals to a thickness of 0.0004 in.,

Jan 1, 1963

By H. P. Stuwe

It is investigated whether a model of sphere packing permits the computation of the lattice parameters and their variation with concentration in various a phases; it does, except in phases containing silicon. Several observations indicate that a particular role is played by the "E-atoms" in the lattice. THE phase is formed by many alloys of the transition elements.&apos;-&apos; Its structure has been described by Bergman and Shoemaker,5,8 Dickins, Douglas, and Taylor,7 and Kasper, Decker, and elanger.&apos; There are five kinds of crystallographically equivalent positions for the 30 atoms in the unit cell which have been designated by these authors I, 11, 111, IV, and V or, in the same sequence, A, B, C,D, and E, as shown in Fig, 1. For the purposes of the present paper two elements of the lattice will be discussed separately: The lattice planes perpendicular to the c-axis which are made up by the atoms A, B, C, and D and which determine the a-parameter; and the close-packed "vertical rows" of E-atoms which determine the c-parameter. asper,9 and Frank and kasper10 described the a phase and other complex structures as sphere packings. In these papers particular attention is paid to the coordination polyhedra. Lattice positions with an equal number of nearest neighbors are considered equivalent regardless of their crystallo-graphic description. This seems to work very well for the explanation of ordering phenomena in the a phase as far as information on this is available." If the description of the phase as a sphere packing is valid it should be possible to draw two further conclusions: 1) It should be possible to calculate the lattice parameters from the atomic radii of the components, and 2) It should be possible to predict the variation of the lattice parameters with concentration within one phase field. Both points are discussed in this paper. PROCEDURE Measurements of lattice constants were done with a symmetrical back-reflection X-ray camera. With Cr-radiation, about a dozen a1 lines could be observed with most samples. From the (550/710)- and 720-lines the a-value was computed. With this value the c-parameter was calculated from the (333)- and (413)-lines. Finally, the position of all the other lines was computed from the lattice parameters thus obtained. In this way the validity of the indexing could be checked and, at the same time, this proce- dure gave an estimate of the accuracy of the lattice parameters. A correction was made for film shrinkage. The atomic radii used in the computations are the Goldschmidt radii for coordination number 12. Actually, some of the atoms have a coordination number of 14 or 15, but the error thereby introduced is small. * For the computation of lattice constants from the atomic radii an "average atomic radius" of all components was used, i.e., complete disorder was assumed. Although this seems to be true for some systems like OS-Cr,4 it is not true in others.11,13,14 The ordering model, however, which has been suggested by Kasper on the basis of extensive research lets both kinds of atoms contribute to both the &apos;&apos;vertical rows" and the lattice planes perpendicular to them, so that even with ordering of this kind the assumption of complete disorder cannot introduce a very serious error in calculating lattice parameters. RESULTS Bergman and shoemaker6 gave a table of the interatomic distances in the Fe-Cr phase. If we disre-

Jan 1, 1960

By R. F. Rolsten

The preparation of pure metals by the thermal decomposition of volatile halides was developed byde boer&apos; and van Arkel.2 This has proved to be a useful technique for the refining of columbium,the mechanical and physical properties of which are seriously altered by such interstitial contaminants as oxygen, nitrogen, hydrogen, and carbon. The conventional deposition surface of a resis-tively heated metal filament was replaced, Fig. I, with an indirectly heated clear quartz surface. The impure columbium (feed) was contained in the annular space between the Vycor deposition bulb and the perforated molybdenum retaining cylinder. The head assembly was attached to the deposition bulb by fusing the Vycor at "A". This unit was then connected to a vacuum system composed of a three-turn glass coil to eliminate fracture due to vibration or expansion, a McLeod gage, a cold trap, a two-stage mercury diffusion pump, and a mechanical pump. The entire system was cold outgassed, and the 10 to 20 g of iodine in the auxiliary bulb (reservoir) was held in dry ice to prevent excessive loss of iodine. The columbium was hot outgassed at 800" to 825oC with the pressure maintained at l0-5 to 10-6 mm Hg for 48 hr and cooled under vacuum to 200° to 240oC. Iodine was sublimed from the reservoir and collected in the deposition bulb. The feed material temperature was adjusted to give the desired iodide vapor pressure and the finger temperature of 1000" C was independently attained and measured with a platinum-platinum 10 pct Rh thermocouple located in the center of the nichrome wire helical heating element, "F." The system was attended until temperature equilibrium was established, whence it operated without attention. This is a definite ad- vantage over the conventional "hot-wire" technique which required either constant attention or automatic control. In operation, the iodine vapor originally introduced into the vessel reacted with the crude metal to form a columbium iodide. This iodide vapor diffused to the hot deposition surface where it thermally decomposed, depositing columbium and releasing iodine. This liberated iodine then diffused back to the crude columbium where iodide was formed again. Thus, by repetition of this process, the iodine carried pure columbium from the crude material to the filament, where it deposited. The success of the "iodide" process in producing ductile metals depends on preventing the transfer of oxygen, nitrogen, hydrogen, and carbon from the feed material to the deposit. In many cases relatively minor quantities of these interstitial contaminants impair the ductility of metals. The material produced by the iodide process has such exceptional ductility that it is frequently and mistakenly assumed that iodide metal is always of very high purity. This assumption is not necessarily true since a number of metallic impurities can also be transferred. In the present investigation, the temperature of the feed material was adjusted to iodinate selectively and the deposition temperature adjusted to decompose the iodides selectively. By this technique, some metals present in the feed material were prevented from codepositing with the columbium. It is quite important to have the feed material symmetrically located around the deposition surface in order to obtain a uniform deposit at the maximum rate. For example, Runnalls and pidgeon3 observed that, when the base of a titanium filament was about 1 cm from the feed material, deposition was rapid on that section close to the feed, but that the rate decreased with distance from the feed so that a

Jan 1, 1960

By R. Schuhmann, G. W. Powell

The concept of local equilibrium is examined and, in particular, the applicability of the concept to two-phase binary diffusion couples is discussed. It is concluded that, if binary solid solutions are treated thermodynamically as ternary solutions, vacancies being the third component, then recent experimental observations of nonequilibrium interface compositions (i.e., compositions different from the phase diagram values) are consistent with a state of complete dynamic equilibrium at the interface. CLASSICAL approaches to diffusion and to other irreversible processes often utilize a postulate of "local equilibrium". That is, one assumes even in the presence of gradients, fluxes, and time changes of properties, a) that the Second Law can be formulated locally, and b) that thermodynamic equations of state among intensive properties are identical with those for equilibrium phases. A direct consequence of this assumption, for example, is that one would calculate the total entropy of a nonequilibrium phase simply by integrating the entropy density (eu per unit volume) over the volume of the phase, assuming that entropy density is the same function of energy and particle density as for equilibrium phases. In the study of diffusion and of diffusion-controlled phase transformations, metallurgists take this concept of local equilibrium for granted and utilize it automatically. In fact, experimental evidence of all kinds is overwhelming that "local equilibrium" is a valid limiting law, approached more and more closely as gradients and rates of change approach zero. On the other hand, we have neither experimental nor statistical thermodynamic proof that this principle of "local equilibrium" should remain rigorously valid at high gradients and for rapid processes. Theoretically, in fact, one would expect deviations because interactions between adjacent volume increments of material affect the evaluation of entropy. Also one observes experimentally structural and small energetic differences between equilibrium and nonequilibrium phases. The purpose of this paper is to present a thermodynamic framework for the study and correlation of deviations from strict local equilibrium in isothermal diffusion and related processes. This procedure, in effect, preserves the principle of local equilibrium by thermodynamically assigning the deviations from static equilibrium properties to a defect component. Although this theory is entirely phenomenological, the possibility does exist of correlating the theory with experimental observations of defects under favorable circumstances. The most immediate application appears to be to the equilibria at the interfaces in two-phase binary diffusion couples. A recent pertinent work by Eifert et al.1 has shown that the concentrations of aluminum in the a(fcc) and B(bcc) phases at the 0-8 interface in Cu:Cu-12.5 wt pct A1 diffusion couples are greater than those given by the phase diagram. It is postulated by these authors that the supersaturation of the o phase at the interface provides the driving force for the a-to-B lattice transformation; the excess free energy of the a phase, ?Gaxs, arising from the greater-than-equilibrium aluminum concentration is about 3 cal per g-atom. On the basis of the schematic free energy-composition diagram shown in Fig. 1, Eifert et al. concluded that the thermodynamic situation at the interface is one of the partial equilibrium; that is, = for the more mobile aluminum whereas will be reexamined in the latter part of this paper. THERMODYNAMIC PROPERTIES OF DEFECTS A rigorous analysis of the thermodynamics of defects requires an added degree of freedom in the Phase Rule sense, because the standard formulations of solution thermodynamics (partial molal properties of components, free energy minimization, equality of chemical potentials, and so forth) involve the implicit

Jan 1, 1970