Search Documents

By I. Javandel, P. A. Witherspoon

The finite element method was originally developed in the aircraft industry to handle problems of stress distribution in complex airframe configurations. This paper describes how the method can be extended to problems of transient flow in porous media. In this approach, the continuum is replaced by a system of finite elements. By employing the variational principle, one can obtain time dependent solutions for the potential at every point in the system by minimizing a potential energy functional. The theory of the method is reviewed. To demonstrate its validity, nonsteady-state results obtained by the finite clement method are compared with those of typical boundary value problems for which rigorous analytical solutions are available. To demonstrate the power of this approach, solutions for the more complex problem of transient flow in layered systems with crossflow are also presented. The generality of this approach with respect to arbitrary boundary conditions and changes in rock properties provides a new method of handling problems of fluid flow in complex systems. INTRODUCTION Problems of transient flow in porous media often can be handled by the methods of analytical mathematics as long as the geometry or properties of the flow system do not become too complex. When the analytical approach becomes intractable, it is customary to resort to numerical methods, and a great variety of problems have been handled in this manner. One such method relies on the finite difference approach wherein the system is divided into a network of elements, and a finite difference equation for the flow into and out of each element is developed. The solution of the resulting set of equations usually requires a high speed computer. When heterogeneous systems of arbitrary geometry must be considered, however, this approach is sometimes difficult to apply and may require large amounts of computer time. The finite element method is a new approach that avoids these difficulties. It was developed originally in the aircraft industry to provide a refined solution for stress distributions in extremely complex airframe configurations. 27 Clough has recently reviewed the application of the finite element method in the field of structural mechanics.6 The technique has been applied successfully in the stress analysis of many complex structures.l, 27, 28 Recognition that this procedure can be interpreted in terms of variational procedures involving minimizing a potential energy functional7 leads naturally to its extension to other boundary value problems. In the field of heat flow, there recently have been introduced several approximate methods of solution that are based on variational principles.2-4, 17 By employing the variational principle in conjunction with the finite element idealization, a powerful solution technique is now available for determining the potential distribution within complex bodies of arbitrary geometry. In the finite element approximation of solids, the continuum is replaced by a system of elements. An approximate solution for the potential field within each element is assumed, and flux equilibrium equations are developed at a discrete number of points within the network of finite elements. For the case of steady-state heat flow, the technique is completely described by Zienkiewicz and Cheung.33 Since the flow of fluids in porous media is analogous to the flow of heat, Zienkiewicz et al. have employed the finite element method in obtaining steady-state solutions to heterogeneous and anisotropic seepage problems. 34 Taylor and Brown have used this method to investigate steady-state flow problems involving a free surface.25 The work of Gurtin has been instrumental in laying the groundwork for the application of finite element methods to linear initial-value problems.12 As a result, Wilson and Nickell have recently

Jan 1, 1969

By G. A. Stamboltzis, N. Arbiter, C. C. Harris

Fracture under low-velocity free-fall and double impact and under slow compression have been investigated. The pattern of breakage and the size distribution of resulting fragments of sand-cement and glass spheres have been determined. Photoelasticity methods were used to simulate the stress distributions in free-fall impact in order to explain the observed patterns of breakage. Oblique fracture planes, occurring only in free-fall impact, develop along lines which coincide with the trajectories of maximum compression as determined by mathematical analysis. Breakage efficiencies for different modes of fracture were compared for both types of spheres. For the same specimen and loading system, static loading and low-velocity dynamic loading induce geometrically similar stress fields resulting in reasonably similar fracture patterns and shapes of fragments. This work had its origin in a desire to understand three fundamental processes involved in autogenous comminution: namely; free-fall impact, double impact, and slow compression. To simplify the investigation, it is preferable to reduce the entire multi-stage fracture process to its most elementary form; that is, to study products resulting from a single-stage operation, such as the breakage of a specimen under single fracture conditions. The goal of this research is to study the energy utilization in the single fracture of brittle specimens and to relate the pattern of breakage and the resultant fragment size distribution with the nature of the material, the specimen size, the manner of load application, and the rate of loading. Published works on single fracture have mainly been concerned with tests on large (> 1-in.) irregular mineral pieces, especially coal and coke. These were mainly friability tests employing a single blow on closely sized pieces. More recently, small mineral specimens in the sieve range and glass spheres have been investigated.5-7 Kick8 has reported free-fall and double impact tests with large large cast-iron, cement, and clay spheres. The establishment of a minimum breaking height independent of size in the free-fall impact tests, and the direct proportionality between the minimum work required for fracture and the volume of the specimen in double impact tests were used by Kick as an experimental proof of his classical law of "Proportional Resistances." In the present work, spherical shapes were chosen because of their simple body geometry and consequent impact and stress field symmetry. Because the study involves several physical principles in connection with brittle-fracture, it may be of interest in fields where the strength of materials is of importance. APPARATUS The equipment for the free-fall impact testing is illustrated in Figs. 1 and 2. This consisted of a massive (17 x 17 x 3-in.) hard steel plate onto which specimens were dropped and a specimen release mechanism which could be set at any desired height up to 10 ft above the plate so that the impact velocity could be varied up to about 25 ft per sec. Double impact breakage is characterized by two points of loading situated at opposite poles. For this mode of breakage the apparatus was modified (Figs. 3, 4 and 5) so that a falling mass provided an impact of predetermined magnitude on a specimen resting on the plate. A spring-loaded device operated immediately after fracture to arrest the falling mass and to record its residual kinetic energy. Precautions were taken to avoid secondary breakage. Slow compression tests were performed on a conventional hydraulic testing machine, the specimen being held between two hardened parallel bearing blocks. In the case of glass spheres, carbide bits were used. Photoelasticity studies were performed on two-dimensional models held in a loading frame equipped with a force gauge (Fig. 6). The models were viewed in a polariscope and isochromatic fringe patterns were recorded photographically. SPECIMENS Sand-Cement Spheres: A series of sand-cement spheres were prepared by molding in spherical glass flasks.

Jan 1, 1970

By C. Laird, C. E. Feltner

An investigation into the effect of a material's slip character on its high strain fatigue properties has been carried out using copper and Cu-7.5 wt pct A1 as representative wavy and planar slip mode materials, respectively. In copper, the constant term and exponent in the Coffin-Manson equation, ?EpNf1/2 = const, are unaffected by the temperature of testing or the history of the metal. On the other hand, in Cu-7.5 wl pcl Al, while the exponent is invariant, the constant is increased with respect to copper and increased further with decrease in temperature. The life in Cu-7.5 wt pct A1 is double that of copper at room temperature and about four times at 78°K. Fractoflaphic studies indicate that these diffevences in life are associated with the early stages of cracking. In copper, severe surface rumpling provides crack nuclei at pain boundaries independent of' testing temperature or the history of the metal. In Cu-7.5 wt pct Al, the planar slip mode of deformation prevents grain boundary folding; crack nucleation is therefore mainly transgranular and the increase in life with respect to copper is associated with a prolonged linking-up stage of such transgranulay cracks. At the lower temperature the increased difficulty of cross slip inhibits multiplication of trans-granular cracks and retards the linking-up process, giving a further enhaancement in life. It is expected that such behavior will be found general in materials of wavy and planar slip mode, respectively. STAGE II crack propagation in high strain fatigue takes place by a repetitive process of crack tip blunting in the tension part of a fatigue cycle followed by resharpening of the crack in the compression part.1-3 This mechanism has come to be called the "plastic relaxation process".3 At very high strains (lives to failure less than 103 cycles), most of the life of a specimen is spent in crack propagation by the plastic relaxation process.4,5 Since the differences in the cyclic strain hardening behavior of metals and alloys are comparatively small in this strain range,6,7 the amount of blunting at a crack tip depends primarily on the applied strain amplitude. The rate of crack propagation in this life range is therefore nearly constant with respect to material and temperature of testing. In the longer life portion of the high strain fatigue range (l03 to 103 cycles to failure), a considerable fraction of the lives is spent in crack initiation and Stage I growth. Grain boundaries have frequently4,5,8 been found to be the sites for this growth which is also associated with a surface folding process.4'5 Such folding behavior may also be general in any material having a grain structure and capable of plastic deformation. It has therefore been suggested475 that the universality of the Coff in-Manson law,9'10 constant,* where ?Ep is the plastic strain range and Nf the number of cycles to failure, is reflected in the approximate invariance of both the crack initiation and propagation behavior. It happens that, in most investigations of the phenomena associated with failure,4,5,8 materials with a wavy slip mode have been used as vehicles for the studies, or else "normally" planar slip mode materials have been tested under conditions giving rise to wavy glide deformation, e.g., stainless steels cycled at high temperature.'' The purpose of the present investigation was to determine whether or not different cracking phenomena would occur when the slip mode of a material is planar and if the associated life to failure would be different than that of a wavy slip

