Search Documents

By E. G. Woods, A. G. Comer

A mathematical model is presented and solved for determinution of the saturation distribution and pressure in a radial gas-storage reservoir. The model consists basically of two parts: (I) the growing gas-bubble core, and (2) the surrounding aquifer. Since the total pressure at the injection well is a function of the two-phase flow in the gas bubble and the unsteady single-phase flow in the aquifer, the resistance to flow in both zones was taken into consideration. The assumptions involved for both the radial equivalent of the Buckley-Leverett two-phase-flow equation and the injection-pressure equation are as follows: (I) the geometry is radial, (2) the gas bubble is free to expand or contract, (3) compression or expansion of the gas within the bubble may occur at the beginning of a time step, (4) the fluids are immiscible, (5) water is incompressible within the gas-storage region whereas it is compressible outside of this region, (6) a stable gas-water interface exists and (7) gas injection occurs at a constant rate or a series of constant rates. The mathematical model was solved numerically using an IBM 650 computer. A comparison is presented between the predicted results of the model, the results assuming steady-state flow and the actual initial injection-pressure history of an operating reservoir. Using the initial field pressure for a basis, the average deviation between the predicted pressure and the actual field pressure was less than 4.3 per cent. INTRODUCTION Underground storage of natural gas in abandoned oil fields, in abandoned coal mines, in caverns and in aquifers has had varying degrees of success. The lack of these first three facilities in the vicinity of most major gas marketing areas is leading more and more to storage in virgin aquifers. For the most part, this storage in aquifers has led to the study upon which this paper is based. Several excellent papers have been published discussing the problems involved in the underground storage of natural gas; however, to these writers' knowledge, none of them considered the two-phase flow of fluids in the gas bubble and the fact that the gas bubble will grow or shrink in size depending upon the injection and withdrawal history of the reservoir. Not only will both of these factors influence the required injection pressure, but also they will have a definite bearing upon the amount of water production occurring upon withdrawal from the storage area. To solve this problem on an IBM 650 computer, several simplifying assumptions were made. Only the drainage portion of the relative-permeability curve was considered, thereby neglecting any hysteresis effects introduced during an imbibition cycle. The quantity of gas withdrawn during the period of the field study was only a small percentage of the total gas in place; therefore, this assumption is justified. For the problem studied in this case, the combined effects of capillarity and gravity were assumed to be negligible. This assumption becomes less valid as the formation thickness or the average pore diameter increases markedly. The justification of these assumptions should be studied in each case. THEORY This mathematical model for the prediction of saturation distribution and pressure for a radial gas-storage aquifer is based upon the equations for radial two-phase fluid flow'.' and for radial, unsteady-state, single-phase fluid flow.' The two-phase flow is considered to take place in a "core" whose radius is equal to the maximum radius that the gas zone will attain, as indicated in Fig. 1. The injection-withdrawal history is approximated by a series of constant flow rates. Within the "core", gas is assumed to behave as a semi-compressible fluid; that is, the gas is assumed to have a constant density based on the average gas-zone pressure for the flow period. Between each constant-rate time increment, the gas density is allowed to change. This involves an iterative solution of the pressure equation. The liquid phase in this zone is considered to be incompressible, whereas it behaves as a compressible liquid outside of this region. The equations of incompressible flow in the "core" may be written for the gas phase

By R. Stickler, A. Vinckier

An optical and electron microscope study was made on a commercial 302 stainless steel heat treated for massive carbide precipitation. Convincirg evidence was obtained that the grain boundary precipitate does not grow away from the grain boundary into the matrix, but grows as a network of large thin dendrites in between the grains. IN austenitic stainless steels the carbon in supersaturated solution precipitates when the steel is heated in the temperature range 450" to 1000°C. During such heat treatments the chromium-rich carbide occurs primarily at the grain boundaries and has adverse effects on the corrosion resistance and low-temperature ductility of the material. It has been shown that the morphology of the precipitated particles differs widely depending on the time and temperature of the heat treatment, the nature of the boundary, the misfit between the grains, and so forth.1"5 It was assumed by Mahla et al.,' Streicher,' and Mc Nutt' that the carbide particles which nucleate in the grain boundary grow in various geometric or dendritic forms away from the boundary into the matrix. However, Kinzel,' Hatwell et a, and Plateau et aL4 found that the particles nucleate in the grain boundary and grow as a network of thin dendrites in the grain boundary interface. We have made a comprehensive studyg of the precipitation occurring in a commercial 302 austenitic stainless steel heat treated for times ranging from 0.15 to 1500 hr at various temperatures from 480" to 1065"~. This study was made to elucidate the influence of heat treatments on both the mechanical and the corrosion properties of this steel. Part of this investigation which is pertinent to the morphology of the grain boundary carbides referred to above is separately reported here. We found that the particles which nucleated in the boundaries grew in the interface of these boundaries, thus supporting the viewpoints of Kinzel,' Hatwell et al.,3 and Plateau et al.4 Furthermore, when carbide particles nucleated in the matrix very close to the boundaries, they possessed a morphology distinctively different from the particles growing in the grain boundaries. MATERIAL A study of metallographic samples and of fracture surfaces of a 302 stainless steel is reported in this paper. The composition and heat treatments are listed in Table I. Polished and etched samples were examined und813r the optical microscope, and carbon extraction replicas of etched micrographic samples and various fracture surfaces3'8'9 were examined under the electron microscope. RESULTS AND DISCUSSION A representative optical micrograph of the 302 steel, condition A, is shown in Fig. l(a). Particles can be seen on the grain boundaries but no traces can be found of particles which grow from the grain boundaries into the matrix. A carbon extraction replica of this sample shows numerous large thin dendrites, Fig. l(b). The lateral dimension of such carbide particles vary up to 100 p, and their thickness estimated from the amoun! of electron penetration i:; between 500 and 1000A. Such large dendrites are extracted from almost all grain boundaries and, if they had grown into the grains, they would definitely be revealed by etching traces in the matrix, Fig. l(a). The reason why these large dendrites on the extraction replica appear to be grown into the grain is due to the mechanism of the extraction replica process. Thus, these dendrites are supported by the carbon replica film only at the edges originally along the grain boundary trace A-A, Fig. l(b), in the polished metal section and have subsequently fallen over onto the carbon replica film during handling. Although such dendrites generally fall only to one particular side of the grain boundary trace, one finds occasionally that parts of a dendrite have fallen to both sides, as shown in Fig. l(b). Further evidence for the growth of carbide in the grain boundary can be obtained from the appearance of fracture surfaces. A small impact specimen of the 302 steel, condition A, was broken at liquid nitrogen temperature. The fracture path follows the

Jan 1, 1962

By Josiah Royce

To evaluate chemically the sample of rock obtained by diamond drilling, it has long been recognized that the analyses of the two components of the sample, core and sludge, must be given appropriate influence in computing the average analysis of any unit of depth. The purpose of this investigation is to set forth what means are available for apportioning the effect of core and sludge on the final analysis, what variables affect the problem, and what combination of applied mathematics will closest approximate the truth under each condition as these variables proceed within their limitations. A drill hole is bored in iron ore exploration principally to test variations in rock composition with depth and is usually directed as nearly normal to the bedding of a horizon to be tested as possible. This practice has a tendency to minimize variation in composition laterally which in any event is not likely to be great. It is obvious that the opportunity for change in analysis of a particular rock is not statistically as great radially in a diamond drill hole where the distance in which such a change may occur is from 0.719 in. (EX bit) to 1.469 in. (NX bit) as there would be longitudinally even in a' run as short as 5 ft. Variations in composition of bedded or layered rocks are usually greater normal to the bedding than parallel thereto. Even in massive rocks, like porphyries, variations are functions of distance. Hence. in either case, variations along the hole are of greater effect than across it. Let us examine Fig 1 briefly to observe a cross section of a unit of depth of a typical diamond drill hole. If we neglect radial change in chemical composition, which we have seen is small compared to lengthwise variation, we can see that if core recovery was 100 pct, the core and its surrounding area, which would be recovered as cuttings (assuming 100 pct sludge recovery), must analyze the same if the sampling is perfect. This is the foundation of this paper, namely, that if core recovery is complete we may assume that the core constitutes as nearly as possible a perfect sample of the ground drilled. Now let us pass to the methods of weighting core and sludge in the final analysis and then develop each in turn. One might employ core analyses only, sludge analyses only, or apportion the influence of each by one of several methods. The most common apportionment is by direct proportion to relative theoretical vol- ' ume occupied by each component in the cylinder hollowed out by the drill, known as the Longyear formula. Another logical treatment would be to consider core and sludge in the final analysis according to the weight of each recovered. This procedure was elaborately set forth and refined by R. S. Moeh1man.l Then to be sure that we have embraced all the possibilities, we must admit that some other method or methods may be devised empirically or by mathematical maneuvers. Should a drill machine be capable of recovering 100 pct of the core in all types of rock, certainly core alone need be analyzed and no complicated mathematics are required; but although some companies have had fair success in obtaining high percentages of core in certain homogeneous ores on the Mar-quette Range, their same methods have proved disappointing on the Menominee, Cuyuna, and Gogebic Ranges where the iron formation often varies widely from flinty chert beds interstratified with soft hematite to wholly leached ore material composed of an unpredictable mixture of hematite and limonite with occasional zones of sugary, recrystallized chert. In spite of numerous mechanical improvements in core barrels, few drilling programs can count on complete core recovery. The sludge alone might be analyzed but for practical reasons this is inadvisable when core is obtained due to the numerous opportunities for contaminating or losing portions of the cuttings. The Longyear formula involves a weighting of the core and sludge in proportion to the theoretical volume of each, that is, where the volume of the core recovered is one-third of the total cubic volume of the cylindrical hole made by the drill, the analysis of the core is given one-third of the weight in the final analysis. This approach to

