Search Documents

By K. G. Davis

Dilute alloys of silver and of thallium in tin have been solidijzed unidirectionally under controlled conditions, to study the segregation associated with a cellular interface under conditions where both thermal and solute convection are present. Autoradiography and radioactive tracer counting techniques were combined with electron-probe microanalysis to study both macro- and microsegregation. It was found that, for concentrations giving only small amounts of constitutional supercooling, cell formation had little effect on the macroscopic distribution of solute along the specimen. At higher concentrations the effective distribution coefficient was higher than that expected for a smooth interface. Node spacing was independent of initial solute content at lower concentrations, becoming greater as keff increased. Silver content at the segregation nodes of silver in tin alloys was independent of initial concentration and considerably in excess of the eutectic composition. SINCE the investigation of cell formation at advancing solid-liquid interfaces by Rutter and Chalmers,' a large volume of work has been dedicated to the determination of solidification conditions under which a planar interface will break down into cellular form. Early experiments were explained satisfactorily by the concept of constitutional supercooling,2 but, due to poor measurement of temperature gradients in the liquid, lack of accurate data on liquid diffusion and equilibrium distribution coefficients, and uncertainty about the effects of thermal and solute convection, these experiments cannot be used as proof for the theory. More recent work, however, has shown that under conditions where convection is eliminated or can be ignored good correlation is observed.3,4 Investigations into segregation at cell caps5 and at cell nodes6-'' have been made, but no measurements appear to have been done on the overall, macroscopic segregation down a unidirectionally solidified rod of material which has solidified with a cellular substructure. This has practical importance in casting, where regions of material with cellular substructure are often encountered, and also in zone refining where the thermal conditions necessary for a planar interface are unattainable. Further, as will be shown, the macroscopic segregation can give information on the following question. Granted that a cellular solid-liquid interface develops from a planar one when the conditions for constitutional supercooling are exceeded, how much supercooling is present after the cells have formed? EXPERIMENTAL PROCEDURE AND RESULTS Specimen Preparation. Specimens 25 cm long with a square cross section 0.6 by 0.6 cm were grown in graphite boats by solidification from one end. Alloy compositions are given in Table I. Two specimens of each composition were grown. The tin was 5-9 grade and the silver and thallium both 4-9 grade. Ag110 and Tl204 were used as tracers. Each composition had the same quantity of tracer so that auto radiographs of specimens containing different concentrations of the same element could be easily compared. Thermocouples inserted through the lid of the boat into a dummy specimen showed that, over the first 10 cm of growth, thermal conditions were quite steady, with a rate of interface advance of 5.8 cm per hr and a temperature gradient in the melt ahead of the interface of 3.0°C per cm. The specimens were seeded from tin crystals of a common orientation to eliminate orientation effects. Dilution of the specimen by seed material was minimized by the provision of a narrow neck between specimen and seed crystal. Macrosegregation. After growth, the specimens were sectioned with a spark cutter. The rods of silver alloy were cut into 1-cm lengths and analyzed for Ag110 using a y -ray counter with fixed geometry. The specimens containing thallium were cut into 2-cm lengths and analyzed for T1 204 by taking 13 counts from each end of the cut lengths through an aperture in lead sheet approximately 0.4 cm square. The results are summarized in Figs. 1 and 2. To find the effective distribution coefficient for the silver in tin alloys under smooth interface conditions, the region of substructure at the bottom surface of one of the 10 ppm specimens, see Fig. 3, was removed by spark machining before counting. Autoradiography. For both alloy systems the samples were polished on sections taken alternately parallel and perpendicular to the growth direction, and autoradiographed by placing the polished surfaces in contact with Kodak "Process Ortho" film. Figs. 3 and 4 show the structures revealed. The alloy containing 10 ppm Ag showed substructure only after a few centimeters of growth, and then substructure was limited to a narrow layer at the base. The "speckled" substructure reported previously in this system4 is here clearly seen to be an intermediate stage between planar and cellular interface conditions. The other samples show a remarkable similarity considering

Jan 1, 1969

By J. E. Zeilinger

