Search Documents

By Leland L. Mink, Roy E. Williams, Lewis M. McNay, Alfred T. Wallace

In 1968 the metal-mining industry of northern Idaho installed tailings ponds as a means of handling and treating mine wastes, metallurgical process wastes and domestic wastes. The installation of the ponds constituted a crash program to minimize pollution of the Coeur d'Alene River. The ponds provided an opportunity to investigate which combinations of wastes tailings ponds can be treated. The research data pertaining to tailings pond discharge lends insight into compliance with the recent EPA guideline for effluent quality.

Jan 7, 1973

By M. L. Jeremic

Abstract-Observations of ground behavior in mining the Choka Dulkan ore body have been described. Ore recovery was limited to approximately 60% because unconsolidated hydraulic fill did not permit either the widening of stope spans nor secondary ore recovery from pillars. Stope undercutting induced high stress concentrations. However, as extraction progressed, ground conditions improved but before extraction was complete they deteriorated again. After completion of a stope, the remnant rib and crown pillars continued to be deformed by other stresses for some time.

Jan 4, 1978

By J. R. Low, A. M. Turkalo

ManY qualitative observations have been made concerning the relationship of the carbide dispersion and the mechanical properties of steel. With reference to the strength and hardness we know that for a given composition a steel having a fine pearlitic structure is stronger than that having a coarse pearlite. The martensitic structure tempered at a relatively low temperature and having finely dispersed carbides gives greater strength than that given by the martensitic structure tempered at a relatively higher temperature and having a coarser

Jan 1, 1959

By M. M. Austin

PURE metallic tantalum, from a practical standpoint, is one of the newer developments in the metal field. Although it was used as filament in incandescent lamps in 1906, only within the last five years has the metal been available in sufficient quantity and suitable form for more general application. In addition to corrosion-resisting properties and high melting point, the metal has an ability to withstand unusual amounts of cold working. It affords some striking demonstrations of block movement during plastic flow.

Jan 1, 1930

By A. Gemperle

Electron diffraction and transmission electron microscopy on foils at room temperature were used to investigate the ordering of Fe-Si alloys containing 10 to 23 at. pct Si. A certain degree of DO3 order was found in all alloys. With the exception of the lowest silicon concentration for which the antiphase domains could not be clearly resolved, the alloys have a domain structure of two-domain type with boundaries having 1/4a01<111> displacement vectors for less than 12.3 at. pct Si and with boundaries having 1/2 a0<100> displacement vectors for more than 12.3 at. pct Si. The alloys with 12.3 at. pct Si have a domain structure consisting of fine domains with 1/2a&apos;o<100> boundaries within much larger domains with 1/4a&apos;o<111> boundaries. The development of these structures can be explained by transition of the alloy from the disordered state into the B2-type order and then into the D03-type order by the mechanism proposed previously for the FeSAl alloys. The existence of the B2 structure in the lower part of the investigated concentration range reported in some articles can be explained by fine domains with 1/2a&apos;o<100> boundaries formed by several disordered planes within large domains with 1/4a&apos;o<111> boundaries. The ordered structure predicted by the theory —with practically no domain boundaries —is found in the alloys having 12.3 at. pct Si where it develops in the B2 structure region. ORDERING in Fe-Si alloys was first studied by phragmenl who found that beginning with 13 at. pct Si the DO3 (Fe3Si) superlattice reflections appear in the diffraction patterns. The equilibrium diagrams constructed later by corson2 and Haughton3 from various measurements proposed the existence of a homogeneous solid solution (a phase) in the range from 0 to 25 at. pct Si. Jette and Greiner4 and Farquhar et al.5 measured the relation between lattice parameter and composition and they considered the break in the curve at 9 to 10 at. pct Si to be caused by the ordered solution a"(Fe3Si). Glaser and Ivanick6 Determined critical ordering temperatures of the alloys containing from 10.9 to 27.9 at. pct Si from the measurement of the electric resistivity of quenched samples. In all cases the critical temperature was lower than the melting point and it was highest for 25 at. pct Si. Lihl and Ebel7 measured the lattice parameter curves at various temperatures up to 1000°C. The region between two breaks on these curves, corresponding to 10 to 12.5 at. pct Si at room temperature, was considered by them to be two-phase (a + a"). They concluded by extrapolation of the measured values that a" in the alloy having 25 at. pct Si is stable up to the melting point. Davies8 studied superlattice reflections in the X-ray diffraction patterns of an alloy containing 8.7 at. pct Si. He found the B2 structure and short-range order in the slowly cooled samples and the DO3 A. GEMPERLE is Research Scientist, Institute of PhysicS, Czechoslovak Academy of Sciences, Prague, Czechoslovakia. __Manuscript submitted January 2, 1968. IMD structure in the quenched and annealed samples. This investigation first reports the presence of the B2 structure, phase a&apos;: in Fe-Si alloys. Meinhardt and krisement9,10 also found its existence in Fe-Si alloys by neutron diffraction. No order was detected by them in the alloy containing 8 at. pct Si. The onset of B2 order was observed at a composition of 9.2 at. pct Si. They found almost perfect B2-type order with partial DO3-type order at room temperature in the 10 to 12.5 at. pct Si range and almost perfect DO3-type order in the 12.5 to 25 at. pct Si range. They established the critical temperatures Tc of both the structures through measurement at higher temperatures. The critical temperature for the B2 structure was found to be always higher than the critical temperature for the DO3 structure of the same alloy. They extrapolated the curves of the critical temperatures and concluded that the alloys with more than 17 at. pct Si have the B2 structure up to the melting point and the alloys with more than 23 at. pct Si have the DO3 structure up to the melting point. The results of Meinhardt and krisement9,10 were confirmed by Dokken&apos;s measurement" of the temperature dependence of the electrical resistivity in the 10.8 to 15 at. pct Si range. On the other hand chessin12 detected by X-ray diffraction a considerably lower degree of order in the 12.7 at. pct Si alloy. ANTIPHASE DOMAIN STRUCTURE IN ALLOYS WITH B2- AND DO$-TYPE ORDER The ordered structures B2 and DO3 can be described in terms of a subdivision of the bcc lattice into four fcc sublattices with a parameter double that of the bcc lattice. Following Marcinkowski13 we will label them I. 11, 111, IV. The B2 structure in the AB alloy is formed by placing A atoms on sublattices I and II and B atoms on sublattices III and IV. In a non-stoichiometric perfectly ordered alloy having concentration (A) > (B), A atoms occupy sublattices I and 11. and A and B atoms distributed at random occupy sublattices III and IV. The DO3 structure in an A3B alloy is formed by placing A atoms on sublattices I, 11, and 111, and B atoms on sublattice IV. In a nonstoichiometric, perfectly ordered alloy having concentration (A)/3 > (B), A atoms occupy sublattices I, 11, and 111, and A and B atoms distributed at random occupy sublattice IV. As further shown in Ref. 13, two types of domains are possible in the B2 superlattice and the antiphase domain structure has associated with it boundaries with displacement vectors 1/4 a&apos;o<ll1> only. Four types of domains are possible in the DO3 superlattice and the antiphase domain structure has associated with it boundaries with displacement vectors 1/4a&apos;o<111> and 1/2a&apos;o<l00>. Bethe14 suggested on the basis of theoretical considerations that in a structure with two sublattices at low temperatures only one domain should be present in the whole crystal at equilibrium. Similarly Bragg15 concluded that at low temperatures the domain struc-