Jan 1, 1968

By D. H. Polonis, Gary A. Dreyer

This investigation involved a study of the reactions occurring during aging of supersaturated a phase in a CIL-Si alloy. The aging processes at temperatures below 552°C were studied by means of metallography, microhardness, and electrical resistivity measurements. During the initial stage of aging in the temperature range 275o to 525oC an anisotropic phase, which is visible in the form of striations in the micro-structure formed prior to the appearance of the -phase. This transition phase has been designated k' due to its structural similarity to the k phase. The rate of K' formation during aging increased with aging temperature. At 450o and 525°C the reaction appeared to reach completion before significant amounts of precipitation were observed. The overall precipitation reaction in this system involved four stages: formation of formation of at a grain boundaries, growth of particles in 7-egions, and coagulation of into large particles. The relative role of each of these processes is dependent on the temperature of aging. The formation of phase from supersaturated a during aging increases the electrical resistivity and slightly increases the hardness of a 4.9 pct silicon alloy. The formation of during aging decreases the electrical resistivity and increases the hardness until coalescence of particles begins during the final stages of aging. The coalescence of results in a decrease of both electvical resistivity and hardness. THE available literature does not record a detailed analysis of the aging processes in Cu-Si alloys. Such an analysis has been difficult due to both the limited data available and misleading microstructural observations which have been reported. Previous work suggests that a rather complicated series of processes may be involved in the aging of supersaturated a. The present investigation is concerned with an analysis of the sequence of these aging processes in a 4.9 pct Si alloy. The currently accepted version of the Cu-Si equilibrium diagram, Fig. 1, was published by C. S. Smith in 1940.' This diagram shows a eutectoid at 5.2 pct Si with the high temperature phase, having its composition range entirely to the right of this value. By quenching alloys in the composition range 4.6 to 5.2 pct Si from the region, Smith was able to retain the phase. However, this phase exhibited a striated microstructure which was anisotropic to polarized light. These striations were suggested by Smith to be duplex twins. Using an oscillating crystal X-ray method Barrett2 studied alloys containing 4.0 to 5.4 pct Si quenched from the and regions. The striations observed in quenched specimens were attributed by Barrett to quenching stresses. He found that striations could also be formed by deformation after quenching and identified them as stacking faults in alloys containing 4.0 to 5.0 pct Si. Barrett related the striations in alloys containing 5.0 to 5.4 pct Si both to stacking faults and to a metastable phase which he called The striated structure reported by smith' was also obtained by Hoffmann, et al. 3 By deeply etching their specimens the anisotropic phase on the surface was removed, leaving a structure which was microscopically similar to pure copper. By means of X-ray analysis before and after deformation at room temperature it was concluded that the anisotropic phase was metastable This conclusion has been questioned by Barrett who attributes the discrepancy to superposition of k and lines on the powder patterns of Hoffmann, et alm3 The decomposition of supersaturated according to smith' and Anderson,4 is very sluggish. However, Nony and Dreyer,5 and Hoffmann, et al.' found that deformation increases the rate of decomposition. X-ray studies by Hoffmann, et al on

Jan 1, 1962

By N. A. Warner

Dense cylindrical specimens of artificial hematite were reduced in hydrogen over a range 0-f total pressures between 0.1 and 1.0 atm and temperatures between 650" and 950°C. Hydrogen reduction at a total pressure of 1.0 atm in the presence of a chemically inert diluent and reduction in hydrogen/water mixtures were also studied Znter.face movement during reduction was followed using metallographic techniques. Substantiu1 partial-pressure gradients are shown to exist in both the gas film and the reduced iron layer. The results are consistent with a mixed-control mechanism involving an interaction of gaseous-difusion effects with a first-order reversible chemical reaction at the iron/wüstite interface. Over the hydrogen-pressure range of 0.1 to 1.0 atm, the rate is not directly proportional to the hydrogen pressure in the bulk gas stream. Nitrogen dilution of the reducing gas masks the true effect of hydrogen partial pressure by influencing the effective molecular dif-fusivity and gives an apparent first-order relationship. The transport of hydrogen and water vapor across the iron layer is consistent with molecular gaseous diffusion rather than a Knudsen diffusion mechanism. ALTHOUGH the literature on the reduction of iron oxides is very extensive, it is only comparatively recently that a definite quantitative model for the reduction of hematite has been proposed. This development has been largely due to the work of McKewan. An extensive literature review of previous work has recently been compiled by Themelis and Gauvin. Exclusive control by gaseous diffusion in the reduction of hematite has been recently advocated by Kawasaki and coworkers.' Although the results of the latter workers are of considerable practical significance, the extreme porosity (approximately 30 pct) of their specimens precludes their application to a kinetic evaluation of the reduction of dense polycrystalline hematite. The salient features of the kinetic model used by McKewan for the gaseous reduction of dense hematite at temperatures in excess of 560°C (when FeO is stable) are: 1) The reduction is a topochemical reaction taking place through 'the series FeO/FesO4/FeO/ Fe; 2) The gas-solid type of reaction takes place only at the FeO/Fe interface and the internal reduction: FeO proceeds by solid state diffusion; 3) The iron layer formed is quite porous and offers negligible impedance to gaseous diffusion; 4) The rate-controlling step in the reduction is a chemical process occurring at the Fe/FeO interface. In view of the relatively fast reaction rates involved when hematite is reduced in a stream of hydrogen, it seemed possible that the simplifying assumptions of negligible diffusional resistance in both the bulk gas phase and the reduced iron layer could introduce serious errors in the interpretation of experimental reduction data. The present work was initiated to elucidate the effect of gaseous diffusion on the reduction rate and to determine its possible influence in the formulation of the reaction mechanism at the metal/oxide interface. EXPERIMENTAL Procedure and Apparatus. Dense cylindrical compacts of artificial hematite were prepared by sintering pressed compacts of reagent-grade ferric oxide, containing 99.6 pct FezO, in air at 1200°C for 80 hr. The resulting compacts had a final porosity of less than 2 pct and a height and diameter of approximately 1 cm. Variations in sintering atmosphere from air to pure oxygen, in sintering time from 20 to 80 hr, and temperature from 1150" to 1250°C had no effect on the subsequent reduction behavior. The reduction of the samples took place in a vertical tube furnace in which was placed a 1-in. transparent silica tube. The 24-in. heated length

Jan 1, 1964

By M. H. Caron