This Bureau of Mines report describes the preparation characteristics of the significant coalbeds of Clay, Gilmer, and Lewis Counties, W. Va. In 1966, almost 1.3 million tons of coal was mined in these counties, of which more than one-half was from the Pittsburgh coalbed. Washability analyses were made of 16 coalbed samples. Of these, only 11 samples could be upgraded to metallurgical quality. For steam or power generation five of the coals were of low sulfur content as received. The remaining samples were generally of medium sulfur content even after upgrading.

Jan 1, 1969

By P. J. Carroll, H. Colijn

Feeders must be considered an integral part of the over-all bin and feeder system. Improper design of either one of these parts will affect the performance of the whole system. Feeder load, material properties, feed rate, physical layout, and environment as well as first cost determine the selection of feeder type. This paper provides a new method of calculating feeder loads with some experimental data as verification. General design criteria for belt and apron feeder, rotary table and plow feeders are discussed and a special section is devoted to vibratory feeders. Feeders are used to provide a means of control for the withdrawal of bulk materials from storage units such as bins, bunkers, silos and hoppers. This control function can only be performed properly as long as the bulk materials flow by gravity to the feeder in a uniform and uninterrupted fashion. During the past few years various papers have been written and presented on feeders and their application.1-9 This paper will be confined to some specific bin-feeder design criteria such as feeder load and power requirements. Special attention is given to the basic design criteria for vibratory feeders. HOPPER AND BIN DESIGN Feeders must be considered an integral part of the over-all bin and feeder system. Improper design of either one of these parts will affect the performance of the whole system. A lot of effort has been devoted in the past decade to the development of functional design criteria for hoppers and bins.1,10 The integrated concept of bin and feeder design requires quantitative analysis of the bulk material characteristics before the design and the selection of components can be attempted. The design of a feeder system must begin with the proper dimensioning of the hopper outlet to prevent arching or ratholing. The hopper opening should be large enough to allow passage of the bulk solid at the required maximum discharge rate." A feeder can only throttle the flow. FEEDER LOAD DETERMINATION Over the past ten years, there have been various attempts to determine the exact load or pressure on feeders mounted directly under hopper openings. Many designers, however, still consider this pressure equal to the hydrostatic head of material above the opening, and because such pressure would be very

Jan 1, 1969

By N. Standish

Heat transfer coefficients have been measured in beds of various packings irrigated with mercury and molten fusible alloy countercurrent to hot gases. The measured coefficients for both systems were found to increase with gas velocities and liquid rates. Correlations were determined which show this dependence and also indicate that heat transfer in these systems is influenced by the liquid flow characteristics and the thermal conductivity of the gas and the solid packings. A heat transfer model has beer2 proposed which explains the various features of the experimental results. On the basis of this study, which gives an insight into the heat exchange in the melting zone of the blast furnace, it was concluded that by comparison with the furnace stack heat transfer coefficients are about 1.5 times higher in the melting zone. EACH year large tonnages of metal are produced in operations which, in part, involve liquid metal irrigation of "packings" countercurrent to hot gases. The melting zone in blast furnaces and in cupolas is a good example of packings irrigated with a liquid melt countercurrent to gases. In all instances of this kind large amounts of heat are exchanged and it is desirable to have some knowledge of heat transfer phenomena involved in these systems. So far the most common method of analyzing furnace efficiencies, fuel requirements, and the general thermal state of the furnace has been through the use of heat balances. As heat balances are essentially statements of the first law of thermodynamics they give no real indication of the factors which govern heat transfer between phases in the various zones of blast furnaces. Hence, rational improvement in production efficiency and the development of theoretical models is only possible if the heat transfer characteristics are known at every stage of the process and related to the important variables involved. This has been generally recognized for some time but it was only recently that Kitaev et al.' have produced a comprehensive treatment of heat transfer in solid-gas countercurrent systems such as the blast furnace stack and the packed bed regenerator. Using their treatment it is now possible to predict the effect of particle size, thermal conductivity, bed porosity, and the flow rates of both the gas and the solid material on the heat transfer in the blast furnace stack. However, the stack of a blast furnace is only one part of an integral unit for which the heat transfer analysis cannot be complete without also considering the heat exchange in the melting zone. The complexity of heat transfer processes in this region of the furnace has so far escaped quantitative description. Yet, the melting zone accounts for a greater amount of heat exchange than all the other zones of the furnace put together. Moreover, if the reduction of oxides in the melting zone proceeds in part in the liquid state the importance of heat transfer on furnace productivity and on the metal and slag temperatures is obvious. THEORY Heat transfer for two-phase flow in packed beds is a complex problem involving a number of heat exchange paths for which interphase areas are not known with any degree of certainty. Analytical solution is, therefore, difficult. This difficulty is emphasized by noting that Rabinovich~ and Luck have only recently solved the steady-state heat transfer for simplified two-phase heat exchangers of known area. However, useful progress can be made for the system considered by making a not unreasonable assumption that the usual heat transfer considerations apply and restricting treatment to the steady state. For these conditions the rate of heat transfer dq in a height dz of a packed bed of unit area is: dq = UaATdz [I.] Integration of Eq. [I] then gives the total heat transferred: assuming both U, the overall heat transfer coefficient, and a, the interphase area, to be independent of bed height. Since a, in these systems, is unknown it is convenient to combine this term with U. The group U, then represents the overall heat transfer coefficient on a volumetric basis. If AT is linear with q, then for a bed of unit volume Eq. [Z] can be integrated to give: is the log mean of terminal temperature differences. From Eq. [3] U, can be readily calculated as q and {AT)im are experimentally obtained quantities, but a difficulty arises in interpreting its meaning. Two approaches are possible depending on whether the effect of packing in the transfer of heat is neglected or not. If the packing is thermally decoupled then the resistance concept gives the relationship: which states that the overall resistance is the sum of the gas phase and the liquid phase resistances (assuming areas are equal throughout). Because the resistance to heat transfer in liquid metals is negligible by comparison with that of the gas,4 Eq. [4] can be simplified, i.e.:

Jan 1, 1969

Two structurally and mineralogically distinct sulphiderich veins cross-cutting the disseminated sulphides of the Nairne Pyrite Member appear to be products of local remobilization of sulphides under fracture-induced free energy gradients during regional metamorphism and tectonic deformation.The contrasting mineralogy of the two vein types may have arisen as a result of differences in the structural orientation of the veins and!or the relative ages of the veins.

Jan 1, 1969

By M. I. Jacobson, M. L. Hammond

The solid solubilities of high-purity silver, aluminum, chromium, copper, and iron in zone-refined beryllium have been determined. Diffusion couple sandwiches of beryllium with aluminum, copper, and silver foil and diffusion couples of beryllium with two-phase alloys in the Be-Cr. Be-Fe, and Be-Ag systems were analyzed by electron beam micro-probe techniques. This paper describes the material used, the experimental methods for preparing the diffusion couples and the alloys, the analytical techniques, and the results obtained. The results are compared with those of previous investigators who used materials of lower purity. ~NVESTIGATORS of beryllium phase diagrams have often been hampered in their interpretations of certain phase boundaries because of the relatively low purity of the beryllium used. This has sometimes resulted in the appearance of three phases in the regions where only two should exist. Advances in purification methods, specifically zone refining, have made available beryllium that is considerably purer than the material previously used for phase diagram studies. The primary object of this investigation was to reex-amine portions of the beryllium-rich end of the binary systems of beryllium with silver, aluminum, chromium, copper, and iron to determine more accurately the limits of terminal solid solubility. The five elements—silver, aluminum, chromium, copper, and iron—include some which are commonly found as the major impurities in beryllium, some of which are being used as solid-solution alloying elements, and some of which are used in ductile matrix alloys. They are thus of considerable importance to physical metallurgists who must work with beryllium but who must rely on the present inaccurate or incomplete diagrams. The determination of more reliable diagrams would aid researchers concerned with the problems of alloying, heat treatment, and ductility. MATERIAL The beryllium was acquired from the Franklin Institute in the form of a single-crystal, 6 in. long by 1 in. in diam, zone-refined in three passes from Pechiney SR-grade beryllium bar. The principal impurities, as determined by mass spectrography at Battelle Memorial Institute, are shown in Table I. Vacuum-melted copper, silver, and aluminum, and electron-beam zone-refined iron were obtained from the Materials Research Corp. in the form of 1/2-in.- diam bars. Iodide chromium was obtained as flakes from the Chromalloy Corp. Analyses supplied with these metals are shown in Table 11. EXPERIMENTAL PROCEDURE Preparation of Diffusion Couples. Two types of dif-fusion couples were used in this study. For the copper, silver, and aluminum systems, the couples were sand-