Jan 1, 1969

By S. S. Cole, H. Sigurdson

The melting points of mixtures of titanium dioxide and other titanates have been reported to a limited extent as binary systems and some results have been reported in conjunction with silicon dioxide. The limited data indicated that a low melting zone might exist in a ternary or quarternary system of CaO-MgO-TiO2-Al2O3, since the eutectics reported in the binary systems of MgO.TiO2-TiO2, Al2O3TiO2-TiO2 and CaO.TiO2-TiO2 were of about equivalent composition of TiO2. The importance of a low melting region with a high titanium oxide value in the developing of high titanium slags is fully appreciated by metallurgists. A comprehensive study was undertaken to establish basic information on low melting titanate mixtures. Compounds in the System CaO-MgO-TiO2 -Al2O3 Preparation of the titanates of CaO, MgO and Al2O3 confirmed reported data that the following compounds could be formed under oxidizing conditions in solid state reactions: CaO.TiO2, 2MgO.TiO2, MgO.TiO2, MgO-2TiO2 and Al2O3.TiO2. Under oxidizing conditions it was not possible to react MgO with TiO2 in mol ratios higher than 1:2 without having unreacted TiO2 in the product, nor was it possible to form calcium or aluminum titanates with a higher mol ratio of CaO and A12O3 to TiO2 than 1:1. By fusing mixtures of Al2O3-TiO2 and MgO-2TiO2 together a series of solid solution products were obtained, which showed an MgO-2TiO2 X ray diffraction pattern shifted to smaller inter-planar spacings. Range of Investigation Data in the literature1 indicated that these titanates had melting points which ranged from 1645 to 1860°C. Since the purpose of the study was primarily to provide useful data for smelting titaniferous ores, the work was restricted to the zone of the system which would include the crystalline phases present in the slags. The explored limits were bounded by CaO.-TiO2-MgO.TiO2-TiO2in the base plane. The system was extended to a fourth component Al2O3-TiO2 since many titaniferous ores contain appreciable amounts of A12O3. In the quarternary system CaO-TiO2-MgOTiO2-TiO2-Al2O3TiO2, a tetrahedron was used to represent graphically the components, with one component at each point of the tetrahedron. In the base plane, mixtures were pre- pared to represent fairly uniform changes of composition over the desired range expressed on a mol percent basis. Al2O3-TiO2 was brought into the system in increments of 10 mol pct up to 40 pct. Mixtures were then made up for each of the Al2O3-TiO2 planes. The ranges employed were then represented by five planes cutting the tetrahedron at 10 mol pct Al2O3-TiO2 intervals. Choice of Equipment Although it was necessary for all smelting work to be done in a strongly reducing atmosphere, it was decided that all melting point determinations should be done under oxidizing conditions. This decision was prompted by the fact that TiO2 reduced to lower oxides in a reducing atmosphere with consequent changes in melting points. The use of a micropyrometer and a platinum strip furnace for work of this type is adequately described in the literature.2 With modifications, equipment was selected which provided a fairly rapid method of determining melting points of refractory oxides. The final assembly, which was used, consisted of a platinum strip furnace and a Leeds and Northrup disappearing filament optical pyrometer attached to a special telescope which magnified the sample about 20 diam. The furnace assembly consisted of a platinum strip 0.005 X 0.3 X 2.2 in. mounted on brass posts on a refractory base. This was enclosed in a black steel shell which had a 2-in. opening at the top for purposes of sighting the telescope. The platinum strip was in series with two 0.04 ohm nichrome resistors connected to the secondary of a 2 kw transformer which supplied 17 volts. The current to the primary was varied

Jan 1, 1950

By David A. Thomas, David N. French

The lrnrdness of WC crystals has been measured with the Knoop indenter at loads of 100 and 500 g on the (0001) and (1070) planes. The hardness as tneasitred on the basal plane is 2400 kg per sq mm and shows little anisotropy. The hardness on the prism plane, however, shows a marked orientation dependence, with a low value of 1000 kg -per sq mm when the long axis of the Knoop indenter is oriented parallel to the c axis and a high value of 2400 kg per sq mm when the indenter is perpendicular to the c axis. Slip lines (Ire observed surrounding the microhardness indentations and they show slip on (1010) planes, probably in [0001] and (1120) directions. This slip behavior can be explained by the crystal structure of TVC, which is simple hexagonal with a c/a ralio of 0.976. The hardness anisotropy call be explained by [0001]{1010} and (1130) {10l0] slii) and the resolved shear-stress analysis of Daniels and Dunn. HARDNESS anisotropy is a well-known phenomenon and has been reported for many metals, with both cubic and hexagonal structure.1-6 For hexagonal tungsten carbide, WC, a wide range of hardness values is reported in the literature. For example, Schwarzkopf and Kieffer7 give a hardness of 2400 kg per sq mm and report a value of 2500 kg per sq mm by Hinnüber. Foster and coworkerss give the average Knoop microhardness as 1307 kg per sq mm with a maximum value of 2105 kg per sq mm. Although these values and the structure of WC suggest the likelihood of hardness anisotropy, no such measurements have been made. We first detected a large apparent hardness anisotropy in WC crystals about 75 p large, in over-sintered cemented tungsten carbide. Prominent slip lines also occurred around many indentations. This report presents further observations and interpretations of hardness anisotropy and slip in WC crystals obtained from Kennametal, Inc. Both Knoop and diamond pyramid indenters were used on a Wilson microhardness tester with loads of 100 and 500 g. EXPERIMENTAL RESULTS The carbide crystals tended to be triangular plates parallel to the (0001) basal plane of the hexagonal structure. The side faces were parallel to the ( 1010) prism planes. Specimens were mounted approximately parallel to these two types of faces and metallographically polished. Laue back-reflection X-ray patterns were used to orient the specimens, which werethen ground to within ±1 deg of the (0001) and (1010) planes. The Knoop hardness values measured on the basal plane are plotted in Fig. 1. There is only a small anisotropy, with a hardness range of 2240 to 2510 kg per sq mm. The additional points at angles from 52.5 to 67.5 deg confirm the sharp minimum hardness at 60-deg intervals, consistent with the sixfold hexagonal symmetry. The average hardness of all values obtained on the basal plane is 2400 kg per sq mm. While the basal plane shows only slight anisotropy, the (1010) plane exhibits marked hardness anisotropy, from 1000 to 2400 kg per sq mm. Fig. 2 shows the hardness as a function of the angle between the long axis of the indenter and the hexagonal c axis, the [0001] direction. The minimum and maximum occur when the indenter is oriented parallel and perpendicular to the [0001] direction, respectively. The anisotropy of the prism plane is contrary to that reported for hexagonal zinc and hard- However, the basal-plane anisotropy is similar to these two metals.1&apos;2 To check the accuracy and reproducibility of the measurements, a series of fifteen impressions was made at 100-g load in the same orientation in the same area of the specimen surface. The average for all was 2040 kg per sq mm, with a range of 1950 to 2130 kg per sq mm, giving an accuracy of about ± 5 pct. Thus the slight anisotropy on the basal plane is almost within experimental error. Fig. 3 shows slip lines around the Knoop indentations on the basal plane. The slip traces are in directions of the type (1130). The presence of slip steps on the basal plane indicates that the slip direction lies out of the (0001) plane. Because WC has a c/a ratio of 0.976,&apos; the shortest slip vector is [0001], which suggests slip systems of the type [0001] (1010). Fig. 4 shows slip lines around the Knoop intentations on the (1010) plane. These slip lines are inconsistent with [0001] slip but can be

