Search Documents

By R. Schuhmann, G. W. Powell

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

Jan 1, 1970

By AIME AIME

WHEN "Chief" Buehler in mid-September set out on his official 10,000 mile swing-around-the-circle visiting Local Sections he decided not to tell his audiences how to organize and operate a state geological survey. As the director for 28 years, of the Missouri Bureau of Geology and Mines (which, he concedes, is the best in existence) lie could have done this to good advantage; but he preferred to discuss affairs of the Institute. In 25 speeches made to member and student groups, he dwelt on sonic or all of the following phases of activities. He stressed the vital importance of the Local Sections in providing the opportunity for members to associate under pleasant auspices with others in their profession. The benefits of Institute membership he classified as tangible and intangible. Of the two lie felt that the latter were more important. A man's progress in his profession depends in large measure upon a broad acquaintance and friendship with others. The Institute is incomparably the best medium for meeting fellow engineers. The primary, function of the new enlarged Membership Committee is to obtain a comprehensive list of engineers eligible to be members, always hearing in mind that the standards of membership must not be lowered in any degree. These men are to be invited by individual letter to apply for membership. Whether they become members or not is of less importance; but it is vital that no qualified man should be outside the fold for the reason that membership had never been suggested to him. Dr. Buehler emphasized repeatedly his feeling that the younger man that had never had Institute membership brought specifically to his attention had been cheated:

Jan 1, 1935

By S. J. Swainson

PHOSPHATE ROCK is a low- priced commodity. This fact has influenced the choice of mining and beneficiating methods to a greater degree, perhaps, than in most other low-grade mining operations. The fact that many yards of everburden and matrix must be stripped and mined per ton of salable "rock" produced and sold has necessitated he development of relatively simple btneficiating plants having high capacity and low capital and operating costs. The law of diminishing returns must always be reckoned with in phosphate mining and washing, and often recovery must be sacrificed in favor of costs in order to achieve the highest over-all economic efficiency. The phosphate-bearing formation, known as "matrix," is composed of (1) phosphate pebbles, which are somewhat rounded in shape and which vary in size from 1-in, to micron sizes, (2) rounded quartz grains usually smaller than 10 mesh, and (3) clay, of varying "stiffness," which is the binder of the matrix. The phosphate pebbles vary in hardness. Their color, which varies from white to black, is not a reliable indicator of the Bone Phosphate of Lime (B.P.L.) content of the

Jan 1, 1944

By J. M. Eldridge, P. Balk

Phosphosilicate glass films were formed, by reacting gaseous P2O5 with SiO2, over a large range of temperature (800° to 1200°C) and gas phase composition (nearly two orders of magnitude of effective P2Ospressure). The film compositions generally corresponded with the liquidus curve, delineating the maximum solubility of the tridymite Phase of SiO 2 in phosphosilicate liquid solution at the temperature of film formation. It is shown that the P2O5 concentration of the phosphosilicate liquid film tends to decrease by reaction with the underlying SiO 2 layer until the liquidus curve is reached. The validity of the thermodynamic argument used in this explanation is supported by the results of a determination of the composition of borosili-cute films, prepared by reacting gaseous B2O3 with SiO2 at different temperatures. The kinetics of phosphosilicate film formation were described by a model predicated on a steady-state diffusion of P2O5 through the film. UNDERSTANDING of the processes leading to formation of phosphosilicate and borosilicate glasses is of great importance for producing passivating layers on FET devices. Passivating films with optimum characteristics are preferably formed in a separate step, independent of the doping of the semiconductor.' The results of an investigation carried out to gain improved insight into the mechanism of glass formation are presented in this paper. It would be expected that application of the known Pz05-Si02 and B 2 O 3-SiO2 phase diagrams should be useful in extending understanding of the glass-forming processes. However, the question of the propriety of treating thermally grown SiO2 in these binary oxide systems by the methods of equilibrium thermodynamics must be considered when this application is attempted. Although Sah et a1.' and Allen et al. 3 investigated the kinetics of formation of phosphosilicate glass (PSG), they failed to adequately relate their diffusion models to the occurrence of experimentally observed phases in the PSG/SiO 2/Si system. Horuichi and yamaguchi4 investigated the diffusion of boron through an oxide layer and described their results in terms of a model similar to that of Sah and coworkers. More recently, Kooi 5 and Snow and Deal6 reported the compositions of PSG films formed by depositing P2 O 5 onto SiO2. These compositions apparently coincide with those at the liquidus curve delineating the maximum solubility of crystalline SiO2 in phosphosilicate liquid solutions. These authors did not discuss the thermodynamic implications of their results on the structure of thermally grown SiO2 films. The structure of thermally grown Sio2 films and that of vitreous silica are generally thought to be quite similar. Since the solubility of a substance depends on its structure, it is relevant that the solubility of vitreous silica in water7 is highly reproducible, like the solubility of thermally grown SiOz in phosphosilicate liquid. Furthermore, the vitreous silica-water system appears to be in true thermodynamic equilibrium (viz., the same solubility value can be approached from both supersaturated and under-saturated solutions). Sosman7 suggested that a type of two-dimensional lattice may form at the silica/solution interface, resulting in the observed solubility behavior that is characteristic of a crystalline solid. An alternative explanation may be that vitreous silica has a microcrystalline grain structure. Other investigators have suggested that vitreous silica has essentially the structure of B cristobalite,' or is composed of microcrystals of p tridymite or cristobalite, or a mixture of both. Presumably the grain size would be sufficiently large to minimize any appreciable contribution of the grain boundaries to the solubility of the crystalline matrix. The present investigation was carried out to clarify the significance of the boundaries in the Pa,-SiO, and B2O3-SiO2 Systems in determining PSG and BSG (borosilicate) film compositions. Furthermore, the kinetic data for PSG film formation were extended, using a wider range of formation parameters than were previously reported. One model describing the kinetics of film formation will be presented that is compatible with the thermodynamics of the Pa5-Si02 system. EXPERIMENTAL PROCEDURE Glass Film Preparation. SiO2 films (1000 to 8000A thick) were obtained by oxidation of silicon substrates in dry O2 at 1100°C. PSG and BSG films were prepared by exposing these layers to gaseous oxides obtained by reacting high-purity POC13 and BBr3, respectively, with O2. A double-columned saturator was used to ensure complete saturation of the N 2 carrier

