Search Documents

By Noel Jarrett, Allen S. Russell, Bernard M. Starner, Stanley C. Jacobs

A process has been developed to refine high-grade commercial aluminum to 99.99 pct purity. This enzploys precipitating titanium, vanadium, and zirconiu~ as borides. The upgraded liquid is partially crystallized under conditions that prevent massive freezing. The resulting crystals are compacted and washed free of adhernig impure liquid aluminum by melting the top crystals while draining off the liquid. Finally, the bottorn portion of the crystal bed is melted and removed as the pure product. Leakage from the crystallizat~on vessel is prevented by a layer of alumina granules that are not wetted by the aluminum. ALUMINUM of 99.99 pct purity is desired in a variety of applications. These include electrolytic capacitors, catalyst carriers, decorative trim for automobiles, aluminum powder, and high-strength alloys. This market was for many years filled by refining commercial-purity aluminum in a three-liquid-layer electrolytic cell of the type originally proposed by Hoopes. To meet the expanding need for lower-cost high-purity aluminum, a refining process based on chemical precipitation and fractional crystallization was developed. The availability of this super-purity metal was announced in November 1961, and three patents1-3 relating to the process were issued in 1965. An alternate process, comprising continuous fractional crystallization of a portion of a liquid-aluminum feed stream under conditions of strong agitation at the surface of the growing crystals, was published by Dewey4 in 1965. Earlier proposals by Larsen 5 and by Regner6 also emphasized the advantage of rapid agitation of the melt in the region of the crystallization zone. johnson 7 more recently proposed extracting purified aluminum from an alloy by mixing crystals of pure aluminum into the molten alloy, allowing them to grow, and removing the pure aluminum as by casting. Aamot's 8 recent patent describes a freeze refining process that involves controlled partial freezing against cooled moving retort walls. The Alcoa process to produce high-purity aluminum is schematically described in Fig. 1. Relatively pure aluminum is treated with boron in a furnace (a) to precipitate high melting point impurities, the supernatant liquid is transferred to a crystallization vessel (b), the liquid is partially crystallized with mechanized "tamping" (c), the remaining impure liquid is tapped out and the crystal bed partially and directionally re-melted to wash the underlying crystals (d), and the balance of the crystals is melted and withdrawn as high-purity aluminum (e). PRINCIPLES OF FRACTIONAL CRYSTALLIZATION Commonly encountered impurity elements form two types of systems with aluminum. Titanium, vanadium, and some other transition elements in periodic groups IV, V, and VI are the high melting point impurities. Initial solidification of aluminum containing small percentages of these elements will produce solid markedly enriched in them at temperatures slightly above the freezing point of high-purity aluminum. Most other elements, such as silicon and iron, produce eutectic type reactions with aluminum. In this case the first solid formed on initial solidification at temperatures slightly below the freezing point for high-purity aluminum is markedly depleted while the remaining liquid is enriched in these elements. The equilibrium change in composition during solidification of aluminum contaminated with silicon, iron, or titanium is illustrated in Fig. 2. Here, for clarity, the temperature scale is expanded beyond the

Jan 1, 1968

By Tell Ertl

The term oil shale is defined. Foreign oil-shale developments are outlined briefly. The richest and most extensive oil-shale deposits in the United States lie in western Colorado. The Bureau of Mines is opening two adjacent areas for underground mining. in western Colorado. From the first area the oil shale is mined selectively as required by the retorting plant. The second area is being developed as a unit of a full scale oil-shale mine to ascertain the most practical methods of mining oil shale and to determine costs. In both of these mining areas, research will be carried on to develop new and improved practices for mining oil shale. Introduction Oil shale is a sedimentary rock containing a solid mineraloid of indefinite composition known as kerogen. On heating oil shale, the kerogen undergoes thermal decomposition and some of the vapor products of the decomposition can be condensed as shale oil. Oil from shale can be a major source of liquid fuels since extensive deposits of oil shale are known to exist throughout the world. The United States considers the establishment of an oil-shale industry only when doubt exists of the adequacy of the domestic supply of crude petroleum. On the other hand, oil shale is not strange to foreign oil men. The industry is an old one in other parts of the world. The French industry has been in existence since 1838, the Scottish since 1859, the Australian since 1860, the Estonian and Swedish The since World War 1, and the Manchurian industry since 1929. The most highly exploited oil-shale deposits are the Carboniferous shales of Scotland. About 20 seams ranging from 4 to I2 ft in thickness are worked. The yield during recent years in scotland has been 18 to 25 gal of shale oil per ton, from 3&"30'000 tons per year' The French and Swedish oil-shale deposits are of lower grade than the Scottish. The Estonian oil-shale deposits are 7 ft thick, and the oil yield averages 50 gal per ton. Over 1,000,-ooo tons are mined each year, half of which is used directly as fuel to substitute for coal and the other half is used for the production of oil. Immense deposits of oil shale that will yield 10 to 15 gal per ton exist in Germany. Because of the impending loss of the Rumanian and Polish oil fields and the bombing of synthetic liquid fuel plants, the Nazis in 1943 started a large-scale oil-shale program in Germany. By V-E Day, it is said, oil-shale production had reached the rate of 2500 tons per day and rapid expansion of the industry was imminent. Australian oil shale, known as torbanite, occurs in beds 2 to 6 ft thick. It is of a very high grade and on retorting yields 80 to 160 gal per ton. The oil-shale industry of Manchuria was the largest in the world during the past decade. The well-known Fushun coal deposit in Manchuria has an overburden of 500 million tons of oil shale about 450 ft thick that will yield 10 to 15 gal of shale oil per ton. The oil shale is

Jan 1, 1949

By P. Hazebroek, M. Prats, E. E. Allen