Jan 1, 1965

By C. G. Dunn, M. Sharp

IN twinned crystals of the face-centered cubic metals the lattice of one twin is a mirror image of the other in a common twin boundary. When several twins appear within large grain in a sheet specimen, the twin one boundaries form a set of lines at the surface of the specimen which coincide with (111) planes of the large grain. Furthermore, for twins of the same orientation, these lines are parallel. Generally, the presence of identically oriented regions with straight parallel boundaries coinciding with a (111) plane of the surrounding crystal is strong evidence for identifying the island regions as twins of the parent crystal. However, Fig. 1, which shows the macrostructure of a large grain of copper with island regions that satisfy these conditions. is not an illustration of (111) twins. Since the reverse side of the specimen has much the same appearance, it was thought at first that these regions, which appear dark in the macrograph, actually were twins. According to X-ray data, however, these regions are second-order twins of the large crystal. With regard to their formation, these second-order twins formed by secondary recrystallization in a cube texture matrix. Growth occurred in the direction of the arrow (see Fig. 1) as the specimen moved slowly into a gradient temperature furnace as described previously.&apos; Nucleation of the second-order twins occurred, therefore, on the ends facing opposite the arrow. If the origin of the second-order twins were due to repeated twinning, some first-order twin structure should be visible on these ends. This proved to be the case, as very small twins were readily found with the aid of a microscope, and probably could have been seen, in some instances, under ideal lighting conditions without aid of a microscope. Fig. 2 shows a cross-section view taken perpendicular to both the surface and the (111) trace of the parent crystal (visible as a straight boundary in Fig. 1) at the beginning point of growth of a second-order twin and where one first-order twin was relatively thick. In the micrograph, A is the large parent grain; B is the first-order twin of A; and C, which is a first-order twin of B; is a second-order twin of A. Between A and B and between B and C the major straight portions are traces of common (111) twin boundaries. The straight portion of boundary between A and C, however, is not a common crystallographic plane to the two lattices; it is a (111) plane of A and a (115) plane of C. Without considering the mechanism of twinning itself, the origin of the second-order twins may be accounted for in terms of repeated twinning and special growth characteristics. After each nucleation, a selective growth process can be thought of as favoring growth of the first-order twin in local spots only and favoring growth of the second-order twin to an extent comparable with that of the parent grain over relatively large areas in a way similar to that described for twinning in aluminum.&apos; It has already been pointed out that the boundary between the large grain (A) and the second-order twin (C), which is responsible for the straight boundary portions in Fig. 1, involves a (111) plane of A and a (115) plane of C. The same combination of planes is not only possible in first-order twins, but actually appears quite frequently.3 Their prevalence in first-order twins and their presence here in second-order twins, together with the necessary occurrence of a large number of common lattice sites at the boundary, is an indication that this combination produces an "energy cusp"&apos; boundary. (Energy cusp boundaries have been described by Shockley and Read.") The configuration of atoms near a {Ill), (115) boundary in first-order twins is of course different from the configuration near the same type of boundary in second-order twins. References 1 M. Sharp and C. G. Dunn: Secondary Recrystallization Texture in Copper. Journal of Metals (January 1952) Trans. AIME, p. 42. 2W. G. Burgers and W. May: Stimulated Crystals and Twinning in Recrystallized Aluminum. Recueil des travaux chimiques des Pays-Bas (1945) 64, p. 5. aD. Whitwham, M. Mouflard, and P. Lacombe: Discussion of W. C. Ellis and R. G. Treuting, "Atomic Relationships in the Cubic Twinned State." Trans. AIME (1951) 191, p. 1070; Journal of Metals (October 1951). 4 W. Shockley and W. T. Read: Dislocation Models of Crystal Grain Boundaries. Physical Review (1950) 78, p. 275.

Jan 1, 1953

By Herman Dykstra, R. L. Parsons

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.&apos; THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy&apos;s Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. &apos;Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If&apos; .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951

By W. O. Philbrook, K. Natesan

The oxidation kinetics of spherical pellets of zinc sulfide made from Santander concentrates were studied using a thermogravimetric technique. The experiments covered a temperature range-. of 740" to 102O°C, 0-N mixtures varying from 20 to LOO pct O2, and pellet diameters between 0.4 and 1.6 cm. Mathematical models were formulated to Predict the reaction rate on the assumption that a single transport or interface reaction step was rate -controlling. Analysis of the data indicated that the process of oxidation was predominantly controlled by transport through the zinc oxide reaction-Product layer. ROASTING processes, which are reactions between solids and gases, are very important because they are employed in the production of a number of basic metals. These processes are highly complicated, and one needs to consider the transport phenomena of heat and mass between the solids and gases in addition to the kinetics of various chemical reactions involved. Because of such complications there is a lack of knowledge concerning the rate-limiting factors, which may strongly depend on temperature, particle size, gas composition, and solid structure. The oxidation of zinc sulfide, which is of commercial importance in zinc production, falls into this class of reactions. The major goal of this work was to elucidate the roles played by different process variables, such as reaction temperature, gas composition, pellet size, and pellet porosity, on the kinetics of oxidation of single pellets of zinc sulfide. Roasting of zinc sulfide single particles has been a subject of both experimental and theoretical investigations.&apos;-&apos; The reaction is exothermic and may be considered to be irreversible. Such a reaction has been found to proceed in a topochemical manner. In other words, as the reaction proceeds, a progressively thicker outer shell of zinc oxide is formed, while the inner core of unreacted sulfide decreases. It has been found experimentally, both in the present work and in the previous investigations,1-9 that the particle retains its original dimensions and the process requires transport of gaseous oxygen across the porous product layer for continued reaction. The reaction may be represented by ZnS(s) + 3/2 O2(g) = ZnO(s) + SO2(g) [1] The solid product considered here is only zinc oxide, since the diffraction patterns of zinc sulfide pellets oxidized partially at &apos;798" and 960°C showed K. NATESAN, Junior Member AIME, formerly St. Joseph Lead Fellow, Department of Metallurgy and Materials Science, Carnegie-Mellon University, Pittsburgh, Pa., is now at Argonne National Laboratory, Argonne, Ill. W. 0. PHILBROOK, Member AIME, is Professor of Metallurgy and Materials Science, Carnegie-Mellon University. This paper is based on a them submitted by K. NATESAN in partial fulfillment of the requirements for the Ph.D. degree in Metallurgy and Materials Science at Carnegie-Mellon University. Manuscript submitted December 2, 1968. EMD lines corresponding to original zinc sulfide and the newly formed zinc oxide. OXIDATION MODEL The generalized model for gaseous oxidation of zinc sulfide is illustrated in Fig. 1. This depicts a partially oxidized sphere of zinc sulfide in a gas stream surrounded by a laminar film of gas. The spherical sample of zinc sulfide of unchanging external radius r0 is suspended in a flowing gas stream of total pressure PT and composition specified by the partial pressures of the individual components. Partial pressures of the gaseous species in the bulk gas phase, at the exterior surface of the pellet, at the ZnS/ZnO interface, and at equilibrium for Reaction [I] are identified by the superscripts b, o, i, and eq, respectively. The overall reaction involves the following ~te~s:&apos;~&apos;~&apos; Step 1. Transfer of reactant gas (oxygen) from the bulk gas stream across the gas boundary layer to the exterior surface of the pellet and the reverse transfer of the product gas (sulfur dioxide). Step 2. Diffusion and bulk flow of oxygen from the pellet surface through the product shell (ZnO) onto the ZnS/ZnO interface and the reverse transfer of sulfur dioxide. Step 3. Chemical reaction at the interface, which results in consumption of oxygen gas and generation of sulfur dioxide gas and heat; at the same time the _________ PARTICLE SURFACE / x^^^^n^X / MOVING INTERFACE core-----/ sJSNxy " /T^T^02 \ V\V$\ ^NWX/ "*/— GAS BOUNDARY x. \\ZnO SHELLV/ / \^ . / &apos; ^ BULK GAS --------N_______ _____,------p£ w \ P« / (f> \ / a \ <i) / <----------------- _(o) ^^--------------(b) °- pso, pso2 ro ri 0 ri ro RADIAL POSITION Fig. 1—Generalized model for oxidation of a sphere of zinc sulfide.

Jan 1, 1970

By J. M. W. Mackenzie, C. J. Wood

In attempting to describe the size classification performance of a hydrocyclone, most workers have elected to use either an equilibrium orbit theory or an non-equilibrium orbit theory. The equilibrium orbit theory has been used by the majority of workers including Lilge,&apos; Bradley; and Yoshioka and Hotta. In applying this theory, it is argued that particles in the body of a hydrocyclone attain an equilibrium radial position where the drag force on the particle resulting from the inward radial fluid velocity is balanced by the outward centrifugal force caused by the tangential component of fluid flow. When considered over the full height of the hydrocy-clone, attainment of this radial equilibrium orbit results in the particle following a conical equilibrium envelope. It is then argued that if this envelope lies outside the envelope of "zero vertical velocity," the particle will report to the underflow, while if the equilibrium envelope lies inside the envelope of "zero vertical velocity," the particle will report to the overflow or vortex finder product. The d50-sized particle which reports in equal quantities to the underflow and overflow is assumed to correspond to particles whose equilibrium envelope is coincident with the envelope of "zero vertical velocity." In considering the equilibrium orbit theory, it is apparent that the horizontal position of the particles in the feed inlet pipe should have no effect on their ultimate destination on the hydrocyclone. Each particle should attain an equilibrium position which depends on the density, size, and shape of the particle; the density and viscosity of the fluid; and the flow patterns within the hydrocy-clone. The nonequilibrium orbit or unsteady state theory has been largely developed by Rietema4 and Mizrahi.6 Mizrahi has listed four main objections to the equilibrium orbit theory. These objections center on the short residence time in the hydrocyclone, the fact that the experimental classification curve is much less sharp than is theoretically predicted, and the absence of negative efficiency conditions in hydro cyclones operating on a feed material which is much finer in size than d50. Proponents of the nonequilibrium orbit theory argue that for a particle to discharge with the underflow it must have sufficient outward radial velocity to reach the downward-flowing region close to the hydrocyclone wall in which the flow lines are parallel to the wall and the ratio of vertical to radial velocity is constant. It is then postulated that a d50 particle entering the cyclone at the center of the feed inlet will just reach this downward-flowing region as it reaches the apex. Thus for uniform distribution of particles across the feed inlet, half the d50 particles—that is, those injected in the half of the inlet area nearest the cyclone wall —will report to the underflow while those injected in the other half will not reach the downward-flowing region and will be carried inward to the center of the cyclone and thus report in the overflow. The exact thickness of the down-ward-flowing region of fluid adjacent to the outer wall of the hydrocyclone is uncertain but Mizrahi considers it to be equal to the apex radius minus the air core radius. Particles larger than d50 have a greater outward centrifugal force acting on them than the d50 particles and may reach the wall even if injected at a distance from the wall greater than Di/2 (Di is inlet diameter). Conversely, particles smaller than d50 may not reach the wall even if injected at a distance less than Di/2 from the cyclone wall. Since the equations put forward by the proponents of both theories yield approximately the same values of d50, it is not possible to decide between these theories by measurement of d50. It should be possible however to examine the theories by injecting a small stream of solids into the feed inlet of a hydrocyclone running on clear water. If the efficiency or classification curve is measured for various horizontal injection positions, then the curves should be coincident if the equilibrium orbit theory holds. If, however, the unsteady state theory describes the cyclone operation, then the classification curves should show finer d50 sizes for particles injected close to the cyclone wall. Experimental A 6-in.-diam hydrocyclone with geometry as in Figs. 1 and 2 was used. Quartz particles were injected as a 50% by wt pulp via an 1/8-in. steel probe. For each in-

Jan 1, 1971

By R. A. Oriani

The partial molal volume of hydrogen is one of the parameters that describe the elastic interaction between the solute and the stress fields about inclusions, dislocations, and cracks. As such the partial molal volume is probably of importance in the elucidation of phenomena such as hydrogen embrittlement and hydrogen yield point. A knowledge of this quantity would also be helpful in thinking about the state of dissolved hydrogen in iron. However, because of the very low lattice solubility of hydrogen in iron the usual ways of determining the lattice expansion are not practicable. It is therefore of interest to apply the thermodynamics of stressed bodies to two sets of measurements of the effect of elastic stress upon the permeability of hydrogen in order to deduce a value of the partial molal volume, VH, of hydrogen dissolved in a iron. Beck ..&apos; and previously de Kazinczy, observed that a uniaxial tensile stress increases the permeability of hydrogen in iron and in various steels. Beck et 01. employed Armco iron and A.I.S.I. 4340 steel, whereas de Kazinczy used a steel the composition of which was 0.13 C, 0.23 Si, 0.46 Mn, 0.006 P, and 0.038 S. Upon releasing the stress the permeability increment disappeared if the stress was below the elastic limit. Both investigators employed cathodic charging to introduce the hydrogen. de Kazinczy measured the permeation through a thin-walled tube by collecting the gas, whereas Beck et al. measured the permeation through sheets of various thicknesses by a very sensitive electrochemical technique. Both investigators measured the steady-state permeation at constant rate of hydrogen ion discharge, and Beck measured in addition the hydrogen diffusivity by a time-lag technique which is independent of the boundary conditions. Beck et al. and de Kazinczy found a linear relationship between log (J,/Jo) and the stress, where Ju/Jo is the ratio of the flux of hydrogen when the metal is under uniaxial tensile stress, a, to the flux under zero stress, for the same temperature and charging current. Beck et al. found in addition that the diffusivity of hydrogen is not changed by stress. Both investigators concluded that the observed change in permeability is due to an increase in hydrogen concentration, and furthermore that the increase in concentration is due directly to the thermodynamic effect of stress upon concentration. Accepting this assessment of the situation for reasons given below, one may use the equation3j4 in order to evaluate I/H, the partial molal volume of hydrogen. This equation is valid for the domain of a/E «¦ 1 (where E is the Young&apos;s modulus) and under the assumption that hydrogen expands the lattice isotrop-ically. From the data of Beck el al, one calculates V7H - 2.0 cu cm per g-atom, and from those of de Kazinczy one obtains 1.8 cu cm per g-atom. That the increase of concentration with stress is indeed of thermodynamic origin is attested to by the facts that the experimental results conform to the thermodynamic relation, Eq. [I], and that the results are the same whether a pure iron1 or each of two different steels172 is used. Neither of these facts wou1.d necessarily be expected if the effect of stress were, rather, to increase the ratio kn/k, of the kinetic factors of the following competing reactions: H(ads) — H (dissolved) Such a change of kinetics at the surface could be an alternative explanation of the effect of stress on the permeability. Although this writer does not deem this explanation to be the correct one for the reasons given above, it must be admitted that unambiguous proof that the phenomenon has a thermodynamic origin does not yet exist. Two kinds of experiments may be suggested. One is to plate a variety of metals on the input surface of the steel and repeat the stress experiment at a variety of hydrogen charging currents. The other is to employ somewhat thicker specimens in order to be able to apply uniaxial compressive stress. Eq. [I] shows that ln(c,/c,) depends on the sign of the stress, but it is difficult to see a physical basis for which k2/kl would depend on the sign of a. The present value of V^ in a iron agrees with phragmen&apos;s5 estimate, which he based on comparisons with the lattice expansion by hydrogen of titanium, zir-

Jan 1, 1967

By R. L. Parsons, Herman Dykstra

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.&apos; THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy&apos;s Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. &apos;Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If&apos; .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951

By V. M. Karve, K. K. Majundar, K. V. Viswanathan, J. Y. Somnay

The effect of calcium, magnesium, iron (both ferrous and ferric) and aluminum ions, which are commonly encountered in a typical beryl ore, was studied in the flotation of pure beryl, soda-feldspar and quartz. The vacuumatic flotation technique was employed. With sodium oleate as collector and in the absence of any activator, beryl floated in a pH range of 3 to 7.5, whereas feldspar and quartz did not float at any pH up to 11.5. The pH range of flotation increased in the presence of the ions studied. With calcium and magnesium ions beryl floated from 3 to 11.5 pH and beyond, soda-feldspar floated beyond pH 6 and quartz floated beyond pH 8. Ferrous ion activation was found to be similar to that of calcium and magnesium. Activation by ferric and aluminium ions was found to be complex and the lower and upper critical pH for all the three minerals was around 2 and 10 respectively. These studies indicated the possibility of separation of beryl from feldspar and quartz even in the presence of calcium, magnesium and ferrous ions between pH 4 and 6. Flotation tests on a mixed feed of pure minerals in a 10 g cell revealed that beryl can be selectively floated from feldspar and quartz if ferric ion is reduced to ferrous state or if it is complexed. Beryl occurs mostly in pegmatites, and hence is associated with feldspar, quartz and micas and small amounts of other minerals such as apatite and tourmaline. The separation of beryl from these minerals is difficult because all the silicates accompanying beryl have more or less the same physical properties. Specific gravities of beryl, feldspar and quartz are 2.70, 2.56 and 2.66 respectively. Electrostatic separation has been suggested but no work has been reported. &apos; The adsorption of sodium tri-decylate tagged with Cl4 on beryl, feldspar and quartz reveal similarity in surface properties. Much work has been reported on the flotation of beryl from ores, either directly or indirectly as a by-product, but little is known about the fundamental aspects of beryl flotation. Kennedy and O&apos;Meara3 laid emphasis on prior cleaning of the mineral surfaces with HF. Mica is removed first by flotation of beryl with oleic acid, around neutral pH. Runke4 introduced calcium hypochlorite conditioning in a final separation stage for activating beryl in a mixed beryl-feldspar concentrate, and after washing to remove the hypochlorite, floated beryl with petroleum sulphonate. The Snedden and Gibbs5 procedure is somewhat similar to that of Kennedy and O&apos;Meara. Emulsified oleic acid is used as collector. Recently Fuerstenau and Bhappu6 studied the flotation of beryl, feldspar and quartz with petroleum sulfonate in the presence of activators and stressed the importance of iron in the flotation of beryl. From the studies conducted in this laboratory, it was found that feldspar and quartz as such do not float with sodium oleate, but in practice selective flotation of beryl from feldspar and quartz in an ore is found to be impossible with sodium oleate as collector. A glance at the chemical analysis of typical beryl ore indicates the presence of several ions like Ca ++, Mg++, Al + + + and Fe+++ in abundance and Ti++++ and Mn++ in traces. Hence, in an attempt to explain the behaviour of feldspar in the beryl flotation, the effect of Ca++, Mg++, Al+++ and Fe+++, which are known as gangue mineral activators7&apos;8 has been investigated. Materials and Methods: Lumps of beryl ore (hand picked) were boiled with 10% sodium hydroxide and washed with distilled water. They were further boiled many times with 10% hydrochloric acid till no positive test for iron was obtained with ammonium thio cyanate. This was followed by thorough flushing with double distilled water. The lumps were crushed in a porcelain mortar and pestle under water. The minus 65 + 100 mesh fraction was used for testing and was always stored under distilled water. Pure feldspar and quartz were similarly prepared and the minus 65 + 100 mesh fractions collected. Inorganic ions tried as activators were ca++, Mg++ , Fe++, Fe ++ and A1 +++ . Calcium nitrate, magnesium chloride, ferrous ammonium sulfate, ferric ammonium sulfate and aluminum nitrate of G.R.E. Merck grade were used. B.D.H. technical grade sodium oleate was used as a collector. The vacuumatic flotation technique developed by Schuhmann and Prakash was used for studying the effect of pH on flotability. 7 The indications given by this work were confirmed by using 10 g miniature cell.&apos;

Jan 1, 1965

By V. F. Zackay, W. A. Goering

ALLOYS of iron and aluminum up to 35 wt pct aluminum are single-phase solid solutions, and are of potentially wide applicability.1-3 In spite of early and continued interest1-4 little progress has been made until recently in the preparation and evaluation of sound alloys containing more than 6 wt pct aluminum. Vacuum-melting techniques for the production of ductile Fe-A1 alloys have been described recently.1-7 A. procedure for air melting these alloys is presented here. Low-carbon iron is induction melted without a slag in a rammed magnesia crucible. At the beginning of melt-down, aluminum pig (99.95 pct Al), charged in a clay-graphite bottom-pouring crucible is placed in a pot furnace at 1800°F. The primary deoxidation of the molten iron after melt-down is effected by the addition of 0.1 pct aluminum and 0.5 pct manganese. (Hilty and Crafts" have reported a significant increase in the deoxidation efficiency of the aluminum and manganese combination over that of the aluminum alone.) A more drastic deoxidation designed to reduce the oxyen content to the lowest possible level is accomplished by plunging metallic calcium to the bottom of the melt. This is done by wiring small cubes of the metal to a steel rod. A circular shield larger than the diameter of the crucible opening is attached to the rod so that any spa&apos;ttering of the molten metal will not endanger the operator. Since the temperature of the molten metal is above the boiling point of calcium, the bath is vigclrously purged by calcium vapor. It is believed that the calcium-vapor treatment permits a homogeneous distribution of calcium in the melt. Owing to the vigor of the reaction the temperature of the molten metal should be kept below 2900°F prior to the calcium addition. A total of 0.05 pct calcium is added in two stages in this manner. The second calcium deoxidation is made just before charging the molten aluminum into the iron, in order that an excess of calcium be present for the remainder of the melt. The aluminum, which has been removed from the holding furnace, is then hydrogen degassed by bubbling chlorine through a quartz tube immersed in the molten aluminum. The hydrogen-chlorine reaction is an exothermic one preventing the solidification of the aluminum during the 5-min chlorination. Approximately 0.1 pct calcium, based on the amount of aluminum, is then added to the aluminum. A further excess of calcium is introduced into the melt in this manner. The oxide dross is removed, fluorspar is added to the molten iron, and the molten aluminum is poured through the fluorspar slag. The fluorspar should be dried thoroughly prior to its use, as any water present will react with the aluminum. Aluminum oxide formed during the pouring operation reacts with the fluorspar slag to form gaseous aluminum fluoride and calcium oxide. A forced-draft ventilating system is required for this operation as aluminum fluoride is toxic. As soon as the molten aluminum has been added, vigorous manual stirring of the melt is required because the slag-aluminum oxide reaction is highly exothermic and tends to take place near the top of the melt. The combination of high temperature and the slagging action of the fluorspar quickly erodes the crucible at the slag line if the aluminum is not stirred uniformly into the melt. It has been found that at least 4 min of manual stirring combined with induction stirring are necessary to ensure homogeneity. The power is shut off 1 min prior to pouring to allow metal and slag to separate. As much slag as possible is removed from the melt, which is then poured directly into cast-iron molds. A mold wash of aluminum oxide is used to prevent ingot sticking. For slab ingots which are to be rolled into sheet, a carbon-tetrachloride vapor atmosphere or a chlorinated-pitch mold wash is desirable, as the aluminum oxide formed in the pouring operation is subsequently removed by the chlorine in the presence of carbon." As in vacuum melting, a pouring temperature of about 2900°F is recommended. Adequate hot-topping is important as iron-aluminum ingots are subject to very deep piping. Ingots are removed from the molds and buried in vermiculite, where they are allowed to cool slowly to room temperature. The ingots are radiographed,

Jan 1, 1959

By J. S. Bowles

THE martensite transformation in lithium, dis- covered by Barrett,&apos; has been studied extensively by X-ray techniques by Barrett and Trautz,² and Barrett and Clifton.V he present paper reports the results of an investigation into the metallographic characteristics of lithium martensite. Such an investigation has not been carried out before. The spontaneous transformation in lithium consists of a change from a body-centered cubic to a close-packed hexagonal structure with the hexagonal layers in imperfect stacking sequence." As far as is known at present, this transformation can be regarded as being crystallographically equivalent to the body-centered cubic to close-packed hexagonal transformation that occurs in zirconium,5 although stacking errors have not been reported in zirconium. From a study of the orientation relationships in zirconium, Burgers5 as proposed that the martensite transformation, b.c.c. to c.p.h., occurs by a heterogeneous shear on the system (112) [111]. The crystal-lographic principle underlying this proposal is that the configuration of atoms in the (112) plane of a b.c.c. structure is exactly the same as that in the (1010) plane of a close-packed hexagonal structure based on the same atomic radius. The pattern in 2v2 both these planes is a rectangle d X 2v2d where v3 d is the atomic diameter. Thus a close-packed hexagonal structure can be built up from a body-centered cubic structure by displacing the (112) planes relative to each other.* This mechanism leads to orientations that can be described by the relations: (110)b.c.e. // (0001)c,p.h.; [111]b.c.c. // [1120]c.p.h Observations confirm these relations. In zirconium, Burgers&apos; measurements indicated an angle of 0" to 2" between the close-packed directions, while Barrett&apos;s measurements on lithium indicated an angle of 3". According to the Burgers&apos; mechanism, the martensite habit plane for this transformation would be expected to be the (112)b.c.c. plane, for this plane would not be distorted by the transformation. One of the purposes of this investigation was to find out whether the observed lithium habit plane agrees with this prediction of the Burgers&apos; mechanism. Experimental Procedure Materials: The lithium was from the same purified ingot used by Barrett and Trautz.² The Bridgman technique was used to produce single crystals. To maintain a temperature gradient in the melt, during the production of these crystals, it was necessary to use a steel mould with a wall thickness of only 0.015 in. Metallographic Techniques: Lithium specimens could be given an excellent metallographic polish by swabbing them gently with cold methyl or ethyl alcohol.? The best results were obtained with methyl alcohol saturated with the reaction product, lithium alcoholate. With higher alcohols the reaction became progressively slower and the attack became an etch pit attack rather than a polish attack. Butyl and amyl alcohols were used for macroetching. After polishing, it was necessary to remove all traces of alcohol from the specimens; otherwise, on subsequent quenching in liquid nitrogen, the alcohol froze to a glassy film. The alcohol was removed with dry benzene. The benzene in turn had to be removed before quenching, but since it does not react with lithium it could be allowed to evaporate. The specimens could then be quickly quenched before they began to tarnish. This operation could be carried out in air on all but excessively humid days when it was advisable to use an atmosphere of dry nitrogen or argon. For examinations at room temperature, the specimens could be transferred directly from the benzene bath into a bath of mineral oil. In mineral oil the specimens oxidized slowly by the diffusion of oxygen through the oil but the structure remained visible for about an hour. Lithium Martensite: Specimens prepared in the manner described above transformed spontaneously to martensite with an audible click when quenched into liquid nitrogen; i.e., M, was above the boiling point of nitrogen (77°K). The disparity between this result and the M, temperature of 71°K, found by Barrett and Trautz, is probably to be attributed to the large grain size and freedom from mechanical deformation of the specimens used in the present work. The relief effects produced by the transformation did not disappear when specimens were quenched from liquid nitrogen into mineral oil at room temperature. This permitted the microstructures to be studied at room temperature where, of course, the martensitic phase was no longer present. Typical micrographs of lithium "martensite" made at room temperature are reproduced in figs. 1, 2, and 3. As anticipated by Barrett and Trautz, the microstruc-

Jan 1, 1952

By D. Havlena, A. S. Odeh

The material balance equation used by reservoir engineers is arranged algebraically, resulting in an equation of a straight line. The straight line method of analysis imposes an additional necessary condition that a successful solution of the material balance equatiott should meet. In addition, this algebraic arrangement attaches a dynamic ineuning to the otherwise static material balance equation. The straight line method requires the plotting of one variable group vs mother variable group. The sequence of the plotted points as well as the general shape of the resulting plot is of utmost importance. Therefore, one cannot progrm the method entirely on a digital computer ar is usually done in the routine solution of the material balance equation. If this method is applied, then plotting and anaIysis are asential. Only the appropriate equations and the method of analysis and interpremtion with comments and discussion are presented in this paper. Illustrative field examples for the various cases treated are deferred to a subsequent writing. INTRODUCTION One of the fundamental principles utilized in engineering work is the law of conservation of matter. The application of this principle to hydrocarbon reservoirs for the purpose of quantitative deductions and prediction is termed "the material balance method of reservoir analysis". While the construction of the material balance equation (MBE) and the computations that go with its application are not difficult tasks, the criteria that a successful solution of the MBE should fulfill have always been a problem facing the reservoir engineer. True and complete criteria should embody necessary and dcient conditions. The criteria which the reservoir engineer uses possess a few necessary but no sufficient conditions. Because of this, the answers obtained from the MBB are always open to question. However, the degree of their acceptability should increase with the increase in the number of the necessary conditions that they should satisfy. Generally, the necessary conditions commonly used are (1) an unspecified consistency of the results and (2) the agreement between the MBE results and those determined volumetrically. This second criterion is usually overemphasized. Actually, the volumetrically determined results are based on geological and petrophysical data of unknown accuracy. In addition, the oil-in-place obtained by the MBE is that oil which contributes to the pressure-production history,&apos; while the volumetrically calculated oil-in-place refers to the total oil, part of which may not contribute to said history. Because of this difference, the disagreement between the two answers might be of paramount importance, and the concordance between them should not be overemphasized as the measure of correctness of either one. In this paper, a third necessary condition of mathematical as well as physical significance is discussed. It is not subject to any geological or petrophysical interpretation, and as such, it is probably the most important necessary condition. It consists essentially of rearranging the MBE to result in an equation of a straight line. This straight line method of the MBE solution has invalidated a few long time accepted concepts. For instance, it has always been advocated that if a water drive exists, but one neglects to take it into account in the MBE, the calculated oil-in-place should increase with time. The straight line method shows that in some cases, depending on the size of the neglected aquifer, the calculated oil-in-place might decrease with time. The straight line method requires the plotting of a variable group vs another variable group, with the variable group selection depending on the mechanism of production under which the reservoir is producing. The most important aspect of this method of solution is that it attaches a significance to the sequence of the plotted points, the direction in which they plot, and to the shape of the resulting plot. Thus, a dynamic meaning has been introduced into the picture in arriving at the final answer. Since the emphasis of this method is placed on the interpretation of the sequence of the points and the shape of the plot, one cannot completely automate the whole sequence to obtain "the best value" as normally done in the routine application of the MBE. If one uses the straight line method, then plotting and analysis are musts. The straight line method was first recognized by van Everdingen, et al,2 but for some reason it was never fully exploited. The advantages and the elegance of this method can be more appreciated after a few cases are carefully treated and worked out by it. SOLUTION OF THE MATERIAL BALANCE EQUATION SATURATED RESERVOIRS The MBE for saturated reservoirs written in AIME symbols is

By G. L. Chierici

Starting from Muskat&apos;s theory of water and gas coning, maximum permissible oil production rates without water and/or free-gas production have been determined, in a broad range of reservoir and well parameters, using the potentiometric model technique. The main assumptions made are as follows: (1) the reservoir rock is homogeneeous (either isotropic or anisotropic); (2) the volume of the aquifer underlying the oil zone is very small, so that it does not contribute to reservoir energy; and (3) the gas cap expands at a very low rate, so that it can be assumed to be in quasi-static conditions. The results obtained are presented in the form of diagrams which can be used for solving two types of problems: (1) given the reservoir and fluids characteristics, as well as the position and length of the perforated interval, determine the maximum oil production rate without water and/or free-gas production; and (2) given the reservoir and fluids characteristics only, determine the position and length of the perforated interval which optimize the maximum permissible oil production rate, without water and/or free-gas production. INTRODUCTION In oil reservoirs where the oil-bearing formation is underlain by an aquifer which does not participate in the production mechanism, water-coning is a limiting factor to the flow rates of producing wells. Production rates are usually kept to a value that will prevent the water from entering the wells. The entry of water into a well lowers its productivity by increasing the weigbt of the fluid column; moreover, the separation of water from the effluent, at the surface, may constitute a very difficult problem in cases of heavy viscous oils. A similar situation is encountered in oil reservoirs with a gas cap overlying the oil-saturated zone; here a downward gas cone is induced by the flow of oil towards the producing wells. Production rates must be low enough to prevent the gas from being produced; producing gas from the gas cap would be a waste of energy. Of course, water-coning and gas-coning phenomena can occur at the same time in the same reservoir if the oil-producing formation is both overlain by a gas zone and underlain by a water zone. Due to its relevant practical importance, the mechanism of coning was studied by many people.2,3,5-8 Defining the conditions for getting the maximum water-free and/or gas-free oil production rate is a difficult problem, often encountered under one of the following aspects: 1. Predict the maximum flow rate that can be assigned to a completed well without the simultaneous production of water and/or free-gas. 2. Define the optimum length and position of the interval to be perforated in a well, in order to obtain the maximum water and gas-free production rate. A systematic study of these problems was made by means of the electrical analog technique. The results of this study are presented here, under the form of a set of curves providing solutions for the above stated problems. These curves are valid only for homogeneous forrnations, either isotropic or anisotropic. Should the formation be non-homogeneous (by horizontal or vertical variation of permeability, shale diaphragms, fractures, etc.), a specific potentiometric study would be required for each specific case. Especially when shale diaphragms of some radial extension are present, the critical rates observed are much larger than would be expected from the diagrams. STATEMENT OF THE PROBLEM In the present study the aquifer is supposed to be of such limited volume that it does not contribute to the energy of the reservoir. Moreover, the gas cap is supposed to expand at such a low rate that the potential gradient in the gas cap is negligible. Under static conditions water-oil and gas-oil interfaces (T1 and T2) are both horizontal. When the reservoir production starts, below each well these interfaces take a cone-like shape (Fig. 1) having as an axis the axis of the well. This shape results from the equilibrium between potential gradients in the oil zone and gravitational forces due to density differences between oil and water and between oil and gas. Assuming the oil-bearing formation to be homogeneous and the oil to be incompressible, the analysis of the problem (see Appendix) shows that the oil-water and gas-oil interfaces are stable only if the oil production rate of the well is not higher than the following values.

Jan 1, 1965

By J. S. Rinehart, W. C. Maurer

The mechanics of impact crater formation in rock, particularly sandstone, has been sutdied, the velocity range being approximately that normally associated with oilwell gun perforators. The bullets were small steel spheres having diameters of 3/16, 9/32 and 7/16 in; impact velocities ranged from 300 to 7,000 ft/sec. The craters have two distinct parts — a cylindrical hole (or burrow) with a diameter the same as that of the impacting sphere, and a wide-angle cup comprising most of the volume of the crater. The burrow is fornred as material in front of the projectile is crushed and pushed aside, forming a cylindrical hole surrounded by a high-density zone. The clip forms as fractures are initiated in front of the projectile and propagate along logarithmic spirals, approximaling maximum shear trajectories, to the free surface of the rock. A most significant observation (made for the first time) was that, below the base of the cup in one type of sandstone, there are a group of similar fractures, not extending to the surface, which are spaced uniformly a few millimeters apart. Each fracture follows roughly the contour of the base of the cup and appears to require a certain threshold impulse to initiate it. These fractures comprise a relatively high fraction of the total, newly exposed surface area. The volume of the material removed by crushing varies as the first power of the impact velocity and the volume removed by fracturing, as the second power of the impact velocity. Penetration varies linearly with the impact velocity and is inversely proportional to the specific acoustic resistance of the target material, the proportionality constant being dependent upon the shape of the projectile. INTRODUCTION Yield of oil from a producing well is frequently enhanced by firing bullets and shaped charges through the well casing into the oil-bearing rock, forming craters and fractures from which oil can flow more readily. The purpose of this investigation has been to develop a better understanding of the mechanics of impact crater formation in rock, particularly sandstone, the velocity range being approximately that normally associated with oilwell gun perforators. FORCES OPERATIVE DURING IMPACT When a projectile moving at considerable velocity strikes a- massive target such as oil-bearing sandstone, intense and complex transient stress situations develop within both the projectile and the rock or sandstone against which it is striking. Usually the struck rock fails, the missile or projectile penetrating into the rock to some depth where it comes to rest or is forcibly ejected from its burrow by expansion of a plug of target material compressed in front of it. When the impact velocity is very high, the projectile itself may fail, breaking apart or becoming distorted; this situation is not considered here, the discussion being limited to nondeforming projectiles. Many experimental studies&apos;.&apos; have been carried out to determine the nature of the mechanics of crater formation and the salient features of the forces coming into play, some of the earliest studies being the French Army experiments performed at Metz between 1835 and 1845.&apos; The stratagem in most instances has been to make a post-mortem examination of the crater, measuring volume and depth of penetration and deducing force relationships from these observations rather than performing the more difficult (usually almost impossible) feat of measuring stresses during penetration. In many materials, the force acting during penetration of the projectile is found to be the sum of two components—(1) a constant force, independent of the velocity, representing some inherent strength of the target material; and (2) a component, proportional to the square of the velocity, representing inertial forces. For such materials, the average force per unit area acting on the projectile at any instant while it is in motion and being decelerated may be written F/A = a + bv2 . . . (1) where v is the velocity of the projectile at that instant, A is the cross-sectional area of the penetrating projectile taken normal to its trajectory, and a and b are constants which are dependent upon the target material and the shape of the projectile. It follows that the total penetration s is given by .........(2) where v, is the velocity of the projectile when it just strikes the target. Values of a and b for spherical projectiles impacting in a loose sand-gravel mixture and compacted earth were obtained in the Metz experiments. For sand-gravel, a and b are 620 psi and 0.0115 (psi) (ft/sec)&apos;, respectively; and for compacted earthworks, a and b are 432 psi and 0.0008 (psi) (ft/sec)&apos;. Figs 1 and 2

By William H. Breeding

The remarks that follow are based substantially on experience covering 45 years, 80% of which has been in placer work, rather than on a review of available literature. Most commercial placers have been deposited by the action of water. The richer and more- difficult-to-mine placers are those in the headwater areas where gradients are steepest. The most lucrative placers are generally in inter- mediate areas where volumes are greater, fewer boulders are present, and gradients are from 3% to 1-1/2%. The higher volume, lower grade placers are in the lower reaches of river systems where gradients are lower. Where gold-bearing rivers have discharged into the sea, wave action can concentrate values on beaches, past and present. Most of the rich, readily accessible placers were mined by our forefathers. Current opportunities exist: (1) in remote areas where infrastructure has been absent in the past, or development has been prohibited by adverse ownership - political or commercial; (2) in deposits that could not be mined by equipment available to our forefathers; (3) in deposits unidentified by our forefathers; (4) where the-price-of-product/cost ratio is substantially better than in earlier years; or (5) a combination of those factors. When I entered the placer business in the late 1930s, and subsequently, a prevailing opinion believed that glacial deposits should be avoided as irregular in mineral content and composition, and unrewarding to explore and develop; yet an operator has been mining a fluvio-glacial deposit profitably for the past 17 years. Rich buried placer channels, of ten called paleo-channels were worked in the last century, generally by hand methods, and under conditions that would be unacceptable today. Exploration and mining equipment now available make some of these channels attractive targets. Well-known examples are in California and Australia. The formation of a commercial placer requires a source of valuable minerals. Above primary deposits, there may be eluvial deposits formed by the erosion of gangue minerals and the concentration "in situ" of valuable minerals. Down slope from these deposits are the hillside or colluvial deposits, and below them are the alluvial deposits of redeposited material. Most of the great placer fields of the world are the result of several generations of erosion and deposition. Well-known examples are in California and Colombia. Gold is a very resistant and malleable material, and gold placers may extend for 64 or 80 km (40 or 50 miles) along a river system. Platinum is less malleable, but is very resistant to disintegration. Diamonds are extremely hard, and (especially gem diamonds) may be found over great lengths of a river system. Cassiterite is less resistant to disintegration, and tin placers seldom extend over two miles without resupply from an additional source or sources of mineralizaton. Tungsten minerals are generally more friable, and within a few hundred yards of the source disintegrate to the point that they are uneconomical to recover. Rutile, ilmenite and zircon placers generally result from the weathering of massive deposits, and may be encountered over extensive areas; most are fine grained and durable. What does a geologist or mining engineer look for in placer exploration? The old adage to look for a mine near an existing mine is still valid. You need a source of valuable mineral. Then you require conditions for concentration, which means a satisfactory gradient and/or other conditions that will permit heavy minerals to settle. Nicely riffled gravel, often called a shingling of the bars, is conducive to placer formation. Coarser gravel is logically associated with coarser gold. Excessive clay and/or high stream velocities in narrow channels can carry gold far downstream and distribute it uncommercially over a large area. When material is extremely fine, in situ weathering and concentration become more important. Placers frequently occur distant from lode mines, and one must remember that in a larger watershed the exceptional floods that occur once in a hundred or a thousand years can move great quantities of material long distances. The carrying power of water is said to vary with the fifth or sixth power of its velocity. I am not ready to disagree with Waldemar Lindgren and accept that many commercial placers are substantially enriched by the chemical deposition of gold from solutions; however, I have seen crystalline gold in clayey material quite distant from known sources of primary gold that is dif-

Jan 1, 1985