Jan 1, 1969

By D. Mills, G. B. Craig

ZIRCONIUM has been float-zone-refined in an electron-beam furnace and the redistribution of oxygen, iron, and tungsten has been measured. The iodide zirconium used in the present experiments initially contained both oxygen and iron in the range 100 to 200 ppm by weight, and tungsten in amounts less than 10 ppm. Float-zone refining of zirconium, using induction heating, has previously been attempted by Kneip and Betterton1 and Langeron.2 Kneip and Betterton were primarily interested in the removal of iron and nickel. They achieved some purification, with respect to both these elements, in material given up to six passes at 150 mm per hr under an argon atmosphere. Langeron reported purification for a large number of elements after four completed passes in a static vacuum system operating at pressures less than 10-6 mm of Hg and with rates of zone travel between 30 and 40 mm per hr. He did not report oxygen concentrations, but stated that there was inverse segregation of this element. westlake5 has also used a minimum of four zone passes to remove iron from zirconium. No rates or conditions were given. Belk6 has shown that tungsten contamination can occur during electron-beam melting. He reported an increase from 0.01 to 0.05 pct by weight tungsten in molybdenum after four complete zone passes. The vacuum system for the electron-beam unit used in the present investigation consisted of a single-stage rotary pump backing a liquid-air-trapped oil-diffusion pump. A pressure of less than 10-8 mm of Hg was obtained with Dow-Corning 705 fluid in the diffusion pump. In order to avoid contamination of the zirconium by evaporation from the tungsten filament, a special focusing cage4 was employed. Three rates of zone travel, 114, 38, and 4 mm per hr, were investigated, with the liquid zone moving from the bottom to the top of the bar. The bars used were 3 mm diam and the total melted length was 130 mm. Oxygen was analyzed by neutron-activation analysis using a neutron flux of 10' neutrons per sq cm per sec to form the isotope N16 by the 016 (n,p) N16 reaction. The standard deviation is quoted for each oxygen determination. The iron and tungsten analyses were performed spectrographically, and the precision is estimated to be ±6 pct. Analyses for tungsten were all below the detectable limit of 10 ppm, confirming the protection given by the focusing cage. No significant redistribution of oxygen was found at the two higher rates of zone travel. The redistribution of oxygen and iron obtained after five passes at a rate of zone travel of 4 mm per hr (1.1 x 10-4 cm per sec) is recorded in Tables I and 11. Burris, Stockman, and Dillon3 have estimated the distribution of solutes during multipass zone refining. Using their curves for an effective distribution coefficient of Keff = 1.5 and 2 for oxygen and Keff = 0.3 for iron, the expected concentrations were estimated for the present material. These are shown in Tables I and 11, along with the experimentally measured values. The calculated concentrations are based on a molten-zone length of 10 pct of the total melted length, whereas in the present experiments the molten zone was approximately 5 pct of the total melted length. The effect of zone length on solute redistribution is most pronounced after a large number of zone passes. Comparison with Pfann's published data8 for solute redistribution for various Keff's and zone lengths indicates only small differences at low numbers of completed zone passes. It is evident that the expected distribution has not been realized in the case of iron. Scrapings of material deposited on the inner surfaces of the electron-focusing cage were found to be magnetic. It is, therefore, concluded that redistribution of iron is masked by evaporation. Deposition of iron inside the focusing cage was observed at all rates of zone travel. The results of the investigation may be summarized as follows: 1) Segregation of oxygen is typical for a solute with a distribution coefficient, Keff > 1. 2) No redistribution of oxygen takes place at high rates of zone travel. 3) The distribution coefficient for oxygen lies between 1.5 and 2.0. 4) Purification with respect to iron occurs mainly by evaporation. 5) The focusing cage is effective in preventing tungsten contamination.

