Search Documents

By J. E. Korski

Quebec Cartier Mining Co. operates a low-grade iron ore deposit at Lac Jeannine, Que., which is located at the extreme southwestern end of the Quebec-Labrador Trough (Fig. 1). The facilities there include one of the world's largest concentrators, a hydroelectric plant of 60,000 hp, and the modern town of Gagnon. A 191-mile railway connects the mine to the seaport situated at Port Cartier on the St. Lawrence River. The harbor in Port Cartier can accommodate both lake and large ocean- going vessels and is also the site of a major grain terminal. The concentrator in Lac Jeannine is designed to produce some 8 million tpy (tons per year) of specular hematite concentrate, containing 66% Fe, from approximately 20 million tons of crude ore. An annual stripping program of 5 to 8 million cu yd of waste rock uncovers the required ore. Approximately 100,000 tons of ore and waste are removed daily on a year-round, 24-hr-per-day basis. The open pit covers an area with maximum dimensions of 8,000 x 2,600 ft (Fig. 2). Present mining plans call for an ultimate depth of approximately 1,000 ft below original ground level. The rocks in the pit area are typical of this portion of the Grenville geological province. They consist largely of marble, quartz rock, pegmatite, iron formation, and stratigraphically older gneisses. The iron-ore formation consists of medium to coarse-grained specular hematite in a quartz host rock. The formation occurs locally in a plunging, anticlinal form. The topography in the vicinity of the mine is rolling, with maximum variations in elevation of about 400 ft. The hills are thickly wooded with birch, black spruce, and some tamarack. Lakes and muskeg areas are plentiful. Precipitation and temperature records have been kept since 1959. The annual rainfall has ranged between a low of 23 in. and a high of 42 in.

Jan 1, 1969

By S. Razi

Injection luminescence of single crystals of rare-earth-activated CdS was investigated. Crystals of CdS:Nd, CdS:Er, and CdS:Yb, grown by the vapor transport technique, were suitably etched and hetero-junctions of Cu2S-CdS were fabricated on one side of the CdS substrate by chemical reaction. The device was then heated to 200°C in an inert atmosphere for 5 min. Static V-I measurements show that at low voltages (V < 1 v) the current is controlled by the junction, while at higher voltages the bulk Properties of the CdS substrate effectively control the current. Because of the large thicknesses no negative resistance was found. Reverse breakdown voltage varied from 5 v, for junctions that were not heat-treated, to 100 v, for diodes that received the proper heat treatment. The rectification ratio, measured in the range of 2 to 10 v, varied from 300 to 3000. Capacitance- voltage measurements indicate a weak dependence of the junction capacity on applied voltage. The capacity was found to increase appreciably with room light illumination. The variation of capacitance with frequency is indicative of a p+-i-n structure. An open-circuit voltage of 0.4 v was measured under an illumination of 100 mw per sq cm. Under forward bias, the bluish-green edge emission characteristic of CdS was observed. No sharp rare earth emission line was observed, which is in contrast to our results using ZnS:Tb, where, under high-voltage excitation, the rare earth emission was dominat over the lattice emission. The absence of rare earth emission implies that either the nature of the rare earth ion is different at the het-erojunction than in the bulk or the probability of energy transfer to the rare earth ion is small. RARE earth elements are playing an increasingly important role in the Synthesis of materials for a variety of purposes, e.g., solid-state lasers, color television screens, and so forth. Interesting systems result when wide-band-gap photoconducting lattices are activated with rare earths. A study of the lumines-

Jan 1, 1969

By J. J. Frawley, W. J. Childs

PREVIOUS studies have shown that if a supercooled liquid metal was perturbed by vibration the probability of nucleation is is During a recent investi- gation to determine quantitatively the amount of energy needed to cause nucleation in a given supercooled liquid specimen, it was found that at low-energy values in the megacycle range, below the threshold energy value for cavitation, the amount of dynamic supercooling was greater than the amount of static supercooling. That is, if a specimen had a given amount of thermal supercooling under static conditions, its thermal supercooling could be increased by the addition of low-energy sound waves. The material used in this investigation was high-purity bismuth. This material was chosen because of its low melting point and the high degree of supercooling which could be easily obtained. The sound gener- ator employed in this work was a commercially built Sperry ultrasonic testing unit. This unit was equipped with a barium titanate crystal as a piezoelectric transducer. These sound waves were transmitted to the supercooled melt by a Pyrex rod, as shown in Fig. 1. The transducer was coupled to this Pyrex rod wave carrier by motor oil. The Pyrex wave carrier was approximately 1.5 ft long to prevent damage to the transducer by heat. The diameter of the Pyrex wave carrier appears to be critical. In the initial tests a wave carrier of 0.5 in. in diam was employed. However, the sound waves had no effect on the amount of dynamic supercooling. When the diameter of the wave carrier was increased to approximately the inside diameter of the Pyrex crucible, the sound waves were more effective. Approximately 200 g of metal was melted under a stannous chloride flux in a 30-ml beaker by induction heating. The temperature was recorded by a Pyrex-encapsulated chromel-alumel thermocoule which was placed in the beaker after the bismuth melted. The Pyrex wave carrier was immersed slightly below the surface of the liquid bismuth and the specimen was reheated to about 500°C. It was observed that the amount of supercooling depends on the amount of superheating the specimen receives; therefore, for meaningful results the amount of superheating was the same for all heating cycles. The specimen was allowed to cool in air and the change in temperature was continuously monitored with a potentiometric strip-chart recorder.

Jan 1, 1969

By A. S. Odeh

This paper analyzes the effect of limited entry to flow at the wellbore on the steady-state productivity of a well. Wells that have been opened to flow along a fraction of their productive interval are termed wells with limited entry. Previous work treated the cases of a partially penetrating well, a well producing from the central portion of the productive interval and a well in which several intervals equally spaced were open to flow. In this paper the open interval can be located anywhere within the productive interval. Thus, in a sense, it generalizes previous work. The finite cosine transform was used to arrive at a solution for steady-state flow of a slightly compressible fluid. The solution was programmed for a CDC 1604 computer. Numerical vaIues for rd = 660 ft, r, = 1/4 ft, and range of sand thickness of 20 to 200 ft are presented in graphical form. The effect of rd and r, values on the result is shown in a table. The correct calculation of skin and damage ratio in the presence of limited entry to flow is explained and illustrated by examples. Moreover, the paper shows how to calculate the net decrease in productivity due to the combined effect of limited entry and perforations. INTRODUCTION In some wells only a fraction of the productive interval is open to flow. Location of this fraction is usually dictated by formation characteristics and reservoir behavior. For instance, if a gas cap exists, the open interval is located away from the gas-oil contact to prevent any possible gas coning. Wells that intentionally have been opened to flow along a fraction of their productive formation are tened wells with limited entry. Obviously, unintentional completions of this type also exist. Limited entry to flow decreases well productivity. Magnitude of the loss depends on the fraction of the formation open to flow, on the thickness of the sand, on the location of the open interval and on the ratio of rd /r, where r , is well radius and rd is the drainage radius of the well. The use of pressure buildup data on producing wells to calculate the condition of the formation around the wellbore is an accepted practice. van Everdingenl and Hurst2 introduced the concept of she skin factor s considered to be due to a thin layer of different permeability immediately around the wellbore. These authors dealt with the case of a well of complete radial geometry, i.e., a well with open-hole completion that completely penetrates the formation. The presence of a low-permeability skin results in a loss of productivity, as does limited entry. Therefore, if pressure buildup data obtained on a well with limited entry are used to establish the presence or absence of skin (i.e., formation damage), and a correction is not made for this loss of productivity, the calculations would result in an erroneous skin value. They might indicate the presence of formation damage when in reality there is none, or they might indicate a value larger than the true value. This could lead to an incorrect basis for planning remedial measures. Muskat3 studied the problem of partially penetrating wells for the case of incompressible flow. He presented equations and figures which allow estimation of loss in productivity. Brons and Marting,4 using equations based on Nisle&apos;s work,5 studied the loss of productivity for three cases. The first was for a partially penetrating well; the second was for a well producing from only the central portion of a productive interval; and the third was for a well in which several intervals equally spaced were open to flow. Their work was for steady-state depletion-type reservoirs wherein the well radius of drainage is established and the fluid is considered to be slightly compressible. Considered in this paper is the problem of wells with limited entry in which the open intervals are located anywhere within the productive sand. The finite cosine transform is used to arrive at a