Jan 1, 1950

By V. Koump, T. F. Perzak, R. G. Olsson

The rate of dissolution of recrystallized alumina in molten iron oxide in equilibrium with iron at 1450°C was studied by rotating alumina disks in the melt. The samples were rotated from 1 to 12 min at speeds from 54 to 3270 rpm. It is concluded from the experimental results that the process is limited by diffusion in the liquid boundary layer. The interdiffusion coefficient for the dissolution of alumina disks in pure molten iron oxide was estimated to be about 3 xl0-5= sq cm per sec. THE erosion of refractory materials is of considerable importance in iron- and steelmaking processes. Since many refractories contain alumina, this investigation was undertaken to develop a further understanding of the erosion mechanism. Shurygin et al.1-3 have investigated the rate of dissolution of alumina in molten silicates and molten fluorides by rotating alumina disks in various melts. The rates were reported to be limited by diffusion in the liquid boundary layer. In the present investigation the rate of dissolution of alumina in molten iron oxide in equilibrium with iron was studied in a similar manner; alumina disks were rotated at different speeds in molten iron oxide at 1450°C. The relationship between the rate of dissolution of the flat surface of the disk and the speed of rotation was used to determine whether the process is limited by diffusion in the liquid boundary layer or by some chemical reaction mechanism. APPARATUS AND EXPERIMENTAL PROCEDURE A schematic diagram of the experimental apparatus is shown in Fig. 1. The iron oxide melts were contained in 2-in.-diam high-purity iron crucibles that were centered in a cylindrical graphite susceptor. Prior to each sequence of experiments, sufficient reagent-grade ferric oxide was added to a crucible to form approximately 350 g of molten iron oxide; the depth of the melt was about 11/2 in. To prevent excessive erosion of the crucible, about 70 pct of the required iron was added to the crucible in the form of iron powder. The mixture was melted in an atmosphere of purified argon by induction heating with a 250-kc generator. After the oxide was molten, a Pt/Pt-10 pct Rh thermocouple at the exterior of the crucible was calibrated against an iron-sheathed thermocouple immersed in the melt and was thereafter used to regulate the melt temperature. All experiments were conducted at 1450°C with a probable error of 10°C. The alumina disks were constructed from pure re-crystallized alumina having a density of 3.7 g per cu cm (porosity approximately 6 pct). For each disk sample a 1-cm-diam alumina cylinder of known length (0.9 to 1.0 cm) was cemented into a ground seat at the end of a 1.27-cm-OD centerless ground alumina tube. The ends of the tube and cylinder were flush and served as the disk surface. When mounted in the apparatus the alumina tube was guided by tungsten carbide bearings and rotated by a variable-speed motor. The bearing and drive assembly could be raised and lowered in order to allow the sample to be quickly moved in and out of the melt. At the start of an experiment, a sample was preheated for several minutes directly above the melt, and then immersed in the middle of the melt to a depth of a in. and rotated for 1 to 12 min at speeds of 64 to 3270 rpm. The speed of rotation was measured by a hand tachometer. At the start and completion of each experiment, a melt sample was withdrawn with a cold copper rod and the alumina content was determined by chemical analysis. During individual experiments the alumina concentration of the melt changed by the order of 0.25 wt pct. Average melt compositions for individual experiments were used in the subsequent computations and were in the range between 0.12 and 4.3 wt pct alumina. The extent of dissolution was determined by sectioning the end of the sample and measuring the final length of the alumina cylinder with a microscope. Since the exposed end of the cylinder, i.e., the disk surface, was etched in a manner characteristic of this geometry, a number of measurements were made on

Jan 1, 1969

By Byoung Whie Lee

In order to investigate the change in surface energy as a function of cold work ad environment, the contact angles between a liquid and the surface of elec-tropolislzed aluminum with various degrees of cold work were measured. Specinlens cut from aluminum rod (99.99 pct Al) were annealed at 600°C for 24 hr and then cold compressed parallel to the drawing di-rectzon of the rod in increments of approximately 10 pcl from 0 to 70 pct. The sections, taken transversely to the drawing direction, were mechanically and electro chemically polished. Contact angles between the sessile drop of distilled water, glycerin, and glycol and the specimen surface were measured optically. The differences in interfacial energy were calculated from the contact angle data. The effects of chemical washing and surface roughness of the specimen on the contact angle were also studied. The liquid contact angle was a linear function of the degree of cold working (D) ad changed markedly with differed liquids. THE plastic deformation of polycrystalline metals induces a preferred orientation in crystal structure. A preferred orientation changes the atomic density of a metal surface. This change in surface atomic density brings about the changes in surface interactions of the metals."' Also, the plastic deformation of metals raises the internal energy by increasing the stored energy.' Consequently, the total Helmholtz free energy of the metal increases and, in turn, the surface free energy of the metal. In order to study the effect of plastic deformation on solid surface phenomena, the liquid contact angle of a sessile drop made on cold compressed and elec-trochemically polished aluminum was measured. The preferred orientation in aluminum, due to uniaxial compression, is known to be [110] direction parallel to the compression axis. However, this description is only a first approximation of a rather complex orientation. In actuality, half the crystallites are more than 10 deg from this orientation, regardless of the amount of compression." For randomly oriented fcc crystals, the average surface atomic density can be calculated as follows: The number of atoms in a unit volume is 4/a3, where a is a lattice constant. Hence, the number of atoms per unit surface area is about (4/a3)2/3 = 2.52/a2 However, when the preferred orientation took place along [110] direction, the surface atomic density is as follows: The area of (110) plane of a unit cell, considering ridges and valleys, is BYOUNG WHlE LEE, Member AIME, is Research Metallurgist, Atomic Energy Research Institute of the Republic of Korea, Seoul, Korea. Manuscript submitted May 29, 1968. IMD 2 x J2- a x Ga = 6a2 Since this plane contains 4 atoms, the surface atomic density is 4i6a2 = 2.3/a2 Upon changing from random to preferred orientations, the surface atomic density decreases about 10 pct. The increase in the stored energy and the change in surface atomic density due to plastic deformation would certainly alter the surface tension of a metal. EXPERIMENT High purity aluminum rod of 6-mm diam (99.99 pct Al, obtained from Johnson, Matthey & Co., Ltd, London, England) was cut into sections approximately 1 cm long. The transverse cross section was mechanically polished. Then, the specimens were annealed in air at 600°C for 24 hr. Thus, eight annealed specimens were uniaxially compressed at room temperature from 0 to 70 pct of cold working with an increment of 10 pct. The extent of cold working was determined by the percentage of reduction in axial length to the original length. The specimens were further polished mechanically and, then, electrochemically with the Modern Battelle Type 1 electropolishing solution with a current density of 0.3 to 0.4 amp per sq cm and an electrode potential of 8 to 12 v at the solution temperature of 78" to 82°C for approxirr~ately 5 min. The Modern Whlte Diffusing Screen \t/ --------------- re Goniometer for Mercury Lamp ^-^ Sessile Drop ¦=• Condenser c=> Microscope l_*~_i bJ Mlcromster Filer with Protmcter Fig. 1—Schematic diagram of apparatus set-up for liquid contact angle measurement.

Jan 1, 1970

By T. A. Pollard, C. R. Sandberg, E. B. Elfrink