Jan 1, 1967

By Eugene McAuliffe

THE last obtainable figures of the value of the coal mining investment are those contained in the U. S. Census Reports, from data gathered in 1919. The values shown therein and set forth below cover lands and equipment and do not include distribution and marketing facilities. The working capital necessary to meet pay rolls to purchase material and supplies and workmen's compensation insurance are likewise excluded.

Jan 1, 1926

By BLAllEY STEVENS

(Canal Zone Meeting, November, 1910.) I. INTRODUCTION. TEH object of this paper is to show how igneous intrusion is governed by definite mechanical laws. A distinction is made between dikes and fissures, and the various characteristics of intrusions are explained on an exact scientific basis. The modification of rock-stresses by intrusion is then considered and extensive isostress regions are recognized. Finally, a complete classification of intrusions is constructed covering all pos-sible cases. The generalizations of observed data, taken from Geikie's Textbook of Geology, conform to the laws of intrusion herein stated, and the deductions based upon these laws. No authority, however, is quoted for the type which I have named a corrugated intrusion. Its existence is a necessary deduction from the laws of intrusion, and even the entire absence of examples in technical publications cannot disprove its reality. In all probability, corrugated fissures would have remained undiscovered, but for the working of mineral veins. As intrusions of tabular form are seldom worked for commercial purposes, it has been left for theory to anticipate the discovery of corrugated dikes. In the light of these laws of intrusion and the corresponding laws of fissures,1 there seems to be a remarkable uniformity of stress-distribution in rocks at moderate depths. This renders .possible the mathematical calculation which should help us to approximate exact science in this and kindred branches of applied geology.

Jan 1, 1911

By F. W. Jessen, R. W. Durie

The dissolution of salt in the development of salt cavities is controlled by free-convection boundary layer flow along the salt surfaces. It is the purpose of this paper to expand upon results published previously, and primarily to illustrate the influence of surface features on the boundary layer and hence the dissolution process in the development of underground salt cavities. General dissolution equations are presented that apply to a surface inclined to the vertical. A salt surface inclined so as to be overlain by water is dissolved by boundary layer flow, while an inverted surface is dissolved by a combination of boundary layer and cellular flow. Detailed investigation of the boundary layer was limited to laminar free-convection flow. Truly laminar free-convection boundary layer flow may be encountered only in laboratory solution cavities, but the boundary layer will control dissolution even where cavity dimensions are such that free-convection turbulence will occur. The factors found to significantly affect the laminar boundary layer will affect dissolution with a turbulent boundary layer. INTRODUCTION Underground salt cavities have been used extensively for the storage of LPG products, and the use of these cavities is expanding each year. The increasing requirements for such cavities suggest the need for a better understanding of the mechanisms involved, for the purpose of the control of cavity geometry and the most effective use of the injected water. A previous paper outlined a quantitative approach to the determination of salt dissolution rates in the laminar free-convection dissolution of salt. It is the purpose of this paper to expand upon the results previously published and to point out significant factors in the application of these results to the dissolution of field cavities. The dissolution of salt in laboratory scale cavities is controlled by the free-convection boundary layer, which is the region of flow along the salt surface caused by the difference in buoyancy forces resulting from the salt dissolution. The transfer of salt into solution at the salt face causes an increase in salinity, and hence a greater density, in the water at the salt surface. The region of higher density tends to fall through the lighter, less saline water and a flowing boundary layer is the result. The water remains stationary at the contact with the salt surface,l so the resulting velocity profile is such that the influence of gravity is opposed by a viscous shear force at the salt surface. For salt to be transferred into solution from a solid surface to a boundary layer of fluid, it must diffuse through a stationary film of water at the salt contact. According to Fick's first law, the rate of mass transfer by diffusion is the product of the diffusion coefficient and the concentration gradient, with diffusion transfer in the direction of decreasing salinity. Once the salt has passed into the boundary layer, it will be transferred by a combination of molecular diffusion and convection. The convection transfer will in turn affect the rate of dissolution by influencing the gradient of the salt concentration at the surface. In the previous investigation by the authors to determine the basic mechanism of the salt dissolution process, an equation was developed that agreed very well with the experimental rates of free-convection salt dissolution from a smooth vertical surface. The technique was taken from the work of Eckert2 in the analogous heat transfer system and will be referred to as the integral method. The rate of salt transfer from a vertical section of height H was found to be:

Jan 1, 1965

By HILAND BATCHELLER

COMMENT on the economic significance of the special alloy steels seems inevitably to reduce itself to an attempt to peer into the future of the industry in which we are interested. We are all familiar with the historical development of iron and with the almost complete replacement of iron by steel. Those who have been connected with the industry since its .early days have had the satisfaction of witnessing the greatest industrial development the world has seen. As the cartoonists of the muck-raking decade liked to picture steel as an ogre armed with a spiked club, his hands reeking with blood and with golden coins dripping from his pocket, fiercely seeking for more widows and orphans to devour, I choose to call on my imagination to create a gigantic but beneficent robot of such marvellous mechanical skill. that it quickly transforms millions of tons of otherwise worthless material into the infinite variety of useful products which are so largely responsible for our present civilization.

Jan 1, 1931

By N. Pessall, A. D. McQuillan

The phase relalionsllips in an isothernzal section of the Ti-Zr-H system at 750°Chave been investigated by measurement of hydrogen equilibrium pressures. The resulting diagram indicates that the hydride phase has an anomalously wide range of homogencity extending from 66.6 to 30 at. pct H in alloys containing about 15 al. pet Zr. Measurement of the heat of solution of hydrogen as a function of alloy composition in ß-phase alloys indicates no anomalous behavior in the region of 15 at. pet Zr, and suggests that the extensive range of homogeneity in the hydride phase at about this composition must he attributed to an enhanced stability of the, hydride phase. THE manner in which titanium and zirconium react with hydrogen and the constitutional diagrams of the two metal-hydrogen systems are qualitatively very similar,'-6 as can be seen from Fig. 1 which shows the constitutional diagrams of the two systems. Quantitatively, however, there are some differences. For example, the exothermic heat of solution of hydrogen in ß titanium is smaller than that in ß zirconium, and consequently the equilibrium hydrogen pressure of a Ti-H alloy is higher than a Zr-H alloy at the same temperature and hydrogen concentration. Other differences manifest themselves in the hydride phases which exist at high hydrogen concentrations. In the Ti-H system, above room temperature, there exists only a single hydride phase in which the metal atoms are arranged on a fee lattice, but in the zirconium system there are two hydrides, and whereas the hydride of lower hydrogen concentration has a fee structure, the higher hydride has a face-centered tetragonal structure. The tetragonal hydride retains its stability up to temperatures in excess of 450oC.7 From the foregoing, it would be expected that the constitutional diagram of the ternary Ti-Zr-H system would have a very simple form, especially as titanium and zirconium are completely miscible in both the a and ß forms.8 The present work, which consists of a study of an isothermal section of the ternary system, was carried out in order to check the assumption frequently made that the electronic structure of the two metals and the mechanism by which hydrogen is held in solution in the metals is the same for each. Although this proved largely to be true, certain inter esting and unexpected features of the ternary system emerged which must be explained before we can really say that the reactions of titanium and zirconium with hydrogen are fully understood. EXPERIMENTAL PROCEDURE The most convenient method of studying hydrogen-metal systems is by the measurement of the hydrogen pressure in equilibrium with the hydrogen dissolved in the metal as a function of hydrogen concentration and temperature. In binary systems, the results are easily interpreted in terms of alloy constitution, since two-phase regions can be recognized by the fact that the hydrogen pressure at constant temperature does not vary with the hydrogen concentration in the alloy when it is in the two-phase condition. This criterion does not apply in ternary systems, but there is usually a marked change of slope in the isothermal pressure-composition curves when a phase boundary is crossed, and with care, and a knowledge of the limiting binary systems, the phase structure of the alloys can usually be deduced. In the