Jan 1, 1969

By E. R. Russell, N. E. Richards

The current ejyiciency of aluminum cells was derived from the metal produced over a period of time and the theoretical faradaic yield. The difference in the actual amount of aluminum in the cathode at the beginning and end of the period must be determined. The weight of aluminum in the cathode was calculated from the dilution of an added quantity of impurity metal. Use of multiple indicator metals, copper, manganese, and titanium, demonstrated that the weight of aluminum in cells can be determined to within 1 pct with routine but careful chemical analyses. Over intervals of the order of 30 days, current efficiencies reliable to within 1 pct can be obtained. INVESTIGATIONS beginning with those of Pearson and waddington ,&apos; through the most recent published work of Georgievskii,9-11 illustrate the direct relationship between the composition of the anode gas and the applicability of analysis of anode gases to the control of alumina reduction cells. McMinn12 noted the lack of an independent method for measuring cell production efficiency over the short term. There is no doubt that changes in the current efficiency are immediately reflected in the composition of anode gases. However, the accuracy of faradaic yields calculated from gas analyses depends upon the degree of interaction between primary anode gas and Carbon.6 A conventional industrial practice of obtaining long-term current efficiency for production units from mass balances and quantity of electricity is generally insensitive to the impact of planned control of any one or more of the influential reduction cell parameters such as temperature, alumina concentration, and mean interelectrode distance. Consequently, there is a real need in the aluminum industry for a procedure to obtain the absolute cell current efficiency over a short term—10 to 30 days—both for the calibration of values obtained from gas analysis6 and for evaluating the effect of controlling specific parameters in the reduction process. The amount of aluminum produced may be determined by considering the cathode pool as a reduction of an impurity metal in aluminum. Analyses over a period will show a decreasing concentration of the impurity due to the accumulation of aluminum solvent. The increase in aluminum inferred from analyses is the amount produced by the cell during the period. Combining weights of the cell aluminum in the cathode at the beginning and end of a specific period, weights of aluminum tapped and the quantity of electricity passed during the interval will yield the current efficiency. Smart,I3 Lange;4 Rempel,15 Beletskii and Mashovets,16 and winkhaus17 have used dilution techniques to determine aluminum inventory in alumina reduction cells. A technique for determining the weight of aluminum in production cells by addition of small amounts of copper to the aluminum cathode was described by smart.13 The precision in values of the aluminum reservoir through dilution of copper in the cathode ranged from about 1 to 3 pct depending upon the quantity of copper added in the range 0.2 to 0.01 wt pct, respectively. Because the method appears so direct and apparently simple, one would not anticipate difficulties in application to industrial cells. The objective of this study was to resolve this problem associated with the trace metal dilution technique for determining the amount of aluminum in a cell. The approach in evaluating trace metal dilution as a basic factor in determining the weight of aluminum in the cell reservoir, and the absolute current efficiency of the Hall-Heroult cell, was to dilute more than one trace metal in the aluminum cathode so that we could discriminate among complications arising from physical mixing, the possibility of separation of intermetallic compounds, loss of the added elements, and chemical detection. EXPERIMENTAL METHODS These experiments are not complex but require standardized procedures. The technique involves addition of the trace metals to the cathode, knowing when these metals are homogeneously distributed in the liquid cathode, timing of the sampling, employing accurate and precise analytical methods, using reliable procedures for monitoring the amount of electricity passed through the cell, and accurate weighing of aluminum removed from the cell during the particular period. More accurate results might be obtained if the increment in concentration of the added indicator metals were of the order of 0.1 to 0.2 wt pct. The method must be applicable to production units and, hence, the contamination of the aluminum minimized. For this reason, the concentration of trace metals in the cathode was kept below 0.07 wt pct and generally at 0.04 wt pct level. Trace quantities of copper, manganese, titanium, and silicon are already present in virgin aluminum and are suitable as additives from electrochemical and analytical points of view. Concentration of silicon is quite dependent upon the characteristics of the raw materials and was not used extensively in this work. Chemical Analyses. All instrumental analyses require calibration against an absolute technique such as a gravimetric, volumetric, or spectrophotometric method which represents the ultimate in sensitivity, precision, and accuracy. On review, the best methods for copper appeared to be optical absorption without

Jan 1, 1969

By Jerry L. Straalsund, D. Bruce Masson