This paper presents an evaluation of compressibility factor data and a discussion of their application to the estimation of gas reserves. A correlation is presented which provides compressibility factors for use in both two-phase and single-phase hydrocarbon systems. Accuracies comparable to those obtained previously for single-phase systems only can be expected. A simple means of predicting the presence or absence of a liquid phase in a condensate system of known composition is illustrated. The correlation is based on 1,030 compressibility determinations from 21 hydrocarbon samples taken from eight oil fields. Of the data used, 75 per cent were from California, 15 per cent were from the Mid-Continent area, and 10 per cent were from South America. The average numerical deviation of the experimental data from this compressibility chart is 1.22 per cent. Charts and tables are included and discussed which illustrate the errors involved through the misuse or nonuse of compressibility factors in estimates of gas-in-place. INTRODUCTION The compressibility factor is a coefficient which expresses the deviation of a gas of given composition from the Perfect Gas Laws. A factor such as this is necessary in any calculation involving volumes of gaseous mixtures, and finds extensive application in the estimation of natural gas and condensate reserves. It is used in the decline curve, the volumetric, and the material balance methods of gas reserves estimation. The behavior of gases at high pressures has been investigated extensively in recent years and considerable data have been assembled on the compressibility characteristics of a number of gases. Notable among the investigators of high pressure gas relationships are Kvalnes and Caddya, Kay" Sage and Laceys, and Brown, Standing" and Katz4. Other contributors are Smith and Watson, Roland and Kaveler', and Stevens and Vance". The data on compressibility factors assembled by these investigators have been presented and correlated in numerous ways; particularly noteworthy is the method first presented by Kay5 which relates compressibility factor "Z" to pseudo-reduced temperature and Pseudo-reduced pressure. Other investigators have used this empirical relationship successfully and numerous charts covering a large number of gases have been constructed on this basis. It is the purpose of this paper to evaluate compressibility factor data and to illustrate the importance of their proper use in the estimation of gas-in-place for both gas and condensate reservoirs. EVALUATION OF COMPRESSIBILITY FACTOR DATA The compressibility chart presented by Standing, Katz, Brown and Holcomb in 1941 is probably the most recent and reliable chart for natural gases. The reliability of the Katz chart has been investigated by the writers for a fairly large number of gases (single-phase) in the pressure range from 500 to 10,000 psia and temperature range from 80 to 400°F. A summary of this comparison is presented in Table I. The numerical average deviation of the chart compressibility factors from the experimental compressibility factors is 1.05 per cent, and the algebraic deviation is -0.28 per cent. The Katz correlation appears to be quite accurate for single-phase hydrocarbon systems. It was based on single-phase systems and does not necessarily provide a means of accurately predicting the behavior of two-phase systems which might be encountered in condensate work. Data published recently by Sage and Oldsl' have enlarged the knowledge of phase behavior in condensate systems by covering extensively the Compressibilities of systems in the two-phase as well as single-phase regions for a wide range of temperatures, pressures and gas-oil ratios. An evaluation of these data indicates that conditions can arise, especially in condensate systems having fairly low gas-oil ratios, in which two phases occur and the usual gaseous compressibility factors do not accurately apply. The measured PVT data of Sage and Olds covering two-phase systems have been correlated with additional such data from the Mid-Continent area and South America; these are presented in Fig. 1. This chart is the result of 1,030 compressibility* determinations from 21 hydrocarbon samples taken from eight oil and condensate fields. Of the data used, 74.75 per cent were from California, 15.45 per cent were from the Mid-Continent area, and 9.80 per cent were from South America. The chart construction is based on the original

Jan 1, 1949

By R. J. Wagner, G. H. Gronhovd, A. J. Wittmier

The paper describes the results of a series of full-scale boiler tests run to determine the ash fouling characteristics of a North Dakota lignite as a function of sodium content of the coal. Four levels of sodium were investigated and rate of fouling was determined both by boiler performance and by specially constructed probe equipment. Discussed is information on chemical and physical properties of the ash deposits. The vast lignite deposits in North Dakota are becoming increasingly important as a source of low-cost fuel for power generation. A new 200-megawatt lignite-burning plant located near Stanton, N.D., was placed in operation early in 1966. Another 172-mega-watt unit in the same area is scheduled for operation late in 1966, and plans have recently been announced for a 212-megawatt unit to be in operation by 1970. If these plants show good performance, western North Dakota could become a major power production center for a large geographic area in the upper midwest. One of the problems associated with the use of lignite in power boilers is its tendency to form troublesome ash deposits on boiler tube surfaces. This problem varies considerably from one lignite-fired boiler to another and is dependent upon such factors as source of lignite, boiler design, load factor, tube metal temperatures, and soot-blowing practices. This report presents data on sodium content* variation in a North Dakota mine, describes the results of a series of tests performed to determine the relationship between ash fouling tendencies and the sodium content, and tells how the mining company and electric utility have coordinated their operations to reduce the ash fouling problem. GENERAL BACKGROUND Until recent years it was generally assumed that North Dakota lignite from a given mine had qualities peculiar to that mine, but that within the mine, the lignite was quite uniform. In 1959, the Otter Tail Power Co. installed a 53.5-megawatt lignite-burning unit at its Hoot Lake plant, Fergus Falls, Minn. Shortly after the unit was put into operation, it was noted that at certain times under high load conditions, ash fouling problems were quite severe, but not consistent or predictable. Each time fouling occurred in the unit the position of the loading shovel at the mine was noted. Within two years it became apparent that there were certain areas in the Knife River Coal Mining Co. mine at Beulah, N.D., where the lignite apparently had a high fouling tendency for the Hoot Lake unit. Results from analyses of the coal did not show large differences in moisture percentage, ash percentage, heating value, or ash fusion temperature between the highly troublesome and less troublesome coals. After three years of operation, the troublesome areas had been quite well defined and the mining company started a procedure of bypassing these areas when load conditions at the plant were high and then returning to the troublesome areas to load the lignite during the summer months when peak loads were lower. It was noted that the lignite from the areas which caused severe fouling under high load conditions produced much less fouling at the lower summer loads. However, bypassing certain areas in the strip pit caused the mining company to disrupt their continuous loading schedule, which in turn upset their continuous stripping schedule and made for less efficient mining operation. Meanwhile, at the Hoot Lake plant, the number of soot blowers on No. 2 boiler had been increased from 30 to 55, and a very rigorous blowing schedule had been adopted resulting in soot blowers operating at least some part of 21 hrs per day. Fouling problems continued, with a resulting reduction in plant economy and reliability. During this period the Federal Bureau of Mines Coal Research Laboratory at Grand Forks, N.D., had initiated various research programs relating to the study of lignite ash and ash fouling problems. One project was a survey in which full seam samples were obtained from various locations in all the major

Jan 1, 1968

By D. O. Hobson, P. L. Rittenhouse

We have recently found that many lots of Zircaloy-2 and -4 tubing tested in tension deform to cross section shapes that are "polyhedral" rather than circular. The noncircular cross sections seem to result from a circumferential variation in plastic proper -ties and, therefore, in texture. We examined the diffracted X-ray intensity from (OOOl), {1010), {ll20}, and {1011} planes parallel to the tubing surface around the entire circumference. Tubing that deformed to a circular cross section showed <5 pct variation in diffraction intensity around the circumference from any of these planes. Greater intensity variations were found for all the tubing materials that developed non-circular cross sections during deformation. In each instance the number of intensity maxima and minima was consistent with the number of sides of the &apos;polyhedral" cross section. THE effects of crystallographic texture on the mechanical properties of zirconium, zirconium alloys, and other anisotropic hexagonal materials such as titanium and hafnium have been intensively studied. If the natural crystallographic anisotropy of these materials could be controlled by selected fabrication operations, a material could be produced with a texture, or preferred orientation, in which the maximum strength directions are in the stress directions of the structure. The permissible stress levels of the structure could then be increased. Rittenhouse and Pickle-simerl have shown that highly textured Zircaloy-2 plate can vary from 52,000 to 70,000 psi in tensile yield and from 62,700 to 122,000 psi in compressive yield strengths as a function of specimen orientation. Similar texture effects have been observed in Zircaloy tubing, whose mechanical properties differ depending upon whether the tubing has basal poles predominantly radial or tangential and upon the mode of testing. These same texture effects might be troublesome if the texture in the tubing is not uniform around the tube walls. This paper is concerned with the effects of nonuniform texture on the strain behavior of commercial Zircaloy tubing. BACKGROUND Previously, various lots of commercial Zircaloy tubing were tested under selected stress conditions as part of a program to determine the four quadrants of the yield locus. One of the testing procedures was uniaxial tension and it is the influence of texture variations on the tensile deformation behavior of the tubing that is of interest here. The effect of texture on the strain behavior in tubing is very apparent from a uni-axial tensile test.2 This behavior is illustrated schematically in Fig. 1 for three different tubing textures: 1) basal poles predominantly in the radial direction; 2) basal poles predominantly in the tangential direction, and 3) basal poles equally distributed between the radial and tangential directions. In the first example most of the grains are oriented for prism slip to operate to allow axial lengthening and diametral reduction of the tube. Reduction in wall thickness requires {1122} twinning, a more difficult operation than slip.* In the second example, the texture is such that *The singleCrystal analogy will be used quite often in this paper when strong texturesare being discussed. The textured material is assumed to react as if it were a single crystal aligned in the same orientation and under the same stress system as the predominant texture._____________________________________________ prism slip operates to give axial lengthening and wall thinning. Diametral reduction requires twinning and again should be more difficult to accomplish. A specimen with basal poles equally distributed between the radial and tangential directions should exhibit about equal diametral and wall strains.* Table I shows the *Diametral strain is accomplished by a combination of circumferential and longitudinal strains. But, for the purposes of this paper, it will be considered as a strain causing a change in the length of a specified diameter.__________ strain behavior of two lots of tubing with different textures. The first tube has a texture with basal poles equally distributed between the radial and the tangential directions. The ratio of true thickness strain to true diametral strain is almost unity. The second tube has a predominantly radial basal pole texture and it is apparent that most of the total strain has occurred by diameter reduction, just as one would predict. There is over three times as much diametral strain as thickness strain. RESULTS We found that specimens from certain lots of commercially fabricated tubing did not retain round cross sections during uniaxial testing. Instead, the speci-