THE most outstanding property of ammonia liquors, used in the ammonia leaching process is their very limited ability to dissolve all compounds present in reduced ore except nickel and cobalt. Although they do have such a high degree of purity, even the smallest amounts of impurities contaminate the precipitate of basic nickel carbonate that may be obtained from pregnant liquors by straight distillation. Figures have been given in a separate paper regarding the extent of such contamination. The product of straight distillation is likely to contain small amounts of Fe, Mn, MgO and SiO2 as well as some sulphate and sulphide. The sulphur can be readily eliminated by calcination of the product at the proper temperature and MgO and SiO2 by reduction and melting of the metal with a suitable flux. After such treatment the metal is relatively pure, and the iron and manganese content will be very small. The cobalt content, however, may vary considerably, depending upon the ore treated. The metal obtained from pure true garni-erite ores usually has only a small cobalt content of less than 1 pct. Metal obtained from medium grade iron-rich nickel ores may show cobalt values from 2 to 3 pct, and the alloy obtained from laterite is likely to contain as much as 10 pct cobalt. By proper treatment a good deal of cobalt may be recovered from lateritic nickeliferous iron ores, and cobalt will become the major constituent of the pregnant liquor if asbolan ores or similar products are treated. The pregnant liquor of "run of mine" nickel ores may have a pinkish shade, indicating the presence of cobalt. It is evident that from such liquors the metals should be recovered separately. Even if cobalt is present in only very small amounts, it would probably be economically justifiable to recover both metals as pure as possible if a suitable process could be devised. The author has made investigations along this line and has patented some processes adapted to the treatment of this type of solutions. Selective Dissolution of Nickel from Basic Nickel Carbonate Containing Some Cobalt: Basic nickel carbonate obtained from true garnierite ores by straight distillation of the pregnant liquors, always contains some cobalt as well as the normal impurities. In. order to improve the quality of the nickel product and to recover most of the cobalt, the following process has been devised. By using liquors with a low CO2 to NH3 ratio, nickel can be redis-solved preferentially to cobalt. It is not a sharp and clean separation but a selective one, for the bulk of the nickel may be redissolved with a small amount of cobalt, and a small residue may hold the bulk of the cobalt. The greatly increased cobalt content of the residue is indicated by a pronounced olive green color, in contrast with the light green color of pure basic nickel carbonate. The following example is given from Dutch patent No. 52788, Klasse 40a. 43 granted June 17, 1942: The original precipitate contained about 0.55 pct cobalt calculated on total nickel + cobalt. Eight hundred and fifty grams of the basic nickel carbonate were redissolved in leach liquor containing 7.12 pct NH3 and 2.12 pct CO2. The weight of the dry residue was about 3.5 pct of the original weight. By analysis it was found that the residue contained: Fe 1.93 pct, Co 6.01 pct, Ni 34.87 pct and also a certain amount of silica. The clear filtered blue liquor was distilled and the precipitate of basic nickel carbonate from this analyzed: Fe 0.008 pct, Co 0.0146 pct, Ni 42.37 pct. It is obvious that this product is considerably more pure, since calculation shows that about 94 pct of the original cobalt remained in the residue. Cobalt and nickel can be recovered by further

Jan 1, 1951

By David E. Nicholas

INTRODUCTION In this paper, a numerical process for selecting a mining method, with the emphasis on underground mass mining techniques, such as caving, induced caving, and stoping, is proposed. In the past, selection of a mining method for a new property was based primarily on operating experience at similar type deposits and on methods already in use in the district of the deposit. Then, the chosen method was modified during the early years of mining as ground conditions and ore character were bet- ter understood. Today, however, the large capital investment required to open a new mine or change an existing mining system make it imperative that the mining methods examined during the feasibility studies and the method actually selected have a high probability of attaining the projected production rates. Although experience and engineering judgment still provide major input into the selection of a mining method, subtle differences in the characteristics of each deposit, which may affect the method chosen or the mine design, can usually be perceived only through analysis of measured characteristics. The parameters that must be examined when choosing a mining method include: 1) geometry and grade distribution of the deposit; 2) rock mass strength for the ore zone, the hanging wall, and the footwall; 3) mining costs and capitalization requirements; 4) mining rate; 5) type and availability of labor; 6) environmental concerns; and 7) other site-specific considerations. This paper encompasses a detailed look at the first two parameters since they, plus mining costs, have the greatest impact on the selection of a mining method. The proposed method selection process is for a project where drilling has defined sufficient geologic reserves, but little or no underground development has been done. Since each deposit has its own characteristic geanetry/grade distribution, and rock mechanics properties, mining method selection should be at least a two-stage process. In Stage 1, the deposit is described in terms of geometry, grade distribution, and rock mechanics properties. Using these parameters, the mining methods can be ranked to determine which are most applicable; they can then be considered in general terms of mining and capitalization cost, mining rate, type and availability of personnel, environmental concerns, and other site-specific considerations. In Stage 2, the most likely mining methods are costed out, based on a general mine plan. Mining and capitalization costs are used to determine a cut-off grade fran which a minable reserve can be calculated; economic comparisons can then be made to determine the optimum mining method and economic feasibility. During the mine planning phase of Stage 2, rock mechanics information would be used to provide realistic estimates of underground opening size, amount of support, orientation of openings, and caving characteristics, and open pit slope angles. If ground control or operational problems should be encountered with the methods being considered, modifications could be made. Although planning on paper extends start-up time, it is cheaper to err on paper than to find the error after mining has begun. METHOD SELECTION - STAGE 1 The main purpose of Stage 1 is to select those mining methods which should be considered in greater detail. The simplest way to do this is by defining those characteristics required for each mining method and then determining whether the characteristics of the deposit are suitable. However, no one mining method is so restrictive that it can be used for only one set of characteristics, as indicated by the classification system proposed by Boshkov and Wright (1973). In the mining method selection proposed, geometry, grade distribution, and rock mechanics characteristics are ranked according to their acceptability for ten general mining methods.

Jan 1, 1981

By C. van der Poel

When oil is displaced from a horizontal formation by another fluid of lower density, the latter tends to override the former in the shape of a tongue owing to gravity segregation. This gravity tongue has an adverse influence on the oil recovery. If the fluids are miscible, diffusion (mixing) takes place at the interfacial boundary of the gravity tongue. This mixing should have a favorable effect on oil recovery. The report describes a laboratory study of the magnitude of the mixing zone under various conditions, so as to assess the effect of diffusion on oil recovery both in laboratory experiments and under actual field conditions. The technique used enables visual observation and measurement of the size of the mixing zone in transparent glass-powder packs. The results show that in experiments in small models and cores the width of the mixing zone may well be of the same order of magnitude as the height of the model. In such cases oil recovery is favorably affected by mixing. It can further be concluded that, under conditions prevailing in the field, mixing of the injected fluid with reservoir oil is equal to that caused by molecular diffusion alone, eddy-mixing not taking place to any appreciable extent. A simple calculation, then, shows that molecular diffusion is too small for a beneficial effect to be expected from the injection of miscible fluids in horizontal or nearly horizontal reservoirs unless pay zones are thin. INTRODUCTION This paper gives results of experiments made on the mixing which occurs when a miscible displacement is carried out in a horizontal reservoir. This mixing takes place at the interfacial boundary of the gravity tongue formed when the lighter injected fluid overrides the oil present in the reservoir. The object of the experiments was to simulate the field case where, for example, propane with a viscosity of 0.075 cp under reservoir conditions displaces an oil of viscosity 0.6 cp. In the experiments the lighter fluid had a viscosity of about 1 cp and the heavier one a viscosity of about 8 cp, so as to obtain the same viscosity ratio. In order to enable the results to be compared with those published in the literature, a set of experiments with viscosity ratio equal to one was also performed. EXPERIMENTAL TECHNIQUE The technique developed for the purpose enables the width of the mixing zone to be studied as a function of time and place. A glass-powder pack is saturated with a water-glycerine mixture of suitable viscosity, which represents the reservoir fluid. The pack is then rendered transparent by dissolving sufficient ammonium thiocyanate in the mixture to obtain a solution with a refractive index matching that of the glass powder. Alkaline water to which phenol-phthalein has been added is used as the lighter and less-viscous displacing fluid. In those places where mixing or diffusion of the two liquids occurs, the slightly acidic ammonium-thiocyanate solution neutralizes the alkali in the water, and the glass pack shows up white against the red-colored invading water. If a black cloth is hung over the back of the apparatus, the transparent part of the glass pack appears black. In this way the position of the two phases and the transition zone between them is clearly visible, as shown in Fig. 1 (where the red-colored invading water is the grey-shaded zone lying uppermost). In all experiments the amount of alkali in the water was chosen such that the upper contour corresponded to a concentration of 5 per cent of the dense liquid. The lower contour is determined by the size of the glass grains and the thickness of the pack and, consequently, varies from experiment to experiment. Corresponding concentrations of the dense liquid range from about 97 per cent for packs of small grain size (permeability = 1 darcy) to about 95 per cent for those with large grains

