Search Documents

By A. T. Corey, C. H. Rathjens

INTRODUCTION Although the oil industry has been aware of the directional variability of permeability in porous rock, the directional variability of relative permeability has been largely ignored. Yet it is apparent that such an effect must be present in a system in which the distribution of oil and gas within the porous matrix is controlled by capillary forces. It is easy to visualize a rock composed of layers of fine and coarse material such that gas flow across the bedding planes would take place only after the average oil saturation had been reduced to a very low value. The fine layers, because of their greater capillarity, would remain saturated and act as barriers to the flow of gas after the coarse layers had been desaturated. Flow of gas parallel to the bedding planes would obviously take place at a much greater liquid saturation. Without more complete information concerning the geology of a reservoir than is generally available, it is not possible to predict exactly how such phenomena would affect the over-all performance of an oil field. It is possible, however, to predict qualitatively the effect of stratification on relative permeability measurements made on laboratory cores. In this investigation the effect of stratification was studied analytically by assuming that two porous materials with different capillary pressure-desaturation curves (but identical relative permeability curves) were in contact and in capillary equilibrium. As a qualitative check on the analytical results, cores having various degrees of visible stratification were used for relative permeability measurements made with fluids flowing both parallel and perpendicular to the bedding planes. A quantitative check was considered impractical because of the difficulty of devising models in which two materials of predetermined properties could be joined without the plane of contact becoming a discontinuity. THEORETICAL CONSIDERATIONS AND ASSUMPTIONS The assumption of capillary equilibrium in an oil-gas system implies that the difference in pressure between oil and gas is everywhere the same. This means that the curvature of the interfaces must be everywhere the same in order to satisfy the equation where PC is the pressure difference between phases, y the interfacial tension and r, and r2 are the major and minor radii of curvature. Depending on the pore size distribution of coarse and fine layers, the volumetric percentages of oil and gas in these layers will differ when equilibrium exists. The exact relationship can only be determined by obtaining the complete capillary pressure-desaturation curves for each of the porous materials in contact. It has been pointed out elsewhere' that the capillary pressure-desaturation curves of sedimentary porous materials can often be approximated by the relation where C is a constant and Soe is the effective saturation to oil based on a percentage of the pore volume effective to flow. In the same paper it was indicated that, as a first approximation, the values of oil relative permeability are given by and the values of gas relative permeability by For this analysis Eqs. 2, 3, and 4 were assumed to apply to each of two components of a hypothetical porous rock in capillary equilibrium. It was also assumed that each of the components had a residual wetting phase saturation of 20 per cent so that 80 per cent of the total pore volume was effective to flow. The permeability of the coarse stratum was taken as 100, and its displacement pressure was such that C in Eq. 2 had the numerical value of 1. The corresponding values for the fine stratum were 10 for the permeability and 10 for C. Units are not specified because they do not enter into the final results. The choice of the permeabilities and displacement pressure ratios was made to expedite the calculations. Any reasonable rock properties could have been chosen without changing the results qualitatively. Several arrangements of the two components were studied. Table 1 summarizes the resultant permeabilities obtained for four types of arrangement. RELATIVE PERMEABILITY CALCULATIONS The first step in the computation of relative permeability for the composite cores was the plotting of the capillary pressure-desaturation curves and the relative permeability curves for the individual components according to Eqs. 2, 3, and 4. At arbitrary

Jan 1, 1957

By J. W. Freeman, E. E. Reynolds, A. E. White

Fram a study of 63 systematic alloy modifications it was found that molybdenum, tungsten, and columbium, added individually or simultaneously, and increases in chromium cause major improvements in 1200°F rupture strengths of Cr-Ni-Co-Fe base alloys. Rupture strengths were a function of the effect of composition modifications on both the inherent creep resistance and the amount of deformation the alloy would tolerate before fracture. THIS paper describes the results of an investigation of a series of alloys with systematic variations of the chemical composition of the following basic alloy: C, 0.15: Mn, 1.7; Si, 0.5; Cr, 20.0; NI, 20.0: Co, 20.0: Mo, 3.0; W, 2.0; Cb, 1.0; N, 0.12; Fe, 32.0 pct. The 62 modifications of this alloy were produced under conditions which minimized all factors influencing properties at high temperatures except composition. Melting, fabrication, and heat treatment were carefully maintained constant. Stress-rupture properties at 1200°F were used as the primary criteria of evaluation of the alloy. The objective of the study was to obtain data for determining the fundamental role of the influence of alloying elements on properties of heat-resistant alloys at high temperatures. In addition the results should be useful in determining optimum chemical compositions, the sensitivity of properties to variations in composition, and the degree to which alloy content could be reduced while retaining worthwhile properties. It is difficult or impossible to develop correlations between properties at high temperatures and systematic variations in chemical composition from published data for wrought heat-resistant alloys developed for gas turbines.' ' The main reason for this is the extreme dependence of the properties on conditions of treatment of the alloys." In most cases variation in final treatments between alloys so influences the properties that the influence of chemical composition is obscured. In addition it is recognized Table I. Basic Alloy and Some Modifications Used Basic AllOy, Variations in Element Pct Composition, Pct C 0.15 0.08. 0.40. 0.60 Mn 1.1 0.03,0.25.0.50,1.0,2.5 S1 0.50 1.2, 1.6 Cr 20.0 10, 30 Ni 20.0 0, 10,30 Co 20.0 0. 10, 30 MO 3.0 0. 1.2.3, 5, 7 W 2.0 0, 1, 5, 1 Cb 1.0 0.2,4,6 N 0.12 0.004, 0.08, 0.18 Fe 32.0 that variations exist between heats of the same alloy which are related to melting practice and that there is a strong possibility that conditions of hot working influence response to final treatments. The development of heat-resistant alloys has been based on the gradual accumulation of data roughly related to composition from extensive testing programs. There is every reason to believe that in most cases the optimum compositions have been achieved by this procedure in the alloys commercially available. There are, however, very little data showing the influence of systematic variations of composition free from the influence of other factors, particularly for alloys of the type investigated. Several investigators of cast alloys have demonstrated compositional effects, notably Grant,1-6 Epremian,t Guy,8 and Harder and Gow.9 Sykes10 eviewed the work on the wrought alloy Rex 78 and the systematic variations of carbon, copper, molybdenum, and cobalt leading to the development of the stronger Rex 337A alloy. From the papers by Wilson11 and Henry12 it is possible to deduce the beneficial effect of substituting cobalt for iron in 0.45 pct C-20 pct Cr-20 pct Ni-4 pct Mo-4 pct W- 4 pct Cb alloys. Wilson mentioned but did not present the extensive compositional studies involved in developing these alloys. Binder" showed optimum properties for 3, 2, and I pct, respectively, for molybdenum, tungsten, and columbium in 20 pct Cr-20 pct Ni-20 pct Co-30 pct Fe alloys for limited systematic variations of these

Jan 1, 1953

By C. E. Ells, A. Sawatzky

