Search Documents

By Dennis E. Stover

Grade control for in situ uranium leaching is maintaining, at desired levels, the uranium concentration in the pregnant lixiviant which feeds the extraction (ion exchange) circuit. This differs from grade control in conventional mining in that in situ grade control is imposed subsequent to rather than prior to leaching without the direct observation and control possible in conventional mining and milling. Grade control is a necessary engineering task through the entire project life, beginning with preliminary evaluation and design and ending only when the final wells are shut in. More specifically, this task is associated with designing the initial well field, maintaining productivity during well field operations, and planning well field expansions to compensate for the depletion of existing wells. Our aim during each of these phases is to maximize uranium productivity while ensuring that we are developing information necessary for future phases. Throughout this task, we remain cognizant of the fact that uranium productivity is a function of both the complex geologic and operational parameters of the project. A minimum list of significant geological factors must include the following: . Host formation mineralogy . Uranium mineralogy . Permeability . Ore quality (grade and thickness) Our concerns are not limited to average or general descriptions of these factors. The interactions among the factors as well as localized variations among them are far more important in obtaining high productivity from a deposit. For example, an ore body may display an average permeability of one darcy but have its high grade ore confined to regions with permeabilities less than 500 millidarcies. This relationship between ore grade and permeability must be understood early in the project life if the operation is to be successful. Similarly, the fundamental operating parameters for the well field include those items which enable us to distribute the lixiviant within the ore as well as the chemistry which will dissolve the uranium minerals. A minimum list of such factors will include: . Well design . Well efficiency . Well bore damage . Pattern size and shape . Oxidant type and concentration . Complexing agent concentration - pH The interactions between these variables and the complex geological variations create a situation which often defies quantitative analysis. Hence, we have adopted a statistical approach when correlating production history with these parameters and the quality of the correlations improves as the number of wells operated increases. However, the correlations are often specific to an ore body. We find that we must tailor the process to meet the unique geologic parameters of each ore body. PHASE 1: INITIAL WELL FIELD DESIGN Designing the initial well field is perhaps the most difficult phase because we have the smallest data base from which to predict performance. We will have conducted numerous laboratory tests and, perhaps, field pilot tests to select the lixiviant system and to characterize the geologic parameters. However, we will learn a great deal more about the deposit and its productivity during actual commercial scale operation. Fully acknowledging its limitations, we will, as a starting point, assume that uranium productivity will be directly proportional to the quality and quantity of recoverable reserves. Future deviations from this assumption will be plentiful and result from our limited understanding of the aforementioned interactions between and among the operational and geological factors. However, in planning the initial field, estimation of the quantity of leachable uranium present at levels above the economic grade and grade-thickness limitations is a priority task. (The concentration of uranium in the ore, the ore grade, is expressed in units of weight percent. The grade-thickness of the ore is the product of the average grade and ore thickness expressed in units of percent-feet.) These parameters provide a direct indication of the quality and quantity of ore. They are determined from gamma logging of closely spaced drill holes (50 to 100 foot centers being typical) accompanied by uranium core analysis or direct uranium assay logging. Gamma- logging measures gross gamma radiation which emanates from uranium decay products and indicates only the possible presence of uranium. Direct uranium assays are necessary to confirm this presence and the true location and quality of the reserve. From this information, an estimate of the in-place uranium reserve is formulated. Utilizing the above data, we begin to assess the quantity and quality of the leachable reserve. Analysis of self-potential and resistivity logs in conjunction with uranium logs will identify those zones which contain uranium and display apparent high permeabilities. Similarly, zones with extensive calcite cementation or dense clay sedimentation are located and uranium within such zones is excluded. (Interpretation of these logs should be confirmed

Jan 1, 1980

By T. L. Chu, G. A. Gruber

Silica films hare been rleposited 011 silicon substmtes at 400° to 600°C by a chemical-transport technique using hydrogen fluoride as the transport agent ill a closed system. This transport takes place from a source materia1 1071: temperature to substrates at higher temperatures, as indicated by the thermochemistry of the transport reaction. The experimental variables of- the transport process, such as the substrate temperature, the pressure pi the transport agent, and so forth, have been studied. The rate -determining step of the transport process appears to he the ),ale of chemical reaction in the source region. The transported films are similar to thermally grown silica films in physical proper-ties with the exception of 'some what higher dissolrrtion rates. SILICA films deposited on suitable substrates serve many purposes in electronic devices. They are used for the fabrication of tunneling devices, the surface passivation of devices, and the shielding of devices from nuclear radiation: and as selective masks against the diffusion of specific impurities into semiconductors. Doped silica films can also be used as sources for the diffusion of impurities into semiconductors. Several oxidation and deposition techniques for the preparation of silica films have been developed to meet the requirements of these applications. The therma1 oxidation of silicon by oxygen or steam at temperatures above 900 C is commonly used in silicon technology. The deposition techniques are perhaps more advantageous since they usually require lower temperatures and are not limited to silicon substrates. Silica films have been deposited on silicon and other substrates by reactive sputtering and chemical reactions. The sputtering of silicon in an oxygen atmosphere is capable of depositing good-quality silica films on silicon' and gallium arenide. Many chemical reactions are known to yield silica at room temperature or higher. These reactions may involve intermediate steps. However, the final step yielding silica should take place predominately on the substrate surface in order to produce adherent films. When silica is formed in the gas phase by volume reactions, no adherent deposit can be obtained. Generally, the experimental conditions of a reaction can be varied so that the surface reaction predominates over the volume reaction. The chemical reactions which have been used successfully for the deposition of silica films are briefly as follows. The pyrolysis of alkoxysilanes in an inert atmosphere or under reduced pressure has been employed to deposit silica films on germanium3 and silicon4 at 650" to 750°C in a flow system. The deposition of silica films from alkoxysilanes has also been achieved at nearly room temperature by a low-pressure plasma. Device quality silica films have been deposited on germanium and gallium arsenide by the deposition of an amorphous thin silicon film followed by oxidation at 600" to 700" . Silica films for high-temperature capacitors have been produced by the hydrolysis of silicon tetrabromide at 950°C in argon and hydrogen atmospheres.7 We have developed a chemical-transport technique for the deposition of silica films on semiconductor substrates at relatively low temperatures. The thermochemistry of the transport reaction, the experimental variables of the transport process, and the properties of the transported silica films are described in this paper. THERMOCHEMICAL CONSDERATIONS The transport of solid substances by chemical reactions in the presence of a temperature gradient has been used for the preparation of films and crystals of many electronic materials. In this technique, a gaseous reagent is chosen so that it reacts reversibly with the solid substance under consideration to form volatile products. Since the equilibrium constants of most reactions are temperature-dependent, the transport of these products to regions of suitable temperature in the reaction system would cause the reverse reaction to take place. depositing the original solid. When the equilibrium is shifted toward the formation of the solid as the temperature is decreased, the solid is transported from a high-temperature zone to a lower-temperature region, and vice versa. This chemical-transport technique can be carried out in a closed or gas-flow system. In a closed system, chemical equilibrium is presumably established in the different temperature regions of the system, and the transport agent regenerated in the deposition region repeats the transport process in a cyclic manner. The local chemical equilibrium may not be approached in a flow system: however, this system offers a greater degree of flexibility. Silica reacts reversibly with hydrogen fluoride and this reaction was chosen for the transport process. The over-all reaction between silica and hydrogen fluoride may be written as: SiO2(s) + 4HF(g-) = SiF4Ur) + 2H2O(^)