By J. C. Wilhoit, E. M. Galle

A three-dimensional photoelastic study was made to determine the stress state around the wall and bottom of a wellbore due to fluid pressure within the wellbore and unequal principal geostatic stresses. Models simulating the wellbore were made from a phthalic anhydride-cured epoxy resin. The frozen stress technique of three-dimensional photoelastic-ity was used to determine the stresses within the model. Two separate problems were solved: (1) hydrostatic pressure on all sides of the model with the exception of the inside of the hole, and (2) uniaxial compression perpendicular to the center line of the hole. The results are presented in the form of contour curves for each stress component. The curves cover a region which extends several hole radii from the bottom and wall of the hole. By proper use of these results, the state of stress may be calculated for any point on the wall or bottom of the hole, as well as for any point within the material that surrounds the hole. A systematic method is given for calculating these stress components for any combination of fluid pressure within the wellbore and system of unequal principal geostatic stresses, provided one of the principal geostatic stresses is parallel with the wellbore. The results show that stresses around the wall and bottom of a wellbore, induced by an uneyual system of principal geostatic stresses, are appreciably different from those induced by geostatic stresses that are hydrostatic. For uneyual principal geostatic stresses, the experimentally determined stresses on the wall of the hole several radii from the bottom are in good agreement with the stresses calculated by elastic theory. INTRODUCTION Interest in the problem of determining the stresses around and near the bottom of a wellbore, due to geostatic loading and fluid pressure within the wellbore, originated in the petroleum industry. The stresses on the wall of the wellbore are of interest due to their relation to lost circulation during drilling, the fracturing of formations during squeeze cementing and hydraulic fracturing of producing formations. The stresses on the bottom of the hole are of interest to the drilling segment of the industry because to produce a hole in the earth it is necessary to remove material from the bottom of the hole. This removal of material is produced by structural failure of the rock, whatever the criterion for failure may be. Cunningham and Eeninkl have shown that for certain conditions the effect on drilling rate of the bottom-hole stresses caused by a hydrostatic overburden pressure is small compared to other effects. This question, however, needs further experimental and analytical investigation. Much information concerning the stresses on the walls of wellbores has been published considering both equal and unequal principal geostatic stresses. The problem of determining the stresses around the bottom of the wellbore is considerably more complicated. However, previous investigators have obtained numerical and experimental solutions to the problem for equal principal geostatic stresses. whitworth2 and woods3 have obtained numerical solutions. word415 obtained a three-dimensional photoelastic solution, Durelli and Deily 6 obtained surface stresses only by photoelastic means, and Cheatham and Wilhoit7 determined the stresses around the bottom of a short cylindrical cavity. While Miles and McGuire, Harrison and Kieschnick and Hubbert and Willis 10 obtained the stresses on the wall of a wellbore for unequal principal geostatic stresses, they did not concern themselves with the stresses near the bottom of the hole. The estimates made by Hubbert and Willis of the differences between the principal geostatic stresses indicate that, to obtain a satisfactory solution for the stresses around the bottom of a wellbore, unequal principal geostatic stresses should be considered. This work is directed toward obtaining such a solution. The materials in the earth's crust are nonhomo-geneous, permeable and anisotropic. Previous

By K. E. Gray, G. M. Myers

In previous publications,5 results of single-blow bit tooth impacts on saturated rocks at various stress states were reported. This paper extends these earlier works to include study of bit impact tests on salt water-saturated Berea and Bandera sandstone samples under conditions of elevated confining and pore pressure. During the tests dynamic filtration and deposition of a mud cake were occurring due to the presence of drilling mud in the borehole and a borebole-to-formation pressure differential. Results indicate that saturation of these sandstones with salt water tends to make them weaker than when saturated with nonreactive fluids. Plastic failure often occurs even when extremely high fluid-loss muds are present in the borehole. The failure mode tends toward plasticity with decreasing fluid loss. Brittle failure of sandstones in mud-filled holes is apparently relatively rare. Depending on the stress level, and the associated failure mode, withdrawal of the bit tooth induces a tensile force that seems rather important relative to the quantity of rock removed by vertical tooth impact. INTRODUCTION A vast majority of the oil wells drilled today involve the use of colloidal muds having a measurable mud filtrate loss. It is known from field experience that reduction of the water loss of a mud generally results in a reduction of the penetration rate. This paper describes an investigation of crater formation at simulated bottom-hole pressure conditions for drilling fluids having different water losses. The literature on single bit tooth impact (cratering) at atmospheric conditions is extensive, but only a limited amount of work has been performed at the stress conditions similar to that found in oil wells. Spherical penetrator cratering tests on rocks at hydrostatic pressure were performed by Payne and Chippenda1e.1 Chisel impact studies on limestone with independently varying overburden and borehole pressure and atmospheric pore pressure were reported by Garner, et al.2 Gnirk and Cheatham3 performed "static" penetrator tests on dry rocks at equal overburden and borehole pressures. Podio and Gray4 studied the effect of pore fluid viscosity at atmospheric pore and borehole pressures and varying overburden pressures. Their results illustrated the importance of fluid saturation of the rock. Yang and Gray5 reported single tooth impact tests on saturated rocks at elevated borehole, pore and confining pressures, but with equal pore and borehole pressures. For the work reported in Refs. 4 and 5, filtration across the hole bottom purposely was avoided. Maurer6 has investigated the effect of independently varying pore and borehole pressures at elevated overburden pressures by using a zero water-loss mud in the borehole. His tests were carried out under conditions that allowed independent control of overburden, pore and borehole pressures, but without control of filtration at the hole bottom. EXPERIMENTAL APPARATUS AND PROCEDURE APPARATUS The single blow chisel impact apparatus used by Garner,2 as modified for this study, is shown in Fig. 1. A cross section of the pressure vessel is shown in Fig. 2. Several pieces of auxiliary equipment were added to the original apparatus in order that the borehole and pore pressure could be applied in a manner analogous to the actual bottom-hole situation. The pore-pressure system included a filtrate volume-measuring transducer and had sufficient surge capacity to allow filtrate to be injected into the sample without significant changes in pore pressure at points remote to the borehole. The borehole-pressure system contained a cross connection to the pore-pressure system so that both systems could be pressured simultaneously. The borehole-pore pressure differential was applied rapidly by the use of a quick and full-opening valve. The operation of this valve triggered two industrial

Jan 1, 1969

By P. H. Turnock, R. D. Pehlke