The effect of neutron irradiation on tensile properties of the Zr-2.5 wt pet Nb-0.5 wt pet Cu alloy has been evaluated for integrated neutron fluxes up to 3 x 1020 n per sq cm (E > 1 Mev). Specimen temperature during irradiation was normally 300°C. The material was studied both in annealed and in quenched and aged conditions. When the alloy (fab-ricatedfrom sponge zirconium) is quenched from a temperature =40°C below the a + ß/ß transformation and aged 6 hr at 535oC, neutron irradiation up to the maximum studied has little effect 071 reduction in area, although there is a drop in uniform elongation and the yield strength is increased by = 25 pct. Conversely, irradiation of the alloy in an annealed condition can result in much greater changes in properties, with yield-strength increases of up to 200 pet. The dose dependence of irradiation hardening obeys the saturation equntion proposed by Makin and Minter to explain hardening in copper and nickel: Act = C[1 - exp(-Døt)/1/2 A rnodel is proposed to explain the effect of niobium content and metallurgical condition on irradiatiotz behavior. THE Zr-Nb alloys form a class of material susceptible to marked strengthening by quench and age heat treatment. With niobium concentrations as low as 2.5 wt pet, strengths can be developed which are double those of annealed zircaloy-2.1 Since the neutron-capture cross section of the Zr-2.5 wt pet Nb alloy is nearly identical to that of Zircaloy-2, significant gains in power-reactor neutron economy could be obtained by replacing stressed in-reactor components of Zircaloy-2 with the heat-treated Zr-Nb alloy. When the solution heat-treatment temperature is in the high (a + ß) phase, then the tensile properties of the quenched (and aged at 500°C) material have been shown to be relatively insensitive to neutron irradiation.1-3 Addition of small quantities of copper to the binary Zr-2.5 wt pet Nb alloy gives a useful reduction to the corrosion rate in air and carbon di- oxide,&apos; and a ternary alloy of composition Zr-2.5 wt pet Nb-0.5 wt pet Cu was chosen for one component in the Douglas Point Reactor. This application utilizes both the corrosion resistance to a moist air-carbon dioxide environment and the strength developed by a quench and age heat treatment. By suitable choice of solution heat treatment and aging temperatures, tensile properties of the ternary alloy can be made nearly identical to those considered as optimum for the binary alloy. One effect of the copper addition, however, is to give marked changes in the aging kinetics of the alloy.&apos; For this reason, it appeared that neutron irradiation might promote overaging in the ternary alloy, particularly at the elevated service temperature (=300°C). This paper describes the effect of neutron irradiation on the tensile behavior of the Zr-2.5 wt pet Nb-0.5 wt pet Cu alloy in several metallurgical conditions pertinent to its use as a reactor material. The metallurgical conditions included those for which aging occurs in the unirradiated material, and fully annealed conditions for which a higher degree of thermal stability could be expected. 1) EXPERIMENTAL Two batches of alloy were used, both made from sponge zirconium. The fabricator&apos;s analysis of principle constituents showed no significant difference between the two batches, Table I. The rod material was rolled to 1.5 in. diam, and then swaged to final sizes of 0.5 to 0.375 in. diam at nominal temperatures of 785" and 625°C for the AM and AS batches, respectively. Subsequent hydrogen analyses at CRNL confirmed that the hydrogen concentration of the as-received rod was indeed <20 ppm. However, it was found that the (a + ß)/ß transus for the AM material was approximately 25°C lower than that of the AS material, indicating rather more difference in oxygen concentration than given in Table I. Except for a small difference in yield-point behavior, the tensile prop-

Jan 1, 1965

By A. G. Guy

Thermodynamic activity rather than chemical composition is basic to the analysis of diffusion. This is the essential conclusion reached by Darken1-3 and by Birchenall and Mehl.4 If so, it is reasonable to expect that diffusion calculations should be carried out using activities. This paper will illustrate a method of employing activities in analyzing interstitial diffusion. For interstitial diffusion (of carbon in iron and probably generally) the diffusion of the solvent atoms, volume changes, and other such complexities not explicitly provided for in the usual diffusion equations cause a negligible error of the order of 1 pct. Therefore it is relatively easy to test the validity of an analysis based on activities in the instance of interstitial diffusion. The mathematical treatment of the more complex substitutional diffusion could be simplified once the adequacy of the activity method had been established for substitutional diffusion. Thus the present work has the purpose not only of providing a superior analysis of interstitial diffusion but also of indicating a possible approach to the more important question of an adequate treatment of substitutional diffusion. Development of the Diffusion Equation Although activities are convenient quantities to determine experimentally, they are derived from the more fundamental fugacities.9 Possible confusion in using activities can be avoided by basing equations involving their use on the relationship, dP = -D1d/ax dt [I] where, dP = the number of grams of solute crossing one cm2 in the time dt sec. D, = the diffusion constant for use with fugacities; it is assumed to be constant at a given temperature, and has units of secs. f = the fugacity of the solute in the solid solution; the units are those of pressure, dynes per cm2. x = the distance in cm. Eq 1 can be considered to be the basic form of the first Fick Law for one-dimensional diffusion. In order to convert Eq 1 into an equation involving activities, a choice of a standard or reference fugacity, f°, must be made: then. a = f [2] where a is the activity corresponding to the fugacity f. However, the value of a is also determined by the relation, a = a\c [3] where a is the activity coefficient and c is the concentration. For a given value of a (such as a = 1) it is evident that the units used for expressing concentration will affect the numerical value of a corresponding to a given amount of carbon dissolved in y-iron. Since weight per cent concentration is so widely used, this unit will be adopted as the standard in this paper. When the value off given by Eq 2 is substituted in Eq 1, the first Fick Law becomes, dP= -D,f°dl [4] Since D,f° is constant (at a given temperature), the second form of Fick&apos;s Law can be obtained in the usual manner, -dt D&apos;f° dx [5] where c&apos; is the concentration in g per cm3. In terms of weight per cent concentration, c, Eq 5 becomes, dc _ 100 d2a Tt = T /f w [6] where p is the density of the solid solution and will be assumed to be constant. In order to simplify Eq6, it is necessary to consider the standard state, f°. Standard states are chosen for convenience. There appears to be one especially convenient choice for the analysis of interstitial diffusion. If the activity coefficient, a, is set equal to unity at infinite dilution of the solute, 100 then — 100D1fl° is approximately equal P to the usual diffusion constant in dilute solutions. This is true since the activity, a, in Eq 6 can be replaced by the concentration, c, with little error in an infinitely dilute solution.* Here f1° is the standard fugacity necessary to achieve this standard state. It is proposed that a diffusion constant Da&apos; be defined, Da1=100/pD1f1° [7] where A is the metal whose diffusion is * The second derivative of a is equal to the second derivative of c in dilute solution only if -da/dr- is. zero. Although theory does not require this condition, in the systems C in Fe, and Zn in Cu, dc is in fact found to be essentially zero.

Jan 1, 1950

By W. O. Philbrook, H. E. Burner, F. S. Manning