Well locations in water floods frequently do not coincide with the ideal well locations associated with the flood pattern which is used. A question of importance is. "How much oil is lost because of this?". This is the problem treated in this paper. It is assumed that each pattern in a five-spot flood has the proper number of production and injection wells but that the wells may be off-pattern. The distance by which a well is off-patern is given in terms of the pattern spacing, defined as the length of one side of a five-spot. The influence of off-pattern wells on the performance of a five-spot flood has been investigated mathematically for unit mobility ratio. All production and injection rates are assumed to remain equal. The loss in recovery resulting from pattern irregularitie: has been calculated at a cumulative water injection corresponding to a 98 per cent water-cut in a regular five-spot. Displacing a single production or injection well by about one-third the pattern distance reduces the total field recovery by 4.4/N per cent of the recovery obtainable with regular spacing, where N is the number of five-spots in the flood (assumed large). When off-pattern wells are separated by at least one regular five-spot and are displaced by less than one-third the pattern distance, the total loss in recovery in a large-scale flood is, at most, 1 per cent of the recovery obtainable from that field with a regular five-spot pattern. However, much higher losses may occur if the off-pattern wells are not separated by a regular five-spot, if the rates are imbalanced, or if the mobility of the injected water is higher than that of the crude. Although the loss in oil recovery due to pattern irregularity often may be low, the off-pattern case will always require injection of more water to reach the same final water-cut. The volume of water injected per unit volume of oil displaced has been calculated for each five-spot adjacent to a five-spot containing a single 08-pattern well. The cumulative water injection at the time of 98 per cent water-cut may be increased by as much as 37 per cent when a production well is displaced one-third the pattern distance toward an injection well. The existence of differences between actual subsurface locations and the idealized well patterns used in secondary recovery operations is well recognized. However, the methods commonly used for predicting waterflood behavior ignore such differences, even when subsurface locations are known. The reason for this is simple. The idealized well patterns are highly symmetrical and, as such, the production history of the pattern (and even of the field) can be determined from that of a few production wells—usually from only one. Consideration of the many irregular well spacings found in a specific field would entail a large amount of work, and then the results would be applicable only to that specific field. In an attempt to provide an answer to the more general problem, this paper considers a water flood having an idealized well pattern, except for one well which is not at its regular location. When several off-pattern wells are present, their contributions are assumed to be additive unless the off-pattern wells are close to each other. Only four positions of the off-center well (at x, y. of Fig. 1) are considered in this paper. These are shown

By R. E. Carter, C. Wagner, F. D. Richardson

By means of inert markers of radio-platinum, it has been shown that cobalt metal oxidizes by out-ward diffusion of cobalt atoms through the oxide. Oxidation rates have been measured at various temperatures and oxygen pressures and have been found to agree with the rates calculated from the Wagner equation and the authors' values for the diffusion coefficient of cobalt in the oxide. The distribution of radio-cobalt in growing oxide loyers has been accurately measured and found to be different from that predicted from the Wagner oxidation theory. Attempts have been made to measure the change of lattice parameter of the oxide with composition. SEVERAL investigators have shown that cobalt metal oxidizes according to the parabolic law,1-3 but no detailed study has been made of the oxidation mechanism. Cobaltous oxide has a metal-deficient lattice, but the defect concentration which is dependent on the oxygen pressure is only approximately known.' From the similarities in the structures of "CoO" and "FeO" and the rate laws governing the oxidation of cobalt and iron, Gulbransen- as suggested that cobalt, like iron," oxidizes by diffusion of cobaltous ions through the oxide layer. To substantiate this hypothesis, he has attempted to measure the oxidation rate of cobalt as a function of the oxygen pressure but has been able to report only that the rate depends on a "small" power of the pressure. On the other hand, Preece," Valensi,' and Arkharov' all have suggested from indirect evidence that oxygen diffuses through the COO to the metal-oxide interface. The present authors, in a previous paper,' have shown that the rate of diffusion of cobalt in cobaltous oxide is dependent on the oxygen pressure. In that investigation, they also found that cobaltous oxide disks made by complete oxidation of cobalt metal were porous at their centers (Fig. 5L, ref. 4). Both these facts strongly suggest that oxidation of the metal proceeds via outward diffusion of cations. In order to establish the mechanism unequivocally, oxidation experiments were made with inert radioactive platinum markers, and the dependence of oxidation rate on oxygen pressure was measured. In addition, an attempt was made to determine the form of the cation vacancy gradient, which the marker experiment showed must exist across a growing oxide layer. This was done by depositing radioactive Co on the surface of a cobalt disk and then measuring the tracer distribution across the oxide layer produced by partial oxidation of the metal. Finally, an X-ray examination was made in an attempt to find the expected change of lattice parameter of COO with oxygen pressure. Experimental The apparatus and materials used were substantially the same as those previously described by the authors.4 Marker Studies: The use of inert markers to follow oxidation reactions was first attempted by Pfeil" in 1929, and radioactive markers have recently been used by Davies et al.' In the investigation reported here, the markers were applied as an (NH4)2PtCI6 slurry in several thin streaks across both faces of four disks of cobalt metal 1 mm thick and 14 mm in diameter. The disks were then oxidized almost to completion at 1200°C in pure oxygen. After oxidation, they were mounted and sectioned and auto-radiographs were taken to locate the markers. Oxidation Rate: The oxidation rate of cobalt was measured at 1148°C at oxygen pressures ranging from 0.0055 to 1.0 atm and at 1000°, 1148", and 1350°C under 1 atm of oxygen. The rate was deter-

Jan 1, 1956

By A. L. Jr. Pozzi, A. O. Garder, D. W. Peaceman

A new numerical method is proposed for the solution of multidimensional miscible displacement problems. Besides the usual stationary grid associated with numerical procedures, the method uses moving points for which positions and concentrations are computed each time step. Based on the method of characteristics for treating combined transport and dispersion, the method applies equally well to any number of dimensions and, in contrast with other numerical procedures, properly takes into account the physical dispersion, no matter how small. Extensive tests show that accurate solutions of one-dimensional problems can be obtained over a wide range of values of the dispersion coefficient, including zero. They also show that the moving points do not need to be uniformly spaced, and that increasing the number of moving points beyond 2/grid interval does not significantly improve the accuracy of the solution. Two-dimmsional calculations using the new method were carried out to simulate laboratory model displacements in which gravity caused overriding of solvent. Good agreement between observed and computed profiles and measured and computed concentrations of produced fluid was obtained for a viscosity ratio of 1.89; at higher viscosity ratios, fair agreement was obtained Computed results were found to be significantly affected by vertical dispersion, and negligibly affected by dispersion in the direction of flow, in agreement with conclusions based on experimental observations. INTRODUCTION Although miscible displacement processes are potentially of great economic significance, limited progress has been made in developing mathematical methods for predicting the outcome of a solvent flood. Because of the complexity of the partial differential equations which describe multidimensional miscible displacement, numerical methods are a natural approach to their solution. The task of finding suitable numerical approximations to one of the partial differential equations, which involves both dispersion and convective transport, has been particularly difficult. Dispersion in the absence of transport is described by the heat flow equation, a second-order equation of parabolic type. This equation has been very successfully treated by numerical methods. Transport in the absence of dispersion is described by a first-order equation of hyperbolic type. This equation has been treated numerically with some success in one dimension, both by Lagrangian and Eulerian techniques, 1.2 but extension to two or more dimensions has not been satisfactorily accomplished. Usually one of two things happens: either the numerical solution develops oscillations or it becomes smeared by artificial dispersion resulting from the numerical process. When transport and dispersion are considered simultaneously, this numerical dispersion may dominate the low physical dispersivicy that generally characterizes miscible displacement. Stone and Brian2 have developed a new difference equation for treating combined transport and dispersion in one dimension. By suitable choice of certain parameters, they succeeded in markedly reducing both oscillation and numerical dispersion. However, a satisfactory extension to higher dimensions has not yet been found. Peaceman and Rachford3 presented a centered-in-time and centered-in-distance difference equation combined with a 1,transfer of overshoot" procedure which was demonstrated to work well in one dimension. The technique was applied to the two-dimensional miscible displacement problem, but subsequent calculations with zero dispersivity showed little change in results compared with the calculations they presented. This indicates that their method involved a numerical dispersion of the same order as the physical dispersion. The present paper describes an alternate method for solving the combined transport-dispersion

