Search Documents

By P. B. Baxendell

The performance of a water-drive reservoir having a gas cap depends primarily on the movement of the gar-oil and oil-water contacts. The movement of the contacts during production depends in turn on fluid withdrawals and how the reservoir pressure changes as fluids are produced from the reservoir; that is, on how effectively the aquifer maintains pressure by replacing withdrawals. inasmuch as pressure changes, fluid withdrawals, contact positions, and produced gas-oil and oil-water ratios are interdependent, the analysis and prediction of the performance of a reservoir produced in a given way must take into account this interdependence throughout depletion. This paper presents an analysis which, within the limitations of the assumptions made, yields an engineering approach to predicting future performance based on reservoir pressure and production history. The most significant assumption is the method of extrapolation of future gas-oil and water-oil ratios. The extrapolation procedure smoothly increases both ratios to preselected values as the remaining oil column undergoes a specified decrease in thickness. This preselection is made on the busis of previous field experience in depletion of similar reservoirs under similar conditions. For computing future performance a vdumebic balance L combined with the differential equation defining pressure distribution in the aquifer to obtain positions of water-oil and gas-oil contacts. From these positions are extrapolated produced water-oil and gas-oil ratios. Reservoir performance can be investigated when oil production rates are dependent upon various factors including the performance of the reservoir itself. Examples of practical application of the procedure are included. INTRODUCTION To predict performance of water-drive reservoirs with gas caps and thin oil columns, it is necessary to describe the motions of the fluids within the reservoir during the entire production period to depletion. These motions depend on the pressure changes and on the withdrawal of oil, gas and water. Production of oil from the oil zone primarily causes water to move in to take its place; production of gas from the cap tends to cause oil to migrate into its place. In addition, the volumes of oil and gas remaining in the reservoir depend on changes in the pressure, since a decrease in pressure causes fluid expansion, gas liberation and oil shrinkage. A method of relating the future pressure to total withdrawals is used to describe the motions of fluids within the reservoir under conditions arising in possible modes of production. The analysis is based on two relationships and will reduce the problem to one amenable to digital computation. The first of these concerns the dependence of the pressure distribution in the entire aquifer furnishing the water drive upon the total reservoir withdrawals. This dependence is dictated by the permeability distribution and the extent of the aquifer; at present, such information is most readily obtained from the performance history by means of the resistancecapacitance reservoir analyzer. The second relation involves withdrawals, pressure in the reservoir, and movement of oil, gas, and water within the reservoir. The analysis is subject to certain simplifying assumptions that are necessary to permit solution of the problem. A comparison of the methods of this paper with those presently practiced is pertinent. One method of analysis is to use the reservoir analyzer to predict reservoir behavior based on aquifer characteristics determined from production history. Inasmuch as the total withdrawal rate depends upon the gas-oil and water-oil ratios, which depend in turn upon, among other things, the positions of the gas-oil and water-oil contacts, these ratios must be assumed in advance. In this paper these ratios are, instead, related to the computed positions of gas-oil and water-oil contacts. Thus, our method can be applied to problems in which the oil production rate is limited by produced gas-oil ratio, or to problems in which it is desired to determine the variable gas injection rate that will maintain the gas-oil or water-oil contact stationary. These problems cannot be worked satisfactorily on the analyzer. Another method presently used is based on a paper by Hurst.' Since his procedure is dependent upon using the solution of the heat flow equation, which requires constant permeability within the aquifer, and our procedure recognizes variations of permeability in

By H. Leslie Bullock

From aluminum to zirconium, the quantitative preponderance of quartz as a gangue material is well recognized. lf this material is to be efficiently removed, its variations must be understood. Variations with temperature are especially important. Too little attention has been given to the thermal polarization of quartz. Under closely controlled conditions, electrostatic upgrading is very reliable. For efficient separation, contact charges must be fostered and charges due to radiation, surface coatings, or thermal polarization avoided. This paper lists thermal transition points of quartz and shows their effect on actual separations. Simple separation tests with all factors except temperature held constant are recommended for determining transition points. With silica making up more than 27% of the earth's crust, its oxides comprising more than 59% of all igneous rocks, and quartz accounting for most of the main free oxides, the mining engineer is in constant contact with quartz, which may occur as a valuable mineral to be purified or, far more frequently, as a gangue material to be removed as completely and economically as possible. As a means of effecting such purification or removal, dry beneficiation is becoming more and more desirable owing to local water scarcities, wet waste disposal problems, or freezing conditions. One method that has been gaining particularly rapid acceptance is electrostatic beneficiation, or the separation of dry free-flowing materials by means of opposite surface charges, differences in potential, or differences in conductivity. Electrostatic beneficiation dates back to the 1870's, but only in recent years have newly developed methods and apparatus and a growing knowledge of solid-state physics widened the field for its economical application. Because the term "electrostatic beneficiation" has been rather loosely used in the literature, it has come to include both electrostatic and electrodynamic procedures. Attracting type separators, however, in which oper- ation is based on differences of surface charge or potential, are truly electrostatic, because the separation occurs in a substantially static field set up between oppositely charged surfaces. Separations with this type of equipment may occur at potentials as low as 1000 v and seldom require potentials as high as 30,000 v. In the new contact charge dielectric separators1 the variation in charge is produced by continuous contact and separation of the particles in the moving feed stream and these charges are fostered by the use of non-conducting support and feed surfaces and by handling the feed in the form of streams of appreciable depth. This favors uniformity of feed and allows higher production rates. The distinctive surface charge differences are set up on the separation of the particles according to Coehn's Law,2 which states that equal and opposite charges are generated on the separation of any two materials in contact and that the substance having the highest dielectric constant will be positively charged. The basic contact charge concept is reliable, but all electrical charges are transient and modified by the electrical conditions of the surroundings. Pyro-electric, photoelectric and radiant effects may modify or totally destroy the contact charges necessary for efficient separation, or contact with conducting surfaces may neutralize them. Such hostile conditions must be carefully guarded against, since they have led to many costly failures in the past. The most consistent difference in surface charges to insure good separation is produced by repeated uniform contact and separation of particles in the moving stream. The thickness of the feed stream possible with this method reduces the effect of contact with the supporting surfaces, but as some contact is inevitable, the best results may be assured by having the dielectric constants of the supporting surfaces between the dielectric constants of the substances to be separated. In general, the hard smooth surface of quartz makes it an ideal substance for electrostatic separation from most minerals. For instance, with calcite, starting with a feed containing 1.9% acid insolubles, one can produce a concentrate containing 0.30% acid insolubles with a tailings containing 21.7% acid insolubles and a yield of 92.8%. The color can be held at 92 or above and the tint at 1.7 or lower. Working with specularite iron ore, laboratory work has given a concentrate of 68.8% Fe., with an iron unit recovery of 96.7%. Excellent results are also in pros-

Jan 1, 1969

By E. A. Farrell

A comprehensive survey is made of the status of the asbestos industry as it relates to marketing the product. Included are descriptions of the various types of asbestos and the grading and classification systems used. The uses of asbestos, distribution practices, and types of ore bodies are all related to marketing. World production, the producers and their capacities and world consumption for 1966-67 are summarized and statistical data are included. Asbestos is a general term describing a family of fibrous minerals of the serpentine and amphibole mineral groups. Asbestos has a long history going back to the time of the Egyptians, when it was used as a lamp wick. The commercial mining of asbestos started in Canada, Russia, and Africa in the 1800's, and the first asbestos products were made in Italy and Russia. The five main types of asbestos are: chrysotile, accounting for 95% of the total mined, amosite, crocido-lite, anthophyllite, and tremolite. Canada, Russia, and Africa are the major producers of asbestos. The commercial utility of asbestos was at first based on the heat resistance of the fibrous mineral in the form of packing, at the start of the industrial revolution. Its current utility is based more on its ability to reinforce binders such as portland cement, rubber, and plastics. Its inertness to the chemical nature of most binders is unique. Most important is its ability to maintain its reinforcing utility when the product is exposed to weather and soil conditions as in asbestos cement boards and pipe, and heat, pressure, and chemical exposure as in brake linings and gaskets and packings. The mineral asbestos is also unique because its fibrous form permits it to be spun and woven to cloth or formed into paper. Many asbestos applications are critical to national defense and at the present time, there are no satisfactory substitutes. Grading and Classification Canadian chrysotile asbestos fiber is graded and priced by length since basically the longer the fiber the higher the utility. The Canadian asbestos industry does not, however, classify the fiber by direct length measurement, but by a dry screening test. The method is called the Quebec Standard Screen (QS) test. One pound of fiber is mechanically shaken in four vertically stacked sieve boxes. The relative proportions remaining on the sieves defines the grade. The longer the fiber, the larger is the amount that stays on the top coarse screens and the less on the lower, finer mesh screens. Other tests can be used to further define the length distribution of fiber such as the wet screen Bauer McNett and the Suter Webb Comb (3 group only). The Canadian grading system divides the milled fibers into 5 main groups: group 3, 4, 5, 6, and 7, with 3 group being the longest, and 7 the shortest. Each group is further divided to subgrades, identified in each group by the letters A to Z, with "A" the longest and "Z" the shortest. See Appendix 1 for the Canadian QS classification system. The Russians also use the QS test on chrysotile. The Africans classify their chrysotile into grades similar to Canadian. The African crocidolite and amosite, however, are classified into actual length groups such as l to 2 in. and 2 to 3 in. Amosite and crocidolite are generally longer than chrysotile but also more brittle. Milled asbestos is not composed of staple length fibers like fiber glass or cotton, but of a mixture or blend of fibers ranging from long to short. Milled asbestos has a fiber length distribution similar to the particle size distribution of a powder. For example, the longest group 3 chrysotile grades have a high percentage of the longest fibers (1/2 to 3/4 in.) and low percentage of short fibers (0.003 in.). The figures in Table 1 give the approximate length distribution of the longest, middle, and shortest groups. A second and further method of classifying fiber is the degree to which the fiber bundles are separated to form a larger number of smaller diameter bundles. This property is normally described as the degree of fiberization, openness or surface area. Air permeability tests are used to measure surface area. Asbestos pro-