The effects of mutual interactions between elements on the solubility of nitrogen in liquid iron alloys containing chrormium, columbium, molybdenurn, nickel, or silicon have been determined. The equilibrium nitrogen solubilities in the liquid iron alloys were measured by the Sieverts method. The second-order effects caused by the presence of chrormium, columbium , or silicon in the melt were found to be significant. The solubility of nitrogen in liquid iron alloys containing several alloying elements has also been measured as a function of melt temperature and nitrogen pressure. The heat of solution of nitrogen in microcomponent iron alloys has been found to be a function of the logarithnz of the activity coefficient of nitrogen, irrvespectilje of the composition of the melt. The heat of solution of nitrogen in pure liquid iron was determined to be 1265 ± 450 cal per g-atom of N over the temperatxue range 1600° to 1800°. Sieverts ' law was obeyed for all melt compositions studied in the pressure interval 0.4 to 1.0 atm. The solubility of nitrogen in liquid pure iron and in many dilute binary liquid-iron alloys has been the subject of many investigations. Pehlke and Elliott' made a comprehensive investigation on this subject in which they summarized and compared previous researches. However, the practical application of these solubility studies requires prediction of the nitrogen solubility in the multicomponent alloys encountered in operating practice. The interaction parameter as defined by wagner2 has been used to predict the behavior of a solute in a complex alloy, and a modification of this approach has met with reasonable success in several instances.'1"5 The necessity for extending the nitrogen solubility data on binary alloys is caused by the lack of direct measurements of the solubility of nitrogen in molten alloy steels. A direct comparison between measured solubilities and those predicted for multicomponent alloys based on data in binary systems has been made in only a few instances, and most of these were reviewed by Langenberg,6 who developed a graphical technique for predicting nitrogen solubilities in multicomponent iron alloys. The solubility of nitrogen has also been studied over the entire Fe-Cr-Ni ternary system by Humbert and Elliott.? Both studies compared the solubilities computed from binary iron alloy data with those measured in ternary alloys in the liquid state at 1600°C. Most of the nitrogen solubility measurements referred to above were for iron alloys at 1600°C; thus, the ability to predict nitrogen solubilities in molten iron alloys has been restricted primarily to alloys at 1600°C. Recently, Nelsson8 devised an approximate method for calculating the change in nitrogen solubility in molten iron alloys as a function of temperature, but an assumption of regular solution behavior limits the applicability of the method. Chipman and Corrigan9,10 have developed an empirical correlation. based on the activity coefficient of nitrogen in the iron alloy at 1600°C, which facilitates the estimation of nitrogen solubilities at higher temperatures. This method also has limitations, due to the necessity of calculating the activity coefficient for nitrogen in a multicomponent iron alloy at 1600°C from the nitrogen solubility data for binary iron alloys. The shortcoming of such a straightforward approach to the problem of predicting nitrogen solubilities is the disregard of mutual interactions between elements in a complex iron melt. This investigation was conducted in an attempt to provide a more reliable method of predicting nitrogen solubilities in multicomponent iron alloys. The effects of the mutual interactions between chromium, nickel, silicon, molybdenum, or columbium in iron were studied, using the Sieverts' technique1,7 to measure the effects of alloy composition, temperature, and pressure on the solubility of nitrogen in ternary liquid iron alloys.

Jan 1, 1967

By Peter M. Fozzard

The San Bartolomt silver deposit was discovered by a United Nations team of geologists working with a counter- part team of staff of the Ecuadorian Government's General Directorate of Geology and Mines (DGGM). The Phase I project over five areas of operation during which this discovery was made was financed through a United Nations Development Program contribution of $860,000 and an Ecuadorian Government participation of about $680,000 equivalent. Existence of mineralization at San Bartolomé was first recognized in 1966 when a reconnaissance geochemical stream-sediment program was carried out over a 10 500 km2 (4,054 sq mile) area centered around the city of Cuenca in southern Ecuador. Over the course of this project more than 10,000 samples of stream-sediment, soil, and rock materials were collected. Average costs per sample were $10.00, as follows: collection, $5; sample preparation, $1.60; assay work, $3.40. Initially in the program (1965-1966) sediments were collected along all motorable roads and tracks. Routine analyses, originally by colorimetry and later by AAS, were for Cu, Pb, Zn and occasionally As, Ni, and Co. Discovery geochemical samples from the San Bartolomé drainage area were among the first to be taken in 1965. As luck would have it, the road along which this traverse was taken passes only some 100 m (328 ft) downstream from the mineralized suboutcrop. The relevant two streams sampled are of inconsequential size, some 1 m (3.3 ft) in width, with a limited flow volume throughout the year. Original anomalies were recognized through cold-extractable Cu/Zn analyses. The pH of all river waters was also recorded at sample sites. Resampling the discovery sediments when AAS became available showed total extractable values of between 150 ppm-300 ppm Zn, 550 ppm Pb, and 10 ppm Ag. In 1967 the first followup was by rough-and-ready means. Dithizone field-kits (THM) were used to rapidly test soils from the surrounding area. The area has no outcrop. This work allowed the definition of an area of interest of some 2 km (1.2 miles) x 500 m (1640 ft). Six approximately 1 km (0.6 mile) to 1.5 km (0.9 miles) lines spaced at approximate 250 m (820 ft) intervals were established and used for IP surveys and soil sampling at 50 m (164 ft) intervals. Where one of the lines crossed the westernmost geochemical discovery creek a first-priority geophysical anomaly coincided with a major Zn/Pb (and later Ag) soil anomaly. A pit was sunk and encountered mineralized vein material at a depth of 1.5 m (5 ft). This lucky glory hole confirmed the geochemical discovery and prompted systematic followup. Immediately, costeaning of the west bank of the discovery creek, the Quebrada Ocashuaico, was carried out and exposed weathered and leached mineralization. Nearly 100 mineralized veinlets ranging in thickness from 0.5 cm (0.2 in.) to 8 cm (3 in.) were exposed across the 180 m (59 1 ft) costean. Sixty mineralized veinlets were determined over a 90 m (295 ft) wide zone. Sampling delineated minable widths averaging between 0.16 and 0.31 kg/t (5 and 10 oz per st) Ag with associated copper (very sporadic) lead and zinc. Over the period 1968- 1972 a considerable amount of detailed work was completed: detailed IP, EM Turam, magnetic and mise-a-la-masse geophysics, very detailed soil- sampling, trenching, drilling and glory-hole drifting, and cross-cutting. Numerous problems were encountered: the area was found to be a most difficult geoelectrical target; zero outcrop; pyrite (nonargentiferous zones) associated with faulting; mineralized groundwater; extensive clay development; loss of diamond drill-core, and no two-dimensional surface information. As the program developed the complexity of the structure was slowly unravelled, using a degree of good guesswork. Apart from photogeological and ground geological mapping at various scales, the following

Jan 1, 1985

By W. A. Young

Polynomial functions of temperature were obtained for the specific heats, thermal diffusivities, and thermal conductivities of zirconium hydrides containing 4 at. pct U. Three hydrides (H/Zr atom ratios of 1.58, 1.65, and 1.70) were studied over the range a" to 900°C and a fourth (H/Zr = 1.81) was studied over the range 0° to 760°C. The specific heats were determined from enthalpy measurements which were obtained using a unique drop calorimeter specifically designed for use with materials in which high temperature phase transitions and/or high dissociation pressures occur. Thermal diffusivities were measured by the flash method using a pulsed laser. The thermal conductiuities were obtained as the product of specific heat, thermal diffusivity, and density. The specific heats agree, within 10 pct, with values derived using a theoretical model in which the hydrogen and zirconium atoms are treated as Einstein and Debye oscillators, respectively. RELIABLE values of the thermophysical properties of the fuel are required to predict the operating temperatures and temperature response of SNAP nuclear reactors. Among the most important of these properties are the thermal conductivity, specific heat, and thermal diffusivity. A considerable number of investigations1-4 have been made of these properties for the Zr-H and Zr-H-U systems.* However, little of the drides, however, this direct method cannot yield meaningful results, since the hydrogen will redistribute under the influence of the thermal gradient, thus forming a concentration gradient; hence, one has a spectrum of compositions, rather than a homogenous alloy. Although the "average" composition of the material may be identical to the initial uniform concentration, the directly measured value of conductivity will be dependent on the thickness of the specimen, due to the highly sensitive dependence of transport properties on hydrogen content. This dependence is strikingly illustrated by the work of Bickel,5 who found that the electrical conduction of zirconium hydrides ranges from primarily hole conduction to primarily electronic conduction, depending upon the hydrogen content. Fortunately, the direct measurement of thermal conductivity is unnecessary, since it can be expressed as the product of the specific heat, thermal diffusivity, and density, all of which can be directly measured with considerable accuracy. EXPERIMENTAL Specimen Preparation. The combined fuel-moderator material used in SNAP reactors is a hydrided zirconium-uranium alloy containing -10 wt pct U. The alloy used in this work was representative of that used in nuclear reactors except that normal uranium was substituted for the enriched uranium required for reactor usage. It was produced by a triple-arc-melt and double-extrusion process. All specimens were prepared from a single cylindrical extrusion which contained 10.30 pet U, 89,35 pct Zr, and 0.35 pct impurities, The specimens for each composition were hydrided simultaneously with ultrapure hydrogen (10 ppm total impurities) using standard fuel production techniques which routinely yield homogeneous, crack-free fuel with negligible increases in the impurity levels. The hydrogen content of each specimen was determined from its weight gain and the density was measured by liquid displacement, Chemical analyses yielded hydrogen concentrations which agreed with the weight gain data within ±0.02 in H/Zr atom ratio) the concentrations of all other elements agreed almost exactly with the initial values after adjustment for the added hydrogen. The specimens used for the determination of specific heat were centerless ground to 2.00 cm diam after hydriding. A thin slice was carefully removed From each end for metallographic examination. In every case, this examination revealed a uniform structure as evidenced by the appearance and distribution of the two phases present in the fuel at the hydrogen concentrations used. TWO specimens (H/Zr = 1.600 and 1.632) appeared to be entirely 6 phase with equi-axed grains; the specimen with H/Zr = 1.756 showed