Jan 1, 1970

By C. R. Sandberg, T. A. Pollard, E. B. Elfrink

This paper presents an evaluation of compressibility factor data and a discussion of their application to the estimation of gas reserves. A correlation is presented which provides compressibility factors for use in both two-phase and single-phase hydrocarbon systems. Accuracies comparable to those obtained previously for single-phase systems only can be expected. A simple means of predicting the presence or absence of a liquid phase in a condensate system of known composition is illustrated. The correlation is based on 1,030 compressibility determinations from 21 hydrocarbon samples taken from eight oil fields. Of the data used, 75 per cent were from California, 15 per cent were from the Mid-Continent area, and 10 per cent were from South America. The average numerical deviation of the experimental data from this compressibility chart is 1.22 per cent. Charts and tables are included and discussed which illustrate the errors involved through the misuse or nonuse of compressibility factors in estimates of gas-in-place. INTRODUCTION The compressibility factor is a coefficient which expresses the deviation of a gas of given composition from the Perfect Gas Laws. A factor such as this is necessary in any calculation involving volumes of gaseous mixtures, and finds extensive application in the estimation of natural gas and condensate reserves. It is used in the decline curve, the volumetric, and the material balance methods of gas reserves estimation. The behavior of gases at high pressures has been investigated extensively in recent years and considerable data have been assembled on the compressibility characteristics of a number of gases. Notable among the investigators of high pressure gas relationships are Kvalnes and Caddya, Kay" Sage and Laceys, and Brown, Standing" and Katz4. Other contributors are Smith and Watson, Roland and Kaveler&apos;, and Stevens and Vance". The data on compressibility factors assembled by these investigators have been presented and correlated in numerous ways; particularly noteworthy is the method first presented by Kay5 which relates compressibility factor "Z" to pseudo-reduced temperature and Pseudo-reduced pressure. Other investigators have used this empirical relationship successfully and numerous charts covering a large number of gases have been constructed on this basis. It is the purpose of this paper to evaluate compressibility factor data and to illustrate the importance of their proper use in the estimation of gas-in-place for both gas and condensate reservoirs. EVALUATION OF COMPRESSIBILITY FACTOR DATA The compressibility chart presented by Standing, Katz, Brown and Holcomb in 1941 is probably the most recent and reliable chart for natural gases. The reliability of the Katz chart has been investigated by the writers for a fairly large number of gases (single-phase) in the pressure range from 500 to 10,000 psia and temperature range from 80 to 400°F. A summary of this comparison is presented in Table I. The numerical average deviation of the chart compressibility factors from the experimental compressibility factors is 1.05 per cent, and the algebraic deviation is -0.28 per cent. The Katz correlation appears to be quite accurate for single-phase hydrocarbon systems. It was based on single-phase systems and does not necessarily provide a means of accurately predicting the behavior of two-phase systems which might be encountered in condensate work. Data published recently by Sage and Oldsl&apos; have enlarged the knowledge of phase behavior in condensate systems by covering extensively the Compressibilities of systems in the two-phase as well as single-phase regions for a wide range of temperatures, pressures and gas-oil ratios. An evaluation of these data indicates that conditions can arise, especially in condensate systems having fairly low gas-oil ratios, in which two phases occur and the usual gaseous compressibility factors do not accurately apply. The measured PVT data of Sage and Olds covering two-phase systems have been correlated with additional such data from the Mid-Continent area and South America; these are presented in Fig. 1. This chart is the result of 1,030 compressibility* determinations from 21 hydrocarbon samples taken from eight oil and condensate fields. Of the data used, 74.75 per cent were from California, 15.45 per cent were from the Mid-Continent area, and 9.80 per cent were from South America. The chart construction is based on the original

Jan 1, 1949

By D. W. Lillie

Property changes produced by irradiation of Al-0.4 wt pet Li alloys at 270°C to a burnup of 0.155 pct of all atoms are described. Metallographic evidence of the formation of internal pores and the concentration of gases at grain boundaries is shown. Some conclusions are drawn concerning desirable property characteristics for resistance to radiation swelling. RADIATION induced gases are of great practical significance in nuclear engineering since they affect the behavior of fissionable fuels and control rod materials. In fuels the gases xenon and krypton are the major offenders, producing serious em-brittlement at low irradiation temperatures and gross swelling at high irradiation temperatures. In control rod materials containing boron, the production of helium from the B 10 (n, 0) Li7 reaction produces similar reactions. Since the trend in reactor design is for operation at increasingly high temperatures, the phenomenon of high-temperature swelling by these induced gases has been the subject of considerable recent study. The objectives of such work have been to obtain engineering information and quantitative data on the effects in specific fuels and control rods and to obtain further understanding of the mechanism of gas damage so that more intelligent corrective steps can be made. Analogue systems are often useful for the latter objective and in this case alloys containing Li are suitable in view of the reaction Li6 + n - He4 + He 4. Both products are gases and no serious radiation problem is produced to complicate post irradiation handling. Considerable caution is, of course, necessary in any operations in which tritium release is possible. The present program was initiated to provide a detailed study of the effects of He + H3 generated by irradiation in Li-bearing alloys. The results on Mg-Li solid solution alloys have been reported elsewhere.&apos; This paper presents information on irradiation of an A1-0.4 pct Li alloy in which the lithium was present as highly enriched Li6 in order to maximize the radiation burnup. PROCEDURE The original intent in the selection of the alloy content was to have Li in solid solution at all temperatures. It became apparent, however, that the solubility of Li in A1 was substantially lower than that reported in the original reference consulted.&apos; This was substantiated by the work of Nowak3 who reported solubility of less than 0.01 pct Li at 100°C. The alloy used (0.4 pct Li) was therefore two phase at room temperature. The irradiation temperature was estimated to be about 270° C, and at that temperature all or almost all of the Li would be soluble under equilibrium conditions. The precise state of the alloy under irradiation is not known. Aluminum used was 99.996 pct purity (max. impurities Si 0.002 pct, Mg 0.001 pct, Fe 0.0005 pct, and Cu 0.0005 pct). The lithium was obtained from Oak Ridge National Laboratory; 96.1 + 0.1 pct of the lithium was ~i~, but chemical impurities were 0.25 pct Ca, 0.05 pct Fe, 0.02 pct Na, 0.02 pct Cu, 0.01 pct Sr, and 0.01 pct Ba. The two metals were melted in an alumina crucible in a resistance heated furnace in an atmosphere of helium. A l-in. ingot was cast in an iron mold, machined to %-in. round, hot swaged (350°C) to 0.080-in. round, cold swaged to 0.060-in. round, and finally cold drawn to 0.040-in. round. Four-inch wire lengths were then cut and annealed at 250°C for 4 hr in helium. Irradiations were carried out in a VT hole in the Argonne CP-5 reactor. The wires were held in an aluminum jig, sixteen wires per capsule, and were maintained in a helium atmosphere for maximum heat transfer. Eight wires of pure aluminum and eight of Al-Li alloy were placed in each capsule. Argome Laboratory statements were accepted as establishing flux and exposure levels. These were based on previous experience and flux map-

Jan 1, 1961

By L. W. Fagg