Jan 1, 1969

The research project on size reduction, sponsored by the mining industry through financial support from individual companies, is now embarking on a second three-year program. The first 3-year program had as its objectives assimilation of the literature on comminution and development of simulation programs encompassing the basic mechanisms of grinding. This involved investigating conditions under which a particle in a mill is captured, the stress to which it is subjected, its break probability, and the resulting particle size distributions. The results were used to obtain the most representative simulation models based on current knowledge of mechanisms. As part of this work a set of computer programs has been developed for the following cases: . Batch mill, by numerical integration for completely general selection and breakage functions. . Continuous open-circuit mill by same method as above. . Continuous open-circuit mill by general numerical integration method, including transient behavior

Jan 1, 1969

By O. G. Paasche, Yuan-Shou Shen

THERE is a disagreement among the various authors about the exact manner of transformation of ZrCr2. Rostokerl and others2 stated that ZrCr2 had a C-14 (MgZn2) type of structure below 1000°C and a C-15 (MgCu2) type of structure at temperatures above 1000°C. Alisova3 and others4 reached the opposite conclusion and stated that the transformation temperature is close to the melting point of ZrCr2. A literature survey shows that various investigators3'= who homogenized the specimens at a temperature higher than 1000°C have concluded that ZrCr2 had the C-15 structure at room temperature. Meanwhile, Jordan et al.4 reached similar conclusions without annealing the specimen. Other investigators1,2,6,7 who X-rayed the specimens in the as-cast condition without annealing reached different conclusions. The investigation reported herein was conducted with the aim of exploring the exact manner of transformation of ZrCr2 by various heat treatment tests. The alloys for this examination were prepared from iodide-reduced zirconium crystal bars, 99.9 pct purity, and electrolytic chromium, 99.9 pct purity. They were melted in a nonconsumable electrode arc furnace with water-cooled copper crucible in a helium atmosphere. The melting loss of each alloy was less than 1.5 pct by weight. Chemical analysis of a randomly selected specimen indicated that there was a very close agreement between calculated and analyzed compositions. Before being heat-treated each specimen was encapsulated in a vycor or quartz tube inside which an argon atmosphere was maintained at a pressure of lower than 1 atm. In determining the crystal structure of each specimen with a Debye-Scherrer camera, the standard procedure8 for X-ray quality analysis (Hanawalt method) was followed. The different series of heat treatment tests in this investigation are tabulated in Tables I and 11. The tests in Series I, specimens from 1-1 to 1-9, which were similar to Rostoker's experiment1 indicated that the transformation temperature seemed to fall between 870° and 900°C and that the crystal structure of ZrCr2 at lower temperature seemed to be of the C-14 type. However, once the compound is transformed to C-15 type, it is impossible to reverse the transformation back to the C-14 type by first heating the specimen above 900°C and then annealing it slowly below 900°C as shown in Experiments II-1 to II-3. Thus, it appears that the specimen of ZrCr2 will transform from C-14 to C-15 structure when heated above 900°C but will not transform from C-15 to C-14 when annealed slowly passing 900° C even after the extremely slow cooling process such as indicated in the experiment of Specimen II-3. As a valid transformation temperature is a temperature at which the transformation is reversible, therefore the temperature 900°C (or other temperature close to 900°C) is not the transformation temperature for ZrCr2 and the C-14 structure is not the stable structure of ZrCr2 at lower temperatures. The C-14 structure is retained at room temperature because the transformation to C-15 structure is very sluggish and the fast cooling after melting does not allow enough time for the transformation to take place. Additional energy is required to alter the metastable condition of the C-14 structure. The sluggishness of this transformation was again demonstrated through another series of experiments. Four specimens with C-14 structure were taken. Then they were annealed at 900°C but each specimen was soaked for a different period of time, Table 11. X-ray diffraction patterns of this group indicated that the C-14 structure gradually disappeared as the soaking period was lengthened. The figures listed under the column "C-14 Structure, pct" were estimated from the intensity of the d = 2.330 line of the diffraction pattern corresponding to the structure. Notice that the intensity of this line became weaker for longer soaking periods. To determine the transformation temperature of ZrCr2, specimens with C-14 structure (as-cast condition) were annealed at 1300°, 1400°, 1500°, 1550°, and 1600°C, respectively. A final specimen was first heat-treated to 1500°C in order to transform it to C-15 structure, then heat-treated at 1600°C again. From the X-ray analyses of this series of tests, Specimen Nos. III-1 to III-6, it is evident that a transition from C-15 structure at lower temperatures to the C-14 structure occurs at some temperature between 1550° and 1600°C.