Jan 1, 1965

By I. R. Kramer, L. J. Demer

T. H. Alden and R. L. Fleischer (General Electric Research Laboratory)— The authors' results indicate clearly and, we believe, significantly that during tensile deformation the surface layers of an aluminum crystal are hardened more severely than the interior of the crystal. A probable explanation of this effect, as the authors indicate, is that dislocations in the primary slip system may be obstructed at the surface or, it should be added, near the surface. The intent of this discussion is to show that the oxide film on aluminum is not likely to be responsible for this effect, but that the results can be understood if it is assumed the secondary slip is more active in the surface layers than in the interior. Prior study has shown that the principal mechanical effect of an oxide film on a single crystal is to raise the yield stress while leaving the rate of strain hardening during the initial deformation relatively unaffected.33 Since the yield stress is unchanged during polishing in the present case, we conclude that continual removal of the oxide film exerts a small effect on the plastic hardening.* It appears that the hardening interactions are occurring not only at the immediate surface, but to an appreciable depth below it, although with decreasing severity. For example, Kramer and Demer found that with removal of 0.004 in. from a specimen, the easy glide region was extended somewhat; but the yield stress did not decrease. The initial yield stress was recovered only after 0.041 in. was removed. Since a very brief polish would permit dislocations trapped behind a surface film to run out,34 extra dislocations must, instead, be trapped to a considerable depth below the surface. The same conclusion is drawn from the observation of decreasing hardening slope with increasing surface removal rates. If the hardening interactions were only at the immediate surface, a full softening effect would be observed at some small removal rate. The view is taken here that strain hardening is principally caused by small amounts of secondary slip.35 The secondary dislocations will interact in various ways with the primaries, interfering with their motion and causing them to accumulate in the crystal. Prior studies of easy glide have shown Diehl's model of hardening to be qualitatively consistent with the effects of impurities,36 of temperature,36 and of crystal size.37 On this basis the enhanced hardening of the surface layers in aluminum arises from increased secondary slip at and to some depth below the surface. Selective removal of this hardened layer is expected to decrease the measurable "bulk" hardening, the effect increasing with the removal rate and decreasing with the applied strain rate. We suggest that the stress on secondary systems is raised by the bending moments arising from interactions with the grips during the deformation. This stress from the grips has been shown to be a maximum37 near the surface, and hence, increased secondary slip should result. Prior investigations of grip effect:; indicate that as the grip stresses are raised by changing the crystal shape, the easy glide slope increases while the extent of easy glide decreases.38-40 It has been shown also that bending moments superimposed during tensile testing may either decrease easy glide, when supporting the moments caused by gripping, or increase it, when cancelling the gripping moments.38 This interpretation of the authors' results, emphasizing the special importance of secondary slip near the surface, is also consistent with the earlier results of Rosi.41 Copper crystals alloyed with silver in the surface layer show greatly increased easy glide compared with pure copper. In addition, the easy glide slope is reduced. The effect of bulk alloying in extending easy glide has been well established and has been interpreted as indicating the relative difficulty of secondary slip in alloy crystals. Since non-basal glide is difficult in zinc crystals, the effects of surface removal during deformation may be less important. Experiments to test this idea are in progress. I. R. Kramer and L. J. Demer (authors' reply)—The authors wish to thank Dr. Alden and Dr. Fleischer for their discussion. Our interpretation of the data in the paper is that dislocation motion is obstructed by "debris" which starts to form at the surface and extends towards the interior of the crystal with further plastic deformation. The fact that we did not find a reversion from Stage II to Stage I by surface removal shows that in Stage II the "debris" fills the entire cross-section of the specimen. Drs. Alden and Fleischer take the view that bending stresses due to the grips are responsible for the

Jan 1, 1962

By J. C. Melrose

The contact angle is one of the boundary conditions for the differential equation specilying the configuration of fluid-fluid interfaces. Hence, applying knowledge concerning the wettability of a solid surface to problems of fluid distribution in porous solids, it is important to consider the complexity of the geometrical shapes of the individual, interconnected pores. As an approach to this problem, the ideal soil model introduced by soil physicists is discussed in detail. This model predicts that the pore structure of typical porous solids will lead to hysteresis effects in capillary pressure, even if a zero value of the contact angle is maintained. The model is generalized to situations in which the contact angle takes on values between zero and 40 degrees. For the imbibition branch of the capillary- pressure function, the model predicts a considerable departme ftom the usually assumed cos 0 relationship. In fact, according to the model, it is possible that a displaced wetting phase will not be able to reimbibe, even when the contact angle does not exhibit hysteresis. INTRODUCTION The subject of capillarity in porous media has long been of interest in many branches of engineering and applied science. The earlier investigators were those concerned with the physics of soils. 1-5 More recently, petroleum engineers and others dealing with the problems of petroleum production from reservoir rock have given much attention to the subject.6-10 Also, important additions to the literature of capillarity have been contributed from the field of chemical engineering.11-12 These attest to the wide range of industrial applications in which capillary phenomena play a role. The present paper is concerned with the role which wettability plays in capillary action in porous media. As is well known, capillary-rise (or capillary pressure) phenomena have frequently been interpreted with the aid of a model employing the concept of a cylindrical capillary tube. This approach has enjoyed a certain degree of success in correlating experimental results.13 The generalization of this model, however, to situations which involve varying wettability, has not been established and, in fact, is likely to be unsuccessful. In this paper another approach to this problem will be discussed. A considerable literature relating to this approach exists in the field of soil science, where it is referred to as the ideal soil model. Certain features of this model have also been discussed by Purcel114 in relation to variable wettability. The application of this model, however, to studying the role of wettability in capillary phenomena has not previously been attempted in detail. In the present paper, additional features of the model are introduced. These features are critical in determining the quantitative behavior of the model. GENERAL FEATURES OF CAPILLARY HYDROSTATICS BASIC PRINCIPLES When the interstices of a typical porous solid are occupied by two or more immiscible fluid phases, the fluids are microscopically commingled. Hence, fluid-fluid interfaces are found within a certain fraction of the pore openings. The fundamental equation of capillarity specifies the configuration of these fluid-fluid interfaces. This is known as the Laplace equation, when derived from mechanics, and as the Gibbs-Kelvin equation when derived thermodynamically.15 Given two fluid phases, a and ß, in hydrostatic equilibrium, the Laplace equation states that the respective fluid pressures in regions close to the interface are related by Here oaß is the surface or interfacial tension and rl and r2 are the principal radii of curvature of the surface or interface. The pressure difference, Pa-Pß, is the capillary pressure, PC. As Buff15 has shown, Eq. 1 states the condition for hydrostatic equilibrium within the two-phase confluent region, which is referred to as the interface. It thus can

Jan 1, 1966

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 J. A. Herbst, D. W. Fuerstenau