The purpose of this paper is to present in a convenient form data and examples necessary in making dynamometer card analyses; also to outline a procedure of well weighing. Many articles and papers have been written delving into the mathematical considerations relative to the shape and characteristics of dynamometer cards. However. it is recognized that there are too many unknown factors involved in such calculations to assure a workable degree of accuracy. For this reason the accepted procedure is to take dynamometer cards on wells in question rather than try to calculate the load curve. The polished rod dynamometer is now recognized as a necessary tool for measuring loads, torque, and horsepower. It is also used to determine pump action and trouble-shoot for any seemingly abnormal pumping condition. The apparently infinite variety of a = maximum load (height x scale constant) b = minimum load. Range of load is difference between maximum and minimum load speed—taken with stop watch stroke —measured at polished rod c = beginning of down stroke in direction of arrow. (End of up stroke) d = beginning of up stroke. (End of down stroke). Polished Rod Horsepower = (Area of card) x Scale const. x Stroke x Length x spm (Length of card) 33000x 12 1.35 1.35------x5300x 24 x 12 2.53 --------------------------= 2.06 23000x12 Appraximate Peak Torque: Upstroke = (2390) (12) (1) = 28,700 in. Ib Downstroke = (1 860) (12) (.866) = 19,300 in. Ib Counterbalance should be increased to make up stroke and down stroke peak torque equal. FIG. 1 dynamometer cards that can be obtained is one reason for the general lack of usage of the dynamometer as a control instrument rather than a means for making routine measurements of leads and horsepower. When it is considered that the dynamometer card is a record of the resultant of all forces acting on the polished rod at any particular instant during the pumping stroke. the problem is then one of breaking down this resultant into its various components. As a means of a quick review we shall consider the examples shown in Figs. 1 to 9: and Tables I and II and then proceed to the interpretation of variously shaped cards caused by some abnormal operating condition. In Table 11, when we were considering the factors involved in calculating the peak polished rod load, it can be seen that the factors involved greatly oversimplify the problem. Certain assumptions are made which may or may not he even close to the actual field conditions, such as the specific gravity of the fluid generally considered as one; that the crank has constant angular velocity; that the down-hole friction is zero; and that the fluid lift is from the pump. In the following examples we shall see what a variation in fluid weight and friction can do to the general shape and magnitude of the dynamometer card. (Figs. 10 to 22.) MAKING THE WELL STUDY It is obvious that it would be impractical to consider in detail all of these factors each time a well study is made, inasmuch as each well study job could conceivably be extended into a research project rather than serve the practical requirements of finding the answer to a specific problem. For this reason it is important that some objective be established previous to the time the well study is made. Load at TU = 3070 Ib Crank angle when polished rod is at position TU = ? ? = 30° Maximum counterbalance effect at polished rod = 2760 Ib Torque at TU = (Load at TU — Max. counterbalance effect-lbs) sin 0 x Length of Stroke 2 = (3070-2760) .5 x 24 = 1860 in. 15 2 Torque at TD = (Load at TD — Max. counterbalance effect-lbs) sin ? x Length of Stroke 2 = (530 - 2760) sin 330&apos; x 24 = 13,400 in. Ib 2 Note: sin 0 from c to d on UPSTROKE will be positive value. sin 0 from c to d on DOWNSTROKE will be negative value Torque at c and d is zero because 0 is zero. FIG. 2 —APPROXIMATE METHOD FOR CALCUlATlNG TORQUE

Jan 1, 1950

By A. E. Dwight

Lattice parameters were determined for eighteen equiatornic alloys of the CsCl-type structure, ten of which were previously un-reported. It was found that fomation of the CsCl-type structure in binary alloys of the transition elements is largely dependent on position of the elements in the periodic table. The relative size of the two elements was not found to be a controlling factor. A recent paper by Beck, Darby, and Arora1 corre-lates the occurrence of CsC1-type ordered structures with the position of the constituent elements in the periodic table for the first long period. It was also suggested that a definite increase in relative bond strength between unlike atoms occurred when, in binary alloys of iron-group elements, the other component is changed from a chromium-group element, to a vanadium-group element, to titanium. A later paper by Philip and Beck2 noted that the lattice contraction increased in the order CrFe, VFe, and TiFe. It was also noted by Philip and Beck2 that the lattice contractions of CsC1-type alloys decreased in the order: TiFe, TiCo, and TiNi, which is an apparent reversal of the contractions expected from the position in the periodic table. It was suggested that the increasing lattice contraction is an indication of increased stability, i.e., greater A-Bbond strength. The present investigation was carried out to determine whether the relation of the position in the periodic table to the formation of the CsC1-type structure was also correct for alloys involving the second and third long-period elements. A systematic search was made for CsC1-type structures among equiatomic alloys and for those found, the lattice contraction was determined. EXPERIMENTAL TECHNIQUE The elements Y, Gd, Ti, Zr, Hf, V, and Cb are designated the A group and the elements Mn, Re, Fe, Ru, Os, Co, Rh, Ir, Ni, Pd, Pt, Cu, Ag, and Au are designated the B group. Equiatomic alloys were prepared for 57 AB combinations. The alloys were arc melted in a multicrucible furnace3 in buttons ranging from 5 to 20 g. Chemical analyses were not made, as the charge weights agreed closely with those of the buttons after melting. The alloy buttons were homogenized at 800° to 900°C. Metal log raphic and X-ray specimens were prepared and heat treated at temperatures from 600° to 1200°C. Specimens for X-ray diffraction were usually ground to a powder in an agate mortar; however, needle-shaped solid specimens were used when the alloy was sufficiently ductile to permit their preparation. Diffraction patterns were taken with a Straumanistype Debye-Scherrer camera using filtered Cu or Co radiation. The lattice parameters were obtained in A by plotting the calculated a0 values against the cos29/sin 0 + cos2?/? function and extrapolating linearly to ?= 900. Metallographic control specimens were polished on cloth wheels with diamond paste and etched with various phosphoric and nitric acid reagents. RESULTS The eighteen equiatomic alloys listed in Table I gave evidence of a cubic structure with two atoms in the unit cell, although two of these cubic structures exist only at elevated temperatures and transform to a tetragonal structure on quenching. Nine of these eighteen alloys gave diffraction patterns with super-lattice lines showing that the structure is of the CsC1-type. The lack of superlattice lines in patterns of the other nine alloys may be attributed to the small difference in atomic scattering power of the components. Metallographic study indicates the occurrence of nine narrow single-phase fields at the AB composition. Any or all of these nine may also have a CsC1-type structure. The VFe alloy was found to have a CsCl-type structure by Philip and Beck2 through the use of CrKa radiation (for which the scattering factor of V is anomalously low), whereas the Cu radiation used

Jan 1, 1960

By P. W. Bennett, E. M. Elkin

A new technique for the recovery of palladium and platinum and sludge from go12 electrolyte eliminates many of the drawbacks of the zinc-dust cementation process. In the electrolytic refining of gold by the Wohwill method, the electrolyte contains about 100 gpl gold as AuCl, and 100 gpl free hydrochloric acid. Any palladium and platinum contained in impure gold dissolve and accumulate in the electrolyte. When their concentration approaches that of gold, it becomes desirable to discard the electrolyte and to recover the gold and the palladium and platinum. As practiced at Canadian Copper Refiners until recently, discarded electrolyte was treated with sulfur dioxide to precipitate the gold. 2AuCl3 + 6H2O + 3SO2 —2Au + 6HC1 + 3H2SO4 The gold precipitate was filtered off, washed with water, and melted with scrap anodes. The wash water was added to the gold-free filtrate and treated with zinc dust to cement the palladium and platinum. PdCl2 + Zn — Pd + ZnCl2 PtCl4 + 2Zn—Pt + 2ZnCl2 The precipitate was separated and sold as Pd-Pt concentrate, and the barren filtrate and wash water were discarded. The receipts of palladium and platinum at this company were too small to warrant their refining. The process suffered from a number of disadvantages. Handling the gold-free solution was disagreeable, as it was laden with sulfur dioxide even after blowing with air overnight. Subsequent treatment with zinc dust was also troublesome. A large excess of zinc dust was required. Spray, acid mist, and vigorous evolution of hydrogen accompanied the reaction. The Pd-Pt precipitate was grossly contaminated with occluded unreacted zinc and with yellow granules of sulfur, thought to be formed by the reduction of sulfur dioxide and its derivatives. Repeated boiling of the precipitate with hydrochloric acid and washing with water decreased somewhat, but did not eliminate, the zinc. After drying overnight, the precipitate hardened to a clinkerlike mass that was difficult to grind. Grinding, handling, and packing created heavy dusting. Even weighing in an analytical balance generated dust, indicating that the powder was easily charged with static electricity. The total palladium and platinum content of the powder varied from 56 to 72 pct, with an average of 65 pct. The powder dissolved in aqua regia with difficulty, leaving a small amount of insoluble residue. However, the main drawback of the process was the strongly hygroscopic Pd-Pt product. The powder showed a marked though erratic increase in weight and a corresponding decrease in assay on exposure to air, even during the course of laboratory analysis. Such changes were the source of disagreement between Canadian Copper Refiners and the buyers of concentrate. The cause of the change in weight is believed to be hydrogen, liberated by the action of zinc dust on hydrochloric acid. Some of this hydrogen was adsorbed so strongly by the precipitate that it could not be driven off during the drying. On exposure to air it oxidized catalytically, forming water. Experiments soon showed that the entire process of treatment of discarded electrolyte then in use would have to be changed. A small quantity of apparatus was moved into the gold room, and experimentation was held to a minimum compatible with the development of a new process. PALLADIUM AND PLATINUM PRECIPITATION It was thought that zinc could be replaced with a reagent other than a hydrogen-evolving metal. Semiquantitative experiments showed that platinum and palladium could be precipitated completely from the gold-free solutions with such well-known reagents as sodium formate, formic acid, formaldehyde, and paraformaldehyde. Sodium formate was the preferred reagent, since the other three are volatile and somewhat toxic. The gold-free solutions used were those saturated with sulfur dioxide, partly freed of sulfur dioxide by aerating overnight, completely freed by boiling, as well as strongly and weakly acidic, neutral, and alkaline solutions. The reactions are complex and can be represented by these over-all equations: PdC1, + HCOONa — Pd + NaCl + HC1 + Co2 PtCl4 + 2HCOONa — Pt + 2NaCl + 2HC1 + 2CCX, Strongly acidic solutions required as much as five times the theoretical amount of sodium formate as did solutions neutralized to pH 4.5. Increasing the pH decreased the consumption of