Jan 1, 1969

By Charles Prasky

The Bureau of Mines investigated the conversion of nonmagnetic iron minerals to magnetic form through reduction Toasting operations. The primary objective was to obtain a concentrate that could be used as a direct shipping material for blast furnace consumption. Crude iron materials from the Steep Rock Lakes, Ontario, area were prepared into a blend of 8 parts siliceous ore, 3 parts pyritic ore, and 1 part tailings. Experimental work included laboratory reduction roasting tests, two continuous 10-day campaigns, and one 12-hour test all in a 34-inch by 36¬-foot rotary kiln. The objective of the project was achieved in producing an iron concentrate containing about 66 percent iron with over 90 percent iron recovery. Desirable operating conditions were developed in the kiln but optimum process variables were not completely established.

Jan 1, 1969

The granulites and related "dry" igneous rocks tend to be passed over during exploration for metals. Preliminary data show that high trace concentrations of Cu, Ni and other metals of economic interest are common in orthopyroxene and some of the other granulite minerals. During retrograde metamorphism of granulites, pyroxenes are converted to amphiboles or phyllosilicates, and significant quantities of metals are released and become available for incorporat- ion into new mineral phases, or may be transported elsewhere. S--rich scapolite is the normal scapolite of granulites. Its retrogression to secondary plagioclase provides much S which scavenges Cu, Ni, etc., from pyroxenes undergoing retrogression, thus producing ore-fluids. A full analysis of a scapolite with 5.10 per cent S03 is given. Beryl is unstable in granulites but the rare Be mineral taaffeite (Al4Pf (Be,Fe)O3) with 5,50 per cent BeO occurs in central Australia as a pale green Tabular mineral easily confused for chloritoid in hand-specimen. Metamorphic differentiation and mobilization of primary ore concentrations are major factors controlling their post-granulite attitude and grade.

Jan 1, 1969

By D. Montagna

Large quantities of zinc-base die-cast scrap containing 3 to 4 percent aluminum are presently refined by mixing with other secondary zinc products such as galvanizer dross and distilling in a retort. The Bureau of Mines investigated alternate aluminum removal processes utilizing low-cost fluxes. Basing the investigation on the metallurgical principle that a reactive impurity can be removed from a less reactive matrix metal by melting under a flux consisting of a salt or oxide of a metal less reactive than the impurity, three separate flux mixtures were studied. In the first flux system, [FeCl3] + NaCI (the active component is underlined), part of the aluminum is' eliminated by the formation of volatile aluminum chloride and the rest by the creation of a high melting point iron-aluminum intermetallic compound that can be separated from the remaining molten zinc by filtration. Aluminum removal in the second system, [Fe203] + CaClz + NaCI, is effected by the formation of the iron-aluminum compound and possibly aluminum oxide. The third system, spent sal skimmings--essentially a mixture of [ZnCl2], NH4Cl, and ZnO, depends solely on aluminum chloride formation. Of the three flux systems investigated, the best results obtained were with system 3. Aluminum was reduced in the molten die-cast alloy from 4.10 percent to 0.01 percent when depleted sal skimmings were used.

Jan 1, 1969

By Mark O. Halverson