A dew-point technique was used to determine the thermodynamic activity of zinc at 430°C in a series of e phase Ag-Zn alloys. The composition of the alloys ranged from 72 to 88 at. pct Zn. This range included the composition at which Massalski and King3 found a reversal in the composition variation of the crystallo-graphic c/a ratio, which they attributed to an overlap of the Fermi surface across the (002) faces of the Brillouin zone. The data is presented in a graph in which RT In yz,, where yz, is the activity coefficient of zinc, is shown as a function of atomic percent zinc. This curve has an unmistakable change in slope at the same composition that Massalski and King observed the beginning of the reversal in the c/a ratio. This change in slope of the thermodynamic data is also attributed to a Brillouin zone overlap. Equations are presented to demonstrate that the thermodynamic activity can be related to the density of states of the conduction electrons, and that the observed phenomena are consistent with this model. It is also demonstrated that the contribution of the density of states can be related to the excess stability, a phenomenological parameter recently shown by Darkeen" to be significant in the interpretation of thermodynamic data of metallic phases. The data seem to indicate that zone overlap has caused a spinodal point, and the resulting misci-bility gap, in the phase diagram. THE problem of developing an adequate thermodynamic model of solid solutions has proven to be difficult, and is still only partially solved. The main approach has been to develop a statistical model, such as that of an ideal solution, regular solution, and so forth, to which can be attached corrections for electronic, vibrational, magnetic, ordering, or other contributions. Such corrective terms are usually derived on an ad hoc basis, and it is difficult to predict in advance what their relative importance will be. This problem has been discussed in general terms by Oriani and Alcock,&apos; who have reviewed several thermodynamic models and a few empirical correlations. The measurements described in the present paper were made to demonstrate in a special case the importance of one such corrective term, the contribu- tion of the energy of the conduction electrons of an alloy. It was our premise that the contribution of the energy of the conduction electrons to the thermodynamic activity of the alloy components could be detected; further, that such an effect would be observed at alloy compositions where other phenomena, also ascribed to the energy and density of states of the conduction electrons, are observed. The idea of the importance of the conduction electrons is hardly new. Hume-Rothery and his adherents have developed the well-known theory of alloy phases in which the sequence of phase fields in binary equilibrium diagrams, especially those involving the noble metals with the IIB, IIIB, and IVB subgroups, can be correlated by replacing the composition variable with the ratio of conduction electrons to atoms, e/a. Jones and others have developed a physical explanation for this correlation, in which they consider the solubility limits of phase fields to be restricted by an intersection between the Fermi surface and a Brillouin zone. The general features of the model are also quite well-known—presumably zone intersection causes the density of states to decrease at critical alloy compositions. The attendant increase in energy of conduction electrons in the original crystal structure allows an alternate structure to become more stable as the concentration of polyvalent solute is increased. In spite of the wide acceptance of these ideas on phase stability, there is only indirect* evidence, such as the variations in lattice parameter recorded extensively by Massalskizy3 and others, that Brillouin zone interactions occur. There are few experimental measurements, other than the correlations of the phase sequence, that substantiate the premise that the energy of conduction electrons affects the solubility limits of alloy phases. Much thermodynamic data of alloys has been found to be consistent with the theory; yet there is a lack of detailed data at compositions where zone intersection and overlap are thought to occur. One would expect that the energy of the conduction electrons would make a measurable contribution to the thermodynamic properties of alloys at compositions near zone intersection and overlap if the theory of Hume-Rothery and Jones is correct. This conclusion cannot be avoided, because the phase boundaries are determined by the requirement that the chemical potential of the components be equal in both phases at equilibrium. An electronic effect large enough to alter the stability of a phase should also affect the thermodynamic activity by a measurable amount.

Jan 1, 1969

By Stanley C. Suboleski, J. Richard Lucas

Simulator I, a program to simulate production in room and pillar mining systems containing up to six face operations, was developed at Virginia Polytechnic Institute under sponsorship of the Office of Coal Research. Since its inception, several hundred analyses have been run, with over thirty companies using the program. This paper describes the simulation technique used in this program, the output achieved from an analysis, and reviews the experience gained in working with Simulator I.

Jan 1, 1969

By Paul A. Flinn, Jean-Marie Robert Genin

The MCssbauer effect was used to study the distribution of carbon in Fe-C alloys of high carbon content during room-temperature aging. In freshly formed rnartensite the carbon distribution is approximately random; carbon clustering occurs dun&apos;ng aging. In the retained austete the carbon atoms are isolated, suggesting an efiective repulsion; and no rearrangement during room-temperature aging is evident. The kinetics of the carbon clustering in martensite has been studied in detail at 0" . 10". and 20°C. The time and temperature dependence observed is sinzilar to that reported for a variety of other aging effects, such as resistance change, hardness increase, and recovery of internal friction. The stabilization of the martensite kransformation also appears to be related to the clustering effect. DESPITE extensive efforts since the early days of Reaumur, the mechanism which provides the strengthening of steel remains controversial in many cases for several reasons: 1) The structure of as-quenched martensite is not well-determined. 2) Experimental methods capable of distinguishing between various hardening mechanisms are not well-developed. 3) The temperature of formation of many low-alloy martensites is so high that the martensite formed during quenching may be already tempered to some extent before any observation. Among all the hardening phenomena, the most puzzling one occurs at room temperature, where no structural changes have been detected by the usual methods of X-ray diffraction or electron microscopy. However, a new and more sensitive technique became available with the discovery of the MBssbauer effect. This effect makes possible the detection of changes in the local environment of iron atoms in solid solution. It has been extensively used in the study of substitutional alloy systems, for example, Fe-Al,1 Fe-Si,2,3 Fe-Rh,4 and Fe-Ni.5 Making use of the possibility of obtaining new structural information, we undertook an investigation of the aging of high-carbon Fe-C martensite to follow the structural changes involved in the aging at room temperature. Preliminary results6 showed that a clustering of the carbon atoms occurs with time. A systematic study of the kinetics of the process showed a striking similarity to the kinetics of stabilization of the marten-sitic transformation, thus suggesting further investigation of the relationship. EXPERIMENTAL TECHNIQUE To eliminate difficulties with the aging of martensite prior to the first experimental measurements, we needed an alloy with an M, temperature not above room temperature. Other investigators7 were misled by results obtained on alloys which transformed above room temperature and underwent structural changes prior to the first Mossbauer measurements. Winchell and cohen8 dealt with this difficulty by using a third element, nickel, to depress the MS below room temperature. The presence of nickel did not seriously interfere with hardness or resistivity measurements, but it would complicate the interpretation of the Mossbauer data by increasing the variety of iron environments. We chose, therefore, to obtain the necessary lowering of the M, simply by increasing the carbon content. After the first quench to room temperature, our specimens are retained austenite, which we later quench into liquid nitrogen. The martensitic transformation takes place suddenly, so that the as-quenched condition of the samples is retained, since the temperature is low enough to prevent any carbon diffusion. A Massbauer spectrum at 77°K can be taken and analyzed; after an appropriate up-quenching to the aging temperature, a new Mossbauer measurement at 77°K can be made. Fig. 1 shows a typical cycle of experiments. The samples were foils 25 p thick of high-purity

Jan 1, 1969

By J. J. Jonas, J. L. Uvira