Jan 1, 1965

By W. F. Hosford, E. W. Kelley

Deformation studies have been conducted at room temperature on single crystals of magnesium and magnesium alloys with thorium and with lithium. Single crystals oriented to suppress shear on the easily activated basal slip systems were deformed by plane-strain compression. Compression along the C axis was accommodated by {1011} banding. Compression perpendicular to the unconstrained c axis activated {1012} twinning, and, after virtually complete twinning, deformation continued by {1011) banding in the twinned material. Compression perpendicular to the constrained c axis was accommodated by the simultaneous operation of (1012) twinning against the constraint and (1011 ) banding. Although this orientation was favorable for {1010)(1210) prism and {1011}(1~10) pyramidal slip, these modes were not observed in pure magnesium or in Mg-0.5 pct Th. However, {10i0)(1~10) prism slip was observed in crystals of Mg-4 pct Li during compression perpendicular to the constrained c axis. Fracture in all materials occurred parallel to (1124) or {l~il) depending on the orientation and composition of the specimen. THE mechanical behavior of the hcp metals is strongly anisotropic. Although several slip systems have been reported the slip is cpmmonly in the directions of closest packing, the (1210),&apos; and this does not produce strains parallel to the c axis. Hence the inherent anisotropy. The deformation mode most easily activated in magnesium at room temperature is (0001)(1210)- basal slip. Also {1010}(1~10) prism slip and {1011)(1210) pyramidal slip have been reported, primarily at elevated temperatures.2"4 However, at room temperature the shear stresses to activate the prism and pyramidal modes are roughly a hundredfold greater than that required for basal slip.&apos;j4 Thus prism and pyramidal slip may be expected only under special conditions of loading. Strains normal to the basal plane can be produced by twinning, however. Many twinning modes have been reported for magnesium,&apos; with (1012) twinning the most common and relatively easy to activate. Magnesium can deform by (1012) twinning when stressed along the c axis jn tension, but not in compression. In contrast, (1011) twinning is activated by compression along the c axis and not by tension. In addition to primary twinning, secondary twinning or slip can occur within the reoriented material of primary twins.&apos; In general at least five independent shear systems must be active to bring about an arbitrary shape change such as that in the individual grains of a deforming polycrystalline material.&apos; Because basal slip can_ provide only two independent shear systems and (1012) twinning can only accommodate an extension of the c axis, other deformation modes must be active in magnesium for an arbitrary shape change to occur. The purpose of this investigation has therefore been to study the various deformation modes in magnesium at room temperature, with special emphasis on those modes that are less easily activated. The effect of the alloying elements, thorium and lithium, has also been investigated. In polycrystalline aggregates, unambiguous identification of deformation modes is extremely difficult and the direct evaluation of the resolved shear stresses to activate them is not feasible. On the other hand, uni-axial tension and compression experiments on single crystals may not activate some of the- deformation modes because basal slip and/or {1012) twinning cannot be suppressed in most orientations. However, it should be possible to activate all possible deformation modes using oriented single crystals and plane-strain compression. Identification of active deformation systems and evaluation of the resolved shear stresses required to activate them should be facilitated. Wonsiewicz and Backofen have recently completed an investigation of the plasticity of pure magnesium crystals at various temperatures utilizing plane-strain compression and selected crystal orientations. This technique has also been used in the present work. The seven orientations selected for study are indicated in Table I. Plane-strain compression along the c axis (orientations A and B) should activate some deformation mode _other than basal, prism, or pyramidal slip, or (1012) twinning. In orientations C and D, prism or pyramidal slip would be expected to take place. When the compressive load is applied perpendicular to an unconstrained c axis (orientations E and F) the three slip modes should be suppressed but not (10i2) twinning. In orientation G, basal slip should occur.

Jan 1, 1969

By F. Weinberg