The equation of continuity for the hydrogen reduction of hematite in a fixed bed of closely-sized particles is solved assuming a flat velocity profile, negligible temperature gradients, md negligible axial diffusion. A kinetic expression from the literature is used which assumes the reduction process is controlled at the oxide-metal interface. The integral fractional conversion is computed, and the importance of particle diameter, flow mte, temperature, bulk axial diffusion, and inter- and intra-particle mass trarnsfer is predicted. Comparison with experimental data suggests one or more additional undefined variable(s) is significant. Equipment modifications are suggested for fidture experimentcll work. WITH the increasing desire to "design" the ore feed of a blast furnace, it has become necessary to define more quantitatively the process variables and to evaluate the extent of their control upon the reduction process. The bulk of the work reported in the literature has centered about the reduction of single particles; however, the blast furnace process more closely resembles reduction of a fixed bed than that of single particles. To simplify this analysis, the isothermal reduction by hydrogen of a natural hematite in a fixed bed of closely-sized particles is examined. With appropriate modifications in the kinetic and flow rate expressions, the approach employed should be extensible to mixtures of carbon monoxide, hydrogen, and inert nitrogen, and to adi-abatic conditions, thus approaching the blast furnace process more closely. The fixed bed variables which are investigated in this study are particle diameter, flow rate, temperature, and bulk axial diffusion, as well as inter-and intra-particle mass transfer. Of particular interest is the effect of particle size upon fractional conversion for a bed operating under constant pressure drop. The constant pressure drop concept bears a close resemblance to the blast furnace process, where regions of larger particles may provide a greater permeability to divert most of the reducing gas flow from the regions of smaller ore particles. Under such conditions one would expect qualitatively that, for sufficiently small particle sizes where the flow rate is small and the total surface area is large, the reduction rate would be limited by the gas flow rate. Initially the gas entering the bed reacts very rapidly and the composition approaches the equilibrium conversion value a short distance from the reactor inlet. This concentration "wave-front" then moves up the bed as the flow of reducing gas is continued, leaving in its wake an ever increasing layer of reduced particles. Providing that the reducing gas leaves the bed at equilibrium concentrations, conversion is expected to increase with an increase in flow rate, and hence with an increase in particle diemater. Conversely, for sufficiently large particles the flow rate will be high but the surface area available for reaction will limit the reduction rate. Since the total surface area decreases with an increase in particle size, conversion, in this case, is expected to decrease as the particle diameter is increased. It is clear then that between these two extremes an optimum particle diameter must exist. THEORETICAL DEVELOPMENT Above the temperature where wiistite is stable (about 560°C) the reduction of hematite proceeds through the series Fe2O3/Fe3O4/FeO/Fe. On the basis of microstructural studies and earlier work, Edström1 has postulated a mechanism for the reduction of an ideal hematite lattice: 1) Iron phase formation at the boundary between wiistite, iron, and gas (for sufficiently porous iron). 2) Diffusion of iron across a wustite layer. 3) Phase boundary reaction of Fe3O4 to FeO. 4) Diffusion across a dense magnetite layer. 5) Phase-boundary reaction of Fe2O3 to Fe3O4. According to the mechanism postulated, gaseous reaction product has to be removed only at the boundaries between wüstite and iron or gas, the only possible exception being enveloped Fe3O4 specks. On this basis McKewan2,3 has proposed that one may consider the hematite reduction process as a simple two-phase system of oxide/metal. If, for FeO and Fe3O4 layers of negligible thickness, the hypothesis is made that the rate of formation of a

Jan 1, 1963

By J. J. Pfeiffer, A. A. Oming, J. W. Myers

The method used to process washery refuse for use as a building material aggregate is described. Results of studies made in investigating this process are summarized. The Bureau of Mines, in cooperation with two commercial coal producers, studied the feasibility of using coal washery refuse in manufacturing lightweight aggregate for building materials. Lightweight aggregate is used in making masonry products of lower bulk density than can be produced from conventional aggregate such as sand, gravel, and crushed rock. Among the advantages of lightweight masonry are reduced deadweight loads which reduce materials cost, lower cost of construction labor, and improved thermal and sound insulation. A process for making lightweight concrete aggregate from coal-washery refuse appeared attractive because of the availability of raw material which contained sufficient fuel to carry out the processing. Vast quantities of such refuse have accumulated throughout the coal-producing areas of the U.S. Generally it is stored in unsightly piles, creating a problem of air pollution and using space that could be devoted to other purposes. Conversion of the refuse to lightweight aggregate could remedy these conditions and effectively utilize a material normally considered waste. The first part of the program was conducted in cooperation with the Truax-Traer Coal Co. to determine the feasibility of utilizing refuse from the company&apos;s coal preparation plant at Ceredo, W. Va., in manufacturing lightweight aggregate. The aggregate was to be used in making cement blocks comparable to those known in the building trades as cinder blocks. This project was completed in 1954. Based on the results of this work, the cooperator began commercial development and formed the Trulite Corp. for this purpose. A plant capable of handling 120 tpd of refuse was constructed for approximately $120,000. It began operating on June 1, 1955.&apos; In 1959 a similar study was made, under a cooperative agreement with the Carbon Fuel Co. of Charleston, W. Va., to rrscertain the suitability of refuse from its plant for conversion to lightweight aggregate. Under this program various modifications of equipment were made, and the investigation of process variables begun in the previous project was extended. PROCESSING METHOD The formation of lightweight aggregate from coal-washery refuse js governed, in part, by the fusion characteristics of the mineral constituents of the refuse when heated above their softening temperatures. The heat required for this purpose is supplied by the combustihle matter in the refuse. Some of this combustible matter is coal with high ash content, but most of it is carbonaceous shale. Whereas the density of aggregate prepared from clay or shale depends largely on expansion, the density of aggregate made from refuse is reduced by burning out the combustible matter and, to a lesser extent, by expansion of the residue in the fuel bed. A chain-grate stoker with an updraft air system was chosen for the conversion process in preference to a downdraft sintering machine or a rotary kiln. In down-draft sintering, which operates on the overfeed burning principle, the volatile matter from the refuse would distill off without burning, resulting in inefficient combustion and interfering with proper operation of the equipment. The rotary kiln is designed primarily for processing material that lacks, or is low in, combustihle matter with the heat of reaction supplied from external sources. Washery refuse, which is relatively high in carbonaceous material, could burn at the wall of the kiln causing clinker formation with resultant decrease in throughput and increase in production cost. With an up-draft air stream on a chain or traveling-grate stoker, the volatile matter would be carried into the ignition and burning zones where the heat could be utilized more efficiently. The combustion procedure in this type of unit may be considered as consisting of four periods: 1) initial ignition of the top of the bed by furnace radiation, 2) ignition travel through the bed until the whole bed is ignited, 3) burnout of the combustible material with high air rates to fuse and expand the residue, 4) a cooling period. Investigation of the process required a study of the variables that affect both

Jan 1, 1962

By Ralph Hultgren, Prodyot Roy