Equipment was constructed which permitted the hot compression of 99.8 pct Armco iron and 2.8 pct Si steel at constant true strain rates of 0.05 to 1 per sec over the temperature range 600" to 1000°C. A maximum true strain of 1.9 could be attained, followed by immediate quenching or by isothermal annealing. True stress-true strain curves were obtained, and the steady-state flout stress, a, strain rate, i, and temperature, T, were found to fit the following expression: The constants A, a, n&apos;, and AH were determined for the two ,materials. The activation energies obtained for Arrnco iron and silicon steel are 66 and 80 kcal per mole, respectively, and approximate those for seu- diffusion and for creep. This indicates that hot compression is a diffusion-controlled thermally activated process. THE hot compression of iron and steel is usually carried out in the austenite temperature range,1"7 as this is the one of greatest interest with respect to industrial hot working. However, a disadvantage of testing in the austenite range is that the hot-worked structure is destroyed by the phase transformations taking place during cooling. The present study was therefore carried out in the ferritic temperature range, and the hot-worked substructures were preserved by rapid quenching after deformation. The substructural results will be described in a separate paper, and the present communication is restricted to the strain rate/stress/temperature relationship observed during hot working. The results of hot-compression testing are usually given in the form of true stress/true strain curves, with a different curve for each temperature/strain rate ~ornbination.&apos;,~,~ Sometimes simple algebraic equations are fitted to the curves, so that certain extrapolations or interpolations can be performed.1&apos;2&apos;5 These equations are, however, devoid of mechanistic interpretations. The present study was undertaken to show that hot compression can be described by an equation of the type used in creep studies. That is, it will be shown that the temperature dependence of the strain rate in hot working is well-represented by an Arrhenius term, with an activation energy similar to that for creep. Also it will be seen that the stress dependence of the strain rate is similar to that for creep EXPERIMENTAL METHOD AND MATERIALS The investigation was carried out with the following three ferrous materials: 1) Armco iron; 2) Armco iron with 0.26 pct Mn; 3) 2.8 pct Si steel. The first two materials were utilized for investigations performed in the temperature range 600" to 800°C; the last one was used for tests at temperatures from 600" to 1000°C. The materials were of commercial purity,* as shown in Table I. Cylinders of five different sizes were prepared by strain annealing. The resulting mean grain diameter was 3 mm in the 1-in. Armco iron and 2 mm in the other sizes and grades. Different ratios of diameter/height were used to enable the determination of the coefficient of friction. The specimen faces were provided with concentric grooves, which acted as lubricant holders during compression, as recommended by English and ~ackofen&apos; and others.2&apos;9 Powdered glass lubricants were applied to the samples before hot compression to decrease friction and barreling, following the practice of previous investigator~.~~~~~~~~ Different mixtures of glass powder were used for each temperature, after the determination of suitable compositions by experiment.&apos;&apos; Equipment for Hot Compression. Hot-compression tests were carried out on an adapted 150-ton vertical hydraulic extrusion press and also on a modified high-temperature creep machine. The modified extrusion press was used for temperatures up to 800°C. The following items were added to the press: a) a high-temperature furnace with built-in upper and lower platens; and b) a transmission assembly, which permitted the attainment of logarithmically decreasing crosshead speeds. These accessories are described in detail in Ref. 11.

Jan 1, 1969

By Richard R. Zupp, David A. Stevenson

A theoretical treatment of austenitic Fe-C-Z solutions (where Z represents a substitutional solute) was developed by extending a solution model which was previously applied to austenitic Fe-C alloys. The model assumes that each carbon atom in the solution excludes H interstitial sites from occupancy by other carbon atoms rmd thai otherwise the assembly of carbon atoms is equivalent to a modified Langmuir lattice gas. The results of the treatment provide theoretical justification for a phenomenological expression previously developed to represent the influence of substitutional solutes on the thernlodynamics of carbon in austenitic Fe-C-Z alloys. The expression is: where Kz is a parameter which characterizes the Fe-C-Z system at a given temperature and Xz and xj are the concentrations of substitutional solute and interstitial carbon expressed as atom ratios. The parameter, Kz, was associated with the differ- A number of investigators have measured the influence of various substitutional solutes on the carbon activity in austenitic Fe-C-Z alloys, where Z represents any substitutional solute. Experimental information is available for alloys containing the following substitutional solutes: Mn,&apos; Si,&apos; Ni,2-5 CO,&apos;-&apos;, Cr,5, 9-11 CU,4 Al,= v,11, 12 and Ti.13 This information is essential for studying many phenomena in Fe-C-Z alloys, particularly diffusion and diffusion-controlled phase transformations For example the pffprts nf various ence between the interaction energies between carbon atorns and host-lattice atoms in the Fe-C-Z solution and the corresponding energies for the Fe-C solution. The carbon activity in a ternary Fe-C-Z alloy, ac, was obtained in terms of the composition of the alloy and the single parameter, Kz: where y° is the Henry&apos;s law constant for the binary Fe-C solution. Both y° and H crm be determined from the data for Fe-C solutions, without recourse to any experimental data for Fe-C-Z solutions. The theoretical treatment was extended to alloys containing several substitutional solutes. By neglecting interactions between the substitutional solutes, an expression for the carbon activity in multicomponent alloys was obtained in terms of the Kz values characterizing ternary Fe-C-Z systems. substitutional solutes on the rate of pearlite growth is explained14 on the basis of their effects on carbon diffusion, which result from their influences on the carbon activity in the solid phases. In addition, recent studies have emphasized the dependence of the pro-eutectoid ferrite transformation15 and other diffusional processes in austenite16 on the influence of substitutional solutes on the carbon activity in austenite. They have indicated the need for suitable analytical expressions to represent such influences. The present authors previously developed general phenomenological expressions to represent the equilibrium partitioning of carbon between Fe-C alloys and Fe-C-Z alloys.&apos;2 The expressions are applicable to all the austenitic Fe-C-Z alloys for which experimental data are available, and are sufficiently simple to recommend their general application. For alloys in

Jan 1, 1969

By J. Brett, L. Seigle, L. Castleman, T. Montelbano