Measurements have been made of solute segregation during dendrilic growth by using radioactive solute elements and ,measuring the activity of den(12-ites cut from decanted specimens. This has been done for both lead awl tin based binary alloys contaitzing the following solute additions: Ag, T1204, was dependet on ko, the equilibrium distribution coefficient in the following way Fay k &apos;c 0.1, C/C 0.6; for k0 >0.1. 0.6 <c,/c,< I. Qualitative obse?-vations were madc of dendritic segregation, by using autoradiographic techniques, for the Sn + Ag110 and Sn + Tlo4 systems. The observation were found to he in general agreement with the measurements ofCA/Co. Autoradiographic were also obtained of scctiolccl delzr11-iie stalks. These indicated that the stalks had a substructure, dclileated by solute corzetlt?atio?zs nlolg the substructure walls. A new dendrite growth direction <JI2> is reported for tila. SOLUTE segregation in dilute binary alloys has been investigated by Pfann,&apos; Smith, Tiller, and Rut-ter,&apos; and others. They considered the case of a slowly advancing plane solid interface, and derived expressions for the distribution of solute in both solid and liquid during solidification. To determine these expressions, they assumed no diffusion in the solid and either complete mixing in the liquid:&apos; or diffusion controlled solute movements in the liquid without any convective mixing.&apos; The present investigation considers solute segregation during dendritic growth, in which case the solid-liquid interface is not plane, and the growth rates are rapid. Segregation under these growth conditions has not been treated mathematically, because of the relative complexity of the system. It has been suggested by Chalmer, on the basis of preliminary results, that an alternative to the diffusion and heat flow controlled conditions during growth is &apos;diffusionless" dendritic growth in which solid is formed with the same composition as the liquid. He suggests this type of growth may depend upon a solvent-solute relationship that permits some solid solubility without excessive increase in internal energy, as is the case for solutions of tin in lead. On the other hand, Montariol,4 and others, have shown experimentally that some segregation does occur during dendritic growth in metals using etching and radioactive tracer techniques to indicate the concentrations of the solute. The present investigation was undertaken to determine, both qualitatively and quantitatively, the extent of solute segregation associated with dendritic growth in a series of binary alloys, as a function of solute concentration. PROCEDURE The solvent materials used were Vulcan Electrolytic tin (99.997 pct purity) and Tadanac lead (99.998 pct purity). The solute materials were Zn, Sn, and Sb (better than 99.998 pct purity), Ag and Co (99.5 pct purity), and T1 (Fisher "purified" metal sticks). Activation of the solute metals was carried out in the reactor at Chalk River, Canada. Master alloys were prepared by induction heating from the radioactive solute metal and the pure solvent, under argon, in graphite crucibles. Pieces of these alloys were then added to the solvent to give the required solute concentration. Dendrites were grown in essentially the same manner as that described by Weinberg and Chalmer, , in which controlled orientation single crystals were grown dendritically in horizontal graphite boats, and the liquid decanted. The crystals were grown and decanted in an atmosphere of tank argon. Before decanting, a sample of the liquid was drawn up in a glass tube and allowed to solidify rapidly. The orientations of the single crystals were such that <loo> was parallel to the growth direction, and (100) in the horizontal plane for lead, and [1101 and (110) respectively for tin. With these orientations long dendrite stalks formed along the bottom of the boat in the dendrite direction (<100> for lead and [I101 for tin) from which secondary branches grew. Only these secondary branches, which grew freely in the liquid from the dendrite stalk to the liquid surface, were used in the measurements. Accordingly, effects due to substrates and oxides on the surface of the liquid need not be considered. In order to measure the solute concentration C, of the dendrites, individual dendrite stalks were cut from the decanted specimens, remelted, and formed

Jan 1, 1962

By D. Muraoka

Introduction Soda ash, anhydrous sodium carbonate, is produced from underground trona deposits occurring in the Green River Basin of southwestern Wyoming. Stauffer Chemical Co. of Wyoming, a jointly owned subsidiary of Stauffer Chemical Co. and Rocky Mountain Energy Co., mines the trona and processes it into high-purity soda ash. Stauffer Wyoming has been in operation since 1962 with a current production capacity of about 1.8 Mt (2 million st) per year of soda ash. The refinery has five parallel operating trains using the monohydrate process which is a series of chemical processing unit operations described below. Screening and Crushing The mined trona is brought to the surface in two counterbalancing ore skips. The ore is initially screened with refinery feed (less than peagravel size material) entering the process. Oversize is conveyed to an outdoor stockpile. Ore is recycled from under the stockpile via cone feeders as the refinery needs require. The ore is conveyed to hammer type crushers and screened. The refinery feed enters the process and the oversize is recycled to the crushers. Calcining The first step in refining trona to soda ash is calcining the trona ore. This is accomplished in concurrent flow, rotary kilns that are direct fired with natural gas. [2 (Na CO NaHCO 2H O) Heat 3Na CO + 5H20 + CO 2 3~3~ 2 ~~2 32] The trona ore is fed to the kiln at the burner end and is intimately mixed with the hot exhaust gases from the burner. Direct contact with the flame is avoided. The ore moves down the kiln by gravity and rotation aided by the kiln&apos;s internal lifters. The exhaust gases are cleaned to 99 + % particulate removal and discharged to the atmosphere. Dissolving the Calcined Trona The calcined trona is conveyed from the kiln discharge to a rotary drum dissolver where it is mixed with water and weak solution (unsaturated solution of calcined trona in water). The calcined trona is soluble whereas the impurities, primarily shale, are insoluble. The discharge from the dissolver is a slurry containing saturated solution (about 30% dissolved calcined trona in water), some undissolved calcined trona, insoluble impurities, and organic impurities. The dissolver discharge enters a rake classifier that makes the first solid/liquid separation in the process by removing the heavy undissolved material from the saturated solution. This undissolved material is conveyed to a second rotary drum dissolver and mixed with water to dissolve any remaining calcined trona. The secondary dissolver discharges to a second rake classifier that separates the solids from the liquid. The solids are now primarily insoluble impurities with very little unrecovered calcined trona. These wastes are conveyed to evaporation ponds for disposal. The liquid from the secondary rake is recycled to the primary dissolvers&apos; weak solution system. Returning to the primary rake classifier, the saturated solution which contains suspended insoluble muds is pumped to a thickener. Flocculant is added to aid settling of the suspended muds. The flocculated muds are collected at the bottom of the thickener and pumped to vacuum cloth filters. The adhering saturated solution is separated from the muds that remain on the filter cloth. The recovered solution is returned to the dissolvers&apos; weak solution system. The muds are reslurried and pumped to the evaporation ponds for disposal. Saturated solution, containing some suspended insoluble muds along with the organic impurities, overflows the thickener. Filtering the Saturated Solution The saturated solution overflowing the thickener enters a holding tank where it is pumped to a bank of pressure leaf filters. The suspended muds are captured on the filter leaves while the cleaned saturated solution is pumped to polishing filters. The polishing filters are also leaf pressure type, intended to remove any suspended solids that may pass the primary filters. Organic impurities are removed with activated carbon. The purified solution is now free of insoluble muds and has most of the organic impurities removed. Crystallizing and Drying the Monohydrate The purified saturated solution is pumped to evaporative crystallizers. The solution is circulated from the crystallizer vessel through a heat exchanger. Indirect contact with steam in the heat exchanger causes the saturated solution to boil, thus crystallizing the soda ash in its monohydrate form. Slurry crystal density is maintained in the crystallizer by simulta¬neous feed of fresh saturated solution and drawoff of the slurry.

Jan 1, 1986

By M. R. Tek

General equations relating the pressure drop necessary to sustain the flow of a fluid through a porous matrix at a given rate have been developed. The results indicate that at high values of flow rate the pressure-flow behavior may not necessarily satisfy the usual Darcy equation. The mathematical analysis, carried through the micro-pore geometry and extended through the macro-reservoir scale, indicate that Darcy&apos;s law, of limited applicability to certain ranges of Reynolds numbers, can be generalized through the inclusion of some additional parameters. The "generalized Darcy equation" has also been formulated in dimen-sionless form permitting the evaluation of its predictive accuracy with regard to literature data. A comparison between predicted and experimental values indicates that the generalized Darcy equation predicts the pressure drops with good agreement over all possible ranges of Reynolds numbers. INTRODUCTION The limits and the nature of validity of Darcy&apos;s law&apos; has been a subject of every-day interest to the industry for many years. It is well known that as the Reynolds number, characteristic of the fluid flow through porous media, becomes large, Darcy&apos;s law gradually loses its predictive accuracy and ultimately becomes completely void. For the last 20 years much has been said and written on this subject. Unfortunately little has been accomplished to bring about a satisfactory agreement, at least on the nature of the threshold of validity of Darcy.&apos;s law. Fluid dynamists, geo-physicists, and engineers all had their individual views, explanations, interpretations and concepts on the subject. To some, a mechanistic analogy with pipe-flow proved a satisfactory explanation.&apos; To others,&apos; turbulence, in its random character, was incompatible with the geometric structure of consolidated porous systems. To some,4 turbulence merely represented a factor influencing the permeability measurements and again to others5,6,7 em-pirical or semi-empirical correlations proved satisfactory from an engineering viewpoint. Deviations from Darcy&apos;s law at high flow rates have been studied by systematic experiments by Fancher, Lewis, and Barnes.&apos; In an article on the flow of gases through porous metals, Green and Duwezs conclude that the onset of turbulence within the pores appears unsatisfactory to explain deviations from Darcy&apos;s law. This view is held by many others. While the subject remained controversial for many years, the development of vast natural gas reserves throughout recent years further justified considerable interest on this problem from the standpoint of gas reservoir behavior. As large amounts of field data became available from the operation of many gas fields, it became evident that the steady-state behavior of gas wells was not, in general, in agreement or compatible with Darcy&apos;s law. This suggested a careful reconsideration of all mechanisms which may account for pressure drops in addition to viscous shear. In a series of articles9,10 . Hou-peurt indicated that deviations from Darcy&apos;s law may be explained on the basis of kinetic energy variations and jetting effects without resorting to assumptions on turbulent flow conditions. Another article by Schneebeli11 indicates that special experiments by Lindquist clearly demonstrated that the onset of turbulence does not necessarily coincide with conditions of deviation from Darcy&apos;s law. This view is also held by M. King Hubbert.12 Starting with the basic pressure-flow relations suggested by Houpeurt, the derivation, development and extension of analytical expressions to -supplement and generalize Darcy&apos;s law has been the objective of this work. MATHEMATICAL ANALYSIS Derivation of Dimensionless Pressure-drop, Flow-rate Relations In considering the flow of a fluid through a porous matrix geometrically represented by a succession of capillary passages in the shape of truncated cones,810 an approximate expression may be derived relating viscous and inertial, i.e., total pressure drop to the physical properties of the fluid, geometric properties of the rock matrix and the rate of flow: ?P/?r = µ/k V [ 1 + c(m4 - 1) p V/16n" mµ w] ..........(1) Let us formally set: c (m4 - 1) / 16n" m = a d ......(2) Such a representation is equivalent to assert that the term [c(m4 — 1)/ 16n"m], variable with various porous media and probably highly variable within a given porous medium, may be macroscopically defined as equal to a lithology factor times the aver-age grain diameter d. In view of the usual grain and pore size distribu-

Jan 1, 1958

By T. H. Alden

A correlation has been found in rock-salt structure single crystals between the latent hardening, measured by the direct stress activation of oblique slip systems, and the stress-strain behavior in simple compression. Materials with high latent hardening, like LiF, strain-harden at a low rate (Stage I) even when severely constrained. KCl, in contrast, shows low latent hardening and a tendency to strain-harden at a high rate (Stage 11). This correlation suggests that oblique slip is essential for Stage II hardening of these materials. THE role of secondary slip in the strain hardening of metal single crystals is a topic of lively controversy. On the one hand are theories in which secondary dislocations participate directly in Stage II hardening in the fcc metals, for example through the production of sessile dislocations,1 forest intersections, or jog formation. Much of the experimental evidence on which these theories are based has been reviewed by Clarebrough and Hargreaves.~ On the other hand, a recent theory5 denies that secondary slip has an essential role in Stage II hardening or in the transition from Stage I to Stage 11. In the latter view, Stage I is not terminated by an increase in the activity of secondary systems but by the exhaustion of undeformed material.6 The experiments reported in this paper will not resolve this controversy since the materials being studied are cubic ionic crystals rather than metals. However, the results do show in an unusual way a direct connection between "secondary" (nonortho-gonal or oblique) slip and strain hardening in these materials. Specifically, a correlation has been found between two independently measured properties, first the latent hardening of oblique (110)(1i0) slip systems as measured by direct stress activation of these systems, and second the stress-strain behavior at small strains. From prior work, it was known that in most cubi-cally oriented rock-salt structure crystals, two orthogonal slip systems operate and exclude the other equally stressed pair, oblique to the first pair.7&apos;8 This observation apparently indicates that a significant interference exists between slip on oblique (110) planes and a relatively small interference between orthogonal (110) planes. The present experiments were begun with the intent of obtaining a quantitative measure of this difference by means of a study of latent hardening in these crystals. I) EXPERIMENTAL METHODS The experiments were basically of two types, first the determination of stress-strain curves by compression along a single (100) axis, and second the measurement of latent hardening by compressive prestrain along one cube axis followed by the determination of the yield stress in a second (latent) cube direction. All tests were done at room temperature in an Instron machine by compression between lapped, parallel steel faces. Three specimen shapes were used. For long crystals (nominal dimensions, 1/8 in. square by 1/2 in. long) the ends were lubricated with an oil-graphite mixture. Superior results with short crystals (about 1/8 in. cube) were obtained using 0.003-in. teflon film.&apos; Thin crystals (1/4 in. square by 3/32 in.) were used for latent-hardening measurements and similar results were obtained with either lubricant. Test specimens were cleaved from Harshaw single crystals which had been irradiated to a dose of lo8 roentgen using a cobalt-60 source. The irradiation raises the yield stress and tends to prevent plastic deformation during sample preparation.10 Prior to testing, the irradiation hardening was removed by an anneal at 400°C. Ideal compression specimens have flat, parallel faces. In short specimens particularly, satisfactory results demand a close approximation to this ideal. In the present work, cleaved surfaces were often used directly as compression faces. The degree of success using this method depends on two factors: 1) the smoothness of cleavage faces, and 2) the extent to which the crystal is deformed or crushed by the chisel during the cutting of a long crystal into short pieces. Unfortunately, only one material of those studied, LiF, was completely satisfactory. NaF, NaC1, KC1, and KBr behaved less well so that