Jan 1, 1970

By Lawrence A. Shepard, Richard H. Krock

A series of ductile, two phase W-Ni-Fe composites, sintered in the presence of a liquid phase, were tested in tension. Identical room temperature stress-strain curves were obtained for specimens containing from 80 to 92 wt pct W (58 to 75 vol pct W particles). The composites exhibited a maximum elongation of 29pct at room temperature, and 10.7 pct at 77 °K. The tungsten particles in the composite elongated by the same amount at these temperatures. Single-phase alloy specimens matching the composition of the composite matrix showed about one half the flow stress of the composites. The test results demonstrated that the mechanical properties of W-Ni-Fe composites are determined by the tungsten particles alone and are independent of matrix volume fraction or mean free path over the composition range studied. A general prediction of the plastic behavior of two-phase alloys or composites from a knowledge of the properties, size, shape, and dispersion of the individual components is not, at present, possible. The excellent theoretical and experimental exposition of the dispersion-strengthening problem applies only to cases where the stronger phase is present as fine particles and in volume fractions of a few percent. Materials in which the stronger phase constitutes a significant fraction of the volume, although of broad engineering interest, have received comparatively little analytical attention. This important group of materials is the concern of the present study. The degree of deformation of each phase during the plastic working of a ductile two-phase alloy has been established by Boas and Honeycombe,1 Clare-brough,2 and Clarebrough and berger.3 These workers empirically correlated the recrystallization temperature of the individual phases with the degree of cold work sustained by each. It was shown that the deformations of both phases are equal and therefore the same as that sustained by the composite whole when the volume fraction of the stronger phase exceeds 35 pct. At yielding, however, deformation begins in the weaker phase. A direct correlation of composite strength, as a function of deformation, with the strength of the individual components cannot be made. In the presence of stronger phase particles, the flow stress of the weaker matrix phase is enhanced by an additional hydrostatic component. The development and potential magnitude of matrix strength promotion upon straining has not even been experimentally established. Quantitative work in this field has been confined to ceramic-metal composites. These materials fracture at strains considerably below 1 pct by cracking of the ceramic particles.4-6 An indication of the potential constrained matrix strength capacity is given by the observation by Nishimatsu and Gur-land4 that the maximum tensile strength of a WC-Co composite is almost double the bulk strength of the cobalt matrix. The relationship between composite strength and geometry is ordinarily expressed in terms of the mean-free matrix path (referred to subsequently as MFP) between the stronger particles.* This factor takes into account both the volume fraction of matrix phase and the average size of the particles. unke17 concluded that composite strength varied directly with the decreasing log MFP, following the Gen-samer8 relationship. However, Gurland and ~radzil~ reported a strength maximum at an MFP of the order of 1 p in WC-Co composites. A ductile two phase metallic composite was chosen for the present investigation of composite mechanical properties with the view of avoiding the difficulties and ambiguities noted above. The W-Ni-Fe heavy alloy of Green, Jones, and pitkin9 contains 80 to 94 wt pct W, and nickel and iron in a weight ratio

Jan 1, 1963

By H. J. McQueen, G. C. Kuczynski

ThE useful properties of k6-Alfenol and the thermenol alloys (summarized by J. E. Nachman and W. J. Buehler') have encouraged research into those alloys which, having compositions near 25 at. pct Al, 75 at. pct Fe, undergo an ordering transformation between 500° and 550°C. The equilibrium diagram of the iron-rich iron-aluminum alloys was first determined by Bradley and Jay.2,3 They distinguished three structures: the disordered, which is bcc; FeAl, which is simple cubic with one Fe atom and one Al atom per lattice site; and Fe3Al, which is fcc with three Fe atoms and one Al atom per lattice site. This last structure may be conceived as eight body-centered cubes with Al atoms at the centers of alternate ones, so that Al atoms have only Fe as first and second nearest neighbors. Later, Osawa and Murata4 determined the boundaries of the different superlattice regions in the iron-rich Fe-Al solid solution by dilatometric and X-ray methods. Two of the six phases, which they found between 0 and 50 at. pct Al, were never identified by other workers. A. Taylor and R. M. ones,' concurrently with and independently of this research, redetermined this diagram by measuring lattice parameters at temperature. Further reference will be made to their results. EQUILIBRIUM DIAGRAM A) Experimental Procedure—The alloys were prepared from NRC gas-free iron (99.9 pct Fe) and 99.99 pct Al. They were melted in alumina' crucubles under vacuum in an induction furnace. If the ingots were unsound or showed more than a 3 pct weight loss, they were rejected. The ingots were homogenized at 1100°C for 24 hr before hot-rolling at 700" to 900°C from 1 in. to 0.25 in. thick. It is interesting to note that alloys containing more than 28 at. pct Al were able to take an 80 pct reduction of area in the ordered FeAl region. The dilatometer specimens and specimens for chemical analysis were cut from these plates. A Leitz dilatometer was used to measure the coefficient of linear expansion as a function of temperature. The specimens were slowly cooled from 900" to 300°C where they were annealed for 10 hr before slowly cooling to room temperature. The expansion was recorded photographically as the alloys were heated and cooled at 1.5oC per min between 100" and 1000°C. The magnetic measurements determining the Curie points were performed on the dilatometer specimens in the dilatometer, by using direct current in the furnace windings both to heat the specimen at 1.5°C per min and to serve as the primary coil. The regular specimen holder was replaced by a thin alundum tube on which was wound, as a secondary, a closely spaced coil of fine copper wire. A reversal of current in the furnace windings produced, in the secondary, an emf which was proportional to the induction in the specimen, and which could be detected on a galvanometer. In this way, it was possible to determine the Curie points of alloys to within * 5°C. B) Results and Discussion—The coefficients of linear expansion, as obtained from the dilatometer curves, are plotted as a function of temperature and composition in Fig. 1. Since in many cases the heating and cooling curves were superimposed, it was decided to measure slopes of heating curves only. Definite high peaks, due to the transformations, were observed on each curve between 400" and 520°C, and, for alloys containing more than 23 at. pct Al, a second set of smaller peaks at temperatures above 660°C. The sharpness of the peaks indicates that ordering to Fe3Al is a first-order reaction. The compositions and the temperatures of these peaks are listed in Table I. When these temperatures are plotted as a function of the composition, the result is the constitution diagram presented in

Jan 1, 1960

By William C. Hagel