Impurity-induced low-temperature recrystallization of cold-worked tungsten was inuestigated with emphasis on the influence of nickel on the reaction. Palladium, nickel, aluminum, manganese, platinum, and iron greatly lower the recrystallization temperature of doped tungsten, which zs normally very high, but the recrystallization temperature of electron-be am zone-refined tungsten wzre is slightly raised by conlacl with nickel. Recrystallization can be induced at low temperature by the presence of solid nickel on the surface of doped tungsten wire, but apparently not by exposure to nickel vapors alone. Approximately 200 ppm of Ni dijjused into 10-mil wires at 1200 from a deposit of nickel on the surface produced total recrystallization, whereas more than 600 ppm of Ni could be absorbed frotn a vapor source without altering the fibrous structure of cold-worked tungsten. Once initiated, nickel-induced recrystallization required a continued source of&apos; nickel for propagation of the recrystallization front. The solubility of nickel in fibrous 10-mil W wire was approximalely 500 ppm at 1150&apos; C, and the activation energy for penetration of the recrystallization front was 52 kcal per mole. In many applications the usefulness of tungsten depends on critical control of its structure. Cold-worked tungsten with the fibrous structure developed by suitable thermo-mechanical treatment has a low, but technologically significant, ductility. It has long been known1 that traces of nickel, and perhaps other metals, are profoundly deleterious in doped tungsten, because they induce recrystallization at low temperature, which produces a brittle, equiaxed grain structure. This effect appears to be an exception to the general observation that recrystallization is impeded and the recrystallization temperature raised by the presence of impurities.2"9 Previous studies7-&apos;&apos; of the annealing and recrystallization of tungsten wire have divided the phenomenon into prior recovery stages, primary recrystallization and secondary recrystallization. The present investigation is concerned principally with primary recrystallization which is defined here as the replacement of the fibrous structure of deformed tungsten by equiaxed grains. The objective of this study was to explore the nickel-induced recrystallization reaction in tungsten and attempt to elucidate its mechanism. As well, an effort to define which other elements give rise to low-temperature-induced recrystallization was carried out. EXPERIMENTAL PROCEDURE The procedure adopted for these experiments was essentially to bring nickel and other elements into diffusive contact with cold-worked tungsten wires. The process of recrystallization was followed as a function of time and temperature by light and electron microscopic observations. First the influence of nickel on the recrystallization temperature of arc-melted, zone-refined, and variously doped tungsten wire was determined by electroplating a deposit of nickel on the surface of the wire and annealing at a variety of temperatures for 3 hr. The chemical analyses of the tungsten wires used in this investigation are given in Table I. The surface of the tungsten wire was etched with Murakami&apos;s slution&apos; and approximately 0.005 in, of Ni was deposited from a Watts-type low pH bath14 for the conditions of these experiments. Variations in plating thickness from about 0.001 to 0.005 in. had no discernible influence on the resulting structures. The wires were then annealed in an atmosphere of dry hydrogen to establish the recrystallization temperature. Concurrently, un-plated specimens were annealed to establish the recrystallization characteristics of nickel-free wire. The criterion of recrystallization was that the fibrous structure be completely replaced by equiaxed grains after a 3-hr treatment of temperature. This provided more reproducible results than use of the first recrystallized grain or a fixed proportion of re-crystallized structure as the critical observation. The structures encountered in longitudinal and transverse sections were examined by both light and electron microscopy at magnifications up to X32,000 using parlo-dion-carbon replicas shadowed with platinum for the latter method. Second, the influence of a variety of metals on the recrystallization temperature of 0.010 in. D alumina-silica doped tungsten wire, AW136-64, was determined. The elements were applied by electroplating whenever possible. Alternatively, they were vapor-plated on the tungsten wire and a greater thickness built up by coating with a dispersion of metal powder in nitrocellulose lacquer. Elements not amenable to either of these procedures were merely slurry coated on the tungsten. The recrystallization temperature was determined as above. Third, the nickel-induced recrystallization process in doped wire was studied more closely by electroplating 8 mils of Ni on 65-mil alumina-silica doped tungsten wire, AW153-NS10, and exposing the wire to temperatures of llOO°, 1200°, or 1300°C for various times in a hydrogen atmosphere. A circular recrystallization front, &apos;Onsisting of equiaxed grains, developed at the periphery Of the coated tungsten wire, and the advance of this front into the fibrous interior was studied. These experiments employed relatively coarse 0.065 in. D wire because the 0.010 in. D wire recrystallized too quickly to permit observation of the pene-

Jan 1, 1969

By J. E. Lawrence

Copper-decorated lattice disorders in silicon have been analyzed by electron transmission microscopy, chemical etching, and P-N junction reverse current measurements before and following different gettering treatments. Lattice disorders acted as the nucleating sites for copper precipitation. Copper precipitates in silicon resist chemical etching, absorb electron beams more than silicon, and cause excess reverse currents. Totally gettered disorders will chemically etch rapidly, be transparent to electron beams at extinction contrast, and not contribute to excess reverse currents. Three different gettering procedures were performed to determine the effectiveness of different EXTRANEOUS metal contamination can alter the minority carrier lifetime and resistivity of bulk silicon as well as cause excessive reverse currents in semiconductor devices. Previous workers have reported such activity for gold,1,2 copper,&apos;-* iron,&apos; and rnanganese.2,5 Gettering is a method employed to render these impurities electrically inactive. However, some property of the junction-degrading metals prevents gettering procedures from being wholly reliable. The comparative effectiveness of gettering with different agents has been the topic of earlier publications.2,3,6-11 Two mechanisms for gettering metallic impurities have been proposed: the distribution coefficient mechanism by Thurmond and Logan7 and the compound formation mechanism by Goetzberger and Shockley.2 Both of these models evolved from data gettering mechanisms. First, phosphorus was allowed to enter a protective SiO2 film but remain away from the Si-SiO2 interface to test the effectiveness of the compound formation gettering mechanism. Second, phosphorus was allowed to enter the backside of silicon wafers to test the effectiveness of the lattice strain gettering mechanism. Finally, nickel was diffused into the backside of silicon wafers to test the effectiveness of the distribution coefficient gettering mechanism. In each case, copper precipitates dissolved as a result of the gettering treatment. Tightly packed precipitates will resist gettering more than loosely packed copper precipitates. acquired by the electrical analyses of gettered semiconductor materials. An electrical analysis, however, cannot identify the fate of the gettered metals. The radiotracer technique has recently been applied to determine the redistribution of metals following different gettering treatments.3,9,10 The spatial resolu tion of the radiotracer technique prevented the analysis of metal precipitates which contribute to the excessive reverse currents of P-N junctions. The gettering processes employed by the earlier workers, almost without exception, allowed the gettering agent to enter the silicon lattice. This procedure makes it quite difficult to differentiate between the distribution coefficient and the compound formation mechanisms. The metallographic analysis of gettered copper precipitates discussed in this work supplements the earlier work by obtaining direct evidence of silicon lattice defects decorated with copper before and following different gettering treatments. The effectiveness of a process which restricts the gettering agent, phosphorus, to a dielectric film and thus away from the silicon lattice is discussed in addition to the ef-

Jan 1, 1969

By P. M. Kelly, E. P. Butler