This paper examines the zero order production phenomenon in the context of the size discretized batch grinding model. A restrictive interrelationship between the selection and breakage parameters of the model, which is mathematically sufficient to ensure zero order behavior, is delineated. Based on this interrelationship, a scheme is developed for predicting the values of the selection and breakage parameters for all size fractions from a minimum of experimental data. To test this scheme, dolomite was ground in a laboratory batch ball mill. Successful simulation of the comminution behavior of this system was achieved using the batch grinding model with the parameter values obtained from the proposed scheme. Historically, a preponderance of comminution research has centered on attempts to relate breakage energy to the performance of comminution machines. Concomitantly, grinding data were interpreted almost exclusively in terms of inherently empirical energy-size reduction relationships1-3 or "laws of comminution"4- 6 which were based on highly oversimplified descriptions of the fracture process. In some instances, these relationships provide a crude basis for the correlation of experimental data, but, invariably, this approach is inadequate for meaningful process simulation. The control and optimum design of comminution circuits require a mathematical model capable of depicting the size reduction behavior of every size fraction for grinding conditions of technological importance. Energy-size relationships do not provide this detailed information. For example, an energy-size analysis for batch ball-milled dolomite has recently been completed in the authors' laboratories.7 This analysis provided a satisfactory correlation between the hypothetical size modulus, which characterizes only the finer portion of the size distribution, and the input energy. However, this relation did not constitute an adequate description of the complete system. In general, one of two distinct viewpoints toward simulation has been adopted in recent comminution research. The first approach focuses on the fracture of single mineral specimens, with the essential aim of representing the over-all process in terms of the breakage characteristics of individual particles and the characteristics of the stress field which the particles experience within the particular size reduction device. As illustrated by Harris's review,' the actual incorporation of single-specimen fracture information into a description of the behavior of a multiparticle comminution system has seldom been attempted. The only model to accomplish this incorporation was the one developed by Schönert.8 The Schonert model treats single-passage grinding machines in terms of a distribution of effective loads acting on a particle, a distribution of required particle breakage energies, and the breakage product distribution. However, the use of single-particle fracture behavior to derive models of the size reduction process is at present limited to machines in which the particle residence time is approximately equal to the time required to apply stress. As a result, Schönert's analysis cannot be used to advantage to simulate the complex environment which prevails within a tumbling mill. For multiple-passage systems a second approach to the description of comminution processes has been more fruitful. This approach entails the formulation of a mathematical model which is phenomenological in nature in that it lumps together the entire spectrum of stress-application events which prevail in a system under a given set of operating conditions. The appropriately defined average of these individual events is then considered to characterize the over-all breakage properties of the device. Thus, to analyze the performance of a tumbling mill, the manner in which the particles of a particular size (or size fraction) are stressed need not be distinguished. Instead, a single parameter is assumed to represent the resistance of that size to fracture, given the average grinding environment which exists in the mill. The isolation of such a parameter and a related set of quantities, which constitute the breakage product size distribution for the average event in this size fraction, allows the formulation of physically meaningful descriptive equations capable of yielding precise and detailed information for simulation.

Jan 1, 1969

By T. G. Nilan

Steels containing 0.4 and 0.8 pet C have been transformed isothermally at pressures up to 34 kbuv. Decomposilion mechanisms are so intimately related to phase equilibvia that, as the equilibria shift under high pressure, the microstruclures of the decomposition producls change, maintaining at pressure the same correspondence between a given phase equilibrium and a given microstructure as a1 1 atm. A bainitic unicvo-structure occurs in these steels at high pressure that IS observed only at 1 aim abole 1.4 pct C. This strrrcl111,e 1s ruliorzalized in terms of the effect of pressure on the transformation to lower bainite. The pressuve depetlderlce of the kinetics of pearlile formation is descrihed by an absolute rule theory amalysis. The activation volume for the transformation is 7 cu cm per mlole, which is indicative of a phase-interface transformation-rate cotntrol mechanism. This fortmulation, when expressed in terms of known solule effects on the energelics of tvansformation, gives promise of explaining the effects of alloying elemenls on hnrden-abilily, particularly that of 'coball. KINETIC processes in solids are functions of the primary thermodynamic variables of state: composition, temperature, and pressure. Until the development of high-pressure techniques' that permitted the generation of sustained high pressure at high temperatures, metallurgical studies of such reactions were largely limited to isobaric conditions at atmospheric pressure. The data so obtained, and theoretical deductions based on these data, were necessarily incomplete For the particular case of austenite decomposition, a knowledge of the phenomenology of the reaction in terms of all the state variables should lead to a clarification of the process. This deeper understanding should facilitate a decision between conflicting theories that have been proposed to rationalize these phe- nomena. An analysis, derived from an understanding of the transformation mechanisms, may aid in the optimization of alloy steel development. The first studies of the effect of pressure on transformations in steel were made by Kulin et al. in 1952.' In a study of the influence of applied stress on mar-tensite formation in a 30 pct Ni steel, they found that the martensite start (M,) was depressed 8°C per kbar* by hydrostatic pressure. Jellinghaus and Friedewold3 were apparently the first to investigate the effect of pressure on isothermal transformations above the Ms. They found that the bainite transformation rate in a 1.2 pct C, 3.8 pct Mn steel was reduced by a factor of 3 by a hydrostatic pressure of 4 kbar. In both of these investigations,2,3 the M, of the steels studied was below room temperature. Austenitizing and quenching were done at atmospheric pressure, followed by the decomposition of the metastable austenite under pressure. The maximum temperature at pressure in the bainite study3 was 350°C. Until the advent of the present high-pressure devices, it was not possible to conduct isobaric high-pressure heat treatments from austenitizing through isothermal decomposition and quenching. The exploitation of high pressure as a variable in metallurgical studies was greatly advanced by the apparatus and techniques developed by Hall4 and co-workers. Utilizing these techniques, Radcliffe et a1.5 and Hilliard6 etermined the pressure dependence of phase equilibria in the Fe-C system. Hilliard and cahn7 examined the pearlite transformation rate in an AISI 1080 steel and also in a high-purity 0.92 C Fe-C alloy at 1 atm and at 34 kbar and found a 700-fold reduction in rate at pressure in the 1080 steel but only a fivefold reduction in the high-purity alloy. In agreement with the shift of phase equilibria under pressure, the microstructures were hypereutectoid at 34 kbar pressure, whereas they are eutectoid at 1 atm. Determination of the effect of pressure on carbon diffu-

Jan 1, 1968

By D. S. Rowley, F. C. Appl