Jan 1, 1971

By T. A. Wilson, S. L. Hoyt

The plastic deformation of slender single crystals of zinc has been described in some detail in the paper by Mark, Polanyi and Schmid,' which has become a classic, and also by one of the present authors in a somewhat shorter account. 2 The study of single-crystal zinc is termed classical because zinc single crystals offer, perhaps, the best material yet available for the study of atomic behavior during deformation and the effect known as "strengthening" in metals. A considerable amount of additional work has also been done on single-crystal zinc, but as it is not closely related to the present paper it will not be considered. Throughout all the former work, the mechanism of plastic deformation first described by Mark, Polanyi and Schmid has been assumed to hold. This may be said to offer the strongest confirmation of its correctness. Even so, the picture of the process has never been as completely portrayed as is desirable. The first object of the present paper is to describe work that has been done in this field, and which is even yet being carried on, in the hope that a more complete picture may be obtained. Mathewson and Phillips3 have recently described a new mechanism of the deformation of zinc based on their study of large rectangular crystals. One of their conclusions was that deformation produced twinning with a rotation of some of the basal planes into positions 94" removed from their original position; a position almost the same as that of the prismatic planes before twinning. A second conclusion was that fracture occurred along these basal planes in their new position, and, therefore that fractures previously regarded as prismatic were in reality basal. Such findings are of great significance and they raise a question as to the generality of such behavior. Professor Mathewson's opinion is that even the slender cylindrical single crystals behave in the same manner as his large rectangular crystals when strained by simple tension.

By H. Udin