Manganese vapor pressures from 1250° to 1500°K were measured by conventional Knudsen and torsion-effusion methods in twelve Fe-Mn alloys with compositions from 9 to 80 at. pct Mn. The Knudsen re-sults agveed approximately with previous measurements found in the literature which indicated the solutions were nearly ideal. However, except for the higher manganese compositions, the initial torsion readings indicated much (up to 50 pct) higher vapor presszcres than the Knudsen method. These high readings decreased steadily with time. The results are interpreted as due to depletion of surface concentration of manganese during evaporation. Thus the initial torsion readings are the most nearly correct and Knudsen methods on alloys are to be regarded with suspicion unless it can be demonstrated that diffusion rates are rapid enough to replenish surface concentrations depleted by evaporation. For the lowest manganese contents and highest temperatures, depletion causes initial torsion readings to be too low. FEW thermodynamic data are available for alloys of the high-melting transition metals due primarily to experimental difficulties at the high temperatures involved. Measurement of vapor pressure is one of the most promising techniques to be applied to this problem. From the vapor pressure of a component of an alloy phase, its activity or its partial molar Gibbs energy can be directly calculated. From measurements over a sufficient range of temperatures and compositions the partial and integral Gibbs energies, enthalpies, and entropies may be determined through Gibbs-Duhem integration. For solid phases of variable composition, which commonly are found in alloys, measured vapor pressures are low because of depletion of the surface concentration of the more volatile component as it is selectively vaporized. Diffusion from the interior tends to restore the depleted concentration. The seriousness of the effect depends on the relative rates of diffusion and vaporization. Although depletion has been recognized as a factor,&apos; direct measurements of its rate are lacking and quantitative estimates of its effect on measurements are uncertain. In the present work the vapor pressures above a series of Fe-Mn alloys have been measured by the torsion-effusion method. This method permits a measurement of the vapor pressure of manganese as soon as the sample comes to temperature and continuously throughout the run. The rate of decline of the apparent vapor pressure measures the rate of depletion and the torsion reading at zero time should be correct. At high manganese concentrations (70 and 80 at. pct Mn), the torsion reading remained nearly constant; surface depletion was negligible. There was a slow drop in measured pressure due to bulk loss of manganese from the sample. However, for lower manganese contents, especially at high temperatures, the depletion effect was considerable and the apparent vapor pressure decreased steadily with time. For these alloys the Knudsen method should give pressures which are too low. The measurement of pressure in the Knudsen method cannot be made until enough manganese has been vaporized to be weighable. The Knudsen result is therefore an integrated average between the initial and final vapor pressures. To verify this, conventional Knudsen measurements were also made with results in approximate agreement with previously published Knudsen measurements for these alloys. As expected, for the lower manganese contents the initial torsion readings of pressure were as much as 50 pct higher than the Knudsen. THE APPARATUS The apparatus, shown in Fig. 1, is capable of operating to temperatures up to 2000°K. A vacuum of 5 x 10"6 mm Hg can be maintained. Temperatures were measured to +3° by a W-Re thermocouple placed in the dummy cell, N. The torsion cell consisted of two alumina crucibles with holes in their covers. They were placed in the molybdenum holder shown suspended from the torsion wire. Rectangular torsion wires (1 by 4 mm) were found to have superior sensitivity and less residual distortion than circular (2-mm-diameter) wires. The precision of angular measurement was found to be *0.025 deg; the variation is no doubt due to temperature fluctuation and vibration. For the Knudsen experiments the alumina cell was set on the support, N. SAMPLE PREPARATION 800-g samples were prepared by melting together electrolytic iron and electrolytic manganese and pouring the melt into a water-cooled copper mold. Melting and pouring were done under a helium at-

Jan 1, 1965

By K. R. Crocker

PRIMARY blasting at the open pit of the Cananea Consolidated Copper Co. has been developed to meet several conditions peculiar to the operation. Mining at Cananea is not confined to the open pit but is in conjunction with the underground mine which has been in operation for many years. Several old stopes included in the pit area were filled with waste obtained in stripping for the pit. The main haulage tunnel for the underground ore passes directly under one section of the pit and at present is but 180 ft below the pit floor. The pit as originally laid out called for an overall final pit limit slope of 60°. However, the upper 275 ft - a fairly soft oxidized capping, was carried at a 50° slope and the lower benches, which constitute 750 ft in combined height, are being mined to an over-all final slope of 60°.

Jan 6, 1950

By B. W. Veal, B. S. Chandrasekhar

THE Laue method of orienting single crystals requires the use of tables of angles between crystal-lographic planes, and between zone axes. Such tables have been published for cubic crystals, and some rhombohedra12 and hexagonal3 crystals. A table of angles between planes in ß-tin has also been published;4 these angles refer to only low-index planes, and are not particularly helpful in the actual indexing of spots on a Laue photograph. Furthermore, it is often useful to know the angles between zones which share a common low-index plane. We have calculated such tables of angles for ß-tin and indium, both between planes in principal zones, and between zones passing through principal planes. The results are given in Tables I and 11. The standard formulas were used, with c/a = 0.54554 for 0-tin, and 1.07586 for indium.

Jan 1, 1962

By V. N. Mackiw, D. E. Weir, D. J. I. Evans

The world&apos;s largest reserves of nickel are found in lateritic ores. Commercial exploitation has been carried out for almost one hundred years; however, none of the plants currently in operation produce pure nickel but rather ferronickel, matte, impure nickel oxide or mixed sulphide precipitate. Over the past decade Sherritt has studied various processes, including the modified reduction roast-ammonia-ammonium carbonate leach process, for the production of pure nickel powder and briquettes from laterites. Nickel can be reduced to the pure elemental form at elevated temperatures and pressures from impure basic nickel carbonate slurries. Several methods were studied in the development of the process for working up the pregnant leach liquor and basic nickel carbonate to produce feed for the aqueous hydrogen reduction step in which pure nickel powder is produced. This paper is concerned only with the laboratory studies using the basic nickel carbonate slurry reduction procedures.

Jan 1, 1973

By L. Adler

Any design of a mine roof in bedded deposits which ignores differential deflections at the supports can quickly lead to dangerous overstressing. As illustrated by the typical case presented on page 1029, if one support deflects only 0.1 in. more than its neighbor, stresses 50 pct higher than anticipated will result. Previous studies1-3 have assumed that the mine roof acts like a beam clamped at both ends, with no deflection or slope&apos; at the supports (Fig. 1). Essentially this assumes:5* ?A = ?, = ?b and dA = dB = 0. However, a host of conditions can cause differential sag or deflection at one support, and even an apparently trivial deviation can have major consequences in the beam: variations in spacing, size, and material of supports; local settlement due to floor heaving or lower mine openings; non-uniformity in wedging or bolt tensioning.8 In one particular instance concrete pillars, poured in place, shrank away from the roof upon setting.&apos; When understood and controlled, differential sag becomes a definite advantage, permitting the transfer of loads from one support to another. Where one support is a bolt or prop and the other a pillar or rib, the load on the former can be reduced while the load on the pillar or rib is increased. This means two things: 1) savings can be effected by using smaller artificial supports such as props and bolts and 2) loads can be increased in the rib to induce eventual failure for longwall caving. In this last respect, another factor underlying coal bumps becomes understandable. The design of a mine roof by beam theory assumes that when a bedded formation is homogeneous, iso-tropic, and elastic, it will act like a uniformly loaded beam, fixed or clamped at both ends,6 s in Fig. 1.&apos;" These are important assumptions. The uniform load is the weight of the beam itself or its body load. For a beam of unit width, this load becomes per foot: w = yhb. w = uniform load y = specific weight, lb per ft3 h = depth or thickness of beam, ft b = width of beam, ft Since b = 1, w - yh [ 1 ] A beam fixed at both supports is statically indeterminate to the second degree, that is, two equations in addition to those for static equilibrium are required to determine the reactions at the supports. These equations result from the geometry. From statistics: ^F,= 0;RA + RB — wL = 0 [2] wU S,M0 = 0 ; MB H------------MA — RBL = 0 [3] it Using superposition in which the left support is kept intact and the right one resolved into an active force and moment, we obtain from Fig. 2:" ?b = ?m + ? b2 + ?B3 = 0 dB = dB1 + dB2 + dB3 / wU Ma RBL \ V d°=(-r+-2—r)-zr [5] Where: Subscripts refer to position on axis of beam R = force, lb M = moment, ft-lb w = uniformly distributed load, lb per ft L = span, ft 0 = slope, radians d = deflection, ft E = Young&apos;s Modulus, lb per ft" 1 = Moment of inertia about the neutral axis, ft4 Solving Eqs. 2 through 5 for the four unknown reactions, we obtain: 2 L3 Fig. 3 is obtained by plotting MA and M, vs d, and then plotting RA and R, vs d, and adjusting their scales to coincide with MA and M, vs d,. It should be noted in Eqs. 6 and 7 that as dB increases, MA increases the same amount that Mn decreases. Eqs. 8 and 9 show a similar relationship between RA and R,. Proof that MA is the maximum absolute moment occurring in the beam is as follows: where M, = moment at x wx-When where M,., = maximum positive moment Eqs. 6 and 8 into above: For small dn: MA> M,,,