Assuming that rock behavior, during cutting with a single diamond, may be approximated by that of a rigid, Coulomb, plastic material, a theory of single diamond cutting action has been developed. Using this theory, the stresses on the diamond cutting surface and the components of the cutting force have been determined. Theoretical results agree reasonably well with available experimental data. The theory of cutting of single diamonds may be used as the basis for subsequent study of the perfonnance of surface set diamond drill bits and other surface set diamond cutting tools. INTRODUCTION For more than twenty years, surface set diamond bits have been used successfully in the mining and petroleum industries to drill and core rock formations that range from medium hard to the hardest known formations. Most extremely hard rock formations cannot be economically drilled or cored by any other means. In view of the growing use and importance of surface set diamond drill bits, it is reasonable to study the way in which these tools cut, and the factors that influence wear of the diamonds. Such knowledge will contribute to better bit design and will improve operating procedures. To analyze the cutting performance of these tools it is appropriate to consider the way in which an individual diamond cuts. Once this is established, the over-all performance of the bits can be determined by considering the total effect of the individual diamonds. In recent years there has been a growing interest in rock mechanics and there now exist several publications dealing with the deformation and failure of rocks. Notable experimental work has been reported by Bredthauer,1 Handin and Hager,2, 3 and Gnirk and Cheatham.4 An extensive collection of data pertaining to rock properties at low confining pressures has been given by Wuerker.5 Cheatham and Gnirk have recently presented an excellent review of the state of knowledge regarding analytical and experimental rock deformation and failure.6 Generally rocks are extremely weak in uniaxial tension, and weak or only moderately strong in uniaxial compression. Both tensile failure and compressive failure are generally brittle in nature. This means that there is little or no deformation or "flow" of the material before shear or tensile fracture occurs. Thus, until fracture occurs, rock behaves very nearly like an elastic material in uniaxial stress. Significant changes in behavior of rock occur under conditions of triaxial compressive stress. There is generally a marked increase in compressive strength and in ductility. It has been demonstrated that the failure becomes plastic (in a macroscopic sense) at sufficiently high values of confining pressure. Sharp wedge chisel indentation studies conducted by Gnirk and Cheatham indicate that many rocks behave in a plastic manner at relatively low triaxial compressive stress states.4 If "failure" under triaxial compressive conditions is considered to be either complete rupture or macroscopic plastic flow, then the general failure characteristics of a rock can be reasonably well represented by the commonly known Mohr envelope. This certainly is not a perfect failure criterion (since the effect of the intermediate principal stress is neglected) but it has been found to be adequate for most engineering purposes. Mohr envelopes usually are constructed graphically using experimental triaxial test data. For purposes of analysis it is convenient to represent the envelope in equation form. Among the several approximate forms appearing in the literature, a parabolic form is suggested by Cheatham,7 and an exponential form is suggested by Paone and Tandanand.8 The exponential form has been adopted for this study, and based on published data, appropriate parameters have been established for 19 rock types. A typical plot of the data compared with the approximate envelope, and a summary of the rock properties and

Jan 1, 1969

By E. E. Underwood, L. L. Marsh

Aluminum single crystals, alloyed with 0.042 atomic pet Cu and 0.11 and 1.1 atomic pct Mg, were subjected to constant stress creep tests, tensile tests, and hot hardness measurements within a temperature range of 300° to 866OK. Calculations based on Dorn's temperature-compensated time parameter, 6, gave a value of DH, = 27,000 cal per mol for the activation energy of early creep in aluminum single crystals. Correlations have been obtained for aluminum alloy single crystals with the parameter E for solid solution strengthening, as well as with the parameter F for solid solution hardening, by using a valence of three for aluminum. Limited measurements on tensile specimens show that the slip band density tends to decrease with increasing temperature and with decreasing solute concentration. INCREASING interest is being shown in the mechanisms of plastic deformation in single crystals during tensile and creep testing. The complexity of deformational processes in polycrystalline materials has led to a search for simpler experimental conditions. Numerous creep and tensile investigations have been conducted with pure, metallic single crystals. To a lesser degree, the effects of alloying elements on the tensile properties of single crystals have been determined. However: the literature dealing with the effects of alloying additions in single crystals under creep conditions is vanishingly small. This paper represents an attempt to narrow this gap in the knowledge of the subject. Previous investigations of creep behavior at the Battelle Memorial Institute1-1 have been of great value in the analysis of the present single crystal data. It is equally desirable, however, to ascertain the extent of correspondence between the behavior of single crystal and polycrystalline materials. For this purpose, correlations, similar to those developed for polycrystalline aluminum alloys by Dorn and co-workers, have been made with the single crystal data. The results from this study have tended to confirm and extend those correlations to the case of alloyed single crystals. Materials and Procedures Three dilute binary aluminum alloys were prepared for this investigation from 99.99 + wt pct Al, 99.8 wt pct Mg, and 99.92 wt pct (electrolytic) Cu. The nominal compositions of the alloys were 0.042 atomic pct Cu, and 0.11 and 1.1 atomic pet Mg. Precautions were taken to avoid contamination of the stock during melting and casting of the alloys. After a heat treatment of 8 hr at 925°F the alloys were extruded, then machined into) threaded tensile specimens with a 3 in. reduced section and a 0.505 in. diam. Spectrographic examination showed less than 47 ppm metallic impurities in each alloy. The single crystals were grown by the strain-anneal method, with a critical strain of about 11/4 pet giving the optimum results. In general, 3 in. crystals were obtained with the Al-Cu alloy, but smaller crystals in the magnesium alloys necessitated the use of 2 in. and 1 in. gage lengths with the low and high magnesium alloy specimens, respectively. After an electrolytic polish, the orientation of that portion of the specimen containing the largest single crystal was determined from Laue back-reflection photographs. Tensile tests were conducted at a constant load rate of about 2 lb per min. Creep runs were made in a constant temperature room, under constant stress at the higher creep temperatures, and constant load at the lower temperatures. The eloneation was measured to within ±5 microin. by a specially designed capacitance extensometer. The ex-tensometer arms were attached to the 3 in. specimens at the shoulders of the test piece or, where the crystals were smaller than 3 in., by knife-edge grips.

Jan 1, 1957

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 L. Valentik

Many attempts have been made during the last 40 years to evaluate the performance of gravity separation equipement, that is, the effectiveness with which light and heavy particles are separated. The most comprehensive treatment of the subject was made by Cerchar at the 1st International Conference on Coal Preparation held in Paris in 1950. The methods suggested by the Conference were accepted and very widely used in the last two decades. This paper discusses an improved method of evaluation in the light of the now-accepted standard presentation. The float-and-sink analysis of the product is presented on a Gaussian distribution curve, resulting in an easier visualization of the inherent difficulties of separation. The ogives of the distribution curve me then plotted, giving a quantitative measure of the deviation from perfect separation as an error distance instead of an error area. Illustrations of the new method are given both for gravel and for coal preparation, but the content is valid and applicable to other types of minerals which are separated by gravity methods. Many attempts had been made during the last forty years to evaluate the performance of heavy-media separation (HMS) equipment, that is, the effectiveness with which floats and sinks are separated.'-' The most comprehensive treatment of the subject was made by Cerchar at the 1st International Conference on Coal Preparation held in Paris. 6 The primary aim was the thorough understanding of the mechanism of separation and the unified presentation of data on gravity separation so that the evaluation and comparison of washery performance could be made from all over the world. No strict overall standardization has been achieved, but after the conference a more or less uniform presentation of performance was accepted, which, during the last two decades, has been very widely used. In this paper, illustration of the old methods and an improved method of evaluation will be given. HEAVY-MEDIA SEPARATION (HMS) PERFORMANCE CRITERIA In the ideal HMS process, all material lower in density than the specific gravity of separation (SGS) would be recovered as floats and all material of higher density would appear as sinks. In order to evaluate the misplaced material, the washery products are tested at the density at which the washing unit is operated. The original type of plot1,7, 8 is shown in Fig 1; this was developed primarily for coal cleaning units. The curve for raw coal represents the cumulative percentages of sink material. The refuse curve is also plotted as a cumulative sink, the percentages being expressed in terms of raw coal. This diagrammatic representation of the results of washing units has the merit of easy visual observance of the degree of separation obtained. The error areas (cross-hatched) are a measure of the amount of misplaced material and therefore they can be used to characterize the quality of separation. The ideal and actual separating performance between floats and sinks can be best seen from the partition curve developed by Tromp,2 where the ordinate is the percentage recovery of the sinks, and the abscissa is the specific gravity (Fig. 2). It can be seen from the shape of the curve that as the SGS is approached, the proportion of material reporting to the improper product increases rapidly. In fact, the SGS can be defined as the density of the material in the feed that is distributed equally between float-and-sink products. When the upper half of the curve is inverted, a shape similar to that of a Gaussian error distribution curve is obtained and therefore the analysis of gravity separation may be carried out by using the law of probability. The shape of the curve in Fig. 2 is determined partly by the density composition of the feed, and partly by the sharpness with which the unit separates floats from the sinks.9, l0