Jan 1, 1965

By T. M. Geffen, H. N. Hall

One-third to one-half of the original oil is left in most reservoirs even after thorough gas driving or waterflooding. In contrast, essentially all of the oil contacted con be recovered by flooding with a solvent (a fluid miscible with the reservoir oil). It is not economical to fill a reservoir with solvent since the value of the solvent would very likely exceed that of the oil. A practical method would be to inject a small bank of solvent to displace the oil. This solvent bank in turn could be driven through the rock by a less valuable scavenging fluid, such as natural gas which is also miscible with the solvent. Such a method of oil recovery was studied in the laboratory using long (up to 95 ft) core systems. It was found that, as the solvent moves through the reservoir, the front end of the solvent bank becomes mixed with the oil and forms a relatively short mixing zone between undiluted crude oil and pure solvent. An engineering correlation was obtained which can be used to predict the length of this zone. Displacement of the solvent by natural gas can be achieved with complete miscibility (above the critical pressure of the gas-solvent system at reservoir temperature) or with partial miscibility at lower pressures. Under conditions of complete miscibility, it was found that the length of the gas-solvent mixing zone could be predicted from the correlation developed for the solvent-oil mixing zone. With partial miscibility, two phases exist and the solvent is displaced by a modified gas drive in which vaporization and condensation occur. Equations are developed to predict the performance of this type of displacement. Two plans are suggested for application of miscible fluid displacement. one for miscible displacement of the solvent and the other applying to lower pressures where partial miscibility occurs during solvent displacement. Calculations based on the laboratory tests indicated that a bank of solvent with a volume of less than 10 per cent of the hydrocarbon volume under consideration will generally be adequate to obtain high oil recovery from the portion of the reservoir contacted in miscible fluid displacement. INTRODUCTION The conventional methods by which oil is displaced from reservoir rock—solution gas drive, conventional gas drive, and waterflooding—are all processes in which the oil is expelled from the rock by direct displacement with a fluid immiscible with the oil. Because of the immiscible type displacement, complete recovery of the reservoir oil is not achieved. In economically successful flooding operations, as much as 50 per cent of the original oil may be abandoned in the reservoir as unrecoverable. The use of a solvent (miscible with oil, such as LPG) as an injection fluid does not have the undesirable displacement features common to the use of gas or water (immiscible with oil). It can be readily appreciated, however, that it would not be economically feasible to displace the oil and entirely replace it with a solvent. A process was therefore devised for use in field operations in which only a small volume of solvent would be required to obtain essentially complete oil recovery in the portion of the reservoir contacted. In this process, the solvent acts as a buffer between the reservoir oil and injected gas, allowing the natural gas to displace indirectly the oil with much greater efficiency than if the gas displaced the oil directly. DESCRIPTION OF OIL RECOVERY PROCESS The solvent oil recovery process will be described in

Jan 1, 1958

By B. Wilhelmi, W. H. Somerton

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

By Oscar Johnson

The industry has been canvassed for information on methods of handling grinding steel at milling plants, an important subject overlooked in most descriptive articles on milling plant arrangements. The data collected are offered in this paper as a guide to designers and planners who have loading and storage problems to solve. WITH capacities of milling plants expanding from a few hundred tons per day to many thousand, replenishment of grinding media is an important phase of mill operation. Information on this aspect of milling procedure, solicited from several mining companies, is summarized here. At International Nickel's Copper Cliff concentrator the normally strenuous job of changing rods is considerably simplified by the use of the Garrow rod charging machine, Fig. 1. The overhead crane, Fig. 2, brings a bundle of rods, weighing four to a ton, to the loading platform of the machine, positioned in front of the mill. The rods roll one by one down the slight incline of the platform onto a conveyor which sends them into the mill. A modified version of the Garrow machine, Fig. 3, is used at the Federal mill of St. Joseph Lead Co. An overhead crane places the portable steel structure in position in front of the discharge end of any one of the five 9x12-ft rod mills. It then raises a bundle of about twenty 3-in. rods to the discharge platform of the charging machine. The rods are 11 ft 7 in. long and weigh about 278 lb. By tripping a lever the operator allows one rod at a time to roll onto the inclined rollers. A slight push by the operator sends the rod smoothly into the mill. The Union Miniére du Haut Katanga rod charging machine, Fig. 4, is mounted on a pair of wheels at the forward end with swiveling castors at the back for steering it to position for use or storage. Wheels are equipped with locking brakes to maintain the machine in position while it delivers rods to the mill. The device is made in two parts. The first consists of the base, which carries the supporting wheels, grooved rollers, and drive. The other part serves as the rod deck, which can carry up to six 3-in. x 12-ft rods. A short boom built integral with the base reaches through the mill trunnion and delivers the rods into the mill barrel. Three grooved guide rollers ,are provided to receive the rods from the sloping charging deck. One rod at a time is fed to the rollers and thus delivered into the mill. Two of the rollers are driven by 11/2-hp motor through a V-belt drive. The third roller is driven off the first by a single V belt in a groove below the rod groove and on the same center line. Rods are delivered to the charger by an overhead crane in bundles of three to six as needed. The power is supplied by cables plugged into power outlets, suitably located near the rod mill discharge end. The Oliver Iron Mining Co. uses a rod charger, Fig. 5, suspended from a monorail by cables. It is constructed .of structural steel and has two decks, each capable of carrying five rods of 4-in. maximum diameter. The rods ride in grooved free-turning rollers. Legs extending down from the mill end of the mechanism act as supports when the unit is ex-