Oxygen-18 exchange between gaseous oxygen, held at a pressure of 125 mm Hg in a PL-IORh chamber, and splzeres of a Fe2O3 containing three or less grains was determined from 9000 to 1250°C. Isotope equilibrum on crystal surfaces appears to be quickly established, and the rate-determining factor is self-difiusion confbvrning to the relation D = 2.04 esp(- 77,900/RT) sq cin sec-'. Changing sphere diameters caused no detectable variation in diffusion coefficients. Oxygen ions diffuse faster than cations at low temperatures, but a crossover may occur at about 1150°C. INVESTIGATIONS by Davies et al.' and Himmel et al2 have suggested that the mechanism of iron oxidation above 600°C involves a predominance of cation migration in the wüstite (FeO) and magnetite (Fe301) layers and oxygen-ion migration in the hematite (a Fe2O3) layer. Evidence for little cation migration during hematite growth consists of: 1) the disappearance of radioactive-silver markers, presumed to vaporize from the oxide-gas interface, and 2) a large discrepancy between parabolic rate constants for magnetite oxidation and those calculated from Lindner's3 data for cation self-diffusion, D:, using Wagner's theory.4 Reversing the calculations gave anticipated self-diffusion coefficients, D:, for oxygen ions of 3.6 x 10-9 sq cm sec-1 at 1100°C and 5.9 x 10-10 sq cm sec-1 at 1000°C. Since then Kingery et a1.5 reported three values for oxygen-ion diffusion in polycrystalline, sintered hematite. These were obtained near 1200°C, at an oxygen pressure of approximately 150 mm Hg, and could be represented by the Arrhenius expression, Dg = 1 x 10 11 exp(-146,000/RT) sq cm sec-1. Extrapolation to temperatures below 1100°C infers that oxygen ions are less mobile than cations. This contradicts the aforementioned oxidation evidence and also coble7s study of the initial sintering of hematite. From the latter, apparent diffusion coefficients for the slowest rate-controlling species could be calculated, and fair agreement was found with an extension of Lindner's data for cation diffusion. Fig. 1 depicts such comparisons graphically. Three points have been included for the determinations of Himmel et al. on cation diffusion in large, natural hematite crystals surrounded by pure oxygen at atmospheric pressure. No dependence on crystal log raphic direction was noted. Since the value of 10-l4 sq cm sec-' at 1000°C actually represents an upper limit, one might interpret that there is a small decrease in D: when single crystals are employed. On the whole, however, cation diffusion appears relatively well-defined,* while the few results tension to lower temperatures. After having gained experience measuring oxygen-ion diffusion in certain glasses7 and in the corundum oxide, a Cr2O3, the purpose of this work was to determine D; values for hematite and thus possibly assist a mechanistic understanding of the kinetics of the oxidation and reduction of solid iron oxides. I) EXPERIMENTAL PROCEDURE A) Material. Hematite powder of moderately high purity was received from the Columbian Carbon Co. An accompanying lot analysis listed 0.03 H2O, 0.008 weight percents, with the remainder Fe2O3. Subsequent treatments were ball milling, drying, prepressing, and hot-pressing for 60 min at 1350°C under 5000 psi. Since graphite dies coated with a slurry of alumina were used for hot-pressing, the surrounding atmosphere at temperature was sufficiently reducing to form some magnetite on the outer surface of each cylinder. After removing any attached alumina, a stabilizing anneal of 50 hr at 1300°C in oxygen was added. The final product was nonmagnetic and exceeded 97 pct of the density (5.24 g cm-3) usually given for FezO3. On polishing cylindrical ends, counts were made of the number of grain boundary intercepts per

Jan 1, 1967

By Frank T. Weems

The basic operating properties of the DorrClone are discussed and certain metallurgical applications which exploit these properties are presented. An effective method of controlling the consistency of the DorrClone underflow is described. THIS paper presents in a general manner some of the applications of the DorrClone (Dutch States Mines cyclone) in the metallurgical and heavy chemical fields. Technical or mathematical analyses of the design factors and of the mechanism whereby the cyclone effects its separation have been capably treated in a series of previous papers.'-' Engineering and operational aspects of a few installations and laboratory tests have been chosen to emphasize the unique classification characteristics of the DorrClone and to illustrate its utility in certain difficult situations. Its use in coal washing and heavy-media operations are not considered. Some 25 or more industrial installations and applications of the DorrClone in the metallurgical fields are operating at present, ranging from single trial units to complete batteries composed of a multiplicity of units. The data were obtained in part from some of these, and in part from operations of units at the Westport Mill, the laboratories and testing plant of the Dorr Co. at Westport, Conn. The DorrClone was developed, in principle, by the research dept. of the State Mines, Limburg, Holland, during World War 11. Its first use was to dewater a —50 micron sand used as a dense medium in coal washing. The metallurgical wet cyclone was first publicized by Driessen and Yancey and Geer.6 ince then the cyclone has been studied as a means of de-watering, heavy media concentration, and fine size classification. The DorrClone is shown in Figs. 1 and 2. pulp to be treated is pumped through the feed nozzle, A, into the feed cylinder, B. The feed nozzle is of rectangular cross-section, the outermost edge being tangent to the feed cylinder walls. Feed shims, C, may be attached to the innermost walls to adjust the open area of the feed nozzle available to entering pulp. The tangential injection of feed under pressure supports the vortex action of the pulp that has been detained in the cyclone; the centrifugal forces of this vortex throw the coarse particles to the walls of the cone, D, where they collect and pass downward and out through the adjustable apex opening, E. The fine particles contained in most of the water that entered with the feed pulp move to the inner spiral of the vortex and flow out through the vortex finder, F, contained in the orifice plate, G, which forms the top of cyclone chamber. The overflow pulp may then spray freely into a spray chamber, H, surmounting the cyclone or pass directly into a 90' elbow, I, of somewhat larger diameter than the vortex finder. Installation Where more than one cyclone is to be installed, best overall efficiency and ease of control are obtained with equal distribution of both pulp volume and weight of solids among the cyclones, and with the same pressure being maintained at each cyclone feed nozzle. Fig. 3 is a suggested arrangement for feeding two or more cyclones in parallel. Feed pressure control: Under normal operating conditions it is not unusual for some surge to exist in the feed that is to be sent to the cyclone. As will be discussed later, surges in the solids content of the pulp are usually not of great operational consequence, but surges in feed volume are a problem. Several systems have been devised for absorbing surges in feed volume, each involving a surge tank. If the surges are small and fairly rhythmical, a tank of sufficient size might be installed to dampen out

Jan 1, 1952

By J. A. Stachura

Continuous mining machinery provides the coal industry with one way to compete for a larger share of the total energy market. Various types of machines are discussed and some of the problems with continuous miners, encountered by operators, are reviewed. Equipment manufacturers are working with mine personnel to provide solutions for problems that arise. While coal production over the past 25 or 30 years has been on a horizontal plane, coal's share of the total energy market has declined. To participate more effectively in this total energy market, it is necessary to produce coal more efficiently. It is the obligation of all management, employes, and mining departments to gear the deep mining industry to the rapid progress and changing of today's modern industry. This can be accomplished in the near future with the selection of the proper type of continuous miner best suited to each operator's individual situation. In most mining operations there is tremendous incentive to undertake the continuous mining program. It can reduce the size of the mine greatly by permitting a minimum of working places; it makes pillar recovery work more efficient from the standpoint of overall cost, amount of coal recovered, and safety. The work force can be reduced materially permitting closer and more efficient supervision. It simplifies maintenance because equipment can be more readily standardized. The trend of the coal market favors the use of continuous mining machines. Although there appears to be a general feeling that continuous mining is still a relatively new program and will be slow in replacing conventional mechanical equipment, the fact is that tremendous strides have been made since the first machines were installed in 1948. This program is advancing at approximately the same rate that mobile loading machines replaced hand loading. From 1948 to 1955 there were approximately 450 continuous mining machines in service. In October 1959, a survey revealed that there were more than 700 continuous mining machines in service. Many operators have expressed a desire to undertake this program, but they feel that they could not do so at this time because of one or more of the following reasons: 1) the thickness of their coal seams, 2) seam characteristics, 3) soft bottoms, 4) bad roof conditions, 5) size consist, 6) insufficient flexibility in machines, 7) difficult ventilation problems, and 8) high maintenance costs. With the realization on coal about the same today as it was in 1948, or slightly less and since coal is still failing to participate to a greater degree in the total energy market, it is not surprising that the coal industry is desperately exploring more economical methods for deep mining. The manufacturers are aware that the coal industry is willing to invest in continuous miners if the equipment is built for maximum flexibility, will produce higher tons per man, and assure long life between overhaul programs. CONTINUOUS MINING MACHINES Before discussing details regarding the selection of a continuous miner, let us have a preview of some of the continuous mining machines which are available to the coal industry today. Jeffrey Manufacturing Co.: The machine shown is the Jeffrey 76 A.M. Colmol. This is their most widely used miner, and has been particularly successful in central Pennsylvania and in high-wall mining in western Kentucky. One of the outstanding features of this auger-type miner is its portability. The entire mining range can be changed from its lowest point to the maximum height without stopping the mining operation. Jeffrey 76 B.M. Colmol: This machine is similar to the 76 A.M. model; however, it is built bigger and stronger for a mining range of 50 1/2 to 72 in. This is the model that is now available (Fig. 2). Jeffrey has added two arms to the top row, omitted the odd arm in the bottom row, thus permitting a 50 pct larger throat opening. This eliminates one