Four Mn-CLL alloys, containing 60, 70, 80, and 90 pct Mn, respectively, have been examined in the quenched and the quenched and aged conditions using electron microscopy and electron, neutron, and X-ray diffraction. After certain heat treatments the alloys transform from fee to fct and in the tetraom1 condition show a domain structure parallel to {101} planes. Neutron diffraction indicates that the domains are antiferrornagnetically ordered. The domain boundary contrast has been examined using bright- and dark-field microscopy, and the contrast effects observed under favorable conditions have been used to deduce the c axis orientation in each domain. The domains are extremely mobile and can be nucleated at precipitate particles and screw dislocations. The domain mobility is responsible for the high damping capacity. In the aged material a Mn precipitates in the Kurdjumov-Sachs orientation and results of both electron microscopy and neutron diffraction indicate that the matrix separates into two components—one rich in manganese and the other rich in copper. ALLOYS of manganese and copper have the unusual combination of a high damping capacity and good mechanical properties and have been the subject of a number of investigations as part of a general interest in high damping capacity alloys for engineering purposes.&apos;,&apos; SO far, however, there has been no reported electron metallographic study of these alloys. The Mn-Cu system has an extensive range of solid solubility at high temperatures, and the equilibrium phases are expected to be y (fee) and a Mn. The high damping capacity is associated with a metastable tetragonal structure of variable c/a ratio, which forms from the high-temperature y phase. This latter phase becomes more difficult to retain as the manganese content increases, and alloys containing more than 82 wt pct Mn undergo a reversible martensitic fcc — fct transformation on quenching. The X-ray work of Basinski and christian3 showed that the Ms temperature for the transformation was below room temperature for alloys in the range 70 to 82 pct Mn and increased linearly with manganese content. When quenched from the y region, alloys in the range 50 to 82 pct Mn are cubic at room temperature, but become tetragonal if aged at temperatures between 400" and 600°C. The martensite transformation occurs on cooling from the aging temperature. Tetragonal alloys have a banded microstructure and the bands analyze to be traces of (110) planes. Similar microstructures have been observed in In-Tl4 and in other manganese-base systems, such as Mn-Au5 and Mn-Ni.6 The mobility of the bands in Mn-Cu alloys can be demonstrated by optical examination of a polished specimen surface subjected to a cyclic stress.7 The bands appear and disappear as the stress is varied, and X-ray measurements of the (200,020) and (002) peak intensities confirm that a reversible reorientation of the tetragonal structure occurs. Meneghetti and sidhu8 investigated the magnetic structure of Mn-Cu alloys and found antiferromagnetic ordering in furnace-cooled alloys of composition >69 at. pct Mn. Magnetic super lattice reflections occurred at the (110) and (201) positions and the proposed structure was fct with the spins along the c axis. A more complete investigation by Bacon et al.9 confirmed this structure. The magnetic ordering temperature Tn was found to increase linearly with manganese content in the same way as the Ms temperature, and at any composition, Tn > Ms. This relationship suggested that the magnetic ordering was responsible for the cubic — tetragonal transformation in the manganese-rich alloys. The purpose of this investigation was to study the mechanism of high damping and the structural changes that occur on aging. The main technique used was transmission electron microscopy, but X-ray and neutron diffraction experiments were also carried out. EXPERIMENTAL Materials and Heat Treatment. The four alloys, provided by the Admiralty Materials Laboratory. were of nominai composition 60, 70, 80, and 90 Mn and all had low impurity levels, <0.05 pct C, <0.2 pct Fe. This material was cold-rolled to 200-µ strip with intermediate annealing and then given a final heat treatment of 24 hr in the range 800° to 900°C followed by water quenching. An identical heat treatment was given a length of 3/4-in.-diam bar of the 70/30 alloy from which the neutron diffraction specimens were machined. It was suspected that the tetragonal structures would be metastable at room temperature, and so the alloys were not aged until required for experiments. After aging in a salt bath the alloys were water-quenched. Thin Foil Preparation. Initial thinning to 50 to 75 µ was possible in a solution consisting of: 50 ml nitric acid 25 ml acetic acid 25 ml water The surface deposit and grain boundary etching was removed by a final electropolish at around 20 V in an electrolyte consisting of:

Jan 1, 1969

By W. B. Gogarty, W. C. Tosch