Jan 1, 1955

By M. H. Gaskell, D. M. Kehn, G. T. Pyndus

A correlation of the bubble point pressure for black oil systems is developed using the standard physical-chemical equations of soiutions. The correlation is based on 158 experimentally measured bubble point pressures of 137 independent systems and is expressed in terins of the usually measured field pnrameters— flash separation gm-oil ratio, rank oil gravity, total gas gravity, and reservoir temperature. The data were obtained on systems produced in Canada, Western and Mid-Continent01 United States, and South America. The average error (algebraic) in the representation is 3.8 per cent, and the maximum error encountered is 14.7 per cent. INTRODUCTION In the absence of experimentally measured properties of reservoir fluids, it is often necessary for the field engineer to make estimates regarding the fluid properties based on the usually measured producing parameters. To aid in these estimations, various correlations have appeared in the literature in recent years. Among the pertinent properties of interest is the bubble point pressure. A correlation for this parameter has been reported by Standing'. However, this correlation was based essentially on California produced crudes and this limitation was pointed out with its presentation. The correlation presented in this paper utilized data on crude oil systems from Canada, Wes- tern and Mid-Continental United States, and South America. CORRELATION DEVELOPMENT The basic assumption used in this development is the same as employed by Standing', There is a wide variety of ways to combine these parameters; however, in this instance the combination was made on the basis of Henry's law2. Accordingly, the relationship proposed is Although Eq. 2 defines an individual system, it is of limited value since H' is a function of gas-phase composition and the system temperature. It was observed that for the systems where the bubble point was measured at several temperatures that the ratio of the bubble point pressures and the ratio of the corresponding absolute temperatures (R) were practically identical. Thus, for correlation purposes the bubble point pressure may be taken as a direct function of the absolute temperature. This relation is valid only for those systems that are not near the critical point. Accordingly, this correlation will be inadequate for systems in the region of the critical point. The solubility of the various hydrocarbons found in the gas phase increases with the molecular weight. Thus, the saturation pressure should be inversely related to the gas gravity. Applying these principles to Eq. 2 and rearranging terms gives: The variables 02 the left side of Eq. 3 were designated as the "hubble point pressure factor". The number of mols of tank oil per barrel is a function of the "mclecular weight" of the tank oil. Although the tank oil is a complex mixture, it was assumed for the purposes of this correlation that a unique molecular weight could be assigned to a given oil. This was designated as the "effective molecular weight", and was related to the oil gravity, This empirical relationship was developed simultaneously with the correlation by assuming values of M,. and working to obtain a smooth curve for both the correlation and the effective molecular weight. The relationship between the oil gravity and the effective molecular weight used in this correlation is shown in Fig. 1. The effective molecular weight is somewhat higher than the molecular weight of the C fraction. The difference between these values is largest for the low-gravity systems. It is noted that this effective molecular

By Nicholas J. Grant, Hanna Ibach, Bill C. Giessen