In this paper, some applications of remote terminal computer time-sharing to mining geophysics are discussed. Examples are presented and evaluations are made. The evaluations are based largely upon one and one-half years experience with the General Electric Mark I and Mark II systems. In both systems, a terminal is placed at the user's location. The terminal consists of a teletypewriter and a paper tape system. Input is on paper tape or keyboard. The output is printed on paper and may also be placed on paper tape if desired. Available languages consist of BASIC, FORTRAN and ALGOL. The time-sharing system is considered very useful. In large part this is because the terminal is available for use as the need arises, much like a slide rule or calculator. Useage is varied. Geophysical field data can be put on paper tape and be computer reduced, for the most part by clerical personnel. Time-sharing is excellent for statistical studies. The system is adequate to generate geophysical parameter values derived from simple models of the subsurface and to make comparisons with field data.

Jan 1, 1969

By L. Michael Kaas

The mineral industries are striving to expand their use of the computer's ability to process exploration, geologic, and mining data. Digital plotting programs and devices are continuing to gain wider use for lowering the cost and the time required to produce many types of maps and cross-sections commonly used by mining companies. Mineral engineers and geologists will soon be using a new tool, computer graphics, to greatly facilitate the use of computers to solve a large number of industry problems. By means of cathode-ray-tube (CRT) displays, an engineer can search large files of geologic data; display elevation contours, geologic structures, and ore grades; construct geologic and mining cross-sections; and compute tonnages for mining plans or ore reserve estimates. Examples of both digital plotting and cathode-ray-tube problem solving approaches are discussed in this paper.

Jan 1, 1969

By Frederick J. Hoppe

The productivity or usefulness of a Front End Loader rises rapidly with its size. Loading is more aggressive. Mother Nature is defeated more easily. Larger tires get a better grip on the ground. Hauling equipment can be loaded more quickly. The productivity or usefulness of a Front End Loader is proportional to its working speed, both automotive and hydraulic. The price of speed is horsepower, which rises rapidly with speed. Horsepower saved shows up as increased productivity. This paper covers the development of the D600, KW-Dart's articulated loader, which combines massive size with new concepts to increase its working speed. Featured is a discussion of the horsepower saving, "Balanced Boom."

Jan 1, 1969

By R. K. Koch

The vapor pressures of liquid titanium and liquid platinum were studied under the Bureau of Mines program for the determination of thermochemical data on the transition elements. Beating was done in a 2-kw electron-beam furnace with a water-cooled copper crucible to contain the melt. Temperature was measured with a disappearing-filament brightness pyrometer. The vapor pressure data for liquid titanium in the range 2,008° to 2,379° K are represented by the least-squares derived equation [log p (atm) = 6.268 - 22, T] Two series of determinations were made on liquid platinum. Results of the first series are expressed by the least-squares derived equation [log p (atm) = 8.240 - 30,460T] for the range 2,0860 to 2,4200 K. Results of the second series are expressed by the least-squares derived equation [log p (atm) = 7.795 - 29,580 T] in the range 2,045° to 2,442° K. Third-law evaluations of each series gave mean standard heats of sublimation of 112.4 kcal/mole, 133.3 kcal/mole, and 133.8 kcal/mole for titanium first series on platinum, and second series on platinum, respectively. Second-law heats of sublimation and equations for heats, free energies, and entropies of vaporization were obtained from sigma plots of each data set.

Jan 1, 1969

By Robert Isaac Williams

The Questa mine is located in northern New Mexico on the western flank of the Sangre de Cristo Mountains about 75 miles north of Santa Fe and ranges in elevation from 6,800 to 10,000 ft. GEOLOGY The rocks of the area are Precambrian metamorphic rocks covered by Tertiary extrusive andesites and rhyolites and intruded by Tertiary granitic rocks. These Tertiary intrusives are localized in an east-west belt about 20 miles long and 2 miles wide, extending from Questa to Red River. The Questa orebody is associated with a north-trending aplite intrusive about 3/4 mile wide and 1 3/4 miles long. In the area of mineralization, the metasomatism of aplite has produced a mottled porphyritic rock. This hybrid rock has caused much of the controversy about the Questa mine intrusive. The orebody is essentially monomineralic and consists of numerous closely spaced molybdenite veins varying in thickness from ' 'paint" to 3 in. and occurring in both the andesites and the intrusive. In the intrusive the veins are often siliceous. Common alteration minerals associated with the molybdenite ore vary with the original composition of the rock. The alteration of the aplite has produced phenocrysts of orthoclase, albite, oligoclase and biotite; veins of quartz, orthoclase, and biotite; and graphic myrmekitic intergrowths of feldspars and quartz. In many places near the contact, the andesite is indistinguishable from the metasomatized aplite. The dominant alteration in andesite, however, is the development of fine-grained biotite with apatite, rutile, and sphene.