Jan 1, 1970

By H. J. Ramey

A general solution is presented for the transient temperature distribution caused by radial movement of a cylindrical heat source through a homogeneous medium of infinite extent. This problem represents a highly simplified model of the movement of a combustion front during the thermal recovery of oil. Numerical solutions are presented for a heat source moving at a constant velocity, and a velocity inversely proportional to the radial location of the heat source. Numerical solutions are presented for both finite and infinite vertical thickness of the heat source, i.e., with and without vertical heat losses. The numerical solutions are used to estimate the fuel concentration needed to maintain a combustion front during the thermal recovery of oil. The sensible heat carried to the combustion front by the injected gas stream is discussed in regard to the fuel concentration. Numerical solutions for a heat source (combustion front) of finite vertical thickness indicate the field conditions that may be necessary to sustain the temperature of a combustion front above the ignition temperature of the fuel over considerable distances from an injection well. The results also have implications in regard to the quantity of heat required during ignition. The results of this computation apply to a highly idealized model of a thermal recovery process. But the results may be used as a guide in engineering consideration of the thermal recovery process. INTRODUCTIO N The thermal recovery of crude oil has received considerable attention since the publication by Kuhn and Kochl of laboratory and field tests of this method of oil recovery. Other publications2,3,4,5,6, since that time have dealt with various features of the thermal recovery of oil and have illustrated the extremely complex nature of this process. Vogel and Krueger7 described an electric analog computer designed to solve a moving heat source problem which was analogous to movement of a combustion front during the thermal recovery of oil. The heat source was maintained at constant temperature and was a cylindrical source of infinite height (no vertical heat loss) moving radially at either constant velocity, or at a velocity inversely proportional to the radial location of the heat source. Jenkins and Ramey2 presented analytical solutions to a similar heat conduction problem and pointed out the possible importance of vertical heat losses. The following presents a general solution to the transient heat conduction problem introduced by Jenkins and Ramey4 as well as numerical solutions for a wide range of conditions possibly similar to those that might exist during field operation of the analogous thermal recovery process. DESCRIPTION OF PROBLEM The problem considered in this paper is determination of transient temperature distribution caused by a cylindrical heat source of infinite or finite vertical height, moving radially through an isotropic medium of infinite extent. It is further assumed that the heat flux generated at the surface of the moving heat source may be a function of time such that the fuel concentration required may be constant, or permitted to vary, and that heat is generated at the surface of the heat source only. In regard to the thermal recovery of oil the assumption that heat is generated at the surface of the front is equivalent to assuming that the thickness of the combustion zone is infinitely small, or that the reaction rate between fuel and oxidant gas is infinite. The last assumption does not appear stringent. See laboratory information published by Martin, et al3 Under these assumptions the temperature distribution throughout an infinite medium caused by a heat source moving radially with an arbitrary velocity may be described by The heat generation function or heat flux at the surface of the source, +(t), is Note that the fuel concentration may be defined as a function of time. Eq. 1 is derived in the Appendix.

By R. L. Schmidt, J. D. Verhoeven

The diffusion coefficients and electrotralzsport mobilities of aluminum, gallium, and arsenic have been determined in molten germanium with the capillary reservoir technique. The diffusion coefficients are slightly higher than those determined by the solidification technique of Brirton et a1.1 It was found that aluminium and gallium migrate to the cathode in molten germanium) whereas arsenic migrates to the anode under an electric field. A possible technique for producing p-n junctions by means of electric field freezing is discussed in view of the experimental data. In the past decade or so the literature has abounded with articles concerning the solid-state properties of the semiconducting elements germanium and silicon. It is surprising by comparison how little has been published concerning the liquid-state properties of these materials. In the production of doped germanium or silicon the diffusion coefficient of the dopant in liquid germanium or silicon is an important quantity. Yet, to the best of our knowledge, no one has ever made a direct measurement of the diffusion coefficient of solutes in either molten germanium or silicon. A number of authors'-4 have determined some diffusion coefficients by the indirect technique developed by Burton et al.1 This technique appears to give numbers of the correct order of magnitude but it does depend upon a number of assumptions in the analysis and also requires prior knowledge of the viscosity of the melt. In the present work a more direct measure of the diffusion coefficients of arsenic, gallium: and aluminum in molten germanium has been made. It has been shown independently by two groups of authors5, 6 that, when an electric current is passed through the interface of a solidifying metal containing an impurity, the resulting electrotransport in the liquid boundary layer at the interface offers a mechanism for control over the redistribution of the impurity accompanying the solidification. These authors show that the degree of control is proportional to the mobility, U of the impurity atoms in the liquid boundary layer, where the mobility is defined as the drift velocity relative to the solvent per unit electric field. In principle the results apply to the solidification of a doped semiconductor and offer an interesting technique for control of the doping of a semiconductor. In the present work the mobilities of aluminum, gallium, and arsenic were determined in molten germanium in order to evaluate the potential usefulness of electro-transport as a technique for solute control when applied to the solidification of doped germanium. EXPERIMENTAL PROCEDURE Both the diffusion coefficients and the mobilities were measured using high-purity alumina capillaries and a resistivity technique for analysis. A schematic diagram of the electrotransport apparatus is shown in Fig. 1 and the details of the cell design are shown in Fig. 2. For the mobility experiments electrical contact with the bottom of the capillary was made with graphite, as shown in Fig. 2. In order to prevent electrical contact with the melt the graphite was coated with a thin layer of alumina using a plasma spray torch. The cell assembly, a graphite electrode, and another assembly containing two blank capillaries were fastened to the lavite disk shown in Fig. 1. To make a run the system was evacuated and the doped germanium was outgassed at the operating temperature of 1030°C. The cell and blank capillaries were then submerged beneath the melt and a slightly positive pressure of helium was introduced to fill them. After passing current through the cell for a given time the cell and blank capillaries were raised from the melt and the furnace was shut off. The temperature of the cells dropped below the freezing point of germanium within 1 to 2 min after the furnace was shut off.

Jan 1, 1968

By C. S. Matthews, M. Prats, R. I. Jewett, J. D. Baker