Jan 1, 1962

US 4,201,749-In the Bayer process production of alumina wherein aluminum hydroxide is precipitated from a sodium aluminate solution for use as seed, sodium oxalate and caustic components are removed from the precipitate by washing it first with a washing liquor which is saturated with sodium oxalate but has little content of caustic, and secondly with hot water. Substantially no sodium oxalate is recycled, yet loss of sodium aluminate is minimal. K. Yamada, T Harato, and H Kato, assigned to Sumitomo Aluminum Smelting Co Ltd May 6, 1980 8 pp. (423-127). Filed. 5-8-79 (37,145) Priority. Japan, 5-8-78 (53/54,803). Same¬British 2,020,642, dated Nov. 21, 1979.

Jan 1, 1982

By M. C. Miller, M. R. Tek, D. L. Katz

Previous studies have shown the influence of geometric variations on the drawdown and buildup behavior of a well completed in a homogeneous, porous medium containing a single-phase fluid. Other papers have shown the influence of leakage across caprocks on drawdown tests on water wells and aquifers. This paper presents two additional cases: the influence of an adjacent, expansible source of fluid, and the effect of caprock leakage on the measured in situ permeability and compressibility. Methods are developed for evaluating the influence of these two deviations from a uniform aquifer upon the in situ permeability and compre.ssibility indicated by pump tests. INTRODUCTION In relatively uniform porous media, saturated with a single-phase fluid, the effective reservoir permeability and compressibility of the water porous bed system may be obtained from drawdown and build-up data on wells, such as the initial wells drilled in aquifer storage reservoirs (Homers and Katz et al.7). In considering the drawdown and build-up behavior of wells adjacent to storage reservoirs, or of wells in an aquifer adjacent to an unknown gas field, questions arise as to the quantitative effect of the presence of compressible fluid on the measured permeability and compressibility of the formation. The solution to this problem is presented, using a constant pressure line source at some distance from two wells, one of which is pumped and the other used as the pressure observation well. The solution is presented in terms of the geometry of the wells and their relationship to the constant pressure line. It is known that caprock leakage can interfere with the determination of an aquifer in situ permeability and compressibility (Mantush4, Hantush and Jacob3 and Witherspoon and Muellerll). This paper presents an added case of the effect of leakage through a caprock of known thickness and permeability on well tests starting with a reservoir of constant pressure above the caprock and contributing to the aquifer as a result of the pump test. Again, both a pumping well and an observation well are postulated. These two cases are presented together since both a compressible fluid nearby and. caprock leakage show similar deviations from expected in situ permeability and compressibility results, and the solutions are similar in form. EFFECT OF A NEIGHBORING CONSTANT PRESSURE SOURCE ON IN SITU MEASUREMENTS Bixel, Larkin and van Poollenl presented a paper showing the influence of a discontinuity of the properties of a fluid formation system on drawdown and build-up performance. Their case is related to the present study but may be distinguished from it in that their numerical work and resulting graphs relate only to the performance of the pumping well itself. Incorrect values for in situ reservoir properties may be obtained from pump test data if the skin effect of the well and after-flow are not considered. Analysis of this problem by the authors indicated that the more reliable method of evaluating the influence of an adjacent gas field or expansible source was to use one well for pumping and another well at some distance from it for observing the aquifer pressure performance. An expansible source of hydrocarbon near a pumping well may be approximated by a line of constant pressure for a pump test. Fig. 1 illustrates a pumping well and an observation well near a line of constant pressure. Pressure behavior of an observation well during a drawdown pump test near a line of constant pressure is given by the diffusivity equation:

Jan 1, 1967

By D. M. Francois-Bongarcon

The idea of parameterizing reserves has only been formalized and scientifically applied for a few years. Its most fully developed application, that of open pits, has been operational only since 1976. Two new versions which have just been issued provide a complete solution to the problem and also a powerful tool for mine planning, and they open up a new field of investigation. This paper describes the basic principles for parameterizing reserves, their application for open-pit mines, the algorithm used, the physical interpretation of the results obtained, their economic use as well as the first prospects offered by a complete parameterization in tm of (I. ti), or of the metal price.

Jan 1, 1984

By R. L. Grayson, Y. J. Wang

This paper presents an application of goal programming (GP) to major mining machinery purchases based on historical operating data on mechanical downtime and production. The GP technique, an extension of integer programming, is used to maximize an objective function comprised of deviations away from management-preferred financial goals. Different capital budgeting techniques are used to evaluate capital proposals, and the results of these analyses are combined and weighted in a maximization of profits objective function. The GP method gives an optimal solution to a multiple goal problem. A sensitivity analysis is performed to determine the effect of a variable production improvement condition or a variable budget control con¬dition on proposal selection.