Jan 1, 1960

By J. E. Hilliard

AS is well-known, the length per unit volume, Lv, of any lineal feature (such as grain edges in a polycrys-talline specimen) can be estimated from a count of the number of point intersections with a random plane of polish. Specifically:&apos; Lv = 2NA in which NA is the expected number of intersections per unit area of the plane. As with all other measurements in quantitative metallography, the estimation of Lv is subject to a sampling error, the magnitude of which will depend on the number of observations. The analysis is usually performed by counting the number of intersections falling within a closed test figure of known area applied to a succession of different areas on the plane of polish. The test figure can be drawn directly on the screen of a projection microscope or else superimposed on the image of the structure by means of a reticule in a focusing eyepiece of an ordinary microscope. The use of micrographs is very inefficient unless, as in the present study, these are also needed for some other purpose. Consider an analysis in which n applications are made of a test figure within which the expected number of intersections is N. The expected total number of intersections counted is then given by NT = nN If (N) is the standard deviation in N for one application of the figure and if the latter is applied to non-overlapping regions so that the N&apos;s are independent in the statistical sense, then the standard deviation in Nt is given by: a(NT) = a(N)vn Eliminating n from Eqs. [2] and [3] we obtain: Thus, k will be a constant for a given specimen and test figure. From Eqs. [ I] and [4] we find As is to be expected, therefore, the fractional error in the estimation of Lv is inversely proportional to the square root of the total number of intersections counted. It is very desirable to have at least a rough estimate of k before starting an analysis, since it will then be possible to decide in advance by use of Eq. [5J how many counts should be made in order to achieve the desired accuracy in Lv. Without such information the investigator, unless lucky, will either do needless (and tedious) work or else obtain a result with too large an error. It is very easy to calculate k for the limiting case of N — 0, i.e., when the test figure is vanishingly small. In this case there is a negligible probability of two or more intersections falling within the figure. There are therefore only two possible outcomes: a count of one with a probability of, say, p and a zero count with a probability of (I-p). Thus and Therefore: We might expect k to depend on N and therefore on the size of the test figure. Such a dependence was observed2 in the analogous case involving the estimation of grain-boundary area per unit volume by a count of the number of intersections of a test circle with traces of the boundaries on the plane of polish. In this case the theoretical value of k at the limit N — 0 was 1.41. It was also predicted that k would display an initial sharp decrease with increasing N. Experimentally, k was observed to decrease to a value of 0.62 at N - 6 and thereafter remain approximately constant. Furthermore, the value of k for N > 6 was found to be insensitive to the form of the structure provided no pronounced anisotropy or grain-size variation was present. In order to determinewhether k would exhibit a similar dependence on N in the estimation of grain-

Jan 1, 1967

By W. W. Mullins

The development, by the mechanism of volume diffusion, of a grain boundary groove on an interface separating a solid phase and a saturated fluid phase is calculated under the following assumptions: 1) isotropy of the interfacial free energy &apos;s, 2) applicability of the Gibbs-Thompson formula relating curvature and chemical potential, 3) a nearly planar groove surface. The groove profile is found to have a .fixed shape and linear dimensions that are proportional to (time)1/3. The constant of proportionality is evaluated and involves ?s in such a way that the latter my be evaluated from grooving kinetics. A grain boundary that intersects an interphase interface will tend to produce a groove along the line of intersection (Fig. 1). The ultimate motivation for the formation of the groove is the reduction in interfacial free energ that occurs as the grain boundary contracts. Smith has discussed the way in which the free energy of the interfaces determines their dihedral angles at the groove root and has given a detailed analysis of the effect of the dihedral angles upon the nature the microstructures developed in polyphase alloys. Herring2 has extended the considerations to the case in which the free energy of an interface depends on its orientation and has thus added torque terms to the condition for interface equilibrium. Although the grain boundary and the interphase interface may quickly achieve the proper equilibrium angles where they intersect, as determined by the balance of equivalent tensions and torques,"2 the groove will generally continue to grow (though with decreasing speed). Thermodynamically, this continued growth reflects the fact that the lowest free energy cannot be achieved until the entire grain boundary has been eliminated; mechanistically, the growth occurs in response to inhomogeneities of chemical potential that are caused by the constraint of fixed dihedral angles at the groove root. The latter point is illustrated in Fig. 1, which shows that atoms* near the groove are on the convex side of a cylindrically curved surface and therefore have a higher chemical potential than atoms on the flat interface beyond. Thus atoms tend to leave the curved groove shoulders and force the groove to deepen. When both phases are solid additional chemical potential gradient may arise from whatever strain fields are necessary to insure continuity of material across the interface. In the discussion that follows, this complication due to strain will not occur since we always assume one phase to be a fluid and the other a solid. The transport processes that may operate to permit growth of the groove are surface (interface)* diffusion, volume diffusion in the phase on either side of the interface, and evaporation-condensation when one phase is gaseous. (Also a melting-freezing process is responsible for groove enlargement on a solid-melt interface as discussed in Sec. VI). The theory of grain boundary grooving has been discussed in detail for the processes of surface diffusion and of one type of evaporation-condensation.3&apos;4 The theoretical predictions for the mechanism of surface diffusion have been confirmed in an experimental study of grooving kinetics in copper.5 In this paper, the theory of groove development by volume diffusion will be presented. Volume diffusion can occur either by atomic migration in the solid S or by diffusion of S atoms in the fluid phase, or solvent, adjacent to the groove surface. As discussed in Sec. V, volume diffusion in the solid should be less important than surface diffusion in forming a groove. Apart from verifying this statement then, the main interest of the present analysis centers on diffusion in an external solvent. The solvent may be either a liquid (e.g., water, or-

Jan 1, 1961

By M. R. Geer, H. F. Yancey