By mathematical analysis it was found that injectivity history of a uniform five-spot pattern can be calculated by rather simple formulas. These calculated injectivities were found to agree rather well with injec-tivities measured in a potentiometric analog. With this as a basis, a simple method was finally developed which allows prediction not only of rate of injection but also of rate and kind of production at the production well, for a uniform five-spot. Using this method formulas for the eflective injectivity and production behavior of a waterflooded reservoir having a wide range in permeability can be calculated by considering that the reservoir consists of several layers, each having a uniform but different permeability. Rather good agreement is shown by a comparison of the observed production history of a field in the Illinois basin and that calculated by the methods just described. INTRODUCTION The importance of being able to predict the injection rates for a water flood is well known, for upon injection rate depends the life of the flood, the size of pumping and treating facilities, and the rate of oil recovery. One method of determining rates of injection is through a pilot flood. However, it appears that even if a pilot flood is resorted to, a theoretical method for calculating injection rate will be valuable in extending field results to times after inter- ference and to locations of other five-spots. It is also important to be able to predict the production history of a flood. Such prediction must generally be done by theoretical means, since in general a small pilot flood will not furnish much quantitative information in this regard because of the distortion after oil-bank interference. Methods for predicting the behavior of five-spot water floods have been proposed by Yuster and Cal-houn' and by Hurst.? However, these methods do not consider the effect of the mobility of the different fluids in the reservoir. Other investigators3,4,5,6 have determined the effect of the water-to-oil mobility ratio on the production history by means of different experimental techniques. Most of the published work applies to homogeneous sand bodies, does not provide information for determining the injection requirements of a flood, and seldom considers the rate of build-up of the oil bank as it fills up the partially depleted reservoir.* The present work was undertaken because of the lack of predictive method which takes into account mobility ratio and the rate of buildup of oil bank resulting from the void space in a partially depleted reservoir. As water is injected into the reservoir, which contains oil, water, and gas in macroscopically homogeneous saturation, it forms several banks ahead of the injection well. In each bank or region there is assumed to be only one mobile phase—water in the water bank, oil in the oil bank and gas in the gas region. The saturation in each region and a plan view of the banks are shown in Figs. 1 and 4, respectively. FACTORS AFFECTING INJECTION AND PRODUCTION RATES Results in this report are determined as functions of the three parameters, F, M, and M, , and of a fourth parameter, r,/L** which relates the size of the well to the flood spacing. The parameter F is the displacement factor defined by Hurst in Ref. 2. It determines the rate of build-up of oil bank as it fills up the partially depleted reservoir. When gravity effects are neglected, the four parameters just given are sufficient to describe the waterflood behavior for horizontal and homogeneous reservoirs having incompressible liquids. Definition of InjectIvity The dinlensionless injectivity (or injection rate) used in this report is defined by

By O. D. Gonzalez, R. A. Oriani

A thermo-osmosis technique has been used to measure the heat of transport, Q* , of hydrogen and of deuterium dissolved in a iron and in nickel, and of hydrogen in Feo.6Nio.4 in the tempevature range 400° to 600°C. For all these systems, Q* is negative and has a large temperature coefficient; an isotope effect can be established for the solutes only in nickel. The magnitude of Q* is considerably larger than the activation energy, E , for migration in the case of these isotopes in a iron, so that all variants of the Wirtz model must be rejected. The phenomenological definition of Q* is developed to show that Q* is related to the mechanisms by which the activation energy is dissipated back into the lattice, and that Q*/E is a function of the ratios of the mean free paths of electrons and of phonons to the distance of jump of the diffusing atom. A correlation is shown to exist between the sign of Q* in thermal diffusion and that of the effective charge in electro migration, and may be understood as due to the large ratio of the mean free path of electrons to that of phonons. THE study of nonisothermal diffusion in the solid state is of interest because certain aspects of the mechanism of diffusion are made manifest which remain concealed though still present in the isothermal case. The thermodynamics of irreversible processes characterizes nonisothermal diffusion, known as thermal diffusion or as the Ludwig-Soret effect, by a quantity, Q*, called the heat of transport, which essentially measures the sign and the magnitude of the steady-state concentration gradient produced by the imposed temperature gradient. The problem then is to understand the quantity Q* at the atomistic level. There does not exist a theory for the atomistic interpretation of the heat of transport. There are, however, two classes of simple models which seem to afford some physical insight into Q*. Wirtz's model1, z and its subsequent generalizations partition the activation energy, E, for the migration of the atom spatially about the region of the migrating atom, so that Q* is the difference between the portion of the activation energy centered about the original site of the jumping atom and that portion about the arrival site. Clearly, the magnitude of Q* cannot on this model be larger than the activation energy for migration. The model of Oriani3 also considers the spatial distribution of E not only before but also after the atomic jump, and Q* is related to the variation of the spatial distribution of E that is produced by the jump. The simplest kind of system with which to compare the consequences of these ideas is the interstitial solid solution, since one may easily choose a temperature at which only the interstitial component moves, the solvent lattice serving as a completely satisfactory frame of reference. Thus, one avoids complications arising when both species in a binary system move. In addition, one avoids the necessity of knowing the energy for the formation of a vacancy, something which is needed for the analysis of thermal-diffusion data on metals which diffuse by the vacancy-exchange mechanism. However, the number of well-measured interstitial systems is extremely small, and, in particular, there are almost no data in the published literature for the temperature dependence of the heat of transport. If one works with solutes which are gases in the pure state at ordinary temperatures, one may use the technique of thermo-osmosis4 which permits thermal diffusion to be measured over a small temperature difference so that the temperature dependence of Q* can be more easily determined. The choice of dissolved hydrogen in iron and nickel was based on these considerations as well as on the desire to look for an isotope effect associated with localized

Jan 1, 1965

By C. Gatlin, M. E. Chenevert

The effects of bedding plane orientation on the elastic constants and the yield strengths of three laminated rocks (one sandstone and two shales) and one isotropic rock (a limestone) were studied. The directional dependence of the elastic properties of these rocks was experimentally evaluated using a triaxial compression cell and auxiliary stress-strain measuring equipment. Symmetry of Poisson's ratios within the bedding plane suggested that horizontal isotropy exists, but the bedding planes do give rise to an appreciable difference between properties in the horizontal and vertical directions. For the three bedded rocks studied, Young's modulus was lower normal to bedding than along bedding. Yield strengths were determined at confining pressures from 0 to 12,000 psi in a triaxial compression cell. The rocks studied showed strength reductions as high as 40 per cent when the test specimen was oriented at 20°-30° to the bedding planes. The mechanical behavior of these rocks suggested that the rock properties of shear strength and/or coefficient of internal friction can vary with direction, depending on the particular rock tested. Tensile strengths were also measured and found to be lowest when failure occurred along bedding. This work shows that bedded formations exhibit sizable directional variations in both their elastic constants and yield strengths. It is suggested that these variations may be accounted for by using the "elastic laminate" model and the "variable coefficients" failure model. INTRODUCTION The nature of rock deformation at elevated pressures has been studied by many workers; papers by Handin,l,2 and Robinson3 illustrate the present state of knowledge. Most investigators have either chosen rocks which were as isotropic as possible (in order to avoid complications of data interpretation and analysis) or they have oriented their samples so that the effects of anisotropies (such as bedding planes) have been avoided. Of the studies performed, few were concerned specifically with mechanical anisotrop-ies. Griggs4 has presented limited data for specimens cut parallel, normal, and 45° to the bedding plane; he relates the strength anisotropy observed to the fabric (bedding) anisotropy. His tests were primarily concerned with large deformations (20 per cent strain); thus, no directional values of the elastic constants were reported. Handin has reported the results of similar experiments. Bott5has discussed rock strength anisotropies due to faults, cleavage, or bedding. He was concerned primari ly with determining the shear stress on such planes and did not mention the effect of friction. Jaeger6 later generalized Bott's work by taking friction into account and presented a limited theory for the failure of rocks having a "single plane of weakness", and also for rocks having a constant coefficient of friction, and a shear strength which varies with bedding plane orientation. Donath and Cohen,7 and Donathg have evaluated rock strengths from shale and slate specimens cut normal and parallel to bedding. A dependence of cohesive strength (7,) on the specimen orientation was also shown. Adler9 has also studied this problem and lists similar results. He assumes that all bedded rocks behave according to Jaeger's single-plane-of-weakness theory. Kalinin and Belorussov10 list results for strength tests parallel and perpendicular to bedding and use this information as a basis for hole deviation analysis. From the literature it is apparent that sedimentary rocks have been tested under widely varying conditions of stress; however, the assumption of isotropy is generally, but not always, made. Since geologic sedimentation often deposits sediments in very definite lavers.,. it seems that more svstema tic attention should be given to the possible effects of this natural bedding. Bedding, as used here, refers to visible regularities of grain size or orientation resulting from depositional processes. The mechanical strength anisotropies, which bedding may present, is not the only anisotronic property to be considered, since the elastic