G. KUCZYNSKI* and B. H. ALEXANDER*—This paper represents a most noteworthy attempt to evaluate experimentally the surface tension of a solid metal. Because of the great importance of such measurements, any proposed method should receive the closest scrutiny before the results can be considered reliable. In regard to the experimental method, we think that the marking of the gauge length by means of tieing knots in the wire may be the cause of some of the spread in the results. Such a knot may be expected to tighten slightly, and thus increase the gauge length, when placed under stress at high temperature. Although this effect would be very small, amounting at most to only a few times the wire diameter. A fairly tight knot in a wire will decrease the wire length by about ten times the wire diameter, thus only a slight tightening of the knot would cause considerable spread in the results. Upon plotting the stress strain curves from the authors' data, the writers found that there was a fairly consistent tendency towards an S-shaped curve, instead of a straight line. Such an effect could be caused by the tightening of the knots. The writers think, however, that the experimental results are fairly reliable, but that there may be other methods of interpreting them depending upon what mechanism is assumed to be responsible for the shrinkage of the wires. The authors have assumed that the stress due to surface tension results in viscous flow. It should be made clear that it has never been demonstrated that viscous flow can occur in metal crystals even at very high temperatures. The experiments of Chalmers13 on tin, which are so frequently quoted as giving evidence of viscous flow at low stresses are by no means satisfactory. In his experiments, Chalmers found that only the initial rate of flow was approximately proportional to stress. He also found that the rate of flow varied markedly with time which, in his experiments, was less than 2 hr. Inasmuch as there is no proof of viscous flow in metals, and the authors have brought forth no conclusive evidence on this point, it may be worth while to investigate other possible mechanisms of material transport which would account for the shrinkage of the wires. The writers wish to point out that in these experiments the shrinkage of the wires can be adequately explained, according to a self diffusion mechanism. Thus, if we assume a concentration gradient for self diffusion which is a function of the radius of curvature of the wires, and assume that diffusion will occur so that the total surface area is decreased, we find the following expression for the self diffusion coefficient: where k = Boltzmann constant r0 = initial radius of the wire T = absolute temperature ? = surface energy 8 = interatomic spacing t = time e = strain at zero applied stress Eq 19 may be used to evaluate the self diffusion coefficient of copper, using the strain measurements obtained by the authors for zero stress as obtained by extrapolating their curves for 5 rail wires. By inserting a reasonable value for the surface energy (1500 ergs per cm2) we find: -66,000 D = 5 X 10e RT [20] The activation energy is of the correct order of magnitude, but the frequency coefficient is much too high, indicating that surface diffusion may be playing an important role. This discrepancy in the action constant is much smaller than the corresponding discrepancy obtained by the authors for the viscosity coefficient. The writers by no means propose that this proves that the shrinkage of the wires is due to self diffusion but we merely wish to point out that there are explanations other than that given by the authors. In this, as in any kinetic phenomena, it is necessary to study the rate of the process before anything can be said about the mechanism. The determination of surface tension given by the authors is based upon an interpretation of the data which embody the concept of viscous flow. The final proof of this concept will be obtained only after the time relationships confirming the authors' Eq 15 have been conclusively established. The rough linearity of the stress strain curves obtained by the authors for experiments run the same length of time should not be considered as proving that viscous flow is occurring. H. UDIN (authors' reply)—All of the test specimens were annealed at 1000°C for an hour or more before preliminary measurements were made. During this anneal the wires recrystallize, and the greatest part of grain growth takes place. Also, the knots sinter at the cross-over points. This does not in itself eliminate the possibility of end errors, although it greatly decreases their probable magnitude. It is still possible that some extension occurs due to creep in shear at the sintered points. If so, this effect would be quite independent of and superimposed on the normal shrinkage or extension of the wire itself. Within the precision of the experimental results, straight lines satisfy the data as well as do any other simple curves. Until data of greater precision are obtained, it is futile to discuss any possible trends away from linearity. The disagreement between Kuczynski and Alexander's Eq 19 and our Eq 18 is one of semantics and mathematics, not mechanism of flow, since Eq 18 is based on the self-diffusion concept of viscous flow. It would be interesting to learn how the mathematics leading to Eq 19 deviates from that of Eyring and of

Jan 1, 1950

By N. S. Stoloff, R. G. Davies

A study has been made of the efject oj different dislocation-precipitate interactions upon the temperature dependence of the flow stress of aged Ni-14 at. pct A1 alloy. It is observed that when the dislocations bow between widely spaced (-20004 coherent Ni3Al particles the flow stress decreases with increasing temperature in the normal way. However, when the dislocations cut closely spaced (-5004 particles the flow stress is independent of temperature from -100 to 600°C, due to a balance between softening of the matrix and an increase in strength of the particles with increasing temperature. The retention of strength at high tempera-tures of commercial nickel-base alloys, which are strengthened by the precipitation of a phase based upon Ni3Al, is thought to be due to the unusual strength properties of Ni3Al. The flow stress of Ni3Al increases continuous1y from -196"C to a maximum at -600"C. It is concluded from a series of thermal-mechanical tests that the sevenfold increase in flow stress over this temperature interval is due to a lattice effect and is not diffusion-controlled. The flow stress of precipitation- or dispersion-hardened materials depends on the resistance to dislocation motion within the matrix and the extra energy required for dislocations to bow between or to cut particles. If the dislocations bow between the particles or if the strength of the cut particles is constant with temperature, then the flow stress of the precipitation-hardened alloy must decrease with increasing temperature due at least to the decrease in elastic modulus of the material. There will be softening also from thermally activated cross-slip or climb, offering an additional degree of freedom for dislocations to avoid particles. For example, in the case of nickel containing a dispersion of thoria,' which most probably deforms by dislocations bowing between particles, the flow stress decreases by about 50 pct between 25" and 650°C. In A1-Cu alloys2 aged to produce the 8" precipitate, dislocations cut the particles, and the flow stress decreases by about 20 pct between -269" and 25°C. However, many commercial high-temperature nickel-base alloys, for example Inconel-X and Udimet-700, exhibit little or no decrease in flow stress with increasing temperature up to about 700°C. A characteristic feature of these alloys is that they are strengthened by the precipitation of a phase based upon Ni3A1. Guard and westbrook4 and flinn' have shown that Ni3Al (and alloys in which a third element such as molybdenum or iron is substituted for part of the aluminum) is unusual in that the hardness and flow stress increase with temperature to a maximum at about 600°C. For the flow stress of a precipitation-hardened alloy to be independent of temperature we propose that the particles must be cut by dislocations moving through the matrix and that the strength of the particle must increase with increasing temperature. Theories of precipitation hardening do not take into account the flow stress of the dispersed particles that are cut during deformation; the only dissipative process usually considered7 is the creation of interface within the particle and between the precipitate and matrix. The purpose of the present investigation has been to study in detail the temperature dependence of the flow stress of a nickel-base alloy strengthened by the precipitation of Ni3Al in two structural conditions such that when deformation occurs it does so by dislocations a) bowing between the particles and b) cutting the particles, respectively. A simple binary Ni-14 at. pct A1 alloy was chosen because considerable information is already available for this system concerning phase equilibria and precipitation reactions and rates.' Dislocation-precipitate interactions in the binary alloy should be similar to those in the more complex commercial alloys. In addition, the mechanical and physical properties of NisAl were studied in detail in the hope of elucidating the mechanism by which the strength increases with increasing temperature up to 600°C. EXPERIMENTAL PROCEDURE For the study of the effect of precipitation of Ni3A1 upon the temperature dependence of the flow stress, an alloy containing 14 at. pct A1 was utilized; a Ni-8 at. pct A1 solid-solution alloy was employed as a comparison material. Vacuum-cast ingots were hot-rolled at 1000°C and cylindrical compression samples, 0.20 in. diam by 0.40 in. high, were prepared from the 1/4-in.-diam rod. Specimens were recrystallized and solution-treated at 1000°C for 1/2 hr and then water-quenched. A preliminary study revealed that, when the Ni-14 at. pct A1 alloy was aged for 1 hr at 700°C, significant precipitation hardening was obtained, and that the structure was free from grain boundary discontinuous precipitation; an overaged condition was produced by annealing the aged specimens at 850°C for 1 hr. To circumvent the difficulties involved in the hot rolling and swaging of Ni3A1, compression samples,

Jan 1, 1965

By Che-Yu Li, J. L. Gregg, W. F. Brehm

Oriented, oxygen-doped niobium bicrystals were tested in liquid lithium. The grain boundaries were attacked preferentially. The depth of the penetrated zone varies as (time)2. The penetration was aniso-tropic, had a high activation energy, and increased with the increased oxygen doping level. A possible model was proposed to account for the experimental observations. 1 HE grain boundary penetration of a metallic system by liquid metal has been studied by several investigators. Their results are summarized by Bishop.' Most of these works show that the penetration by liquid metal corresponds to the phenomenon of liquid metal wetting. In the case of a grain boundary, wetting will occur when twice the solid-liquid interfacial tension is smaller than the grain boundary tension resulting in the replacement of the grain boundary by two new solid-liquid interfaces. Other possibilities exist; for example, the atoms of the liquid metal may diffuse into the grain boundary region due to chemical potential gradient. The gradient can be produced by impurity segregation or simply be due to the increase in solubility in the grain boundary region. The penetrated grain boundary in these cases may remain solid at the test temperature. The Nb-Li system has been of considerable interest because of its possible technological applications. For fundamental interest it provides a possibility of studying the grain boundary penetration process which is not controlled by the wetting mechanism. The pure niobium is not attacked by the liquid lithium, but if niobium containing more than 300 to 500 ppm oxygen by weight is exposed to liquid lithium, corrosion occurs at the solid-liquid interface and preferentially at grain boundaries. Previous investigators2-' have proposed that this preferential corrosion at grain boundaries is caused by oxygen segregation there, with subsequent inward diffusion of lithium to form a Li-Nb-0 compound. These investigators also found that the corrosion could be retarded by adding 1 pct Zr to the niobium to precipitate the oxygen as ZrO2 upon proper heat treatment. However, there are no quantitative data on the kinetics of the grain boundary penetration process to test the validity of the proposed corrosion mechanism. In this work an investigation of this penetration process in oriented bicrystals was made as a function of the oxygen doping level in the bulk niobium and the grain boundary orientation. A possible model for the penetration process based on the experimental results was proposed. EXPERIMENTS Oriented niobium bicrystals were grown by arc-zone melting oriented single-crystal seeds.7 These bicrystals contained simple tilt boundary. The [001] directions in the two grains were tilted about a common [110]. The bicrystals were 31/2 in. long and 5 by 4 in. in cross section with the straight, symmetric, planar grain boundary longitudinally bisecting the crystal rod. The bicrystals were doped with oxygen by anodically depositing a layer of Nb2O on the surface in a 70 pct HNO solution at 100 v, using a stainless-steel cathode. The specimens were homogenized by annealing in evacuated quartz tubes at 127 5°C. Oxygen content of the niobium was measured from microhardness values, after DiStefano and Litmman.' Supplementary checks were made with vacuum-fusion analysis.7 Individual test specimens cut from the doped bi-crystal rods, about by by % in. in size, were tested inside double jacket sealed capsules. The inner jacket was niobium, the outer was stainless steel. The niobium inner jacket eliminated the problem of dissimilar-metal mass transfer.' The lithium (99.8 pct pure, obtained from Lithium Corp. of America) was handled only in a purified argon atmosphere in a Blickman stainless-steel glove box. After introduction of lithium, the capsules were sealed by welding. Further detailed experimental procedures are given in Ref. 7. The capsules were heat-treated in vertical Marshall resistance furnaces. Temperatures were controlled to When heating above 1100°C, it was necessary to seal the furnace work tube and flow argon through to prevent failure of the stainless-steel outer jacket of the capsule. Tests were made on 6" 2", 16" 2, and 33" i2" bicrystals at oxygen levels up to 2600 ppm by weight in the 6' and 16" crystals and with 1300 ppm oxygen in the 33' crystals. The oxygen levels were controlled to 100 ppm. Most of the quantitative data were obtained from 16" bicrystals between 800" and 1050°C. The capsules were quenched into water after the test and cut open with a water-cooled abrasive wheel. The capsules were then submerged in water, which dissolved the lithium and freed the specimen. Measurement of the depth of the penetrated zone in the grain boundary was done either on metallographically prepared surfaces or directly on the grain boundary plane after the specimen was fractured in tension in the grain boundary plane. The depth of penetration measured by both methods agreed well. Further details describing these techniques have been reported elsewhere.'p7

Jan 1, 1969

By Rudolph G. Wuerker

Charles T. Holland (Head, Dept. of Mining Engineeri*, Virginia Polytechnical Inst., Blacksburg, Va.) Mr. Wuerker has presented a very interesting discussion of the use of triaxial test methods for investigating the strength properties of rocks. Such methods, no doubt, eventually will develop considerable information of interest to those concerned with the design of mine layouts, particularly in the field of pillar design. From his discussion of my recent article, "Cause and Occurrence of Coal Mine Bumps" (Holland Mining Engineering 1958, p. 933-1002), it is evident that in one place at least I did not make my meaning clear to him and perhaps others. To clear the matter up I think it best to quote from the article, somewhat more fully than did Mr. Wuerker, as follows: "4) In actual operations — because rocksare not perfectly elastic, homogeneous, nor isotropic and because local yield does occur — the maximum stress as demonstrated by Phillips (Ref. 22, pp. 64, 65) and indicated by much experience in mining, does not occur at the walls of the opening but at a short distance inside the pillar. Furthermore, the maximum stress does not reach as great a value as theoretical considerations and laboratory experimental methods indicate.* Actual distance inside the pillar, measured from the wall, at which the maximum stress exists, has not been determined. Observations in many mines, however, indicate that this distance could have a mini-value of one to six or eight times the bed thickness and that it is probably affected by width and height of the opening, depth of cover, and relative values of the elasticity and plasticity of materials comprising the roof, floor, and coal seam. The actual value of the stress produced probably lies between the theoretical maximum and the average stress concentration that would be produced if the weight of the strata above the unsupported opening were evenly distributed over the pillars for a distance equal to the opening width." The footnote reference in the above quotation referred to the following: "*For example, the Pocahontas No. 4 coal bed in southern West Virginia is mined under cover up to 1800 ft thick. Development openings are driven 18 to 20 ft wide, and the bed is about 6 ft thick. According to the work of Panek, the tangential wall stress at mid-bed height under these conditions would reach values between 4000 and 5000 psi. Actual tests of 3-in. cubes of this coal show its compressive strength would be much less than this, perhaps as low as 400 psi. Yet the pillars usually show no evidence of failure in these headings. In this same bed at a depth of 800 ft, the author has seen an opening 225 ft between supports lying between two old groves approximately 1100 ft apart. According to the theoretical considerations, the stress in the pillar walls would have been about 18,000 psi, yet the pillar showed little or no evidence of weight. In view of these observations, it is clear that the wall stress does not attain the maximum values indicated by theory." (Underlining added to original wording.) By referring to Fig. 2A of my paper it will be noted that theoretically the maximum pillar stress would occur at the pillar wall, i.e., at the passageway surface of the pillar. Obviously this cannot be correct in the cases of stress ranging from 4000 to 18000 psi since the coal at the surface of the pillar is under no constraint and cannot have a strength much greater than 400 or 500 psi. Hence, my conclusion that the maximum stress does not occur at the wall but back in the pillar some distance from the wall. Since these stresses are pushed back in the pillar from the wall, it is also obvious that the loads transferred to the pillar from the opening will be spread over a greater area and hence Pillar stresses will not rise to the values postulated by theory and photoelastic experiment. Further since to visual inspection the coal along the pillar wall did not appear to be failed the conclusion was reached that the stress shift was caused by local elastic or plastic yield and by difference in the elastic modulus of the rocks composing the mine floor, mine roof, and coal bed. Later on under the heading "Strength of Mine Pillars" (pages 1000-1002) the effects of constraint is briefly described. Also a formula taking into account constraint is developed relating pillar strength to the uniaxial strength of coal and the L/T ratio of the pillar. Since my paper was written, reports of experiments conducted in South Africa (Denkhaus, et. al., 1959), in Sweden (Hast 19581, and in Canada (McInnes, et.al., 1959) reveal that the conclusion expressed relative to the existence of a low stress area existing around the edges of pillars and solid faces as described above is generally correct. But it seems possible that where the wall stress developed is less than the unconfined strength of the rock composing the pillar and where the roof, floor, and pillar

Jan 1, 1961

By J. R. Cahoon, H. W. Paxton

The possibility of using unidirectionally solidified, two-phase alloys as an approximation to fiber composite materials is investigated. The short-term me.chanical properties and failure modes of unidirectionully solidified A1 (rich)-Cu alloys containing ap -Proximately 0, 17.5, and 27.7 vol pct of 0 phase 'fibers" are determined at temperatures from 25" to 500" and compared with those obtained for conventionul SAP alloys. In a previous publication,' hereafter referred to as I, the possibility of understanding some of the room-temperature mechanical properties of unidirectionally solidified castings was explored. For Al(rich)-Cu and Al(rich)-Mg two-phase alloys over a substantial range of compositions, the yield and ultimate strengths and common ductility measures were very adequately predicted from the principles of fiber strengthening4 and the analysis of ductility outlined by Gurland and Plateau." The results obtained in I suggest the possibility of using unidirectionally solidified, two-phase alloys to simulate fiber composite materials where the inter-dendritic second phase or constituent acts as the reinforcing material. Recent attempts concerning the fabrication of fiber conlposites have concentrated on producing composites with a good bond between fiber and matrix and with very long fibers so that their maximum contribution to the strength of the composite may be realized. However, these objectives are difficult to attain in practice and present fabrication processes are either extremely laborious or costly.13 The slow, unidirectional solidification of eutectics has received considerable attention as a method for producing composite materials. 5,6 This method can fulfill both of the above objectives but it is currently laborious, expensive, and has the additional disadvantage that the volume fraction of reinforcing phase cannot be easily varied. On the other hand, unidirectionally solidified, two-phase alloys, also with a good bond between the phases, are relatively easy to make and the volume fraction of reinforcing "fibers" can be easily varied by changing the average composition of the alloy. The disadvantage of the cast alloys is that the mechanical effectiveness of the "elongated interdendritic reinforcements" (EIR)* may be reduced due to their rela- tively short lengths, the w factor in Eq. [2] of I. However, if the EIR have a high strength their contribution can be considerable. For composite materials containing discontinuous cylindrical fibers of various lengths the ultimate strength is given by1 where it is assumed that the composite fractures when the fibers fail. In Eq. [I], a, is the stress in the matrix just prior to failure of the composite, Vf is the total volume fraction of fiber reinforcing constituent, Vf(l+) is the volume fraction of fibers whose lengths exceed the critical length, I,, which is defined as the shortest length of fiber in which the stress can build up sufficiently to break the fiber. af is the fracture strength of the fiber material, w is a factor accounting for the discontinuity of those fibers whose lengths exceed I,, 1-/d is the average aspect ratio of those fibers whose lengths are shorter than I,, and t is the shear stress in the matrix at the fiber-matrix interface. The factor w is dependent on the length of the fibers and also on whether deformation of the matrix occurs plastically or elastically. However, for a given length of fiber, w is smaller when elastic deformation of the matrix is assumed.' It is of interest to consider the properties of simple unidirectionally solidified, two-phase alloys at elevated temperatures in view of the possibility of using suitable modifications for high temperature service. Knowledge of the creep behavior of these materials is still rudimentary (although under active investigation) and the present paper concerns itself with short time tensile properties of some alloys similar to those investigated in I (i.e., unidirectionally solidified Al(rich)-Cu alloys). Unidirectionally solidified alloys containing 5.6, 17, and 23 wt pct Cu were tested parallel to the direction of solidification at temperatures from 25" to 500°C. In the present investigation, the alloys were homogenized for 2 days at 535°C giving a matrix of homogeneous a phase (5.2 wt pct Cu) and an interdendritic constituent (EIR) which was completely Q phase (53 wt pct Cu). EXPERIMENTAL Alloys of nominal composition 5.6, 17, and 23 wt pct Cu (containing approximately 0, 17.5, and 27.7 vol pct 8 phase, respectively, after homogenization at 535°C) were prepared by melting 1200 g of A1 (99.99 pct) in a high purity graphite crucible and adding the appropriate amount of freshly cleaned copper chips (99.9 pct). The molten alloy (at 700°C) was poured into a preheated graphite mold (also at 700°C) and the ingot unidirectionally solidified by impinging water on the steel baseplate of the mold. The alloy was degassed immediately

Jan 1, 1970

By D. E. Coates, G. R. Purdy, S. V. Subramanian

The morphological stability of the planar solid-liquid interface in dilute ternary alloys, undergoing steady-state unidirectional solidification, is analyzed in terms of both the constitutional supercooling principle and the perturbation methods recently developed by Mullins and Sekerka. First, various steady-state solutions for the two solute distributions ahead of a planar interface are examined. The nature of the solutions depends on the size and concentration dependence of the off-diagonal diffusion coefficients. W~thin the framework of the constitutional supercooling principle, a cumulative contribution to instability frorn the two solutes is found to exist in the absence of diffusional interaction. It is shown that the latter can produce a further enhancement of instability or can have a stabilizing influence, depending on the form of the liquidus surface and on the sign of the solute-solute interaction. A perturbation analysis, which ignores diffusional interaction, verifies the cumulative influence of lhe solute fields and demonstrates that the Mullins-Sekerka stability criterion for binary systems (with capillarity accounted for) can be readily extended for application to ternary systems. SOME time ago, Tiller et al.' calculated the solute concentration distribution ahead of the planar solid-liquid interface of binary alloys undergoing steady-state unidirectional solidification. An earlier qualitative proposal that the transition from planar to nonplanar growth morphologies is associated solely with the onset of constitutional supercooling in the liquid layer ahead of the moving interface2 was used in conjunction with this calculation to put the now well-known constitutional supercooling (C-S) stability criterion into quantitative terms. Mullins and Sekerka,3 in a recent and very elegant analysis, established a more complete criterion (hereafter referred to as the M-S criterion). Interfacial stability was investigated by determining the time derivative of the amplitude of a sinusoidal perturbation of infinitesimal amplitude which had been introduced into the originally planar shape of the moving interface. Of particular importance is the fact that capillarity was included in the boundary conditions of their calculation. The purpose of the present paper is to extend all of this earlier work on dilute binary systems for application to dilute ternary alloy solidification. The analysis is divided into three sections. In the first the two solute distributions ahead of a moving planar interface are considered. Mathematical solutions are de- termined for situations in which: a) diffusional interaction is negligible, 6) diffusional interaction must be considered but circumstances permit use of constant diffusion coefficients, and c) the concentration dependence of off-diagonal diffusion coefficients can be described by first-order dilute solution approximations. In the next section, a stability criterion analogous to the C-S criterion is developed and the influence of diffusional interaction on interface stability is analyzed. Finally, the perturbation formalism of Mullins and Sekerka, with capillarity included in the boundary conditions, is extended for analysis of ternary systems in which diffusional interaction is negligible. The study of interface stability in binary systems usually commences with the assumption that the equilibrium distribution coefficient and the slope of the liquidus line are constant at values corresponding to infinite dilution. Similar assumptions have not been introduced into the present treatment; that is, we do not assume planar solidus and liquidus surfaces joined by tie lines which yield constant distribution coefficients. The latter involves the assumption of no ther-modynamic interaction between solute species in both the solid and liquid. We consider a ternary phase diagram for which the solidus and liquidus surfaces are, in general, nonplanar and of course pass through the corresponding binary solidus and liquidus lines. These lines are not assumed to have constant slope. In the dilute regions we are concerned with, the following assumptions are made: i) The solidus and liquidus surfaces are of a form such that both the solidus and liquidus temperatures are monotonically varying functions of each solute concentration. ii) The tie lines are such that the equilibrium distribution coefficient of a given solute is greater than unity for every point on the solidus (or liquidus) surface or it is less than unity for every point. STEADY-STATE SOLUTE DISTRIBUTIONS IN THE LIQUID As will be demonstrated in the next section, a knowledge of the steady-state solute profiles is not a necessary prerequisite for the formulation of a ternary C-S stability criterion. However, in that details, such as the complete description of the equilibrium liquidus temperature profile, require an evaluation of the solute distributions, the overall treatment is enhanced if these distributions are determined. Consider a ternary system (solvent plus solutes 1 and 2) for which a planar solid-liquid interface is in unidirectional motion at constant velocity V. At this stage it is unnecessary to limit ourselves to dilute solutions. For a stationary frame of reference the generalized forms of Fick's equations are:

Jan 1, 1969

By A. N. Holden

A DISLOCATION mechanism has been described by Cottrell' by which metals can yield locally, I. form Liiders bands, giving rise to a characteristic stress-strain curve with a sharp yield point and appreciable strain at constant or decreasing stress. It is undoubtedly the best mechanism that has been suggested to date." In its present development, however, the dislocation mechanism provides a more satisfying explanation for the sharp yield point than for the extensive localized flow occurring at the lower yield stress. The primary objective in this paper is to extend the dislocation mechanism to account for localized cataclysmic flow by a dislocation collision process and to give experimental evidence to support such a process. Only the yielding of iron containing carbon -will be discussed, although other metal-solute systems are known to behave similarly. Cottrell Mechanism In brief, Cottrell explains the yield point in the following way: The dislocations in iron which must propagate to produce slip usually lie at the center of local concentrations of carbon atoms, since segregation about these dislocatlons relieves some of the local stress resulting from them. A dislocation surrounded by a "cloud" of carbon atoms is thus anchored, and a higher stress is required to set it in motion than to move a free dislocation. Considering all available dislocatlons to be anchored in this fashion, the iron exhibits a yield point when the first dialocations break free and move through the lattice causing slip. This first breaking away of a dislocation enables other dislocations to break loose by "interaction" and the process becomes a cataclysm producing local deformation or Luders bands. The yield point in the stress-strain diagram for iron is absent in freshly deformed material, but returns gradually with time; the phenomenon is one aspect of what is called strain aging. The rate at which the yield point returns following straining depends on the temperature of aging. According to Cottrell the rate of return of the yield point in strained iron is limited by the rate of diffusion of carbon at the aging temperature, the mechanism is onr: of reforming the solute atmospheres around carbon-free dislocations that had stopped moving coincident with the removal of stress. If the specimen is retested immediately after straining and unloading, carbon will not have had time to diffuse to, and re-anchor, dislocations and the yield point will not occur. The carbon diffusion limitation for the rate of strain aging apparently applies if the criterion for strain aging is either the change in hardness" or the change in electrical resistance" of the strained speci- men with aging time. The possibility exists, however, that the yield point actually returns to strained iron at some rate other than that deduced from hardness or electrical resistance data. Therefore, as a preliminary experiment, the rate of yield point return in a rimmed sheet steel strained 6 pct in tension was measured at 27°, 77°, and 100°C. A plot of yield-point elongation for each of these temperatures against aging time appears in Fig. 1. The aging process is described by curves which rise to a plateau value of elongation that seems independent of temperature, but at a rate that depends on temperature. Very long times lead to a further rise in the yield-point elongation above the plateau value. However, if the later increase in yield-point elongation is ignored and the log of the time to reach half the plateau value of elongation is plotted against 1/T, a straight line results for which an activation energy of about 25 kcal pel- mol may be assigned. Within the accuracy of this sort of experiment this is approximately the activation energy for the diffusion of carbon in iron (20 kcal per mol), and the carbon diffusion limitation suggested for the yield-point return on strain aging is valid. The Cottrell mechanism thus explains in a qualitative manner the occurrence of a yield point in iron and its return with strain aging. It fails, however, to explain some of the other experimental observations that have been made of the yielding behavior of iron. For example, it is known that the yield point in iron becomes less pronounced with increasing grain size. Annealed single crystals of iron have very small yield-point elongations .if indeed they have any,' compared to a polycrystalline steel. If the only requirement for a yield point is that the dislocations in the lattice of the annealed. material be anchored by carbon atoms, the difference in the behavior of single crystals and polycrystals is not explained. That a dislocation mechanism may be entirely consistent with little or no yield point in an annealed single crystal will become apparent later when dislocation interaction is discussed. Strain aging produces a definite yield point even in single crystals. This accentuation of the yield-point phenomenon in single crystals after strain

Jan 1, 1953

By Pradeep K. Rohatgi, Clyde M. Adams

Structures produced upon solidification of the eu-tectic composition (33 wt pct Cu) aluminum copper alloy have been examined as a function of freezing rate dfs /d? , the rate of change of fraction solid (fs) with time (8). Slow (dfs/d? = 0.0016 sec-1), intermediate (dfs/d? = 0.02 sec-1) and rapid (dfs/d? = 0.4 to 7.30 sec-1) freezing rates were used. The lamellar Al-Cual2 eutectic is arranged in the form of rod-shaped colonies at rapid freezing rates. The colonies are aligned parallel to the direction of heat flow, whereas the lamellae within the colonies are aligned at various angles, as high as 90 deg, to the direction of heat flow. The colony spacing (C) is proportional to the square root of inverse freezihg rate. The relationship is C = 15.5(dfs/d?)-1/2 where C is in µ and 8 is in sec. The ratio of colony spacing to lamellar spacing is greater than 20.0 and increases with a decrease in the freezing rate. A duplex dendritic structure is produced at intermediate freezing rates. A fine lamellar eutectic is arranged within the dendrites (exhibiting side branches at an angle close to 60 deg from the main stem) and a coarse irregular eutectic appears in the interdendritic regions. The duplex eutectic structure is also produced at slow freezing rates. However, at slow freezing rates there is a Platelat of CuAl2, along the center of the main stem of each dendrite and the other lamellae are arranged perpendicular to the central platelet. THE eutectic between CuA12 and a! aluminum has been reported to freeze in a lamellar form by several workers.'-3 chadwick4 has measured the interlamel-lar spacing as a function of growth rate. Kraft and Albright2 have reported on irregularities in the lamellar structures, and have proposed growth models which account for the formation of faults during solidification. In certain instances the lamellar eutectic has been found to exist in colonies. The colony formation315 has been attributed to the breakdown of a planar liquid-solid interface due to rejection of impurities. The aim of the present work is to study the structures produced from the eutectic aluminum-copper alloy under relatively fast solidification rates, such as encountered in casting and welding operations. The solid-liquid interface presumably remains planar under conditions of slow unidirectional freezing which produce lamellae aligned parallel to the direction of heat flow. The local growth velocities are the same over the entire interface and are equal to the rate of growth of the all-solid region. The spacing between the eutectic lamellae is inversely proportional to the square root of the growth rate of the all-solid region. Under the freezing conditions used in the present study, the solid-liquid interface is cellular or dendritic and the local growth velocities are different in the different regions of the interface. The relationship between the growth rate of the all solid region and the local growth velocities varies with the location and the shape of the interface. The growth rate of the all-solid region is, therefore, an inadequate parameter to describe the eutectic micro-structures which depend upon the local growth velocities. For this reason the structures have been examined as a function of freezing rate, dfs/d?, where fs is the fraction solidified at time 0. The freezing rate was varied by a factor of 4000. The relationship between the freezing rate, dfs/d?, and the growth velocit of the all solid region depends upon the specimen geometry and the shape of the interface. EXPERIMENTAL PROCEDURES The A1-33 pct Cu alloy used throughout this study was made in an induction furnace, using electrolytic copper and aluminum of commercial purity (99.7 pct), the primary impurities being silicon (0.12 pct), iron (0.14 pct), and zinc (0.02 pct). Three ranges of freezing rates were investigated: 1) A spectrum of rapid freezing rates (ranging from 0.40 to 7.30 sec-1) was obtained in arc deposits made on 2-in. thick cast plates of the eutectic alloy. The arc was operated at constant power and was made to travel at constant velocity on the surface of the plate that was in contact with the chill surface during solidification. The pool of liquid metal formed under the moving tungsten arc solidified rapidly by heat extraction through the unmelted plate. Conditions of unidirectional heat flow were achieved near the fusion zone interface, especially in the center of the arc deposits. The great advantage of the arc technique is that rapid cooling and freezing rates can be varied in a qualitative way. The correlation between the arc parameters and the solidification rate is given by the following relationship:6-8

Jan 1, 1970

By R. C. McFarlane, T. D. Mueller, J. E. Walstrom

In evaluating uncertainty, experiments are usually performed repeatedly and then conclusions are drawn from the distribution of results. With the advent of high-speed electronic computers, it is possible to perform experiments using mathematical models constructed to simulate complex experiments or operations. Statistical methods are then applied to the results of the simulated experiments. This procedure forms the busis of this paper. Demonstrated is the need for properly accounting for uncertainty in petroleum engineering problems. How uncertainty affects solutions is evaluated in three example illustrations. The method used to evaluate uncertainty in petroleum engineering studies is the Monte Carlo simulation procedure.'-" INTRODUCTION The solution to most technical problems may be derived from interrelationships among several quantities called variables or parameters. There may be only a few variables or several hundred. Interrelationships among parameters may be explicit or implicit, well established or only approximate. Some variables that fully or partially depend on the magnitude of others are called dependent variables. Input variables for most practical problems are not precisely known; there is usually an uncertainty in their value. The degree of uncertainty may vary from one variable to another. Variables that are known accurately are called determinates.' For instance, the gravity of crude obtained from a particular pool may be known precisely, and therefore is a determinate. The degree of precision with which a quantity can be determined increases as data describing the pool are accumulated during the development of the field and the producing life of the pool. The uncertainty of a parameter may result from difficulty in directly and accurately measuring the quantity. This is particularly true of the physical reservoir parameters which, at best, can only be sampled at various points, and which are subject to errors caused by presence of the borehole and borehole fluid or by changes that occur during the transfer of rock and its fluids to laboratory temperature and pressure conditions. Uncertainty may also result in attempting to predict future parameter values. This type of uncertainty is particularly evident in investment analyses involving future costs, prices, sales volumes and product demand. Uncertainty in the solution to investment problems is often called risk, and its study is called risk analysis.' Uncertainty also enters into biological and sociological analyses in which indeterminate factors are often important due to limited control of the experimental material. It is customary, in evaluating uncertainty, to perform repeated experiments and to draw conclusions from the distribution of the results of these experiments. With the advent of the high-speed electronic computer, it is possible to construct mathematical models which simulate complex experiments or operations and to perform the experiments repeatedly, utilizing the models. Statistical methods are then applied to the results of the simulated experiments This method forms the basis of the investigation reported here. PROBABILITY DISTRIBUTIONS FOR VARIABLES The uncertainty in the value of a variable may be indicated by a probabilistic description accomplished by expressing the quantity by a probability distribution. Many recognized probability distributions can be used to describe physical quantities. Recent studies used various types of distributions to describe core analysis data.',' However, for the examples in this paper, the uniform and triangular distributions are believed to reasonably approximate the data used (Fig. 1). The uniform distribution confines the variable between an upper and a lower limit. The variable may lie anywhere between the two limits. This distribution is used when no one range of values for a variable is more probable than any other, but information or intuitive reasoning indicates the variable will lie somewhere between the chosen limits. The triangular distribution is used for a variable when more data are available to indicate a central tendency of distribution. This allows postulating a "most likely" value to the distribution and upper and lower limits. In this case, as for the uniform distribution, the variable is not expected to assume a value less than the lower limit or greater than the upper limit. However, with improved quality of data it can be postulated that the variable will tend to assume a value close to the most likely value, and that there will be a decreasing probability for values away from the most likely value. The area under either of these probability distributions is equal to unity since it is assumed that there is a 100 percent probability that the variable will lie somewhere under the curve. An ordinate erected at any particular value of the variable divides the area under the curve into two parts: the area to the left of the ordinate represents the probability that the value of the variable will be equal to or less than the value of the variable at the position of the ordinate, and vice versa. The probability is zero that the variable will have any specific deterministic value. If two ordinates are drawn for any two values of the variable, the probability that the variables will have a value lying between these ordinates is equal to the area under the curve lying between the ordinates.

By R. E. Cech, T. D. Tiemann

It has been shown that hydrogen may be made to serve as a rapid and eflicient reducing agent for Cu, Ni, Co, and Fe sulfides if a scavenging agent for hydrogen sulfide is intimately mixed with the sulfide particles being reduced. Accelerated reduction kinetics are demonstrated for nickel sulfide. Copper, nickel, and cobalt sulfides, when treated at certain temperatures in a combined reducing agent-scavenging agent system, are converted to voluminous masses of fibrous metal product. Studies have been carried out to determine the conditions which lead, on the one hand, to irregular poly crystalline fibers and, on the other, to long single crystal filaments a few microns in diameter. A mechanism is proposed to account for the formation of single crystal filuments. The sulfide minerals of Cu, Ni, Co, and Fe are an important source of these metals yet there has been comparatively little scientific effort devoted towards understanding reduction mechanisms of these minerals. This may be, in part, due to the fact that the most convenient reducing agents for carrying out such studies, viz., hydrogen and carbon, do not react appreciably with sulfides. We have found that the reaction of hydrogen with metal sulfides can be markedly accelerated by placing a scavenging agent for hydrogen sulfide in close proximity to the metal sulfide. A brief series of experiments demonstrating relative reduction rates is reported in this paper to illustrate the effect. With the reduction process thus accelerated we have observed an unusual type of reduction behavior on some of the sulfides investigated. Under certain conditions the metallic product of the reduction reaction takes the form of filaments growing outward from the sulfide particles. The present paper deals largely with efforts to classify the various types of growth forms observed. This study has shown that filamentary growths from sulfides take a much greater variety of forms than has heretofore been reported by Ercker,1 Hardy,2 and Nabarro and Jackson3 in their reviews of metallic growths from copper and silver sulfides. THERMODYNAMIC CONSIDERATIONS The thermodynamics for hydrogen reduction of metal sulfides is quite unfavorable. For the sulfides considered here equilibrium constants typically range from 10-3 to 10-5. These low equilibrium constants impose severe kinetic limitations on reduction since hydrogen sulfide must be transported out of the system at concentrations of only a few hundred ppm. Unless extremely high gas flow rates are employed the atmosphere surrounding any sulfide particle will always be essentially in equilibrium with the sulfide. If, however, one places an efficient scavenging agent for hydrogen sulfide in close proximity to the metal sulfide particles the concentration of H2S near the metal sulfide will be held to a very low value. This would permit the reduction reaction to proceed with little or no inhibition from a buildup of reaction product gas. It is well known that calcium oxide is capable of removing hydrogen sulfide from a hydrogen gas stream of low dew point.4 If a sufficient quantity of calcium oxide is mixed with the metal sulfide particles the reaction: CaO+H2S=CaS+ H2O [l] will substitute moisture in place of hydrogen sulfide in the gas stream and this will not affect, in a direct manner, the reaction: MeS +H2=Me + H2S [2] A convenient method of considering the thermodynamics of the combined reducing agent-scavenging agent system is to consider the atmosphere when the partial pressure of hydrogen sulfide is the same over both the metal sulfide and the scavenging agent, i.e., pH2S (1) =pH2S (2). As a consequence: pH2O (1) pH2(2) =K1K2 The chemical driving force for reduction will depend inversely upon the moisture content of the gas and will be 0 when, in the system, pH2O = pH2.K1K2. Table I lists values of the equilibrium constants for reduction and H2S scavenging reactions for a number of sulfides at several temperatures. Data are taken from Rosenqvist4,5 and Kelly.6 The equilibrium constant products calculated from this data show that the limiting level of gaseous reaction product has been increased by a factor of 10' to l04 as a result of substituting a reducing agent-scavenging agent system for a simple reducing agent system. One possible side effect which must be considered is the possibility that the moisture evolved in the scavenging reaction might cause the atmosphere in the system to be sufficiently oxidizing to favor the formation of oxide rather than metal. This possibility was examined by comparing the equilibrium constant products listed in Table I with equilibrium constants for hydrogen reduction of the respective metal oxides. It was found that for copper, nickel, and cobalt the combined reduction-scavenging reactions could not develop a sufficiently high oxidizing potential in the

Jan 1, 1970

By W. L. McMorris, A. H. Brisse

IN the last several years the coal industry has intensified its effort to solve the growing problem of cleaning and recovering fine mesh coals. On one hand these has been increasing civic pressure for cleaner streams, and on the other hand there has been increasing production of fine mesh coal, resulting directly from adoption of the modern mining methods so essential to the economy of the coal mining industry. Cleaning fine coal with the same precision possible with coarser coals is a difficult task, and for coals finer than 200 mesh it has been impractical. Furthermore, the inclusion of —200 mesh material in the final product markedly increases costs of de-watering and thermal drying, which are necessary steps if coal is to meet market requirements. Consequently these extreme fines have generally been wasted. As a result, problems have been created in many districts because there has not been enough area for adequate settling basins. Wasting of coal in the -200 mesh slimes may account for a loss in washer yield equivalent to 2.0 to 2.5 pct of the raw coal input. With rising mining costs the value of such a loss is constantly increasing and a need for a better solution to the fines problem becomes more pressing every day. From an operating viewpoint, also, continuous removal of extreme fines from the washing plant circuit permits good water clarification practice, improving significantly the overall cleaning efficiency. The obvious desirability of recovering a commercially acceptable coal from washery slimes prompted U. S. Steel Corp. to investigate the merits of the Convertol process developed in Germany." Although this process has been used commercially in Europe for some time, little if any consideration has been given to its possible adoption in the U. S. until very recently. Fundamentals of the Convertol Process: In the Convertol process, droplets of dispersed oil are brought into intimate contact with the solids suspended in the coal slurry to be treated. This contact causes oil to displace the water on the surface of the coal by preferential wetting, or phase inversion, after which the coal particles are allowed to agglomerate in a manner permitting their re- moval from the slurry by centrifugal filtration. The clay and other particles of mineral matter suspended in the slurry do not have the affinity for oil the coal particles have. Consequently the oil treatment is preferential to coal to the extent that more than 95 pct of the oil used reports with the clean coal recovered. Figs. 1 through 3 will clarify the steps involved in the process. Fig. 1 shows the suspended material in the slurry to be treated, which is a thickened product containing 40 to 45 pct solids. Oil is now injected into the slurry under vigorous agitation to produce good oil to coal contact conditions, which result in preferential oiling of the coal particles. These coal particles are then permitted to agglomerate by gentle stirring in a conditioner to form flocs, as shown in Fig. 2. At this point in the process the agglomerated oiled coal can be washed and partially dewatered on a vibrating screen, as shown in Fig. 3. Finally, the washed flocculate can be further dewatered in a high-speed screen basket centrifuge or in a solid bowl centrifuge. Commercial Application of the Convertol Process in Germany: The original Convertol process was developed by Bergwerksverband zur Verwertung von Schutzrechten der Kohlentechnik, G.m.b.H., a German research organization controlled by the Coal Operators Assn. of the Ruhr Valley. The process as reduced to commercial practice in Germany' is shown in Fig. 4. In this process a thickened slurry (40 to 45 pct solids) mixed with a predetermined percentage of oil is fed from a surge tank to the phase inversion mill. After the phase inversion step, the slurry is usually discharged directly to a highspeed screen centrifuge. From 3 to 10 pct oil is used, depending on type of oil, size consist of coal to be recovered, and operating temperature. The top size of fine coal cleaned in Germany by the Convertol process is limited by the size of the openings in the centrifuge screen basket. Any mineral matter coarser than the basket opening, which is generally 60 to 80 mesh, must remain with the oiled coal. If the coal fines have been effectively cleaned down to about 80 mesh, the cleaning performance of the process is practically unaffected by the presence of coarse coal particles. However, since recovery of coal much coarser than 80 mesh is mow economical by conventional methods, it normally becomes more costly to allow substantial percentages of this coarse coal in Convertol process feed. Where the general plant layout does not permit effective cleaning of coal sizes down to 80 mesh or lower. there is some justification for a coarser Con-

Jan 1, 1959

By A. T. Chatas

A description of the geometrical characteristics of spherical reservoir systems, a discussion of unsteady-state flow of such systems and examples of engineering applications are presented as backgmund material. The fundamental differential equation, a description of average spherical permeability and the introduction of the Laplace transformation serve as theoretical foundations. Engineering concepts are investigated to indicate particular solutions of interest, which are analytically obtained with the aid of the Laplace transform. These are numerically evaluated by computer, and presented in tabular form. INTRODUCTION A tractable mathematical analysis of unsteady fluid flow through porous media generally requires incorporation of a geometrical symmetry. The simplest forms include the linear, cylindrical (radial) and spherical. Most analytical endeavors have concentrated on cylindrical symmetry because it occurs more often in petroleum reservoirs. Nevertheless, some reservoir systems do exist that are better approximated by spherical geometry. Review of technical literature revealed but a single reference to unsteady spherical flow in petroleum reservoirs.l The motive and purpose of the present work was to remove this gap in technical information, and to provide the practicing engineer with some useful analytical tools. The mathematical details associated with the particular solutions of interest involved use of the Laplace transformation. Hurst and van Everdingen previously demonstrated the efficacy of this operational technique, and in many respects the present treatment was patterned after their earlier work.2 PRELIMINARY CONSIDERATIONS GEOMETRICAL CHARACTERISTICS Geometrically, a spherical reservoir system is defined at any instant of time by two concentric hemispheres whose physical properties of interest vary only with the radial distance. Every physical property is thus restricted to be a space function of only one variable: the distance along a radius vector emanating from the center. Such a system is composed of an outer region and an inner region, separated by a defined internal boundary. The inner region simply extends inward from this boundary, whereas the outer region extends outward from it to an external boundary. The position of the internal boundary is presumed fixed, so that the size of the inner region remains constant. On the other hand, the position of the external boundary at any given instant of time is determined by the distance into the system that a sensible pressure reaction has occurred, Thus, the external boundary may change position with time. It initially emerges from the inner region and advances outward to its ultimate position. When this ultimate position coincides with a geometric limit, the reservoir system is said to be limited. When it coincides with points subject to pressure gradients furthest removed from the internal boundary, yet short of a geometric limit, the system is said to be unlimited. In this investigation two different boundary conditions are imposed at the ultimate boundaries of limited systems. The first requires that no fluid flow occur across this boundary; the second that the pressure remain fixed at this boundary.3-5 UNSTEADY-STATE FLOW In a strict sense virtually all flow phenomena associated with a reservoir system are unsteady-state. The transient behavior of these phenomena requires accounting, however, only when time must be introduced as an explicit variable. Otherwise, steady-state mechanics may be used. Analytically, steady-state conditions prevail in a reservoir system only over that portion of its history when this relation is satisfied:

Jan 1, 1967

By Ralph B. Stevens

INTRODUCTION Underground fires continue to be one of the most serious hazards to life and property in the mining industry. Although underground mines are analogous to high-rise buildings where persons are isolated from immediate escape or rescue, application of technology to locate and control fire hazards while still in their controllable state is slow to be implemented in underground mines. Even in large surface structures such as hotels, often only fire protection systems which meet minimal laws are implemented due to the high cost of adding extensive extinguishing systems, isolation barriers, alternate ventilation, escape routes and alarm systems. Incomplete and ineffective protection occasionally is evidenced where costs would not seem to be a factor, such as the $211 million MGM Grand Hotel fire November 21, 19801. Paramount in increasing fire safety and decreasing the threat of serious fire is early warning followed by proper decision analysis to perform the correct action. However, very complex fire situations can be produced in structures such as high-rise buildings and underground mines simply because of the distances between the numerous fire-potential locations and fire safe areas. Other complexities arise when normal activities occur that emit products of combustion signaling a fire condition to a sensitive fire/smoke sensor. For example, the operation of diesel equipment or the performance of regular blasting can produce combustion products that reach the sensitive alarm points of many sensors2. Smoke detectors for surface installations provide fire warning when occupants are at a distant location or when sleeping, thus greatly reducing injuries and property damage. However, when installed in the harsh environments of underground mines, fire and smoke detection equipment soon becomes inoperative, unreliable, or requires excessive maintenance. The U.S. Bureau of Mines has performed many studies and tests to improve fire and smoke protection for underground mine workers3. This paper describes several USBM safety programs which included in-mine testing with mine fire and smoke sensors, telemetry and instrumentation to develop recommendations for improving mine fire safety. It is hoped that the technology developed during these programs can be added to other programs to provide the mining industry with the necessary fire safety facts. By recognizing fire potentials and being provided with cost-effective, proven components that will perform reliably under the poor environmental conditions of mining, mine operators can provide protection for their working life and property equal to that which they provide for themselves and their families at home. The basis of this report is two USBM programs for fire protection in metal and nonmetal mines4,5 and one coal program6. The data was collected beginning in May 1974 and continuing through the present with underground tests of a South African fire system installed at Magma Mine in Superior, Arizona, and a computer-assisted, experimental system at Peabody Coal Mine in Pawnee, Illinois. The conduct of each program was as follows: • Define the problem and its magnitude in the industry • Develop concepts to solve or diminish the problem • Review available hardware or systems approaches to fit the concepts • Install and demonstrate the performance of a prototype system through fire tests in an operating mine. MINE FIRE FACTS Whether in coal or metal and nonmetal mines, the potential severity of fire hazard is directly related to location. As shown in Figure 1, fire in intake air at zones A, B, C or D can cause contamined air to route throughout the mine quickly if not detected, isolated or rerouted. Causes and location of former metal and nonmetal fires are represented in Table 1; the cause and location of fatalities and injuries is shown in Table 2. Coal-related fires and their impact on deaths and injuries are graphed in Figure 2; their locations are described in Table 37. Significantly the table shows that the hazard to personnel was three times greater for fires occurring in shaft or slope areas, and the percentage of deaths and injuries was four times that of other areas. Number of Persons Affected A 129-mine sample indicated that from 8 to 479 employees per shift work in underground metal and nonmetal mines, and that deeper mines have larger populations, as shown in Figure 3. Coal mining relates similar employment, and a 16-state sample of 670 mines employing at least 25 persons shows the distribution in Figure 4. Drift mines accounted for 58 percent of the sample but employ only 45 percent of the underground workers.

Jan 1, 1982

By C. S. Matthews, H. C. Lefkovits

Drilling progress is often delayed by sticking of the drill string. The development of preventive and remedial methods has been hampered by incomplete understanding of the sticking mechanism. A recent lahorntory investigation hns indicated that one type of sticking may be attributed to the difference in pressure between the borehole and formation. This paper shows, by means of soil mechanics, that the primary cause for differential pressure sticking is cessation of pipe movement, whereas diflerential pressre and stanrtding time determine the severity of the sticking. The analysis stresses the importance of using low-weight muds with low solids content and low water loss to alleviate diflerential pressure sticking and describes why packed hole drilling, long strings of drill collars, and a large deviation from the vertical are conducive to sticking. Finally, preventrve and remedial methods ore evaluated, and a theory is presented on the release of stuck pipe by spotting oil. INTRODUCTION Since drilling with long strings of oversize drill collars has become standard practice in many areas, the incidence and severity of the stuck pipe problem has increased. It has been noticed that in the majority of these cases the sticking could not possibly be attributed to key seating or caving of shales. It appeared that, due to the differential pressure between the mud column and the formation fluid, the collars were pressed into the wall and so became "wall stuck". Points to note about differential pressure sticking are: (1) sticking is restricted to the drill collars, (2) the collars become stuck opposite a permeable formation, (3) the sticking occurs after an interruption of pipe movement, (4) circulation, if interrupted, can be restarted after the sticking is noticed, and (5) no large amounts of cuttings are circulated out after restarting circulation. Helmick and Longleyl investigated pipe sticking by differential pressure in the laboratory and found an empirical relationship between the differential pressure, the sticking time and the required pull-out force. In this paper an explanation of the mechanism is given based on Terzaghi's theory of clay consolidation. A qualitative description is given in the following paragraphs while the derivation of fonnulas is given in Appendices. This paper is a first attempt to explain pressure differential sticking and many points will require additional theoretical and practical investigation before the problem can be fully understood. PRESSURE DIFFERENTIAL STICKING AS A CONSOLIDATION PROBLEM In any borehole, where the mud pressure is higher than that exerted by the formation fluids, a mud cake is formed opposite the permeable sections of the hole and a continuous flow of filtrate takes place from the mud, through the cake and into the formation. This radial flow pattern requires a certain distribution of the hydraulic and the effective (grain-to-grain) stresses inside the mud cake. Any quantitative or qualitative change in the external pressure conditions will produce a change in the flow pattern and, consequently, also in the internal stress distribution inside the cake. In view of the low permeability and the high compressibility of a clay mud cake, the adjustment of the internal stress distribution is slow and is accompanied by a change in volume. Time dependent stresses are thus created which gradually diminish as the new state of equilibrium between internal and external pressures is approached. Some 30 years ago, Terzaghi developed his "Theory of Consolidation" to account for the time-dependent stresses and settling of clay formations under the influence of external loads. He derived a differential equation by which the time-dependent hydraulic stress and the consolidation can be computed for any point inside the layer during the consolidation process. His theory is based on the assumption that the change in stress is solely due to a change in water content and it may only be applied to one-dimensional consolidation phenomena. Other investiga-tors5,10 have expanded his theory to include processes of more than one dimension. The difference between the external pressures on the mud cake before and after sticking is a qualitative one (isolation of part of the cake by the static contact with the drill collars after pipe movement has been stopped)', and the time-dependent stresses thus created may be investigated by means of Terzaghi's theory. By this analysis the changes in the nature of the contact surface between the drill collars and the mud cake during the sticking can be explained; and the friction force between the two may be computed as a function of the sticking time, the borehole dimensions and the mud cake characteristics.

By R. S. Wagner, W. C. Ellis

A new mechanism of crystal growth involving oapor, liquid, crnd solid phases explains many observations of the effect of implurities in crystal growth from the vapor. The role of the impuuitq is to form a liquid Solution with the crystalline tnalerial to be grown from the vapor. Since the solution is n prefevred site for deposition firorti the uapor, the liquid becorrles supersaturated. Crystal growth occurs by precipitatzon from the supersaturated liquid crt tlie solid-liquid zntevfnce. A crystalline defect, such as a screw dislocation, is not essetztial for VLS (vapor -liquid-solid) growth. The concept of the VLS mechanism is discussed in detail with reference to tire controlled growth of silicon crystals using gold, platinum, palladium, nickel, silver, or copper as an implurity agent. RECENTLY a short communication' described a new concept of crystal growth from the vapor, the VLS mechanism. In this paper we present a detailed description of the process and its application to the growth of silicon crystals and we discuss its relevance to existing concepts of .'whisker" crystal growth. Crystal growth from the vapor is usually explained by a theory proposed by Frank2 and developed in detail by Burton, Cabrera, and Frank.3 In this theory a screw dislocation terminating at the growth surface provides a self-perpetuating step. Accommodation of atoms at the step is energetically favorable, and is possible of much lower supersatu-ration than required for two-dimensional nucleation. Crystals of a unique form resulting from aniso-tropic growth from the vapor are "whisker" or filamentary ones. Such crystals have a lengthwise dimension orders of magnitude larger than those of the cross section. For most filamentary crystals both the fast-growth direction and directions of lateral growth have small Miller indices. The special growth form for a whisker crystal implies that the tip surface of the crystal must be a preferred growth site. sears4 proposed that, according to the Frank theory. a whisker contains a screw dislocation emergent at the growing tip. Such an axial defect provides a preferred growth site and accounts for unidirectional growth. The hypothesis was extended by Price. Vermilyea. and Webb," still implying the presence of a dislocation at the whisker tip. They postulated that impurities arriving at the fast-growing tip face become buried while those arriving on the surface of slow-growing lateral faces accumulate and thereby hinder growth. These considerations led to a whisker morphology. There is increasing evidence that most whisker crystals grown from the vapor are dislocation-free. Webb and his coworkers6 searched for an Eshelby twist7 in zinc? cadmium, iron. copper, silver, and palladium whisker crystals. They found unequivocal evidence for an axial screw dislocation in only one element, palladium. However, not every palladium crystal examined contained a dislocation. Observations with the electron microscope have failed to show dislocations in whisker crystals of zinc, silicon.9 and one morphology of AlN.10 Since many whiskers are completely free of dislocations, an axial dislocation does not appear to be required for whisker growth of many substances. A significant advance in understanding whisker growth has been a recognition of the need for impurities. This requirement has been clearly demonstrated for copper,11 iron,13 and silicon9-1 whiskers. For silicon, detailed studies proved conclusively that certain impurities, for example, nickel or gold, are essential. Another pertinent phenomenon which has received little attention is the presence of a liquid layer or droplets on the surface of some crystals growing from the vapor. Crystals in which this has been observed include p-toluidine,14 MoO3,15 ferrites,16 and silicon carbide.'" The liquid layers or globules were considered to be metastable phases, molecular complexes, or intermediate polymers originating from condensation of the vapor phase. The possibility has been suggested that the halide being reduced is condensed at the tip18 or adsorbed on the surface11 of a growing metal whisker, for example copper. The literature on whiskers discloses illustrations of rounded terminations at the tips. These appear. for example, on crystals of A12O3,19,20 sic,21 and BeO.22 For BeO, Edwards and Happel suggested that during growth of the whisker the rounded termination consisted of molten beryllium enclosed in a solid shell of BeO. A recent paper9 on the growth of silicon whiskers contains many observations pertinent to an understanding of the mechanisnl of whisker growth. These observations are summarized as follows. 1) Silicon whiskers are dislocation-free. 2) Certain impurities are essential for whisker growth. Without such impurities the silicon deposit is in the form of a film or consists of discrete polyhedral crystals.

Jan 1, 1965

By David W. James, Howard Jones, George M. Leak

Conditions are described for producing, by primary recrystallization, a matrix suitable for the growth of large grains near (110)[001] orientation in silicon iron strip by secondary recrystallizaliun in a steep temperature gradient. The orientation distribution of these large grains is expressed in terms of rotational deviations about the cross-rolling direction, the rolling direction, and the normal to the sheet, the deviational spread increasing in that order. With the aid of cowplenientary published data on the orientation dependence of growth rate, it is shown that this observation is consistent with the oriented-growth theory of recrystallization lextures. It is conclutled that growth-rate and orientation-distribution data obtained in a steep thermal gradient should be used with caution to account for isothermally Produced recrystallization textures. SEVERAL authors have reported methods of growing large grains by re crystallization of a small-grained matrix in silicon iron 1- B and pure a cr The present study was a preliminary in the growth of single crystals and bicrystals for surface relaxation," grain boundary mobility, and grain boundary diffusion studies. The method was to control the growth of a seed crystal into a suitable primary re crystallized matrix by feeding through a steep temperature gradient. The driving energy for growth derived from the grain boundary energy released as the seed crystals grew into the matrix. Thus, stability of the matrix against normal grain growth was considered to be essential for success. It was known that the manganese sulfide dispersion present in commercial silicon iron performs this function during secondary recrystallization to the (110)[001.] texture.12 Hence commercial, rather than high-purity, material was used throughout. The paper describes the growth conditions for grains large enough to be used as seed crystals for further growth into single crystals. The orientation distribution of the seed crystals is analyzed and its significance for the theory of recrystallization textures is discussed. EXPERIMENTAL PROCEDURE Strip material was supplied by the Steel Co. of Wales, Ltd. The chemical analysis in weight percent was Si, 2.90; C, 0.015; Mn, 0.059; P, 0.011; S, 0.027; Ni, 0.032; 0, 0.009; Fe, balance. A gradient furnace of similar design to one described previously4 was loaned from B.I.S.R.A. It consisted essentially of a vertical water-cooled copper slot projecting downwards into the hot zone of a molybdenum furnace. Hydrogen was passed through the furnace to protect both heating element and specimen from oxidation. Strip specimens up to 8 cm wide and 0.2 cm thick were sealed into the furnace at the mouth of the copper slot. A coating of light oil on the strip surface maintained the seal during translation of a specimen. The maximum temperature gradient in the region just below the copper slot was 500°C per cm over 1 cm, with the hottest point controlled at 1175°C. Several large grains would usually grow by secondary recrystallization from the primary matrix when a specimen was immersed in the hot zone for about 30 min. A back-reflection X-ray camera was constructed to facilitate rapid and accurate orientation determinations of the large grains produced. It was possible to reproduce a standard geometry, with regard to strip and camera, without the tedium of careful alignment on each occasion. Specimens, typically 4 cm wide and 75 cm long, were cut with the longitudinal axis parallel to the rolling direction of the original strip. The surfaces were cleaned by immersion alternately in a hot aqueous solution containing 2 pct hydrofluoric acid plus 10 pct sulfuric acid and in cold 10 pct nitric acid. The nitric acid etch was just sufficient to reveal the grain structure. Rolling and annealing treatments to prepare the matrix (discussed below) were followed by growth of seed crystals in the gradient furnace. The matrix was transformed to a single crystal by growth of a selected seed crystal connected to the matrix by a thin neck. 4,5 Growth was promoted by controlled feeding into the gradient furnace. Several single crystals of controlled orientation were grown successfully from seed crystals by twisting the interconnecting neck in a reorien-tation jig.4 EXPERIMENTAL RESULTS AND DISCUSSION Growth Conditions. A suitable matrix for growth of large grains was prepared starting from primary re-crystallized strip 1.9 mm thick. This was cold-rolled in two stages each being followed by a recrystallization anneal at 800°C for a few minutes. Such treatment gave the required growth matrix only if the two cold-reduction stages were each performed in several passes and in the following ranges: the first, 30 to 70 pct; the second, 10 to 50 pct. Immersion in the temperature gradient otherwise resulted in an equiaxed

Jan 1, 1967