Jan 1, 1984

By K. R. Janowski, A. B. Chase, E. J. Stofel

Ruby crystals of exceptional optical quality have been grown from fluxes composed of Bi2O3 + PbF2 or Bi2O3 + BiF3. Careful control of growth conditions was required to minimize the number of growth defects consisting primarily of dislocations, occluded flux, and growth twins. Slow growth rates coupled with the addition of small amounts of La2O3 as a flux dopant were helpful in producing nearly perfect crystals. Chromium concentrations up to 1 pet were incorporated into these crystals without adversely affecting their perfection. Crystals up to 10 mm in diameter having as few as one to ten dislocations per crystal were typical of those produced under optimum growth conditions. Under these same conditions the occlusions were eliminated, and growth twins were either eliminated or reduced to one to three twin boundaries per crystal. INTEREST in the nature of lattice imperfections produced during the growth of ruby crystals (Al2O3-corundum structure) has been stimulated recently by the need for ruby lasers of good optical uniformity. As a result, ruby crystals grown by the flame-fusion (Verneuil) process have been studied in detail and have been found to contain high densities of dislocations (106 to 107 dislocations per sq cm) and mosaic structures that impair their optical quality.'-3 Attempts to improve the quality of Verneuil-grown crystals by controlling the growth variables have not been entirely successful. On the other hand, rubies grown by the flux process have received little attention even though their more nearly equilibrium growth conditions suggest that they would have fewer growth defects. Nelson and Remeika have shown recently that small flux-grown rubies can indeed be made into very low-loss lasers. However, they did not report information on crystalline imperfections in these crystals.4 Stephens and Alford have reported briefly on defects found in a flux-grown corundum crystal but they neither described these defects in detail nor related them to the method of growth.5 The present investigation was undertaken to establish the relationship between growth conditions and crystalline perfection of flux-grown rubies, sapphires, and other crystals of the corundum family. EXPERIMENTAL PROCEDURE The corundum crystals used in this study were grown by slowly cooling solutions of A12O3 in molten-salt solvents of PbO + PbF2, Bi2O3 + PbF2 or Bi2O3 + BiF3. The materials employed were Linde A (a Al2O3 powder) with reagent-grade chemicals as solvents. Some of the solutions also contained small additions (1 pet or less) of one or more of the transition metal elements chromium, manganese, iron, and titanium as crystal dopants. La2O3 was added to some of the solutions to alter the growth habit of the crystals. The La2O3 promoted the occurrence of growth facets parallel to (0112) and {1011} planes, where the indexing system is the morphological system described in detail by Kron-berg6 and used throughout this paper. The melts from which the crystals were grown were prepared by mechanically mixing the various constituents and placing them in platinum crucibles. The covered 50- and 100-ml platinum crucibles were heated in muffle furnaces at 1250°C for short periods of time (1 to 4 hr) and then program-cooled at either 2° or 4°C per hr to 1000°C. At this lower

Jan 1, 1965

By C. A. Stetefeldt

Jan 1, 1895

By Ray W. Carpenter, Ralph Hultgren, Raymond L. Orr