Jan 1, 1964

By Louis Raymond, John E. Dorn

The object of this investigation was to analyze the recovery that arises when the stress on a specimen undertaking creep is reduced. For this purpose annealed specimens of high-purity aluminum were precrept under a stress of 1000 bsi to a strain of 0.08 following which the stress was reduced for various periods of time to 10, 250, 500, or 700 psi. When the original stress was reapplied the subsequent creep curve lay above that for the unre-covered state and below that for the original annealed state. Analyses on the kinetics of this recovery as a function of the temperature gave a stress-sensitive activation energy that decreased as the reduced stress was increased from a value of 64,000 cal per mole at 10 psi to 37,000 cal per mole at 750 psi. Recovery was also detected and measured during creep under the reduced stress. Following a short initial period, the creep rate under the reduced stress increased monotonically until it reached the secondary-creep rate for the reduced stress. The temperature dependence of this phenomenon was also shown to be correlatable in terms of the previously deduced activation energy for recovery. The activation energies for creep of most pure metals at high temperatures have been shown to agree well with those for self-diffusion.&apos;j2 Since the true secondary stage of creep is usually due to the steady-state balance between the rate of strain hardening and the rate of recovery, it is generally thought that the activation energy for recovery of the creep-induced substructure equals that for creep itself. A shoft time ago, however, Ludemann, Shepard, and Dorn~ found that the activation energy for recovery of the creep-induced substructure in high-purity aluminum under zero stress was almost twice that for self-diffusion, namely about 65,000 cal per mole; obviously recovery under reduced stresses differs in some significant way from the recovery that accompanies the secondary stage of creep. The major purpose of this investigation is to study the effect of stress on the re- covery of the creep-induced substructure in order to provide a better understanding of the recovery mechanism itself. EXPERIMENTAL TECHNIQUE High purity aluminum, containing 0.004 pct Cu, 0.002 pct Fe, and 0.001 pct Si, used in this investigation, was in the form of 0.100-in.-thick sheet which has been cold-rolled to the H-18 temper. Creep specimens were milled from the sheet with their tensile axes in the rolling direction. All specimens were then heated at 686°K for 1 hr followed by air cooling in order to produce an annealed structure which exhibited a uniform equiaxed grain size of about 4 grains per mm. Tests were run in creep machines fitted with Andrade-Chalmers type of lever arms so contoured as to maintain the stress constant to within 0.05 pct of the reported values. Constant temperatures to *O.l°K were obtained by complete immersion of each specimen in a temperature-controlled and agitated bath of molten KN02-KNOs mixture. Where changes in temperature were involved, the change was effected in less than 2 min by manually replacing one bath by another controlled at the second temperature. Displacements over the gage section were sensed by linear differential transformers, the output of which was autographically recorded. The calculated strain measurements were sensitive to 5x EXPERIMENTAL PROCEDURE The following analyses are based on extensions of the previously announced effect of the temperature on the creep strain,2 namely for a = constant, where e = the total true tensile creep strain for a given applied true tensile stress, t = the duration of the test, R = the gas constant, T = the absolute temperature, Q, = the activation energy per mole for creep which is independent of the stress, / = a function of 8, = and of the stress, and a = the stress. The validity of this correlation for high-purity aluminum is demonstrated in Fig. 1 for temperatures in the near vicinity of 600°K; the activation energy for creep, Q,, which is approximately that for self-diffusion, is insensitive to the applied stress

Jan 1, 1964

By G. Venturello, C. Antonione, G. Della Gatta

Results from work on the effect of inferstitials on recovery and recrystallization of&apos; very pure iron (99.995 pet) doped with nilrogen up to 400 ppm are reported. Nitrided specimens were obtained by heating the iron in a static atmosphere of NH3 + H2. The samples were cold-rolled 80 pet, subjected to a series of isochronal and isothermal anneals, and submitted to examination by X-rays, micrographs, and hardness tests. Small additions of nitrogen show a strong effect in reducing recovery of&apos; the mechanical properties of high-purity iron. At the temperatures of. the experiments, this effect proved greater than in the case of carbon. At 400°C pure iron recovers more than 50 pet of the total hardness increase, carburized iron vecovers 25 pet, and nitrided iron only 10 pet. On the other hand. the addition of nitrogen has, similarly to carbon, a very small effect on recrystallization; the effect is, however, slightly higher than that of carbon. When the concentration Of nitrogen is above the solubility limit in a iron at the temperattue of the experiments, an increase in the frequency of nu -cleation is also observed. THE object of the present work is to improve knowledge on the effect of interstitials both on recovery and on primary recrystallization of pure iron. In fact. at present there are no available data on the effect of nitrogen on high-purity iron. The only work which at present can be somewhat related to this problem is a work of Leak et al.1 in which grain boundary internal friction at high temperature in the presence of increasing amounts of nitrogen is studied by a torsion pendulum. The present work is a continuation of a previous one2 on the effects of carbon. EXPERIMENTAL PART Preparation of the Material. Pure iron was prepared following the same methods used in the pre- vious work2 and was nitrided using controlled quantities of ammonia in an apparatus similar to that used for carburizing. The apparatus consists essentially of a gas-tight quartz tube containing the specimens. The tube, maintained at constant pressure by means of a mercury seal, is placed in a resistance furnace. For nitriding, the NH3 + H2 atmosphere with a suitable concentration of NH3 was introduced into the tube, which had been previously carefully evacuated and degassed by heating at 600°C. The hydrogen used was carefully purified by passing it through a palladium filter. The ammonia required was prepared with the following reaction: CaO + 2NH4CI —CaCl2 + H2O + 2NH3 by heating the mixture of CaO and NH4C1 in a glass flask and then collecting the ammonia in a graduated Hempel burette. For each nitriding program a given amount of NH3 prepared by this method was introduced together with the H2 into the tube previously evacuated. Nitriding of the sample was then performed at 570°C for 48 hr. After this, to ensure a homogeneous distribution of nitrogen, the samples were further heated for 72 hr at 800°C. With this method which uses a static nitriding atmosphere, carefully controlled conditions of purity are obtained. Contrarily, with the conventional continuous-flow methods, using large quantities of circulating gas, there is a greater risk of introducing impurities into the samples. Furthermore, the required nitrogen concentration in the specimen is easily predetermined and obtained. A plot of the relation observed between the initial concentration of NH3 in the gas and the quantity of nitrogen introduced into the iron samples is given in Fig. 1. The data refer to samples of 1.2-mm thickness and to the nitriding conditions previously described. Furthermore, as a comparison, the data relative to two 0.5-mm-thick samples are reported. The presence of a reservoir in the cold zone of the reaction bulb made possible the gradual dissociation of ammonia on the sample in the hot zone. Internal-friction methods were used to determine the content of nitrogen dissolved in the iron samples. Some analyses were also carried out with chemical methods. The following table, Table I, shows a comparison of data obtained with chemical and internal-friction methods. The correspondence between chemical methods and internal-friction

Jan 1, 1964