John Grifien (Pittsburgh)—I wish to congratulate the authors on this paper, which, I am sure, will promote a clearer conception of the various criteria which have been advanced as measures of coal-cleaning performance. Their clear descriptions of the more important criteria that have been advanced and explanations of their meaning and application have been sorely needed and, if widely promulgated and followed, should remove much of the confusion which, heretofore, has been evident. I suggest that much of this confusion has also arisen from lack of common agreement as to the exact meaning of the words and terms used in our discussions of this subject. I am sure I have been equally guilty in this respect and suggest it is time that we select a committee to study vocabulary as well as criteria, so that common understanding and agreement can be had. Thus, I suggest, that the title of this paper should have used the words "coal-cleaning performance" and not "coal-washery performance." My understanding is that "washery" refers to a complete plant and it has additional functions to that of coal-cleaning, such as control of size, moisture, etc. Also, I do not agree with the authors that the formulas given in their paper, the Fraser and Yancey formulas, are necessarily the only ones that can be called efficiency formulas. When a dictionary is consulted one realizes that the words "efficiency" and "efficient" have not precise meanings but only very general ones and must be further defined for each particular use, if they are to have any precise meaning. To my mind, the first of the Fraser and Yancey formulas given should never be used as expressing "efficiency" alone but only by the fuller phrase which truly expresses the kind of efficiency it states, which is "coal recovery efficiency." For complete clarity we then need only define "coal" for the particular circumstances under which the formula is used. This would recognize the fact, fully explained by the authors, that the limits in coal qualities imposed by the market have imposed a factor—and a variable factor—which has influenced any answer obtained by this formula. Such a procedure makes it entirely practical and logical to use the values of any of the qualities of coal—ash, sulphur, phosphorus, specific gravity—or a combination of them to define what is coal for that particular use of the formula. Further, it is clear to me that the criteria to be used must be dictated by the objective in view and the use for which the data are obtained. The authors clearly point out that the Fraser and Yancey coal-recovery efficiency formula is influenced by three factors, "one is determined by the coal, one is dictated by market or use conditions, and one is an inherent characteristic of the cleaning unit itself." The formula which uses the inverse of misplaced material as a measure of efficiency is also a coal recovery efficiency formula and is influenced in like manner by these three factors. As the authors point out, both are useless in obtaining data which enables a direct and valid comparison of the performance of a given cleaning unit on two coals where the other two factors vary or of the performance of two or more cleaning units of different types unless the values of the other two factors are identical. Much of the confusion in our discussions of and literature on coal cleaning has come from failure to understand clearly this fact. Many such comparisons have been made from which erroneous conclusions have been drawn. I am, therefore, particularly glad to see the authors&apos; emphasis on 4?-1e criteria which indicate the "sharpness of separation. I feel that these are the essential criteria

Jan 1, 1952

By Albert S. Trube

The purpose of this paper is to clarify the definition of compressibility and to present a uniform basis upon which instantaneous compressibilities of liquids and gases can be compared. The equations gaverning the instantaneous compressibilities of imperfect gases are derived and the concept of pseudo-reduced compressibility is introduced. Part of the data presented by Brown, Katz et a1 on compressibility factors for natural gases has been rearranged. A graph of pseudo-reduced compressibility vs pseudo-reduced pressure for various pseudo-reduced temperatures is presented. The need for additional work in relating the compressibilities of liquids and gases is discussed. This information should be of value to reservoir engineers in making non-steady state performance calculations in gas reservoirs. It should be of further use irz pointing the direction for additional research in the nature of liquid and gas compressibilities. INTRODUCTION With the increasing use of steady and non-steady state well and reservoir data, there is a corresponding increase in the importance of the various factors entering into such calculations. Increasing emphasis is being placed on the necessity for obtaining reasonably accurate estimates of the physical properties of the reservoir fluids well in advance of the more accurate laboratory data. One such factor is the isothermal coefficient of expansion of the media which are transmitting and attenuating the non-steady state pressure waves. The average isothermal coefficient of expansion, or "compressibility" is a complex function controlled by the physical properties of the formation and the fluids contained therein. The isothermal expansion coefficients for reservoir gases are usually quite variable, in many cases being highly-pressure sensitive. The coefficients for reservoir liquids tend to be pressure sensitive, but not nearly so much as reservoir gases. The coefficients for solids, usually expressed in terms of a "modulus of elasticity" are relatively insensitive to pressure variations within their elastic limits. For this reason, and also because many previous applications have been limited to rel- atively small pressure ranges, there has been a tendency to ignore the variable nature of isothermal expansion coefficients and treat them as constants. Also, the term "compressibility" by which these coefficients are generally designated is commonly confused with a similar term, z, used to define the deviation of an imperfect natural gas from the perfect gas laws. A clear distinction should be made at the outset between the term "compressibility", which is an isothermal coefficient of expansion of a substance, and the term "compressibility factor", z, which refers to the deviation of a gas from the perfect gas laws. Although the scope of this paper is limited to the compressibility of single phase natural gases, it is definitely related to the problem of accurately estimating the compressibilities of single phase hydrocarbon reservoir liquids, which will form the basis of a future presentation. BASIC PRINCIPLES The coefficient of isothermal compressibility of a substance, c, is usually determined from pressure-volume or pressure-length -measurements depending upon whether the substance is single phase gas, liquid, or solid. A convenient method for making such estimates for a finite change in pressure and volume at constant temperature is to use the well known equation V1-V2/V1 (p2 - p1) .....(1) Eq. 1 is negative because the volume of a confined substance decreases as the pressure is increased. In this case V1 > V2 and p2 > p1. This equation is useful in approximating the compressibilities of single phase gases and liquids undergoing small pressure changes. It is evident, however, that this equation is almost identical with the determination of Young&apos;s modulus of elasticity for solids. If the assumption is made that change in length is directly proportional to change in volume, as would very nearly be the case for a steel rod in tension within its elastic limit, then E5=-L1 (p2 - p1)/L1 - L2 .......(2) in which E. is the isothermal expansion coefficient, or Young&apos;s modulus of elasticity, for a solid. And further, for this special case L1 (p2 - p1)/L1 - L2 .......(3)

Jan 1, 1958

By M. E. Wadswort, A. S. Peck