Heats of formation of Au-Zn alloys of compositions Au3Zn and AuZn were rneasured at several temperatures by liquid tin solution calorimetry. The data for Au3Zn show that much smaller heat and entropy effects are associated with the a2- a1 transformation than for the a1 - a transformation. This result is consistent with reported X-ray diffraction studies which indicate that the a1 and (a2 phases have ordered structures which only slightly differ, whereas the a phase is disordered. The heat of formation of AuZn(ß') does not change significantly between 322o and 8000 K, confirming that the ß' alloy remains highly ordered to high temperatures. Not many measurements have been made of the energies of order-disorder transformations in alloy phases. In many cases these energies may be conveniently determined by measuring heats of formation as a function of temperature. Many alloy phases approximately obey Kopp's law of additivity of heat capacities; as a consequence their heats of formation do not change with temperature. In ordered alloys, however, a decrease in degree of order with increasing temperature will contribute an endothermic component to the heat of formation. Thus, for ordered phases where other anomalous contributions to heat capacity in the alloy or its pure components are absent, a measured change in the heat of formation may be attributed to the energy of disordering. A portion of the Au-Zn phase diagram1" is shown in Fig. 1. The composition AuZn(ß') has a well-ordered superlattice structure of the CsCl type at room temperature. The alloy is considered to remain ordered to its melting point,1 although, from X-ray diffraction data and electromotive force measurements of zinc activities, Terpilowski3 concluded that considerable disordering began at -673°K. The Au3Zn composition is a disordered fee solid solution (a) at high temperatures. Between 693o and 543oK, the ailoy orders to form the a1 phase, having a structure based on the L12(AuCu3) type, but with a regular distribution of ordering faults in the c direction, making it tetragonal. Below 543oK, the alloy orders further by means of slight atomic shifts, without diffusion, yielding a complex long-period orthorhom-bic superlattice structure, a2.4-7 Heats of formation of Au3Zn(a2) and AuZn(ß') were measured at 363°K by Biltz et al.8 by aqueous solution calorimetry. Heat capacities of several alloys between 14.9 and 38.4 at. pet Zn in the range from 373" to 723°K were measured by Iwasaki ef al.6 Stoichiometric Au3Zn, however, was not included in the alloys measured, the nearest composition being 26.1 at, pet Zn. Also, the heating rate used for the measurements, 2 deg per min, may not have permitted equilibrium to be reached in the neighborhood of the transformations. Heat contents of solid and liquid AuZn were measured between 488o and 1138°K by Kubaschewski.9 It was decided to determine the ordering energies of Au3Zn and to check on the possible disordering of AuZn by measuring the heats of formation of these alloys at several temperatures. The present paper reports measurements by liquid tin solution calorimetry of the heats of formation of the three phases of Au3Zn at temperatures within their regions of stability, and of AuZn at a low (322°K) and a high (800°K) temperature. EXPERIMENTAL Alloy Preparation and X-Ray Examination. Ten-gram ingots of Au3Zn and AuZn were prepared by melting weighed amounts of gold (99.95 pct Au) and zinc (99.99 pet Zn) together in sealed evacuated Vycor tubes at 1060°K, followed by rapid quenching. The resulting ingots were homogenized at 900°K for 10 days. The Au3Zn alloy lost no weight in preparation, the loss from the AuZn alloy was enough to introduce an uncertainty of only 0.1 at. pet in its composition. Filings for X-ray examination were taken from various parts of each ingot and strain-annealed at 723°K for 20 min in vacuo. The Au3Zn filings were quenched rapidly to retain the disordered phase. Diffraction patterns were determined with an X-ray dif-fractometer using copper Ka radiation.

Jan 1, 1968

By G. E. Perry, J. F. Bruskotter, D. G. Russell, J. H. Goodrich

The depletion performance of gas wells has been investigated by mathematical simulation techniques. The gas well model which was studied consisted of a single well located in the center of a bounded, cylindrical, homogeneous reservoir. Dependency of gas compressibility and viscosity on pressure was considered in studies of well performance at both constant mass flow rate and constant flowing pressure conditions. To carry out the investigation, the nonlinear, second-order, partial differential equation which describes Darcy flow of a nonideal gas through porous media was solved numerically. Some of the previous investigations of gas well performance have been of limited general use, because assumptions were introduced to simplify either the gas properties or the basic differential equation. Other studies have been rigorous in these respects but have presented a very limited set of calculated results. The present study was attempted to present a rigorous theoretical model and sufficient numerical results to permit meaningful conclusions to be drawn. It was found that all terms must be retained in the partial differential equation to make accurate predictions. The neglect of higher-order terms, e.g., terms of the order of the "gradient squared", leads to serious material balance errors at large times and to conservative estimates of gas well performance. The higher the gas flow rate and/or the lower the permeability-thickness product of the formation, the more pronounced are these deviations. For example, in a well draining 640 acres in a 25-md-ft formation (8,120 MMcf gas in place) at a constant rate of 993 Mcf/D, the rigorous solution predicts a bottom-hole pressure decline from 4,000 to 1,000 psia in 8.7 years. If higher-order terms are neglected in the differential equation, this decline in pressure is predicted to occur in 5.3 years. With the results of the numerical solution of the differential equation as a basis, simple, easy-to-use approximations for predicting gas well performance for Darcy flow conditions have been developed. These simple approximations are based on the familiar equations for flow of a single, slightly compressible fluid. The approximate methods possess a high degree of accuracy and enable the prediction of long-term gas well performance to be made quickly and accurately without the use of a digital computer. Both transient and stable flow period approximations were developed. INTRODUCTION In recent years income from the sale of natural gas and associated products has represented an ever-incre as ing fraction of the industry's total revenue from operations. To meet the surge in demand for natural gas, the industry has depended heavily upon established reserves and has actively pursued development of new reserves. The search has progressively led to reservoirs which yesterday were too tight and/or deep to yield the desired return on invested capital. More than ever before, evaluation accuracy is now required to forecast the criteria upon which engineering recommendations and management decisions are based. Considerable effort has been expended by both research and operations personnel on the development and application of methods for analyzing and predicting the performance of gas wells. Fundamentally, the problem is the familiar one of extracting data during the drilling, testing and early production life of a well and applying these data within an accurate simulation model to predict long-term behavior. During the past 30 or more years a voluminous literature dealing specifically with gas field problems has been generated. A recent book' lists a comprehensive bibliography of published material through 1959. Over 1,200 references are cited. Since then 39 additional articles on natural gas technology have been published in Transactions volumes of the Society of Petroleum Engineers of AIME. Most existing theory for predicting gas well performance requires that one or more idealizations (e.g., steady-state flow, ideal gas of constant viscosity, small and constant compressibility and constant-viscosity fluid) be applied. Although existing theory may apply directly or be adapted by various artifices to describe specific gas well and reservoir behavior, no widely applicable method is available, and existing methods appear to be subject to appreciable error unless better limits of applicability are defined. The objectives of this paper are (1) to present numerical solutions to the partial differential equation describing gas flow under conditions of general interest in gas well performance prediction work, (2) to present solutions which possess a high order of accuracy for both transient and stabilized flow periods of a well producing at constant rate or constant pressure, and (3) to develop and present simplified relationships which can be used as high-order approximations to the exact numerical results for fast and accurate predictions of gas well performance at the operating level. Combined, these objectives are designed gen-

Jan 1, 1967

Jan 1, 1882

By J. B. Clark, A. S. Yue

The zone-melting technique can be adapted for the de-termination of the eutectic composition in complex metal systerrzs. The application of this method is demonstrated in a simple eutectic system, Mg-AE, in which the eutectic composition is known and in a complex ternary system, Mg-Al-Zn, in which the literature is uncertain as to the composition of the ternary eutectic. The advantages and limitations of this unique approach for the determirution of the eutectic composition are discussed. ThE eutectic composition in phase diagrams is usually determined by thermal analysis of a series of arbitarily selected alloy compositions. The temperatures of the liquidus and eutectic arrests for each alloy are plotted and the eutectic composition is estimated by extrapolation to the entectic temperature of the curve formed by the liquidus arrest temperatures for the series of alloys. Although this approach readily yields the eutectic temperatures, several dekrminations are needed to estimate the eutectic composition, even in a simple system. In a complex ternary system, a great many determinations are often necessary to estimate, even approximately, the ternary eutectic composition. Thus, the conventional method is inherently tedious and time consuming. The application of zone melting provides a relatively efficient method of determining the eutectic composition.' As an illustration, consider the zone melting of an alloy composition C, in the simple eutectic system A-B, shown in Fig. 1, under the ideal conditions of infinitely rapid diffusion in the liquid, no diffusion in the solid, and ease of nucleation of both phases. As the zone moves along the bar, the composition of the zone is enriched by solute rejected during the freezing of the a phase and the zone composition approaches the eutectic composition CE.Once the zone composition attains the lowest melting composition, i.e., the eutectic composition, it remains at this composition as the zone moves along the bar and finally freezes as such at the end of the bar. Chemical analysis of the end of such a zone-melted bar then should yield the eutectic composition CE of the system A-B. A corresponding analysis can be made for a ternary eutectic system. Under the ideal conditions noted above, only one pass of the molten zone is required to yield the eutectic composition. Experimentally, however, a solute-rich layer forms at the solid-liquid interface2 and the composition of the zone does not attain the eutectic composition as rapidly or as uniformly as in the ideal case outlined above. But with agitation of the melt, the ideal conditions can be approximated experimentally. Also, as shown below, zone melting permits a cumulative build-up of the solute at the end of the bar to the eutectic composition. The purpose of this paper is to demonstrate this unique method for determining the eutectic composition. The method is here applied to two systems: first, to the binary Mg-A1 system, in which the composition is well established, and second, to the complex Mg-Al-Zn system in which there is doubt concerning the composition of the ternary eutectic. EXPERIMENTAL PROCEDURE The zone-melting apparatus used in this study is described in detail elsewhere.3 Briefly, the binary alloy rods (5/8 in. in diam and 9 in. long) and the ternary alloy rods (5/8 in. in diam and 6 in. long), with the compositions given in Table I, were zone

Jan 1, 1962