The system Ta-Rh was investigated over the entire comnposition range using metallogvaphic and X-ray techniques as well as thermal analysis. Terminal solubility limits, solidus temperatures, and the crystal structures of two new intermediate phases were determined. There are five intermediate phases: o, tetragonal, isostructural with o(Fe-Cr); a1, orthorhomnbic : a2, orthorhombic, isostructural with CozSi; 0 3, a high -temperature phase of unknozun structure; and a TaRh3, cubic, AuCu3 structure. a, aL, a,, and a melt peritectically; and a TaRh, has a maximum melting point. Five peritectic and one eutectic reactions occur. THE Ta-Rh phase diagram has not been fully treated in the literature; however, two inter metallic phases have been described.&apos;-&apos; This study attempts to work out the diagram covering the details of phase boundaries as well as the crystallography of the intermediate phases. EXPERIMENTAL METHODS Starting materials were tantalum powder of 99.8 pct purity, supplied by National Research Corp., Cambridge, Mass., and rhodium powder of 99.9 pct purity, supplied by Bishop Co., Malvern, Pa. Analyses as given by the manufacturers are presented in Table I. In addition to the spectroscopic analysis, an oxygen analysis has been made of the tantalum powder. After preliminary trials, the standard procedure employed in making the alloys was as follows (see Ref. 4 for further details). The powders were weighed, mixed manually, compacted under 16,000 psi, and arc-melted at not more than 500 amp under thoroughly titanium gettered welding-grade argon of 400 mm Hg pressure. The gettering process was repeated after each melt. It had been shown previously5 that under these conditions no oxygen pick-up exceeding 0.001 pct in a hafnium button could be detected. The buttons, ranging from 5 to 10 g in weight, were inverted or, if possible, crushed and remelted three times. Melting presented no difficulty, except for tantalum-rich al- loys, which showed considerable spashing during melting. The average weight loss of 1 wt pct rose to 2.5 wt pct in these cases. Subsequent melting, however, resulted in losses of less than 0.1 wt pct. The melts obtained frequently showed small in-homogeneities, mostly in the tantalum-rich part of the system below 25 at. pct Rh. This was not due to insufficient blending of the elements, but to crystallization taking place from the cooled bottom. This condition was not improved by further re-melting; therefore the following precautions were taken. Cross sections of as-cast buttons were made. For the investigation only buttons were used where these inhomogeneities did not exceed 1 at. pct (estimated from the amounts of haSeS resent). In addition, only the center portions of such buttdns were utilized to minimize errors. Smaller buttons showed smaller concentration gradients and were used as a check. Finally, all alloys were homogenized to eliminate coring (see below). After this, the error in alloy composition due to inhomogeneity was cl at. pct. In order to confirm the concentrations of the alloys and to provide an impurity check, wet-chemical analyses and oxygen-fusion analyses of four key alloys were carried out. The results, together with the concentrations as calculated with the weight balance method, are presented in Table 11. It was found that for three alloys the agreement is very good (< 0.2

Jan 1, 1964

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

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

By Edwin S. Smith

Current research in the phenomenon of liquefaction has indicated some of the in situ material properties which should be considered in the planning of waste deposits. Practical experience in tailings disposal, however, must be combined with the results of theoretical and laboratory studies to develop the most economical and safe techniques against liquefaction failure. One of the many responsibilities of mining engineers and mill superintendents is providing for the safe and economic disposal of wastes (or tailings) that remain after a mineral has been extracted from an ore. The cost of disposal of these tailings exerts a considerable influence on the minimum grade of ore extracted from the mine. It varies to a great extent with each project and is dependent on many diverse factors, such as the local topographical, meteorological and seismological conditions together with the applicable water pollution laws of the area. The process used to extract a mineral from an ore determines the physical characteristics of the tailings and, thus, is one of the more important factors to be considered in the planning of disposal facilities. Today&apos;s trends in ore-processing operations complicate the problems of waste disposal: recently developed extraction methods result in the production of finer tailings, the use of lower grade ores produces greater volumes of wastes, and the enhanced land values in the vicinity of industrial plants preclude the former practice of spreading wastes laterally and force the expansion of waste deposits upwards to greater heights. The frequent occurrence of stability problems in waste deposits indicates that many mining organizations may not be employing current soil mechanics concepts and practices. One of the most spectacular, and often most disastrous, types of tailings dike failure is that resulting from liquefaction of the fine particles combined with a breach of the main retaining dike. TAILINGS DISPOSAL Since the many variations in site conditions, ore bodies and extraction processes dictate the general disposal procedures and mill superintendents incorporate special facilities to provide the most economical disposal method consistent with accepted safety standards, it is only natural that there are almost as many different disposal techniques as there are disposal areas. Methods used in the construction of retaining dikes include: (1) transporting granular material from an adjacent borrow pit and (2) mechanically segregating the total effluent and allowing the coarse fraction to form a retaining dike by hydraulic deposition. The location of a plant and the surrounding topography exert the greatest influence on the layout of tailings disposal facilities. Most of the above-ground waste deposits can be classified as side-hill, (Fig. 1,) cross-valley, (Fig. 2,) or stockpile deposits. Side-hill disposal areas are probably the most common type of layout because of the natural advantage of gravity flow, without the disadvantage of river diversion usually needed for tailings ponds formed behind cross-valley dikes. Stockpiles, which are placed on relatively level ground, require four retaining dikes and the effluent generally has to be pumped up into the disposal area. Changing production techniques or developments adjacent to disposal areas have produced variations or partial combinations of these three categories. In a series of articles, Lenhart describes many of the problems encountered in the disposal of tailings and certain innovations used to overcome them. He discusses weaknesses in design and layout and shows that disposal techniques are often learned the

Jan 1, 1970

By George D. Graffin.

The first use of graphite is lost in the mists of time. It was used by primitive man to make drawings on the walls of caves and by the Egyptians to decorate pottery. As early as 1400 A.D. graphite crucibles were being made in the Haffnerzell district of Bavaria. Through the Middle Ages graphite was confused with other minerals, especially galena and molybdenite. Two common names that are still used for the mineral are "Plumbago," meaning leadlike, and "Black Lead." The latter name implied that graphite is either composed of lead or at least contains a large percentage of it. Gessner is credited with having recognized it as a separate mineral in 1565, but its composition was not determined until 1779, when Scheele demonstrated that graphite oxidized to carbon dioxide, thus proving its carbon constitution. In 1789, Werner named it "graphite" from the Greek word "graphein," meaning to write. The United States has long relied on foreign countries for its graphite supplies. This reliance has been due more to the lack of suitable domestic graphites than to any particular preference for buying minerals overseas. Ceylon (now Sri Lanka) graphites were imported as early as 1820, and Madagascar graphites have been imported since the properties on that island were opened up about a hundred years later. Mexico has long supplied the United States with most of its requirements for amorphous graphite. Physical Properties and Classification Graphite is found in laminated, flaky aggregates disseminated in schistose rocks. It also occurs in veins and exhibits a foliated or fibrous structure. Graphite is a black lustrous mineral that crystallizes in the hexagonal system, with rhombohedral symmetry. The crystals have tabular form and are six-sided; the faces are commonly striated. The flakes have perfect basal cleavage and are opaque. When well crystallized, the flakes have a black metallic luster while the amorphous material is black and earthy, with a microcrystalline compactness. The flakes feel greasy. Graphite is an excellent conductor of heat and electricity. It melts at approximately 3550°C at a triple point under 1250 psi. It sublimes between 3300-3500°C at 14.7 psi. At 3726°C a second triple point occurs at approximately 100,000 atmospheres. Thermal oxidation in the presence of oxygen begins at 300°C and the rate increases with temperature. The three forms of carbon (charcoal, graphite, and diamond) are distinguished by chemical and physical tests. The specific gravity of charcoal is 1.3 to 1.9, of graphite 2.266 g/ cm3 (crystal density), and of diamond 3.5. Graphite has a hardness of 1 to 2 (Mohs&apos; scale). Natural graphite can be divided into three classes: disseminated flake, crystalline vein (fibrous or columnar), and amorphous. Flake graphite is a lamellar form found in metamorphic rocks, such as marble, gneiss, and schist. Each flake is separate, having crystallized as such in the rock. Crystalline vein graphite is found in the form of well defined veins or as pocket accumulations along the intrusive contacts of pegmatites with limestones. The graphite of such deposits is of two types, foliated and columnar. The important Sri Lanka (Ceylon) deposits are vein type. Amorphous graphite is commonly found as minute microcrystalline particles more or less uniformly distributed in feebly metamorphic rocks, such as slates and shales, or in beds consisting practically entirely of graphite. The latter usually are metamorphosed coal seams and carry as much as 80-85% graphitic carbon, while the former, being altered carbonaceous sediments, commonly range from 25 to 60%

Jan 1, 1983

By D. E. Thomas

Oxidation rate constants were determined for Cu-Pd and Cu-Pt alloys as a function of alloy composition and temperature. Reaction products were identified. Relationship between oxidation rate constants and diffusion constants in the reaction zones is discussed. THE mechanism of oxidation of alloys is consider--&apos;¦ ably more complex and not so well understood as that of pure metals. The oxidation of a one-phase binary alloy of which one component is a noble metal should logically be the simplest case of alloy oxidation, for in this case one is concerned presumably only with the oxidation of the base metal and the diffusion process by which it reaches the reaction site. This had led a number of investigators to study such systems. Probably the earliest work was that of Tammann and Rienacker&apos; on the tarnishing of Cu-Au alloys at temperatures up to 300°C. These alloys oxidize according to the logarithmic law, that is, the film thickness is proportional to the logarithm of the time. The rate of film formation was found to decrease with increasing gold content. Raub and Engel² investigated the oxidation of Au-Cu, Cu-Pd, Au-Ag-Cu, Au-Pd-Cu, and Ag-Pd-Cu alloys at 750°C in air. These alloy systems form a continuous series of solid solutions at the temperature of oxidation with the exception of those containing silver. Au-Cu alloys oxidize parabolically, the rate constant decreasing slowly with increasing gold content up to about 15 atomic pct* Au, and more rapidly thereafter. An alloy with 14 pct Cu forms only a very thin layer of CuO on the surface and probably also a subscale.† Alloys containing more than 14 pct Cu have scales consisting of CuO and subscales of Cu,O embedded in nearly pure gold. The Cu-Pd alloys obey the parabolic law except for slight deviations for alloys containing 53 and 84 pct Cu. The overall oxidation rate decreases rather smoothly with increasing palladium content. A scale consisting of a thin layer of CuO and a thick layer of Cu²O and a subscale of Cu²O embedded in nearly pure palladium was observed for alloys containing 84, 53, and 30 pct Cu. An alloy with 19 pct Cu formed only a subscale of CuO embedded in nearly pure palladium, and an 8 pct Cu alloy formed only a thin film which was presumed to be PdO. The ternary alloys considered by Raub and Engel oxidized in much the same way as did the binaries. Wagner and Grunewald³ found that Ni-Au alloys, oxidized at 900°C in oxygen, formed a scale of NiO and a subscale containing NiO and gold. The overall oxidation rate increases with increasing gold content, passes through a maximum and then decreases to zero for pure gold. The fact that gold increases the oxidation rate was attributed to "pores" in the scale. Kubaschewski&apos; investigated the air oxidation of alloys of platinum or gold with up to 10 pct Cu or Ni at temperatures in the range 300" to 1550°C. Parabolic behavior was noted for all alloys except Au-Ni, with rates decreasing in the order Au-Cu, Au-Ni, Pt-Ni, Pt-Cu. The diffusion coefficients for the latter three systems are known, and they decrease in the same order. However, activation energies for the oxidation of these alloys do not agree with the activation energies for diffusion in the respective alloy systems. Kubaschewski and Goldbeck7 nvestigated the oxidation of Ni-Pt alloys in air at temperatures of 600"

Jan 1, 1952

By C. F. Weinaug, P. M. Blair

Many difference equations used to approximate reservoir flow problems treat the phase pressures implicitly but not the mobility-density coefficients. Such difference equations are neither wholly explicit nor implicit, but might be described as mixed. Mixed equations are relatively easy to apply. But the associated time truncation error is relatively large, and when used to solve problems characterized by high flow rates, these equations may be unstable for practical size time steps. This paper outlines the development of a completely implicit difference analogue for reservoir simulation, along with a Newtonian iterative method for solving the resulting nonlinear set of algebraic equations that arise at each time step. While this implicit equation requires two to three times more work than does the mixed equation, it is shown to markedly decrease the time truncation error and to yield a stable solution for much larger time steps than does the mixed equation. INTRODUCTION Calculation of multiphase, multidimensional flow in porous media is generally accomplished by numerical methods that involve approximating systems of partial differential equations by systems of partial difference equations. The early development of difference equation techniques was directed toward the solution of linear differential systems. However, the equations of multiphase fluid flow through porous media are highly nonlinear in that mobility and density often are strong functions of pressure and saturation. Thus, solution of the flow problems by difference equations involves solving sets of algebraic equations whose coefficients change from step to step of the calculations. Current practice is to evaluate most coefficients at the beginning of a time step and then to apply difference methods well suited for solving linear problems. At least two major difficulties arise from the use of this practice. First, evaluation of coefficients at the old time level and evaluation of pressures at the new time level results in larger time truncation errors near displacement fronts than if all quantities in the distance derivative were evaluated at the same time level. Second, evaluation of mobility coefficients at old time levels results in an unstable difference equation in regions of high flow rate. The first of these difficulties often results in optimistic recovery behind the displacement front, and the second requires the use of small time steps or special techniques in solving coning problems1 or problems of gas percolation.2 Evaluation of all quantities in the distance differences at the new time level results in a completely implicit difference system. Such a difference system has a lower time truncation error than equations in which mobilities are evaluated at the old level and pressures evaluated at the new level. The algebraic equations resulting from the fully implicit difference equation are, in this case, nonlinear and require some iterative method for their solution. This paper develops fully implicit difference equations for two-phase flow in porous media, describes their solution by Newtonian iteration, and gives examples of problems more easily solved by the new method than by previous methods. EVOLUTION OF DIFFERENCE EQUATIONS FOR NONLINEAR PROCESSES The completely implicit difference equation used in these calculations results from evolution beginning with explicit difference equations, proceeding through what we will call "mixed" equations, and arriving now at a completely implicit equation. Let us illustrate the evolution by considering different analogs of the nonlinear equation

Jan 1, 1970

By M. J. Pool, with Appendix by Larry Kaufman, R. W. Krenzer

Thermodynmic properties of the Cu-Mn system have been deterrrzined in the temperature range 973" to 1273°K by measuring the vapor pressure of manganese in equilibrium with alloys of compositions varying from 0.13 to 0.91 rnole fraction manganese using the Knudsen effusion technique. Activities of manganese calculated from the measured vapor pressures show large positive deviations. Calculated heats and entropies of mixing are positive with a maximum occurring at approximately 0. 7 mole fraction manganese. ThE investigation of thermodynamic properties of binary systems involving transition metals is of considerable interest at the present time. The incomplete subvalency electron levels of transition metals lead to unusual behavior in some alloys, and considerable information can be obtained by determining the thermodynamic characteristics of these systems and relating the thermodynamic behavior to observed physical properties. The Cu-Mn system, shown in Fig. 1,&apos; provides an example of an alloy system in which the determination of the equilibrium phase diagram and various physical properties has received considerable attention, but very little work has been done on measurement of the thermodynamic properties. Elemental copper and manganese have comparable electronegatives, atomic sizes, and cohesive energies, but differ in crystal structure. They form a continuous series of solid solutions only at high temperatures and the structure varies in a continuous manner from fcc at the copper-rich end to fct at the manganese-rich end.&apos; The purpose of this investigation was to measure the vapor pressure of manganese over a series of Cu-Mn alloys and from these values determine the thermodynamic properties of the binary system. Vapor pressure measurements were made in the temperature range 973" to 1273"K, remaining within the solid-solution region for each alloy studied. EXPERIMENTAL PROCEDURE The relative thermodynamic properties of solid Cu-Mn alloys were determined from vapor pressure measurements by means of the Knudsen effusion tech- nique. Several discussions of the theory can be found in the literature. In the case of copper and manganese, it is assumed that manganese is the only component subliming, so the vapor pressure of manganese over the various alloys can be determined directly by measuring the weight loss of the effusion cell. This is a valid assumption since the vapor pressure of manganese is about four orders of magnitude higher than the vapor pressure of copper in the temperature region of interest. The vapor pressure of manganese over the alloys at any given temperature is then related to the weight loss of the effusion cell per unit area per unit time. An upper limit of 1.5 x 10"s has been measured for the ratio of polyatomic to monoatomic molecules for manganese vapor in equilibrium with solid manganese at 1180°K, so it was assumed that manganese vapor is monoatomic. The effusion cells, fabricated from high-purity alumina, are 26 mm high, 18 mm OD, and have a wall thickness of 1 mm. Three different sized holes of 2.03, 2.54, and 3.18 mm diam were drilled in the lids from the inside face to minimize chipping, and the lids were then slip-cast onto the cells to provide a gas-tight seal. Alumina was selected for this particular application both because of its resistance to thermal shock and reaction with manganese vapor. In the ideal case, the orifice of the effusion cell should have a knife edge so that effusion through the orifice

Jan 1, 1970

By W. A. Backofen, D. Lee

Titanium of commercial purity (RC-70) and two all-a (hcp) alloys (4Al-1/4O2 and 5Al-2.5Sn) were tested in sheet form under conditions of combined-stress loading. Plane-strain compression and plane -strain tension together with uniaxial compression and tension provided information for constructing the tension-tension and compression-compression quadvants of the plane-stress yield locus. Coordinate values were obtained from yield-strength measurements; the slope at each point was based on a plastic-strain ratio. The materials represented a wide range of plastic anisotropy, which was reflected in correspondingly varied departures of the yield locus from that for iso-tropy. The tension-tension and compression-compression quadrants were not identical. If the uniaxial yield strength is insensitive to direction and the sign of the applied stress, the Mises or perhaps the Tresca criterion is adequate for predicting the onset of yielding under combined-stress loading. The prediction is not so straightforward in plastically anisotropic material. Depending upon details of texture, combined-stress yielding resistance may be substantially different from that indicated by the usual criteria.&apos; However. a continuum theory of yielding under such conditions is available for reasonably direct application.2,3 Whether or not it can be used successfully must, of course, depend upon how closely the conditions of the theory are satisfied in the material in question. In its present form the theory contains no provision for Bauschinger-type effects, nor does it allow for the consequences that might grow out of deformation by more than one mechanism, e.g., by slip and twinning.4 With those limitations, more importance must be attached to experimental information about the yield surface for anisotropic materials. Studies of the yield surface have generally been confined to states of plane stress (tension-tension and tension-compression) established in thin-wall tubes by externally applied forces and internal pressure. That approach would be entirely suitable for anisotropic materials if they were available in tube form or in sufficient bulk to be machined into test specimens. 6-9 It would not serve for the textured and anisotropic sheets of much current engineering interest, unless sheet could be formed into high-pressure tubing without structural change. Since there would seem to be little chance of the latter, practical alternatives are needed. One such alternative is the subject of this paper. It involves the use of a few selected but relatively simple uniaxial loading tests on sheet specimens. The underlying idea is presented here with some details relating to application and a few results bearing on how well it has worked. THE BASIS OF TESTING Conventional Measurements. Yielding for a general case of anisotropy is described by the plane-stress locus in Fig. 1. Reference directions in the sheet specimen are related to the stress coordinates as shown with the insert: x is the rolling direction, y the transverse direction, and z the through-thickness direction. The four reference-direction yield strengths (t, tension; c, compression), Xt, Xc, Yt, and Yc, are all different. Nine paths that can be followed in uniaxial loading experiments are identified. Paths 1, 2, 3, and 4 represent conventional tension and compression testing to establish Xt, Yt,.Xc, and Yc. Useful information also comes from the ratio of the transverse plastic-strain components associated with yielding and flow in conventional tension or compression. A condition to be satisfied in locus construction is that the plastic-strain vector be everywhere normal to the locus.10 The projection of the vector which would

Jan 1, 1967

By M. Mills, M. W. Clegg

This Paper Presents an approximate analytical solution to the problem of the pressure distributions arising from the production of a compressible liquid in a partially penetrating well. The limits of validity of the approximation are discussed and it is shown that the basic approximation is the same as that used in previous numerical of the problem. The variation in pressure at the production well with time for different penetration ratios is investigated and is shown to agree with previous work. The solutiOns are used to assess the effects of partial penetration at large distances from the well and at the wellbore itself. In the latter case the effects of anisotropy are also studied for constant, but unequal, horizontal and vertical permeabilities. It is &apos;uggested that measurements of the Pressure drawdown below the producing interval in a well should give some indication of a thickness-vertical permeability product for the formation; this would chmacterize the effectiveness of the formation below the producing interval. INTRODUCTION In many oil and gas reservoirs the producing wells are completed as partially penetrating wells; that is, only a portion of the pay zone is perforated. This may be done for a variety of reasons, but the most common one is to prevent or delay the intrusion of unwanted fluids into the wellbore. This paper studies the effects of partial penetration on the pressure distribution resulting from the production of a compressible liquid at a constant rate. To simplify the mathematical analysis the case of a single well located in an infinite reservoir of constant thickness is considered, and it is assumed that this well is perforated over an interval adjacent to the upper (or lower) impermeable boundary of the formation. A well producing from an arbitrary interval within the bulk of the formation could be treated in a similar way, but this action would merely add to the complexities of the solution and contribute little to the appreciation of the over-all effects of partial penetration. In a recently published paper by Odeh,10 the more general problem of producing from an arbitrary interval within the pay zone was studied for the quasi-steady state flow of a compressible liquid, This is the case in which there is a a potential at the outer boundary and a constant flux across the same boundary. This contrasts with the problem studied here, which is concerned with the nonsteady-state behavior in an infinite reservoir. The exact solution of the partial penetration problem presents great analytical problems because the boundary conditions that solutions of the partial differential equations must satisfy are mixed; i.e., on one of the boundaries the pressure is specified on one portion and the flux on the other. This difficulty occurs at the wellbore, for the flux over the nonproductive section of the well is zero, and the potential over the perforated interval must be constant. In the case of constant rate production from the well, this uniform potential is time dependent and unknown, and the additional condition that where hl is the thickness of the perforated interval, must also b&apos;e satisfied. This problem may be overcome in the case of constant rate production by making the assumption that the flux into the well is uniform over the entire perforated interval, so that on the wellbore the flux is specified over the total formation thickness. This approximation naturally leads to an error in the solution since the potential (pressure) will not

Jan 1, 1970

By A. Steiner, E. Miller, K. L. Komarek

Equations have been derived to calculate the chemical potentials of the components of liquid binary alloys from liquidus and enthalpy data. The equations are applicable to systems with intermetal-lic compounds of limited solid solubility. The liquidus curve of the Mg-Sn system was accurately redetermined, and the melting point of the compound Mg2Sn was found to be 770.5° * 0.3°C. The thermo-dynamic properties of the liquid alloys were calcu -lated from the liquidus data. The activity of tin shows both positive and negative deviations from Raoult&apos;s law and the relative integral entropy is positive. THE phase boundaries in an equilibrium phase diagram are functions of the thermodynamic properties of the coexisting phases. Specifically, the shape of the liquidus curve near an intermediate compound is determined by the properties of the compound and the liquid phase. Hauffe and wagnerl have derived an equation for the chemical potential differences (ui) of the liquid alloys near the congruent melting point when the heat of fusion of the compound and the liquidus curve are known. In their derivation the temperature dependence of the ui value was neglected, and liquid of stoichio-metric composition was chosen as the reference state. To obtain chemical potentials over the ntire composition range of the phase diagram with the pure components as the reference state, equations have been derived which incorporate the temperature dependence of the chemical potentials. Such a thermodynamic analysis can be carried out when the liquidus curve and the enthalpy values are known. The equations have been applied to the redetermined liquidus curve of the Mg-Sn system to obtain the thermodynamic properties of liquid Mg-Sn alloys. The Mg-Sn phase diagram has been studied by various investigators and the results have been compiled and critically assessed by Hansen.2 Preliminary thermal analyses showed that the basic features of the diagram were correct. The system has one congruent melting compound, Mg2Sn, two eutectics, and some solid solubility of tin in magnesium. However, the published liquidus temperatures show considerable scatter. In order to perform an accurate thermodynamic analysis, the liquidus curve between the magnesium-rich and tin-rich eutectic compositions was therefore redetermined by careful thermal analysis using highest-purity starting materials. EXPERIMENTAL PROCEDURE The magnesium metal (Dominion Magnesium Ltd, Toronto, Canada) had a purity of 99.99+ pct with the following impurities (in ppm): 20 Al, 30 Zn, 10 Si, <1 Ni, <1 Cu, <10 Fe; tin (Cominco Products, Spokane, Wash.) contained 99.999+ pct Sn. The graphite crucibles were machined from graphite rods (United Carbon Products Corp., Bay City, Mich.) which had an average total ash content of less than 30 ppm. All the graphite parts were baked out in vacuum at 950" to 1000°C for a minimum of 24 hr to remove volatile organic materials. The graphite crucibles (3 in. long, 1-5/8 in. OD, 1-3/8 in. ID) were tightly closed with a threaded cover (1/2 in. thick) which had a central thermocouple well (2-1/4 in. long, 5/16 in. OD, 3/16 in. ID). Cover and thermocouple well were machined from one piece of graphite. A thin sheath of quartz was inserted in the well to protect the thermocouple

Jan 1, 1964

By M. Kaufman, A. E. Palty

The phase structure and aging characteristics of two nickel-base alloys, Inconel 718 and 702, were investigated. Wrought and cast Inconel 718 showed ArisCb as the major hardening phase, as well as ture was positively identified. Inconel 702 had Ni3Al as the major hardening phase, as well as Cr7C3-of the phase amounts with temperature was determined. 1 HE value of phase structure and aging rate information in the utilization of alloys has been well established. As part of a continuing program, the results obtained on two nickel-base alloys, hconel 718 and Inconel 702, are presented here. Inconel 718 is a relatively new alloy considered for application in both sheet and cast form. Its novelty lies in the use of high columbium tantalum levels in a Ni-Cr-Fe base. No previous work on this type of alloy has been reported in the literature. Inconel 702 is being used now as an oxidation-resistant sheet alloy. Occasional fabrication problems spurred work on its basic structural behavior. The composition of 702 is essentially nichrome plus aluminum and titanium for strengthening and added oxidation resistance. This puts it in the basic class of other nickel-base alloys which have already been studied. The two alloys will be treated separately. INCONEL 718 Phase Study Methods. The phases present after the various conditions investigated were determined by X-ray diffraction analysis of two different elec-trolytically extracted residues and microscopic examination. The extraction methods, using 10 pct HC1 in alcohol and 10 pct H3PO, in water and the X-ray technique have been described previously.&apos; The effect of temperature on the phase behavior was studied on sheet material using the following heat treatments: 2250°F, 2 hr, ice brine quench 1300 100 hr, water quench 2250°F, 2 hr, ice brine quench 1400°, 100 hr, water quench 2250°F, 2 hr, ice brine quench 1500°F, 100 hr, water quench 2250o 2 hr, ice brine quench 1700°F, 48 hr, water quench 2250°F, 2 hr, ice brine quench 1850°, 24 hr, water quench 2250°, 2 hr, ice brine quench 2000°F, 6 hr, water quench The high solution temperature was an attempt to dissolve as many phases as possible without causing melting. Aging exposures were made longer than normal heat treating times in order to approach equilibrium conditions more closely. Cast material was examined only in the as-cast condition and with the normal heat treatment (1700°F, 1 hr, air cool 1325oF, 16 hr,air cool). Materials. The full phase study was performed on sheet specimens approximately 1/2 in. by 2 in. by 1/16 in. The cast pieces were cut from the gating system of a vacuum-cast test piece. An all weld-metal specimen was made by running many weld beads on the sheet using filler made from the sheared edges of the sheet. All external contamination was removed by belt-sanding followed by an electro-polishing treatment. The typical composition of wrought Inco 718 is as follows: The alloy is air-melted. Cast hco 718 may have higher molybdenum. Due to the expectation of finding a Ni-Cb phase (Ni, Cb) for which no X-ray diffraction pattern was available in the ASTM Index, a special 5 Ib heat to the stoichiometric composition of this intermetallic

Jan 1, 1962