A new recovery process for producing oil under both secondary and tertiary conditions utilizes the unique properties of micellar solutions (also known as microemulsions, swollen micelles, and soluble oils). These solutions, which displace 100 percent of the oil in the reservoir contacted, can be driven through the reservoir with water and are stable in the presence of reservoir water and rock. Basic components of micellar solutions are surfactant, hydrocarbon and water. They may also contain small amounts of electrolytes and co surfactants such as a1cohol.r. The specific reservoir application dictates the type and concentration of each component. A salient feature of [he process is the capability for mobility control. Micellar solution slug mobility, by way of viscosity control, is made equal to or less than the combined oil and water mobility. Mobility control continues with a mobility buffer that prevents drive water from contacting the micellar solution. Laboratory and field flooding have proven that the process is technically feasible and that surfactant losses by adsorption on porous media are small. Introduction projects are under way to recover the maximum amount of oil under the most favorable economic conditions.&apos; : New techniques are being developed to increase oil recovery,3" Polymer solutions are becoming an important means of controlling mobility in a waterflood. Thermal methods such as in-situ combustion and steam injection are being used in reservoirs containing highly viscous crudes. Surfactant flooding is receiving attention as a method of reducing interfacial tension to increase recovery.*&apos;" Exotic recovery processes have been considered primarily for &apos; perations. Economics are unfavorable in most cases for tertiary recovery. studies at the Denver Research Center of the Marathon oil CO. have led to a new oil recovery method.* Micellar solutions (sometimes called microemulsions, swollen micelles, and soluble oils) are used to recover oil by miscible-type waterflooding. Basically, these solutions contain surfactant, hydrocarbon, and water. The method can be used in either secondary or tertiary operations. First, thc concept of thc process is considered in terms of the requirements for an effective miscible waterflood ing operation. Next, micellar solution properties are described including structure, composition, and phase behavior with reservoir fluids. Fluid characteristics are then considered as related to mobility control, and, finally, laboratory and field results are presented to illustrate the efficiency of the process. Concept of the Process Unit displacement efficiency and conformance determine the effectiveness of any oil recovery mechanism. In theory, a miscible waterflood should be capable of a 100-percent unit displacement efficiency with a correspondingly high conformance. Requirements for the slug of a miscible waterflood include (1) 100-percent displacement of oil in the reservoir contacted, (2) controllable mobility, (3) the capability of being driven through the reservoir with water, (4) a low unit cost to enhance economics, and (5) the ability to remain stable in the presence of reservoir water and rock. Micellar solutions satisfy requirements for the slug of a miscible waterflood process. Our discovery that these solutions acted as though they were miscible by displacing all fluids in the reservoir and by being displaced by water solved the miscibility problem. Adequate mobility control is possible by variations in solution viscosity through compositional changes. Economic requirements are met since micellar solution costs below $6/bbl appear possible, Mi cellar solutions stabilize surfactant in the presence of reservoir rock and water, thus reducing the importance of the problem of surfactant loss by adsorption. Fig. 1 illustrates schematically how these solutions are used. Operations start with injection of a micellar solution slug that serves as the oil displacing agent. Next, a mobility buffer of either a water-external emulsion or water solution containing polymer (thickened water) is injected to protect the slug from water invasion. Finally, drive water (water used in a regular waterflood) is injected to propel the slug and mobility buffer through the reservoir. Reservoir oil and water are displaced ahead of the slug, and a stabilized oil and water bank develops as shown in Fig. 1. Stabilized bank saturations are independent of original oil and water saturations. This means that, for a particular type of reservoir, the displacement mechanism is the same under secondary and tertiary recovery conditions. Oil is produced first in a secondary operation. For tertiary conditions, water is produced first. Movement of the slug through the reservoir is stabilized by the mobility buffer. An unfavorable mobility ratio usually exists at the interface between the buffer and drive

Jan 1, 1969

By R. A. Wattenbarger, H. J. Ramey

A systematic study has been made of the application of the real gas pseudo-pressure m(p) to short-time gas well testing. The m(p) function can be used in real gas flow problems to account for the variation of viscosity and z-factor with pressure. A mathematical model was solved numerically to generate solutions of various real gas flow problems. The effects of turbulence, formation damage and wellbore storage were included in the model. The analysis of simulated well tests showed that the interpretation methods normally used for liquid flow are generally accurate when the m(p) is used. For practical rates without turbulence, the solutions are not rate-sensifive, and the flow capacity.. kh and skin effect can be determined accurately from either a buildup or a drawdown test. However, the kh calculated from a drawdown test can be significantly low when turbulence is present. A case was simulated in which this error wa.s 36 percent. Turbulence does not affect the determination of kh from buildup tests. The proper determination of kh from p and pa buildup and drawdown plots is developed by analyzing their relationship to the m(p) method. A simple equation is given that can be used for long-range gas well performance forecasting. This expression is compared with a method pre-sended by Russell et al." Introduction The formation flow capacity and wellbore damage condition can be determined for liquid-producing tests by means of buildup and drawdown tests. These tests make use of short-time pressure transient data rather than stabilized flow tests that are often used for gas well testing. Tracy&apos; presented a method of gas well testing that was based on ideal gas equations and utilized p2 buildup plots. This method was shown to be good for low pressure wells. MatthewsZ suggested plotting p and using an average slope and average gas properties. This method was more successful on high pressure wells. Al-Hussainy, Ramey and Crawfordh howed that variation in gas properties could be simplified by using the real gas pseu do-pressure m(p). Some cases showed that the use of the m(p) function Original manuscript received in Society of Petroleum Engineers office Aug. 8. 1967. Revised manuscript received June 11, 1968. Paper (SPE 1835) was presented at SPE 42nd Annual Fall Meeting held in Houston Tex., Oct. 1-4. 1967, and at SPE 4th Annual Eastern Regional ~eetin$ held in Pittsburgh. Pa., Nov. 2-3. 1967. 0 Copyright 1968 American Institute of Mining, Metallurgical, and Petroleum Engineers, Inc. &apos;References given at end of paper. provided an accurate method of interpreting buildup and drawdown tests. Because the work of Al-Hussainy et al. included only a few actual cases, there was a need to explore the use of the m(p) method for a greater variety of flow conditions. The occurrence of turbulent flow around the wellbore often is an important factor in gas well testing. Swift and Kiel,&apos; and Carter et a1.&apos; treated turbulence in gas well testing, but did not consider the variation of viscosity and z-factor with pressure. This paper presents the results of an investigation of the application of the m(p) method to buildup and drawdown testing. Vhe investigation included the effects of turbulence, formation damage and wellbore storage. Gas Flow Equations To formulate the mathematical model, many of the assumptions usually used in well testing theory are applied. The system has radial geometry with a closed outer boundary and is composed of a horizontal porous formation that has uniform and isotropic rock properties and uniform thickness. Allowance is made, however, for a radial region of reduced permeability near the wellbore. This region represents formation damage. The geometry of the system is shown in Fig. 1. Darcy&apos;s law does not always apply to gas flow. A more general expression is needed for non-Darcy, or "turbulent" REGION OF DAMAGED PERMEABILITY SEALED BOUNDARIES ClRCULAR WELL BORE Fig. 1—Radial flow model. AUGUST, BT7

Jan 1, 1969

By d&apos, C. Antonio, S. H. Reichman

In the course of our studies of the annealing behavior of vapor-deposited silver films we have determined the electrical resistivity as a function of annealing temperature for films in the thickness range of 630 to 3240A. Fig. 1 shows the resistivity of the silver films vs film thickness in as-deposited and annealed conditions. The films were prepared by depositing 99.999 pct pure silver from molybdenum boats onto Vycor substrates. Deposition was conducted at a residual pressure of 5 x 10"6 Torr at rates of 20 to 408, sec-&apos;. Thicknesses were determined using a multiple-beam interferometer. During annealing, resistance was measured using the four-probe method. Measurements were made every 5°C during two consecutive anneals to 300°C under a continuous vacuum of 2 x X6 Torr. The heating and cooling rates during annealing were approximately l° and 2°C min-&apos;, respectively. Typical resistivity vs temperature curves are shown in Figs. 2 through 4. There is an inflection point evident at 120°C on Fig. 2, the first anneal of a 3240A film. This inflection point was found on all films above approximately 1!00A on both first and second anneals. Below 1400A the inflection occurred at approximately 180°C. Figs; 3 and 4 show the first and second anneals of a 1370A film. The inflection point occurs at 180°C on the first anneal and approximately 120°C on the second. Our interpretation of this phenomenon is based on oxygen liberation from the film surface. It is well-documented that the presence of oxygen on silver surfaces causes a reversible increase in resistivity.1"3 The liberation of this oxygen, for example during vacuum annealing, would therefore cause a decrease in resistivity as observed on Figs. 2, 3, and 4 as a change in slope. This is accounted for by the decomposition of Ag20 on the film surface. The free energy of AgzO becomes positive between 125" and 225°C.4 The thickness dependence of the inflection points can be explained as an orientation effect. The films used in this study were preferentially oriented with their (111) planes parallel to the Substrate.&apos;&apos;~ In addition it has been sten that the thinner films, below approximately 1500A, possess a decided "hill and valley"

Jan 1, 1969

By W. M. Robertson

The formation of grain boundary grooves on surfaces of poly crystalline samples of cobalt and iron immersed in liquid lead has been studied. The grooves form by volume diffusion of the solutes cobalt and iron in the liquid. The diffusion coefficients of the solutes in liquid lead are derived from the measured rate of grooving. The diffusion coefficients are described by the relation D = Do exp (-Q/RT), with, for cobalt, Do = 4.6 x 10-4 sq cm per sec and Q = 5300 ± 800 cal per mole, and for iron, DO = 4.9 x 10-3 sq cm per sec and Q = 10,500 ± 1500 cal per mole. LIQUID metal-solid metal interactions occur at solid-liquid interfaces. Interfacial energy provides a driving force to change the morphology of the interface. Mullins1,2 has derived expressions for the kinetics of interface morphology changes driven by capillarity. These expressions can be applied to an isothermal system of a solid in equilibrium with a liquid saturated with the solid. Surface profile changes can occur by volume diffusion of the solute in the liquid, by volume self-diffusion in the solid, and by interfacial diffusion at the liquid-solid interface. A groove will form at the intersection of a grain boundary with a solid-liquid interface, reducing the total interfacial free energy of the system. The solid-liquid interfacial energy ? must be greater than half the grain boundary energy of the solid ?6 for Mullins&apos; calculations to apply. If ? is less than ?b/2, then the liquid penetrates the boundaries, separating the grains rather than forming grooves. Boundary penetration did not occur in the work described here. where CO is the equilibrium volume concentration of the solid in the liquid, Dv the volume diffusion coefficient of the solid in the liquid, ? the interfacial free energy of the solid-liquid interface, O the atomic volume of the solid crystal, k Boltzmann&apos;s constant and T the absolute temperature. Eqs. [1] and [2 ] also apply to grooving by volume self-diffusion in the solid,1 with CoODv = D Self, where DSelf is the volume self-diffusion coefficient of the solid. For a grooving mechanism of interfacial diffusion at the solid-liquid interface, the groove width is given by2 where CS is the interfacial concentration of the diffusing species, and DS is the interfacial diffusion coefficient. Eqs. [1] and [3] can be used to determine the mechanism of groove growth. A t1/3 dependence of the growth indicates volume diffusion and t1/4 indicates interfacial diffusion. In some cases, volume diffusion and interfacial diffusion both can contribute substantially to the grooving process, causing the time dependence to be intermediate between t 1/3 and t1/4.3 For these cases, the relative contributions of the two processes can be separated.4 However, in many cases, one process will be dominant, and the data can be analyzed on the basis of Eq. [1] or Eq. [3] alone. The time dependences for volume diffusion in a liquid and volume self-diffusion in a solid are the same. However, the self-diffusion contribution of the solid is usually negligible compared to volume diffusion in the liquid. After the grooving mechanism has been determined, Eq. [1] or Eq. [3 ] yields the kinetic parameter A or B. The kinetic parameter can be used to calculate values for the unknown quantities in the product CD?. Usually C is known or can be estimated. If ? is known, then D can be calculated. In a measurement of grain boundary grooving of copper in liquid lead,&apos; the time dependence indicated volume diffusion in the liquid. The quantities Co, Dv, and ? were obtained from the literature, giving excellent agreement between the observed values of A and the values calculated from Eq. [2 ].5 In a study of the grooving of several refractory metals in liquid tin and liquid silver, A1len6 educed that grooves formed by volume diffusion in the liquid. In a study of nickel in a nickel sulfide melt, Steidel, Li, and spencer7 found volume diffusion grooving kinetics. Both Dv and ? were unknown, so they could not obtain either one separately, though they did obtain a reasonable value for the temperature dependence of the product Dv ?. Several methods have been used to obtain surface profiles. It can be done by sectioning through the interface7 or by chemically removing the liquid from the solid surface after solidification of the liquid.6 However, if the liquid dewets the solid on removing the solid from the melt, then the interface can be observed directly. This method was used previously&apos; and was utilized also in the present study. EXPERIMENTAL PROCEDURE Lead of 99.999 pct purity was obtained from American Smelting and Refining Co. Cobalt sheet was obtained from Sherritt-Gordon Mines, Ltd., with a nominal purity of 99.9 pct, the principal impurities being nickel, iron, copper, carbon, and sulfur. The sheet was

Jan 1, 1969

By H. T. Kennedy, S. M. Avasthi

An equation developed for gaseous hydrocarbon mixtures predicts molal volumes with an average absolute deviation of 0.73 percent when applied to 264 natural gas and condensate systems including 2,043 PVT points. Another equation developed for liquid hydrocarbon mixtures predicts molal volumes with an average absolute deviation of 1.12 percent when applied to 346 crude oil systems including 1,759 PVT points. Both equations require composition of the mixture to be expressed as mole fraction of methane through heptanes-plus, hydrogen sulfide, nitrogen and carbon dioxide, together with the characteristics of the heptanes-plus fraction in addition to the temperature and pressure. The equations cover wide ranges of the variables involved, and their accuracy is considerably better than that of other available methods. The equations were differentiated to allow calculation of the coefficients of isothermal compressibility and isobaric thermal expansion. (In this paper the coefficient of isothermal com-pressibility and the coefficient of isobaric thermal expansion will be expressed as compressibility and thermal expansion coefficient, respectively.) Equations to calculate these quantities are presented. INTRODUCTION Calculations of reservoir performance for petroleum reservoirs require accurate knowledge of the volumetric behavior of hydrocarbon mixtures, both liquid and gaseous. Compressibilities are required in transient fluid flow problems, and thermal expansion coefficients are important in thermal methods of production. An accurate laboratory investigation of the PVT behavior of each reservoir fluid encountered would be costly and time consuming. For this reason various correlations for predicting fluid properties have been developed and recorded in recent literature. Correlations have been presented in the form of graphs, tables and equations. Since an increasing number of studies are being conducted with the aid of electronic computers, recent efforts have been directed toward development of correla-tions suitable for computer programming. Application of computers permits the use of more complex correlations which otherwise are not feasible. Moreover, methods for predicting reservoir performance, particularly those based on the compositional material balance, depend upon the capability of accurately expressing the molal volumes and other fluid properties as functions of pressure, temperature and composition. The coefficient of isothermal compressibility c is defined by

Jan 1, 1969

By L. G. Austin, R. R. Klimpel

This article demonstrates that the Gilvarry and Klimpel-Austin equations for the random fracture of solids are incorrect by deriving intuitively correct expressions for simple cases and showing that the foregoing equations do not reduce to right answers. It is the purpose of this article to point out serious errors in some existing treatments of the random fracture of solids. Gilvarry &apos; gives the equation where B(y,x) is the cumulative fraction by weight of particles which fall below size y when particles of size x are broken. The constants i, j, k contain shape factors and densities of flaws leading to propagating fracture surfaces. The authors believe this equation to be incorrect. Klimpel and usttin&apos; give where r~ is the number of fracture planes originating from edge flaws (per L where L is the edge length of the original particle), rs is the number of fracture planes originating from flaws on the particle surface (per S where S is the surface of the original particle) and rv is the number of fracture planes originating from flaws within the particle volume. This equation is also incorrect. The failures of these equations will be demonstrated by deriving intuitively correct expressions for simple cases and showing that the Eqs. 1 and 2 do not reduce to the correct answer. In this process, the flaws in statistical reasoning inherent in the derivations of Eqs. 1 and 2 become apparent. THE KLIMPEL-AUSTIN EQUATION Eq. 2 was derived by considering a set number of fracture planes occurring spatially at random within a given particle. This case will be treated here in its simplest form: a line with 1, 2, 3, etc., fractures occurring at random. This is illustrated in Fig. 1, with the original length of the line as Vo. The breakage of a long, thin, brittle fiber is a corresponding physical case, with the volume of a fragment proportional to its length. Assuming independence of the occurrence of the fractures, the probability of any one fracture occurring in length dV is dV/V0. The probability of the first fracture occurring at a length V to V + dV is dV/Vo. For the case of only one fracture per particle, formation of a number-size distribution of the left-

Jan 1, 1969

By G. R. St. Pierre, R. D. Blackburn

Oxygen contents of&apos; liquid-iron alloys saturated with oxides are analyzed. Minimum points observed in particular Fe-X-O systems can be predicled by a simple extension of first-order interaction methods. If minimum oxygen content occurs at N; where r is the ratio of metal to oxygen in the oxide phase. RECENTLY, wisher&apos; reported on solubility measurements of A1203 in liquid iron at 1580°C which support the earlier work of d&apos;Entremont, Guernsey, and chipman2 and Gokcen and Chipman. Swisher&apos;s diagram is reproduced in Fig. l. Minimum points in oxide solubilities have also been found in the Fe-Cr-0 system,4 Fig. 2, and the Fe-Ti-0 system,&apos; Fig. 3. Minima positions can be used to confirm approximately the measured first-order interactions, &, cgr, and E&~, in liquid iron and to deduce compositions of components of constant activity in the oxide phases. Consider the solution of Alz03 in liquid iron: for which an equilibrium constant can be written as: differentiating Eq. [2] with respect to Na~, the atom fraction of &l in the metallic solution, and noting that If the activities of Al and 9 are expanded by a first-order interaction model, then:

Jan 1, 1969