Infrared data disclose that phenacite reacts with oleic acid to form a chemisorbed oleate monolayer on the mineral surface. The absorption band characteristic of the antisymmetrical C = 0 oleate structure appears at a spectral wavelength of 1623 cm-&apos;. Chemisorption also results in the displacement of OH- and silicate from the mineral surface. Infrared absorbance at 1623 cm-1 is correlated with float-ability at different pH for -60-mesh phenacite particles. Phenacite, Be2SiO4, is an ore mineral often associated with quartz, microcline, and chrysoberyl in pegmatites and mica schist. Experimental phenacite beneficiating processes have been confined principally to the flotation of low grade ores occur-ing in several mountain ranges of Utah and Nevada. The U.S. Bureau of Mines &apos; has recently reported on the separation of phenacite from ore deposits of Mount Wheeler, Nev. Infrared techniques were previously applied by the authors to characterize the role of sodium oleate and oleic acid in the flotation of several non-sulfide minerals.23 A parallel study was projected to the flotation of the orthosilicate mineral phenacite with oleic acid. Infrared spectra provided a &apos;molecular fingerprint&apos; of the mineral surface and permitted identification of the reaction mechanisms. Infrared data were complemented with flotation and chemical tests. EXPERIMENTAL PROCEDURES The experimental procedures used in this investigation, for which data are presented in graphical or spectral form, are described in chronological test order. Phenacite crystals from Mt. Antero, Colo., were obtained from Ward&apos;s Co., Monterey, Calif. The crystals were fragmented in a porcelain mortar to -60 mesh and then pretreated with 0.5 molar HC1 for 15 min to remove iron. Then, the mineral particles were washed free of acid with distilled water and dried. Part of the dried material was ground to a fine powder with a Fisher Grinder using an agate mortar and pestle. Next, the powder was placed in a graduated cylinder and mixed with 500 ml of distilled water. The mineral in suspension was decanted from the settled fraction after 5 min, and the colloidal suspension was used for adsorption tests in which a high specific mineral surface area was required for adsorption and infrared analysis. Conditioning was performed in 50 ml glass centrifuge tubes for both infrared and flotation test determinations. Twenty ml of colloidal suspension was used for the infrared adsorption tests. One and a half g of deslimed -60-mesh HCl pretreated phenacite was added to 20 ml of distilled water for flotation conditioning tests. Six mg of oleic acid was used for each infrared adsorption test, and 0.2 mg of oleic acid was used for the flotation conditioning tests. HCl and NaOH were used for pH modification. The pulps were stirred at a speed of about 1500 rpm for 3 min. Then, the solids were separated out and either dried at 50°C or immediately transferred to a float cell. The pH of the residual liquid was determined with a Beckman Zeromatic pH meter. The flotation procedure is described as follows: 1) A 325-1111 glass micro-float cell was filled with distilled water. 2) Conditioned mineral was placed in the cell and the pulp was aerated through the glass fritted bottom at 7 psi pressure. 3) Floated mineral was skimmed from the cell with a glass slide until flotation ceased. (The end point was reached in all tests in less than 2 min.) 4) The flotation products were dried, weighed, and the recovery was calculated. Samples for infrared analysis were made into mulls using either oleic acid or series 13-21 Halocarbon oil obtained from Halocarbon Products Corp., Hacken-sack, N.J. The mulls were placed between two sodium chloride windows. A characteristic phenacite infrared band at 1868 cm-&apos; was used as an internal standard for quantitative determinations in which absorbance at this position was quantified by the baseline-tangent method4 and related to absorbance at 1623 cm-&apos;. The ratio obtained by dividing absorbance at 1623 cm-&apos; by absorbance at 1868 cm-&apos; was defined as an absorbance quotient for use in graphical presentation.

Jan 1, 1968

By D. E. Peacock, B. Harris

Thermal-expansion and electrical-resistivity measurements have been carried out below 400°K on niobium and two niobium alloys containing tungsten. For anonaly in the expansion us temperature curve below aboct 300°K. This reslt is consistent with an earlier obsercation that Young&apos;s modilus joy Cb-132 lms an unusiully strong -tenperati-e dependence below roon temperature. No anonulous behavior was observed in the resistiuity vs temperature curve for either alloy or in the magnetic susceptibility us tenperature curve for Cb-132. These results are disclssed in terms of various electronic and bonding theories of the transition metals and their alloys, but it has not been possible to reach any definite conclusion as to the cause of the results obtained. ALONG with investigations of the low-temperature mechanical behavior of niobium and its alloys, some of the physical properties of these materials have also been studied. In the work described here, thermal expansion and electrical resistivity were measured below 400°K for niobium and two of its alloys, Cb-132 and D-43. In addition, the magnetic susceptibility of Cb-132 was measured over the same temperature range. The alloys are representative of two distinct tion-hardened alloy such as might find service as a turbine blade, while D-43 (Nb, 5.3 at. pct W, 1.1 at. pct Zr. and 0.08 at. pct C) is a solution- and dispersion-hardened alloy, of considerably lower strength than Cb-132. which is a potential vane material. Their significant common feature is the presence of tungsten as the principal source of solution hardening. EXPERIMENTAL TECHNIQUES Thermal-expansion measurements were made with a quartz-tube dilatometer using baths of liquid nitrogen, liquid bromoethane, and a liquid polyalkylene glycol derivative to cover the range 80 to 400 K. Results are recorded as expansion relative to that of quartz, and the probable maximum error in any value of relative expansion. All, is about il pct. The metallurgical condition of the specimens, which were in rod form. will be discussed in the section dealing with experimental results. Resistivity was measured on recrystallized wires 0.030 in. in diam using a potentiometric technique. The same baths were used for temperatures between 80&apos; and 400°K as were used for the expansion measurements. and in addition the resistance of D-43 and Cb-132 was measured at 4.2 K by immersion in liquid helium. The maximum expected experimental error in any resistivity value is about 3 pct and that in any change in resistivity about i0.05 pct. Temperatures in both sets of experiments were measured with liquid-in-glass thermometers with an accuracy usually better than 1°K. EXPERIMENTAL RESULTS 1) Thermal Expansion. In Fig. 1 is shown the relative linear expansion of recrystallized niobium. D-43. and Cb-132. The curve for Cb-132 is composed of results from several runs on three different samples. Two of the samples were annealed for an hour at 1500°C to cause recrystallization and the third was annealed for an hour at 2000&apos;C to remove some mi-crosegregation of tungsten which was known to be present in the wrought alloy received from the manu-

Jan 1, 1967

By B. B. Rath, Hsun Hu

In wedge-shaped bicrystals of zone-refined aluminum it is observed that (111) pure tilt boundaries migrate under the driving force of their own inter-facial free energy. The boundary velocity is a power function of the driving force. The driving force exponent decreases with decreasing angle of misorien-tation. For example, at 64O°C, the exponent decreased from 4.0 for a 40 deg to 3.2 for a 16 deg tilt boundary. An evaluation of the driving force acting on the boundaries during their motion indicates that for low driv-forces, up to about 2 x l03 ergs per cu cm, the velocity is relatively independent of misorientation, whereas at higher driving forces a 40 deg tilt boundary exhibits the highest velocity. The measured activation energy for boundary migration approaches that for bulk self-diffusion at low driving forces, decreasing from 33 to 27 kcal per mole as the driving force is increased from 1 x l0 to 5 x l03 ergs per cu cm. These results are compared with current theories of grain-boundary migration. In previous experimental studies of grain boundary migration the driving force has been limited to a difference in stored energy across the boundary. This stored energy has been introduced into the crystal either by prior deformation1-3 or by grown-in lineage structure. A part of the energy stored in the deformed crystal is released by recovery either prior to or concurrently with grain boundary migration, thus introducing an uncertainty as to the magnitude of the driving force responsible for grain boundary migration. The grown-in lineage structure, though thermally stable during annealing, neither provides conditions under which different levels of energy may be stored in the imperfect crystal nor provides a control of orientation difference across the migrating boundary of a growing grain. Furthermore, because of variation in the lineage structure, it is difficult to determine accurately the energy stored in the imperfect crystal. Several investigations of grain boundary migration during normal grain growth have also suffered from difficulties in estimating the driving force because of uncertainties in the principal radii of curvature.~ In the present investigation the velocity of pure tilt boundaries in zone-refined aluminum bicrystals of selected orientation (40, 30, and 16 deg around the [Ill] tilt axis) has been measured in the absence of a dislocation density difference across the moving boundary, thus eliminating the previous experimental difficulties. The driving force for boundary migration is derived from a gradient of the total interfacial free energy of the migrating boundary in wedge-shaped bicrystals. A similar method was attempted by Bron and Machlin in a study of grain boundary migration in silver. However, they found that one of the crystals was deformed and consequently the motion of the boundary was partly due to a difference of stored energy across the boundary. The observed behavior of boundary velocities as affected by the driving force is examined in the light of the predictions of the current theories of grain boundary migration.7"10 The effect of boundary misorientation on velocity is compared with the theory of " which is based on a dislocation core model for high-angle boundaries. EXPERIMENTAL METHOD Seed-oriented bicrystals of zone-refined aluminum, 2.5 cm wide, 0.5 cm thick, and 12 cm long, containing tilt boundaries with a common (111) axis, were grown from the melt in the direction of this axis. Spectro-graphic analysis, reported earlier,&apos;&apos; indicated the purity of the crystals to be 99.999+pct. Three such bicrystals containing 16, 30, and 40 deg tilt boundaries were used. Wedge-shaped specimens were prepared from these bicrystals by spark cutting followed by electrolytic polishing. The angle of the wedge was usually 40 deg and the specimens were usually 0.25 cm thick. The intercrystalline boundary was located within 0.2 to 0.5 cm from the tip of the wedge. Fig. 1 shows a section of an oriented bicrystal containing an outline of a wedge-shaped specimen. The crystallographic directions shown in Fig. 1 represent the orientation of one of the crystals (the larger section of the bicrys-tal); the orientation of the other crystal differs only by rotation around the common [lil] axis. The parallel faces of the wedge always corresponded to the common (171) planes in both crystals, whereas the orientation of the side faces varied, depending on the misorientation angle. The bicrystal orientations were determined