Jan 1, 1966

By R. I. Jaffee, C. M. Craighead, G. A. Lenning

The effects of increasing hydrogen content, introducing a notch, and changing the strain rate on properties of titanium and one of its alloys were investigated over a range of testing temperatures from —196° to 200°C. Both high purity and commercial purity A-55 titanium were used as representative a materials, while a commercial Ti-8 pct Mn alloy was used for an a-P alloy. It was found possible to analyze the data, using the ductile-to-brittle transition temperature concept. Increasing hydrogen, the presence of a notch, and increasing the testing speed raised the transition temperature for the a materials. The presence of hydrogen and notches raised the transition temperature of the a-B alloy also. However, increasing the testing speed generally decreased the transition temperature of the a-p alloy. IN previous papers, the authors have shown that hydrogen is most detrimental to the notch-bend impact toughness of a-titanium' and to the tensile ductility of a-B titanium alloys.' These general effects have been confirmed by Kessler and co-workers," and by Kotfila and Erbin." Kessler showed that the slower strain rate in notch-rupture tests was a more sensitive measure of hydrogen embrittle-ment of a-B alloys than the tensile test. In the present paper the authors cover, with tensile and notch-bend tests, the effect of testing speed, notches, and testing temperature on hydrogen em-brittlement. High purity and commercial purity titanium, examples of u materials, and a commercial a-,f3 Ti-8 pct Mn alloy were investigated. The purpose of the investigation was to correlate the factors producing embrittlement. Materials The high purity titanium for this investigation was obtained by arc-melting crystal bar from Foote Mineral Co. into 1 Ib ingots. As-cast hardnesses for these ingots ranged from 60 to 90 Brinell. The arc-melted ingots were forged to % in. diam rods at 1470°F. These rods were swaged to & in. diam at room temperature prior to machining of the specimens for hydrogenation and testing. The commercial A-55 titanium and the Ti-8 pct Mn alloy were obtained from Rem-Cru Titanium Inc., as % in. diam rod. The A-55 titanium was hot swaged to 1/4 in. diam at 1300°F prior to machining of the specimens for hydrogenation and testing. The a-p C-11OM alloy was hydrogenated prior to hot swaging to 1/4 in. roc1 at 1300°F. Hydrogenation prior to fabrication permitted testing of the a-P alloy in the equilibrated and stabilized condition. The equilibration and stabilization treatment used was an air anneal for 1 hr at 1300°F, followed by furnace cooling to 1100°F and then air cooling. Hydrogenation was done in a Sieverts apparatus described previously.' The hydrogenation treatment consisted of heating for 1/4 to 1/2 hr at 820 °C in purified hydrogen, slow cooling to 680"C, holding 1 hr for homogenization, and rapid cooling to room temperature. Intended hydrogen levels for the a-tita-nium materials were 0.05, 0.5, 1.0, 1.5, and 2.0 atomic pct—about 15 to 430 ppm. For the Ti-8 pct Mn alloy, the intended hydrogen levels were 0.05, 0.5, 1.0, 2.5, and 5.0 atomic pct—about 15 to 1100 ppm. The material for the 0.5 atomic pct level of hydrogen was vacuum annealed prior to hydrogenation to increase the accuracy of the addition. Vacuum annealing to remove hydrogen consisted of heating for 6 hr at 820°C, and resulted in a hydrogen content of about 10 to 15 ppm. When hydrogen contents were above 0.5 atomic pct—100 ppm—the hydrogen was added without prior vacuum annealing.

Jan 1, 1957

By Gilbert Santoro

In this study tantalum carbide filaments of various compositions in the fcc region were prepared by heating a tantalum wire in a measured amount of hydrocarbon vapor. Such properties as tensile strength, transverse-rupture strength, microhard-ness, lattice parameter, electrical resistivity, and magnetic susceptibility were measured as a function of carbon content. Some of these properties showed a maximum or a minimum between a carbon-to-tantalum molar ratio of 0.80 and 0.85, the composition range at which the color changes from metallic gray to gold. The recent literature has shown that in the large homogeneity range for the fcc phase of the Ta-C system there is a marked variation for such properties as microhardness,' electrical resistivity, ' thermodynamic properties,z'3 lattice parameter and superconductivity. In addition to the technological interest of these variations within a single phase, important theoretical benefits may ensue from a continuous investigation of properties vs stoichiometry in tantalum carbide and other refractory defect compounds, for example the effect of the metal-metal distance and crystal structure upon bond formation and electronic structure. Whether there is a similar variation of the tensile strength is of practical importance. Of theoretical interest would be the knowledge of how the magnetic susceptibility varies with carbon composition. Both of these properties, as well as the transverse-rupture strength, were determined for tantalum carbide filaments in the fcc phase region. The lattice parameter, the electrical resistivity, and the microhardness were also investigated and compared with values appearing in the literature. EXPERIMENTAL PROCEDURE Material Preparation. Tantalum carbide specimens of various carbon content were prepared by statically carburizing high-purity tantalum filaments in measured quantities of pure hydrocarbon gas. The filaments, which were 10 mil in diameter and 20 or 40 cm long, were heated by their own resistance at temperatures from 1800" to 2200°C using alternating current in some cases and direct current in others. Tantalum wire from two sources was used. One had a reported purity of 99.96 wt pct Ta. The ingot from which this wire was drawn had been electron beam-melted as the final step in its purification. The reported purity of the other wire, which was prepared by a powder metallurgical process, was 99.88 wt pct Ta. The impurities of both tantalum sources are listed in Table I. Spec-trographic-grade toluene or research-grade propane was used as the source of carbon. The apparatus for the carburization of the filaments is shown schematically in Fig. 1. The system was evacuated to 1 x l0-' mm of Hg when using the electron beam-melted tantalum and to less than 1 p when using the less-pure tantalum. Neither the difference in the base vacuum, nor in the hydrocarbon gas, nor in the tantalum filaments resulted in any detectable difference in the properties studied in this work. However, only the carbides prepared from the higher-purity tantalum were used in the magnetic-susceptibility determinations. Once the base vacuum had been reached, the system was flushed with the hydrocarbon and re-evacuated. For any desired composition of the carbide, hydrocarbon gas was allowed to leak into a constant-volume flask to a predetermined pressure. These pressures were measured by a Cartesian diver-type mercury gage with a range from 0 to 20 mm and with a scale expanded to nine times. The hydrocarbon gas was allowed into the reaction tube and the filament was quickly raised to the reaction temperature. The slack in the filament due to thermal expansion was taken up by drawing out the movable electrode with a screw mechanism. The filaments were kept at temperature from a few minutes to 200 hr. In the majority of cases the heating time was 6 hr. The reaction tube was cooled by cascading water on its outer surface. By means of an ammeter and a voltmeter the electrical resistance at the reaction