Jan 1, 1969

By J. H. Swisher

When an Fe-Mn alloy is internally oxidized, the inclusions formed are MnO which contains some dissolzled FeO. In the internal oxidation reaction, not all of the manganese is oxidized; some remains in solid solution as a result of the high Mn-0 solubility product in iron. Taking these factors into consideration, the rate of internal oxidation of an Fe-1.0 pct Mn alloy is computed as a function of temperature, using available thermodynanzic data and recently published data for the solubility and diffusivity of oxygen in iron. The predicted and experimentally determined rates for the temperature range from 950 to 1350°C are in good agreement. ThE rates of internal oxidation of austenitic Fe-A1 and Fe-Si alloys have been studied extensively.1"4 Schenck et al. report the results of a few experiments with Fe-Mn alloys at 854" and 956C, and Bradford5 has studied the rate of internal oxidation of commercial alloys containing manganese in the temperature range from 677" to 899°C. When Fe-Mn alloys are internally oxidized, the inclusions formed are solutions of FeO in MnO, the composition depending on the experimental conditions. Since the thermodynamics of the Fe-Mn and FeO-MnO systems have been investigated,6"9 and since the solubility and diffusion coefficient of oxygen in y iron have been determined recently,' it is possible to predict the rate of internal oxidation from known data. The calculations used in predicting the rate of internal oxidation will first be outlined, then the results of the prediction will be compared with the experimental results of this investigation. PREDICTION OF PERMEABILITY FROM THERMODYNAMIC AND DIFFUSIVITY DATA Oxygen is provided for internal oxidation in these experiments by the dissociation of water vapor on the surface of the alloy. The dissociation reaction is: + H2(g) + [O] [1] where [0] denotes oxygen in solution. The equilibrium constant for this reaction is known as a function of temperature:' log As oxygen diffuses into the alloy, oxide inclusions are formed which are MnO with some FeO in solid solution. The reactions occurring are: [Mn] + [0] = (MnO) [31 and [Fe] + [0] = (FeO) [41 where [ Mn] is manganese dissolved in iron and (FeO) is iron oxide dissolved in MnO. The overall reactions may be written as follows: [Mn] + HOte) = (MnO) + H2(£) [5] and [Fe] + H20(g) = (FeO) + Hz(R) [61 The standard free-energy changes and equilibrium constants for Reactions [5] and [6] are known.6 Therefore the equilibrium constants for Reactions [3] and [4] may be obtained by combining known thermodynamic data for Reactions [I], [5], and [6]. For Reactions [3] and [4]: K = and For the present purpose, both the Fe-Mn7,8 and FeO-~n0' systems can be considered to be ideal, i.e., [amn] = [NM~] and (aFeO) = (NM~~) = 1 - (NFeO) where the Ns are mole fractions. These relations, together with Eqs. [I] and [8], permit us to compute both the oxide and metal compositions as a function of temperature and oxygen potential at any point in the specimen. For cases where the oxygen concentration gradient between the surface and the subscale-base metal interface is linear, the kinetics of internal oxidation is an application of Fick's first law: where dn/dt is the instantaneous flux of oxygen into the specimen, g-atom per sq cm sec; 6 is the instantaneous thickness of the subscale, cm; Do is the diffusion coefficient of oxygen in iron, sq cm per sec; p is density of iron, g per cu cm; h[%O] is the oxygen concentration difference between the surface and sub-scale-base metal interface, wt pct. B6hm and ~ahlweit" derived an exact solution to the diffusion equation for systems in which there is a stoichiometric oxide formed. They showed that the oxygen concentration gradient is given by a rather complex error function relation. For the Fe-Mn-0 system and for most other systems that have been studied, however, variations in oxide compositions are small and rates of internal oxidation are sufficiently slow that the deviation from linearity in the concentration gradient of oxygen is negligible. The mass of oxygen transported across a unit area of the specimen for the total time of the experiment is given by the mass balance equation:

Jan 1, 1969

By G. A. Josey, T. P. Floridis

INFORMATION on the effect of alloying elements on the solubility of sulfur in solid iron alloys is limited to the effects of carbon,' manganese,' and Silicon.Refs. 3 to 5 provide information on the combined effects of silicon and manganese also. In the investigation reported in this communication, the effect of chromium, copper, nickel, and vanadium on the solubility of sulfur in solid iron was determined at 1000". The specimens for this investigation were prepared by melting electrolytic iron and reagent-grade alloying additions in an induction furnace under an argon atmosphere. The melts were allowed to solidify in

Jan 1, 1969

By K. H. Coats

This paper describes an approximate technique for handling the problem of percolation of evolved gas upwards through the oil column in computer simulation of natural depletion. This technique has been incorporated into a general model for simulating three-phase flow in one, two or three dimensions. The mathematical model for perfonning these calculations is described in detail in a separate article. 1 While vertical gas percolation occurs during natural depletion in a reservoir of any configuration, it is especially pronounced in the pinnacle reel or bioherm. The reel may have an areal extent less than one section with a thickness of up to 800 It. The upward percolation of evolved gas often limits the time step in calculations as little as one day. Using this small time step, computer expense for a single 30- to 40-year simulation on even a one-dimensional basis has exceeded 92,000. The method described here for handling the percolation allows time steps of 60 days or more (depending on reservoir size and production rate), resulting in a considerable reduction in computing expense. The method also allows calculations in which secondary gas caps build up in tight zones of the oil column below the main gas cap. The validity of the method is indicated in connection with an example pinnacle reef field. Calculated results using a 2-day time step (where the method in question is not needed and not invoked) are compared with results using the method and a 60-day time step. The comparison shows good agreement. Results from one- and two-dimensional simulations of the reef are presented along with the corresponding computing times on the CDC 6600 computer. A three-dimensional simulation was also performed and the required computing time is given. INTRODUCTION A general analysis for simulating three-dimensional, three-phase flow in reservoirs has been developed. 1 In applying this model to reservoirs undergoing natural depletion, a time-step restriction was encountered due to the flow of evolved gas upwards through the oil column and toward the gas cap. A remedy for this problem has been incorporated in the analysis and is described here. The time-step restriction is encountered only in calculations which include flow in the vertical or near vertical direction. The restriction occurs to some extent in reservoirs of any geometrical configuration, but it is especially pronounced in the pinnacle reef or bioherm where the ratio of thickness to areal extent is unusually large. This paper describes the problem, a method of handling it, incorporation of the method in the three-dimensional, three-phase model, and, finally, a test of the method's validity in an application to an example pinnacle reef reservoir. THE PROBLEM During the early stages of reservoir depletion, pressure falls below bubble point in progressively lower regions of the oil column. Continuing production results in evolution of dissoived gas throughout the oil column. This gas then percolates upwards toward the top of the reservoir. If the gas encounters a sufficiently low permeability zone in its travel upward, then it can accumulate and form a secondary gas cap. If no zones are encountered which are tight enough to hold gas against the gravity forces, then the gas travels on until it reaches the top of the sand or the main gas cap. Numerical simulation of multiphase flow in reservoirs is generally performed with time steps such that the flow of a fluid into or out of a block in one time step is a fraction (considerably less than one) of the total amount of that fluid present in the block. Use of a time step where the time step's

Jan 1, 1969

The mineral composition and microtextures of a cubanite-rich sulphide ore are described and cubanite breakdown textures are illustrated. The occurrence of cubanite in macroscopic dimensions is featured.The origin of the ore is tentatively related to a composite ore fluid arising initially from an ultramalic source and extracting some metal components from an older granite mass along the margin of which the serpentine was emplaced. Some alternative possibilities as to origin are also mentioned.

Jan 1, 1969