Jan 1, 1970

By J. C. Martin

A theoretical analysis of steam stimulation is presented for single sands. The analysis includes the effect of steam production and most of the effects of heat conduction. The results show the effects of a number of important variables on the performance of an idealized well under steam stimulation. Calculated responses are presented which indicate the effects of steam production, amount of steam injected, water production, formation thickness and formation damage. Results indicate that steam production can cause large reductions in the heat contained in the formation. This effect can be eliminated by drawdown control. Water production reduces the amount of oil produced during stimulation. The optimum amount of steam to inject depends on economic factors as well as the well response. In many cases, the increased temperature resulting from stimulation reduces oil viscosity near the well sufficiently to overcome the effects of formation damage even if the damage is not removed during steam injection. Calculated responses for thin sands are more favorable than anticipated. INTRODUCTION Little has been published on the theory of steam stimulation&apos;-&apos; despite the interest it has created and the wide variation in well responses. The results of the present analysis provide an insight into steam stimulation, and the methods employed provide a foundation for future work. Analyses presented in Refs. 1 and 2 are very limited and apply to gravity drainage conditions. Ref. 3 contains an analysis similar to the one presented here. The idealized models used and the assumptions made in Ref. 3 are different from those used in this paper. The analysis assumes that after steam injection has heated a small portion of the volume within the radius of drainage of a single uniform sand, a shut-in soaking period is allowed before returning the well to production. The effects of gravity, capillarity, transient pressure and water-sensitive sands are neglected. The injection and soaking times are assumed short compared to the stimulated production time. The initial temperature is assumed uniform; thus, the results apply primarily to first-cycle stimulation. The effects of gas production other than steam are neglected, and the water-oil ratio during production is assumed constant. Steam stimulation involves the simultaneous variation of the temperature, pressure and saturations. General mathematical equations for these variations are complicated and can be very difficult to solve. Simplified equations based on idealized models are used to reduce the mathematics sufficiently to allow approximate solutions to be obtained. DISCUSSION INJECTION An idealized model for which heat conduction is neglected is used to represent the behavior of a well during steam injection. The mathematics for this model is presented in Appendix A which also contains an approximate solution for the behavior of the no-conduction model. SOAKING During soaking, the well is shut in. Only the temperature, pressure and saturation distributions at the end of soaking are needed in the analysis. During soaking the heat is considered to be conducted in a uniform medium from an initially uniformly-heated circular cylinder confined to the producing interval. At the end of soaking the saturations and pressures are assumed to correspond to the cold zone. Analysis of the heat flow during soaking is included in the next section. The radius of the heated cylinder is calculated from the following heat balance (for constant quality steam injection). At the end of the soaking period it is assumed that little or no free gas is present near the well, and that the soaking period has been sufficiently long that the steam zone has had time to expand and the steam has condensed. The condition where there is no soaking is considered in the next section. PRODUCTION In this section, an approximate method is presented for solving the equations of heat and fluid flow associated with the production of oil and water during steam stimulation. Where initial pressure drawdowns are sufficient to cause steam to flow into the wellbore, steam production is assumed to occur within a short initial adjustment period (Appendix B). The production practices followed soon after the well is returned to production can have a large influence on the amount of oil produced during the stimulation cycle. Under most conditions, there is a period of time in which some or all of the water in the heated zone is converted into steam and produced. This flashing of hot water is caused by the pressure in the heated zone dropping below

By J. W. Marx, R. H. Langenheim

The injection of heat-bearing fluids may offer a wider application to secondary and tertiary recovery from conventional oil reservoirs than inderground combustion since the process is more easily controlled and the reservoir requirements, in general, are less critical. In all cares, conduction heat losses to overburden and un-derburden will impose an economic limit upon the size of the area which can be swept out from any one injection point for any given set of conditions. The present report describes a method for c.rtin~ating thermal invasion rates, cumulative heated area, and theoretical economic limits for sustained hot-fluid injection at a constant rate into an idealized reservoir. Full allow-ance is made for non-productive reservoir heat losses. Required solutions can be obtained directly from the numerical tables given here, as soon as all operating conditions are specified, withou extensive mathematica1 ~nanipulation on the part of the user. The injection of a heat-bearing fluid into a reservoir is frequently proposed as a means of secondary or tertiary oil production. This form of thermal recovery may offer a wider application to conventional oil reservoirs than underground combustion since the process is more easily controlled and the reservoir requirements, in general, are less critical. Because of its large gross heat capacity, steam appears to be the most efficient heat injection medium, although mixtures of steam and other gases, hot water, hot oil and hot non-condensable gases also have been employed. Whatever the nature of the injection medium, economic evaluation of such thermal recovery processes will depend upon thermal invasion rates, or rates of productive reservoir heating, at any given time after the start of heat injection. In all cases, conduction heat losses to overburden and underburden will impose an economic limit upon the size of the area which can be swept out from any one injection point, for a given set of reservoir conditions, at any given heat injection rate. The present report describes a method for estimating thermal invasion rates, cumulative heated area and theoretical economic limits for sustained heat injection at a constant rate into an idealized reservoir, making full allowance for non-productive reservoir heat losses. Required solutions can be obtained directly from the nu- merical tables given here, as soon as all operating conditions are specified, without extensive mathematical manipulation. It should be noted that ordinary water flooding may be used to move an annular-heated zone outward, even though the economic areal limit for straight heat injection has been reached. Thus well spacing need not be confined to the areal limits as determined here. However, the present calculations do not cover such supplementary heat scavenging programs. The operating element consists of a radial flow system, concentric about the point of injection, with the idealized step function temperature profile illustrated by the dashed line in Fig. 1. The solid curve shown in Fig. 1 gives a qualitative picture of the true radial temperature distribution. It will be shown later that the expressions given below, which are derived from the simplifying step idealization (dashed line in Fig. I), also can be interpreted in terms of the more realistic temperature distribution shown in Fig. 2. For this case, the thermal invasion radius is simply redefined as the distance from injection well to the mid-point of the temperature distribution.