Jan 1, 1961

By W. R. Kanne, P. R. Strutt, R. A. Dodd

A study has been made of the creep behavior of ß-NiAl of stoichiometric composition in the temperature range 475" to 775°C. Single crystal tensile specimens were deformed under a constant applied load. From slip line studies and tensile axis rotations it was concluded that the crystals deformed by pencil-glide on (100) and (110) planes with a (100) slip direction. Electron microscope thin-film studies showed evidence for profuse cross-slip, and the dislocation Burgers vector was determined to be (100). There was no tendency for the dislocations to lie in definite slip planes, although a strong alignment of segments of edge dislocations was frequently seen. The activation energy for creep in the temperature range 475" and 775°C was judged to be much lower than that for creep at temperatures above 775°C. It is possible that the rate-controlling process is the thermally activated cross-slip of screw dislocations. In recent years considerable interest has developed in the possible use of high-melting-point intermediate phases for high-temperature applications necessitating good creep resistance. Ordered alloys such as Ni3A1, modified by alloying additions, have received particular attention, culminating in the recent development by Pratt and Whitney of single crystal turbine blades.' The intermediate phase NiAl has considerable theoretical interest from a creep standpoint because not only is it ordered up to the melting point (-1650°C at stoichiometry) but it contains a high concentration of structural vacancies on the aluminum-rich side of stoichiometry. These vacancies are known to exert a strengthening effect under some conditions,' but in addition they might be expected to affect high-temperature creep rates through dislocation climb. The present paper deals with the slip line and dislocation structures developed during creep of stoichiometric NiAl between 475" and 775C, and with possible attendant creep mechanisms. MATERIALS AND TEST TECHNIQUES NiAl ingots of stoichiometric composition were prepared by induction melting 99.996 pct A1 and 99.98 pct Ni (kindly provided by the Kaiser Aluminum Co. and the International Nickel Co., respectively) in 100 cc pure alumina crucibles in a helium atmosphere. No chemical melting flux was used in order to avoid con- tamination of the alloy. Each ingot was slowly direc-tionally solidified to avoid porosity and to produce very large grains up to 2 in. in length. In this way, ingots were obtained which were suitable for sectioning into single crystal creep specimens. Chemical analysis indicated that all alloys were within the range 49.6 0.3 at. pct Ni. For the present investigation, identical rectangular single crystal blanks were cut vertically from one particularly large crystal in an ingot using a series of three diamond cut-off wheels with appropriate spacers. Thus, all sliced crystals had the same orientation with the tensile axis parallel to the long side of the rectangle; the initial orientation of the tensile axis is shown as S, in Fig. 1. The blanks were homogenized at 1250 to 1300°C for 48 hr and subsequently cooled slowly to room temperature. A gage length of rectangular cross section and about 0.75 in. in length was prepared on each specimen using an electrolytic lathe and a dilute nitric/phosphoric acid electrolyte. The gage length cross section was 0.250 0.0005 in. by 0.040 ± 0.0005 in. Holes to accommodate tensile grips were electrolytically machined using a rotating acid jet. After creep testing, this same technique was used to cut discs for transmission electron microscopy. Thus, all shaping was done electrolytically, the only exception being that (122) sections were spark-machined from some speci-

Jan 1, 1970

By E. L. Dougherty

A phenomenological theory for a one-dimensional unstable miscible displacement similar in type to the Buckley-Leverett model but including the effects of mixing is proposed An equation giving the fractional flow of pure solvent through an oil phase containing dissolved solvent is derived including the effects of gravity and nonbomogeneities. Numer ical integration of the resulting pair of nonlinear hyperbolic partial differential equations gives volumes of dissolved and undissolved solvent as functions of space and time. Parameters in the model which characterize mixing due to dispersion were determined from experimental data; the parmeters correlated with viscosity ratio. The results indicate that the rate of dispersive mixing is proportional to volumetric flow rate Incorporated into our equations were Koval's heterogeneity factor H. which characterizes mixing due to channeling. This allowed satisfactory predictions of displacement behavior observed in horizontal floods in nonhomo-geneous cores; but in most cases the values of H which we used were considerably larger than those used by Koval. INTRODUCTION It has been shownl,2,3 that for favorable viscosity ratios, the diffusion equation with convection satisfactorily describes the behavior of a miscible displacement in a porous media. Numerous attempts to develop a satisfactory mathematical description for the case of unfavorable viscosity ratio have been reported, but they have met with only partial success. These attempts, which are reviewed by oval! have taken two approaches: (1) development of a flow model akin to the Buckley-Leverett approach for immiscible displacement neglecting the effects of mixing5,6 and (2) simultaneous application of Darcy's Law for flow and the diffusion equation with a convection term for mass transfer.7,8,9 We set out to construct a mathematical model of Type 1 including, though, the effects of mixing. Based upon a set of hypotheses, we derived a pair of nonlinear hyperbolic partial differential equations which describe in a one-dimensional system the combined effects of flow and mixing. To work with these equations it was necessary to develop a fractional flow formula. The combined system was integrated numerically using the method of characteristics. In our equations the mixing process is characterized by four parameters. Three of these, labeled ß, P1 and P2, account for dispersive type mixing. The fourth is the heterogeneity factor H, proposed by Koval to account for channeling due to nonhomogeneities in the porous media. Values of the parameters which would cause agreement between theory and experiment were determined by trial-and-error for miscible floods conducted in horizontal cores of both homogeneous and heterogeneous materials. Calculations were also performed for vertical floods in homogeneous cores. The purpose of this paper is (1) to present the details of the mathematical analysis, (2) to present the results of the calculations, (3) to consider what light the results shed on the mixing process in an unstable miscible displacement, and (4) to provide a firmer foundation for the correlation technique developed by Koval for predicting the behavior of unstable miscible floods. STATEMENT OF PROBLEM The problem is to construct a mathematical model which describes in one dimension the observed behavior in an unstable miscible displacement. The approach is phenomenological in that the equations are based on assumptions which violate certain physical precepts known to apply to the displacement phenomenon. However, the assumptions do allow us to account mathematically for the more essential phenomena. The work is of value if we can quantitatively predict experimental results which heretofore could not be predicted, even though the synthetic nature of the model belies complete explanation of the observations. We assume that the system is comprised of two contiguous flowing phases. One phase, which we call free solvent, has the properties of pure solvent. We call the fraction of the pore volume occupied by this phase the solvent saturation, designated s. The oil phase consists of oil containing dissolved