Jan 1, 1963

By F. P. Bullen

Empirical relationships between fracture stress, orientation angle, and diameter of crystal have been determined at 77°K. Orientation ranges of markedly different behavior were found—a law of constant normal stress&apos; of value a (diameter)-1/2 for the fracture of ductile crystals, and a condition of shear stress (or strain) for more brittle crystals. The observations are not consistent with current theories. An interpretation is advanced which is also applicable to observations on the effect of prestrain at room temperature on the subsequent fracture stress at 77°K and to the effect of cyclic stressing on the cleavage strength.&apos;&apos; The law of constunt normal stress&apos; and the brittle ductile transition are also explained. The interpretation is more consistent with the initiation of cracks by intersecting dislocations than with theories based on stress -concentration by dislocation arrays. ZINC single crystals are particularly suited to the study of cleavage because fracture occurs on the basal plane over a wide range of crystal orientations. Analysis of the conditions of stress and strain at fracture in crystals of different orientations should indicate which parameters control the cleavage process. Unfortunately, controversy has arisen over the correct empirical relationship between tensile fracture stress and orientation. Schmid&apos;s observations1,2 favored a &apos;law of constant normal stress&apos;, as observed in other materials.2 For zinc, however, the observed values are far below the theoretical strength and cannot represent the true limit of cohesion between neighboring atomic planes. Hence, the interpretation of such a &apos;law&apos; is not straightforward. Deruyttere and Greenough3,4 found a complex variation between tensile fracture stress and orientation; this variation did not agree with a &apos;law of constant normal stress&apos;. Two theories have been advanced to account for their observations: a) the propagation of cracks from low-angle boundaries,5 and b) the release of energy from piled-up dislocations during crack-propagation.4 The present work resolves the apparent discrepancy between the observations and shows that neither of the above theories are applicable to the tensile fracture of zinc single crystals. A phenomenological explanation, along the lines suggested by Gilman,&apos; is advanced and successfully applied to previously unexplained effects. EXPERIMENTAL DETAILS &apos;Crown Special&apos; redistilled zinc was used, except for one comparison series op tests using &apos;Tadanac&apos; electrolytic zinc. Crystals of 1 mm diam, subsequently called &apos;1 mm crystals&apos;, were grown from the melt in vacuo in precision-bore Pyrex tubes internally coated with graphite. Several specimens 1 in. long were cut from each crystal and chemically polished. Jigs were used to minimize handling strains, and crystals were mounted in the Polanyi machine the day prior to testing to allow recovery from any such strains. One-mm crystals were chemically polished for long periods to obtain 0.1 mm (approx) crystals. One-mm crystals were cemented into miniature gimbals by &apos;Araldite&apos; casting resin. The Appendix gives the reasons for using gimbals and the results obtained by other methods. More complete details of all techniques are given elsewhere.7 The symbols and terminology used are as follows: X = orientation angle (angle between tensile axis and line of greatest slope in the basal plane). T = tensile stress (on true cross-section) S,N = shear and normal stress (components of T with respect to the basal plane) ? = shear strain D = crystal diameter. The subscript &apos;f&apos; will be used to denote values at fracture. PART I-ANNEALED CRYSTALS EXPERIMENTAL OBSERVATIONS One-mm crystals were used to establish the variation of fracture stress at 77°K with orientation at fracture (Xf), Fig. 1. For 18 deg = Xf = 55 deg, a &apos;law of constant normal stress&apos; was observed. For Xf > 55 deg, the fracture condition approximated to a constant shear stress. At Xf< 18 deg, twinning occurred before fracture so that the results were not typical of homogeneous single crystals,4,8— such specimens will not be considered herein. The dependences of fracture stress upon Xf were of similar type for 6 mm,* 1 mm, and 0.1 mm crystals,

Jan 1, 1963

By A. R. Elsea, E. E. Fletcher

Cast X-40 alloy with the lowest nitrogen content studied had a 100-hr rupture stress at 1500°F about equal to the reported value for the commercial alloy. Increases in nitrogen content progressively decreased the rupture strength, and with about 0.1 to 0.2 pct N a minimum strength was reached, but there was no accompanying change in microstructure. Further increases in nitrogen resulted in the fomzation of progressively increasing amounts of a lamellar, pearlitic-type structure at the grain boundaries, with accompanying increases in rupture strength, and above about 0.25 to 0.33 pct N, the rupture strength exceeded the normally reported values. Two series of melts representing different melting stocks exhibited minimum rupture strengths at different nitrogen contents. This suggests that some unidentwed trace element not removed by vacuum melting has a strong influence on the response of this alloy to variations in nitrogen content. In a series of melts with manganese contents ranging from 0.6 to 3.2 pct, the material with the lowest manganese content had the highest 100-hr rupture strength at 1500°F. The manganese variations had no appreciable effect on the microstructure. CAST cobalt-base alloys such as X-40* exhibit high strengths at elevated temperatures (1500°F and above), but they also show considerable scatter in high-temperature properties. Since the design of parts for critical applications must be based on the minimum expected properties, this broad scatter appreciably lowers the design strength of the alloy. If the relatively poor properties exhibited by a few test parts from any given group could be improved, the spread in properties would be narrowed, and the design stress could more closely approach the average expected properties. The purpose of this research investigation was to learn what were the causes of these erratic properties and, if possible, to find methods for eliminating them. Earlier unpublished work showed that vacuum melting markedly influenced the high-temperature strength of cast X-40 alloy. There were indications that this behavior may have been the result of changes in the nitrogen or manganese content of the castings. Consequently, in the current investigation, the influence of variations in nitrogen and manganese content on the structure and high-temperature properties of cast X-40 alloy was studied. In this program, X-40 melting stock was vacuum melted to produce a base material that was very low in nitrogen content and in volatile trace elements. Additions were then made to this base material to produce alloys of the desired composition. The metal was cast into investment molds producing coupons which were machined into stress-rupture bars. These bars were subjected to stress-rupture tests and me-tallographic study to determine the influence of changes in analyses (nitrogen and manganese content) on the structure and properties. Material—The melting stock used for these experi-ments was in the form of shot from two master heats of commercial X-40 alloy. Only one master heat was used in a given series of melts. Thus, any differences in structure or properties that showed up in the specimens from a given series would be the result only of the intended variations in composition. The certified analyses of the master heats were as follows: C Mn Si Cr Ni Lot No. 2 0.50 0.71 0.66 25.44 10.53 Lot No. 3 0.54 0.73 0.70 25.95 10.00 W Fe Co P S Lot No. 2 7.45 1.34 Bal. 0.011 0.013 Lot No. 3 7.58 1.53 Bal. 0.005 0.016

Jan 1, 1960