Search Documents

By W. M. Robertson

The rate of crack Propagation thvough a solid meta1 in the presence of a liquid metal has been calculated. The crack is assumed to propagate by solution of the solid in the liquid under the influence of an applied stress, with volume diffusion of the dissolved solute through the liquid controlling the propagation. The calculated velocities can be tens of centimeters per second or higher. THE embrittlement of normally ductile solid metals by liquid metals has been the subject of extensive study, as summarized in the book by Rostoker, McCaughey, and arkuus.' In experimental studies of liquid-metal embrittlement it is usually concluded that the observed crack-propagation rates are too high to be accounted for by diffusion processes. Diffusion is considered to be an unimportant factor and various other mechanisms are invoked to account for the experimental observations. The similarities of liquid-metal embrittlement to stress-corrosion cracking, to hydrogen embrittlement of steel, and to low-temperature embrittlement of bcc metals are frequently used as analogies for these mechanisms.1"3 The present paper assesses how rapidly a crack could propagate by a mechanism of volume diffusion through the liquid metal. We assume that the crack propagates by solution of the solid in the liquid at the crack tip under the influence of the applied stress and diffusion of the dissolved solute away from the crack tip. The model uses only macroscopic thermodynamic, kinetic, and elastic concepts. The calculations indicate that this dissolution-diffusion mechanism can give high crack velocities. In some previous treatments4'' of crack propagation in a liquid-metal environment the fundamental step has been the breaking of atomic bonds at the crack tip. The liquid greatly decreases the energy necessary to break these bonds. The present model assumes the crack is filled with liquid metal; the presence of the liquid is essential because diffusion through the liquid is assumed to control the rate of propagation. With proper modification this model could also apply to stress corrosion in a liquid environment. THEORY Consider a crack in a homogeneous elastic continuum, ignoring all the structural details of the solid. This neglects a host of metallurgical details which can be significant in the problem of liquid-metal embrittlement;' but this simple assumption is a useful approach in evaluating the role of diffusion through the liquid. The stress distribution around the crack tip establishes a gradient in the chemical potential at the tip, which causes a diffusion flux through the liquid away from the tip. The solid-liquid interfacial tension works against crack propagation, since propagation requires creation of interfacial area. From these two considerations we derive the diffusion flux, which then yields an expression for the velocity of the crack tip. The concentration of the solute in the liquid, C, is equal to the equilibrium concentration in the presence of an unstressed, flat surface, Co, plus the excess concentration due to the stress, AC(a), plus the excess concentration due to capillarity, AC(y): We calculate AC(o) and AC(y) and then obtain the diffusion flux. Consider (Fig. 1) a long flat crack of length 2L, thickness 2r, and tip radius r, where r remains constant as the crack propagates in the positive x direction. There is no stress ay in the y direction across the flat surfaces of the crack. The magnitude of the stress at the crack tip due to the applied stress, a,, is given by the relation6 ^=^[l-2^] 12] which, for long narrow cracks, becomes The stress increases the chemical potential in proportion to the strain energy.? The excess chemical potential due to the stress, Ap(a), is given by where 52 is the atomic volume of the solid and E the Young's modulus. Inserting Eq. [3] gives Now consider the crack to be filled with a liquid metal. The excess solubility due to stress, hC(a), is given by the excess chemical potential due to capillarity, Ap(y), is given by8 where y is the solid-liquid interfacial energy, and v1 and YZ are the two principal radii of curvature of any surface point. At the crack tip, r1 = -r and vz = m, so TRANSACTIONS OF THE METALLURGICAL SOCIETY OF AIME

Jan 1, 1967

By Craig R. Barreft, Oleg D. Sherby

The steady-state creep characteristics of pure polycrystalline copper were studied in the temperature range 400" to 950°C and in the stress range 400 to 7000 psi. Tests were conducted in dry deoxidized hydrogen and all specimens exhibited high ductility before fracture. Two activation energies for creep were established as a function of temperature. Above about T/T, = 0.65, the activation energy for creep was 48,000 cal per mole whereas between 0.5 and 0.65 T, it was 28,000 cal per mole. The observed temperature dependence of the activation energy for creep follows the same trend as exhibited by the self-diffusion coefficient of copper in copper single crystals and silver in silver single crystals. It is therefore believed that the rate-cmtrolling mechanism in creep of copper in the stress and temperature range studied is diffusion-dependent and is probably associated with a dislocation-climb process. The low activation-energy region is attributed to enhanced volume difFusion from the presence of dislocations which act as short-circuit paths for diffusion. THERE is extensive experimental evidence1,2 that the activation energy for creep, Q,, for pure poly-crystalline metals tested at temperatures above one half the absolute melting temperature (0.5 Tm) is a constant and equal to the activation energy for volume self-diffusion, Qsd This equality is generally interpreted as meaning that creep in this temperature range is controlled by a dislocation-climb process. The existence of high-temperature creep activation energies not equal to the respective self-diffusion values has been suggested both theoretically3"5 and experimentally.6-8 Two examples of apparent nonequality between Q, and Qsd have been reported for the cases of silver7 and copper,' where activation energies for creep much lower than those for self-diffusion persist to rela- tively high temperatures (up to 0.6 T,). These observations suggest that dislocation climb may not be the only important process in the high-temperature deformation of these metals. That dislocation climb does indeed occur relatively slowly in these two metals due to their low stacking-fault energy has been suggested by many investigators.9-10 Although the creep behavior of copper has been studied extensively at low temperatures, the high-temperature activation energy for creep is not well-defined. The data of Feltham and Meakin8 indicated that Q, is temperature-dependent above one half the melting temperature, equaling a value of -30,000 cal per mole at temperatures up to 550°C and increasing to a value of -49,000 cal per mole at 700°C. No data has been reported above 700°C. The value of 49,000 cal per mole is very near the value of Qsd for copper reported by Kuper et al." as 47,100 cal per mole. In view of the lack of high-temperature creep data and the apparent temperature dependence of Q, for copper the present investigation was undertaken to determine the value of Q, over the temperature range of 400" to 950°C. It was hoped that such a study would clearly document the temperature dependence of Q, and yield significant information as to the nature of the rate-controlling creep mechanisms. MATERIALS AND PROCEDURES The material used in this investigation was high-purity OFHC copper. After rolling and recrystal-lization treatments spectrographic and chemical analyses showed only 0.0008 pct Ag, 0.0002 pctCa, 0.0002 pct H2, and 0.002 pct O2 as trace impurities indicating a purity of greater then 99.995 pct. Specimens were machined from cold-rolled sheet and the gage sections polished with 2/0 emery paper prior to annealing treatments. Annealing was done in both vacuum and dry deoxidized hydrogen atmospheres. No variation in creep strength was found as a result of the different annealing atmospheres. The final annealing temperature was either 700°, 800°, or 1000°C depending on the temperature range to be studied. These treatments resulted in average grain diameters of 0.03, 0.4, and 1.0 mm, respectively. The precise mechanical-thermal history of these

Jan 1, 1964

By S. W. Kennedy, M. Cohen, L. D. Calvert

A high-temperature X-ray diffraction method has been used to study the composition and the kinetics of formation of oxide scales at 800 °C on iron and pure iron-nickel alloys containing 25.6, 75, and 84 pct Ni. This technique has been combined with chemical analysis, and metallographic examination, and determination of weight-gain time curves. The 75 and 84 pct Ni alloys form a three-layer scale consisting of porous NiO next to the metal, a spinel layer, and a -Fe203 on the outside. The 25 pct Ni alloy produces a scale of spinel and Fe,O,, in which, however, traces of FeO and NiO are observed near to the metal. On all the nickel alloys the Fe20, is orieqted with the (107) plane parallel to the surfaces. The exposed surface of the Fe20, is in the form of needles whose axis is perpendicular to the specimen surface. Although the alloys show nonuniform oxidation rates as determined from weight-gain measurements, the X-ray data indicate that the total scale and individml scale layers grow uniformly and parabolically. Deep subscale formation is observed in the 25 pct Ni alloy only. This appeared to be associated with cracking of the scale and penetration of the oxygen along the grain boundaries. ALTHOUGH film formation has been investigated by high-temperature electron-diffraction techniques, for example by Gulbransenl and Hickman,' there have apparently been no investigations of scaling by methods permitting study of the scale during formation at the reaction temperature. Methods involving quenching of scales are based on the working assumption that no significant changes take place on cooling. In general, however, compounds may disproportionate, and the different thermal expansions of constituents may cause disruption of the scale on cooling. In particular, some scales cannot be quenched without spalling. These difficulties are avoided by X-ray diffraction examination at the temperature of formation. High-temperature diffraction procedures also provide a method of determining the kinetics of formation of each phase separately in a multicom-ponent scale. The scaling of binary iron-nickel alloys Containing up to 30 pct Ni at a temperature of 800°C has been studied by Benard and Moreau using cooled specimen~,~ whilst Foley, Druck, and Fryxel14 made an electron-diffraction study of cooled oxide films formed on an iron 42 pci Ni alloy. In the present work the scaling of three iron-nickel alloys containing 25.6, 75, and 84 pct by weight of nickel, at 800°C in oxygen at atmospheric pressure was investigated by X-ray diffraction techniques using a Norelco Geiger-counter diffrac-tometer. The oxidation of iron was investigated for comparison. In addition normal metallographic examinations and chemical and X-ray analyses were made, and weight-gain time curves obtained. EXPERIMENTAL Materials Used—The analyses of the alloys were as follows: Percent by Weight Nickel 25.6 pct 75.4 pct 83.2 pct Aluminum 0.004 None found 0.002 Cobalt — 0.05 0.05 Chromium — — — Silicon — 0.01 None found Carbon All less than 0.05 pct The chemical analysis of the iron was: Al— 0.002 i 0.001 pct, C—0.007 pct, Cr less than 0.001 pct, Cu—0.002 ± 0.001 pct, Mn—0.003 ± 0.001 pct, N—0.0002 pct, 0-0.013 pct, P—less than 0.0002 pct, Pb—not detected, S—0.0017 pct, Si—0.0095 pct, Sn—none detected. This iron was the same as that used to make the alloys and was obtained from the same source-Vacuum Metals Corp. One specimen of pure iron from The British Iron and Steel Research Association was used for preliminary investigation. Specimens were machined to shape from cast rod, polished (the final polish being with 600 Sic paper), degreased in trichlorethylene, washed, and etched. The etching and neutralizing procedure adopted by

Jan 1, 1960

By C. W. Binckley, F. R. Burgess, E. R. Haymaker

A method of establishing stabilized back-pressure curves for gas wells producing from formations of extremely low permeability is presented. Actual well performance under many different operating conditions is shown by the stabilized back-pressure curve. By use of the method. it is possible to conduct back-pressure tests with a critical-flow prover on wells that stabilize slowly, and save approximately 88% of the gas ordinarily vented to obtain satisfactory test data, with a great reduction in time required for testing. INTRODUCTION The reasons for establishing dependable back-pressure curves on gas wells have been pointed out by previous publications. The publication most referred to. of course, is the United States Bureau of Mines Monograph 7, titled "Back-Pressure Data on Natural Gas Wells and Their Application to Production Practices". The technique generally established therein has been accepted and used by many engineers; and, when proper tests are conducted, the results can be used for the analysis and solution of several practical problems concerning field operation and development. Even where formations of low specific permeability are encountered, the determination of a well's actual performance by the back-pressure test method permits the engineer to analyze many problems in individual well operation and also to predict necessary future field development. Such problems as the determination of the ability of a well to produce into a pipe line at a predetermined line pressure, the design of gas gathering systems and meter settings, and the determination of the time and the number of wells required to be drilled to meet future market obligations, can be solved, in part., by the use of a reliable back-~ressure curve. In addition, the computed well delivery rates determined by data from backpressure tests ordered by state regtilatory bodies, when compared with the true back-Pressure curve, permit the operator to ascertain whether such data represent unstable or relatively stabilized delivery rates for given pressure conditions of the well. The technique of back-pressure testing, as described in this report, was developed by Phillips Petroleum Company engineers from data obtained during a testing program that started in 1944 and has been continued to date. Three hundred and eleven back-pressure tests were conducted on 299 wells located in the southern part of the Hugoton Field. The gas-bearing zone is composed of several dolomitic formations of the Permian Age; the important ones are the Herington, Upper Krider, Lower Krider, and Winfield. The average bottom-hole temperature is approximately 91 °F.. and the initial wellhead shut-in pressures range from 400 to 440 psig. The spacing pattern is 640 acres per well with each well located near the center of the section. The range of back-pressure potentials on wells tested was from 500 to 23,000 Mcfd. All gas wells were acidized, and the quantity of acid used, expressed in 1574 hydrocloric acid, varied from 12,000 to 22,000 gallons per well. The quantity and concentration of each treatment depended on the stage, the formation being treated, and experience gained from previously completed wells. The gas in the Hugoton Feld is a "dry" gas. It has a gasoline content of approximately 0.25 gallons per thousand cubic feet, as determined by charcoal test, and its specific gravity averages about 0.71 as compared to air (air = 1.00 at 60°F.). Of the wells tested, 71 were completed with 7" casing, 3 with 9 5/8" casing, and 1 with 6%" casing set on top of the upper producing formation with the well bore through the gas bearing formations being open hole. Two hundred and twenty-four were completed with 5 1/2'' O.D. casing set through the gas bearing formations and perforated. For the purpose of establishing reliable back-pressure curves in the area, Phillips Petroleum Company personnel has computed data on the basis of 24-hour flows per point. Early in the program, many tests were actually permitted to flow 24 hours to obtain data for each plotting point, at great expense in man power and time. Presently, however, such tests have been replaced by tests of short duration flows which can be made to effect results that correspond to the tests obtained by flows of much longer duration. METHOD When a gas well producing from a reservoir of low permeability is opened for production through a constant size orifice, both the rate of flow and working pressure decline. first at a high rate and later at a lower rate until after several hours the decline becomes difficult to ascertain. In this paper the rae of flow and working pressure are considered to be stabilized when it becomes difficult to observe changes in working pressure during a period of three hours by the use of a deadweight pressure gage. Stalibization of pressure in the literal sense is never obtained in a producing gas well. In formations of low permeability. such as those in the Hugo-ton Field, most wellhead working pressures approach stabilization closely enough to be used satisfactorily in the determination of a back-Pressure potential curve after flow periods of 24 hours. We shall therefore describe the backpressure curve calculated from ohserved rates of flow and working pressure at the end of 24-hour flow periods.

Jan 1, 1949

By J. J. Gilvarry, B. H. Bergstorm

The first part of this paper describes a new approach to the problem of energy relationships in fracture and comminution. The basic theoretical method used (as contrasted to previous empirical or semi-empirical approaches) is an attempt to verify or disprove the hypothesis of von Rittinger, In the second part of the paper, theoretical conclusions for the three headings of single fracture, plural fracture, and comminution are tested. Agreement between theory and experimental results is demonstrated. The classical hypotheses on the energy relationship in brittle fracture are those of von Rittinger-that the energy required is proportional to the new surface formed— and the alternative one of Kick.' Over the years, repeated attempts have been made to discriminate experimentally between the two assumptions, or to replace both of them, as has been done by Bond? However, no semblance of agreement seems to exist among different investigators. This paper is the first in a series essaying a new approach to the entire problem of energy relationships in fracture and comminution. The purpose of this program of research is to verify (or disprove) the hypothesis of von Rittinger. In contrast to previous empirical or semi-empirical approaches to the problem, the basic method adopted is theoretical. The prime aim is to determine the distribution function for fragment size theoretically and to verify the predicted distribution experimentally. This approach should make possible determination of the total surface of the particles by integration over the corresponding distribution function; one may hope that this direct procedure will yield more definitive results than use of gas adsorption or other methods. The surface area so obtained can then be compared with experimental energies, as obtained by Bergstrom, Sollenberger, and Mitchell,4 for example. It is convenient to distinguish between single fratture, plural fracture, and comminution. Single fracture of a brittle solid is defined as fracture by an external stress system which is removed instantly and permanently when fracture is initiated. Plural fracture consists of single fracture of the original specimen, followed by a sequence of only a few secondary fractures; comminution differs in that the sequence of secondary fractures consists of a large number of repetitive steps. For single fracture, Gilvarrys has given a rigorous derivation of the proper distribution function for fragment size, based on a closely defined physical model and deduced strictly by the laws of probability. Gilvarry and Bergstrome have compared the predictions of the theory with experiment, and have found excellent agreement. The prior work of Gilvarry and of Gilvarry and Bergstrom considers only single and plural fracture. The purpose of the present paper is to extend the discussion to the case of comminution. For prefatory purposes, the theory of Gilvarry for single fracture will be outlined. The considerations will then be applied to the cases of plural fracture and comminution. SINGLE FRACTURE The derivation is based on the Poisson law.7 For Points distributed with mean density y over a domain, this law states that the probability p(t)dt of one point lying in the range t to t + dt with none in the interval 0 to Ms p(t)dt = ye-yt dt. For this distribution to be valid, it is necessary and sufficient that the points be distributed at random, individually and collectively in the sense of ~~~.7 The former qualification requires that the position of One Point be independent of another, and the latter implies that the probability of a region containing a

Jan 1, 1961

By W. van der Knapp

A theoretical and experimental ana1ysis is given of the change in volume of a porous medium due to changes in external and internal pressures. The result enable one to deduce directly the effect of large increlents in stresses. It is shown that the three-dimensional representation of the volume of a nonlinear porous system is a cylindrical rules surface, the generating lines of Which have a slope according to the solid rock compressibility. Analytical expression for such a ,surface is very simple. It follows that the bulk and pore compressihilities depend only on the difference between pore fluid and external hydrostatic pressure, i.e., the effective rock frame pressure. From measurements on a number of sandstones as well as limestones it can be conduded that for practical values of this effective frame pressure the pore compressibility falls in the range between the compressibility of water and that of undersaturated crude oil. In general, pore compressibility is found to be higher the lower the porosity. For a particular limestone reservoir pore compressibility and porosity could he related by means of a simple empirical formula. Results from the described study find a direct applicxtion in material balance calculations and in prob1ems of liquid flow through porous media. Since the elastic constants of the rock bulk material enter into equatiorls for velocity of acoustic waves, the results are also of importance for calculating rates of wave propagation in these media. It is shown for sandstones with 1.5 to -10 prr cent porosity that laboratory measured porosities may differ from those lender reservoir conditions by about I per cent of their value. The difference for low porosity limestotnes can be of the order of 10 per cent. The subsidence caused by elastic deformation of reservoir rock is shown to he very small. General formulas relating pore and bulk volume to external and internal pressure can be derived, as has been shown by Biot1, Gassmann2 and Geertsma.2 Their theories describe the differential change in volume caused by a variation in pore fluid tension and external hydrostatic tension. It has been known for a long timei that the stress-strain relations of a packing of spheres are essentially noninear. Our experiments show that for almost all reservoir rocks which arc elastically deformable the relation between volume and loading is also a nonlinear one, resembling that of a packing of spheres very closely. The behavior of these media is often described by using the linear stress-strain relations for small incremental stresses and strains in the vicinity of a pre-stressed condition." The present study differs from these considerations in that it enables one to deduce the effect of large increments in the stresses directly. In analogy to the stress tensor for a homogeneous body the stress tensor for a porous body can be written in the conventional way: where the tensor components, o, and t, represent the total forces acting on the faces of a unit cube of porous material. This tensor can be split into two parts. a hydrostatic and a deviatoric part, as follows2,3,1 where the component of the hydrostatic stress system, a. equals The hydrostatic stress causes a change in volume whereas the deviatoric stress system gives rise to a change in shape only. One of the difficulties in describing changes of volume occurring in reservoir rocks is that the relations connecting these volume changes with external and internal stresses are essentially non-linear. A complete description requires the bulk volume, V,,, as well as the pore volume,V,,. to be known as a function of both the changes in pore fluid tension (v) and the hydrostatic tension on the outer boundary.

By O. P. Arora, M. Metzger, G. R. Ramagopal

Single-phase aluminum containing 0.0001 to 0.06 pct Cu was studied in strong acid, mainly through observations of hydrogen evolution. The strong influence of copper was exerted almost entirely through the imposition after a certain delay time of an auto-catalytic localized-corrosiott reaction. Additions of cupric ion to the acid produced lower accelerations. The significance of the quantity and distribution of copper was discussed, and the implications for intergranular corrosion and neutral chloride pitting were indicated. AN investigation of intergranular corrosion in single-phase high purity aluminum exposed to hydrochloric acid indicated the copper content of the metal to have an influence on corrosion at lower levels than previously suspected.&apos; The work reported here was a closer examination of the action of copper but dealt with general corrosion to gain the advantage of having a continuous measure of corrosion through the volume of hydrogen evolved, the reduction of hydrogen ion to hydrogen gas being the principal or only cathode reaction in strong hydrochloric acid. Previous work on the hydrochloric acid corrosion of aluminum was sometimes insufficiently structure-conscious and the need for care in evaluating it arises from the low solubility of the iron impurity,&apos; and of some alloying elements, and the known or possible presence in many of the compositions studied of second phases leading to greatly increased corrosion rates.3 These increases are attributed to the presence of low hydrogen-overvoltage cathodes provided by the second phase.3&apos;4 For the present single-phase work, a few studies which used high-purity base material and small copper additions5-&apos; provide the essential information most unambiguously. The corrosion rate was shown to be increased markedly by the introduction into the acid of small quantities of the ions of copper (and of certain other metals) which cement on the aluminum and provide cathodes of low overvoltage.5 When there was sufficient copper in the aluminum, the same result was produced during the course of corrosion leading to a rate which increased with time as the reaction was stimulated by one of its products (autocatalytic reaction). In 2N (7pct) HC1, an accelerating rate was observed at 0.1 pct Cu but not at 0.01 pct.5,7 The present work dealt with corrosion rate and morphology and their correlation with the quantity and distribution of copper catalyst for copper contents from 0.0001 to 0.06 pct. PROCEDURE A lot of high-purity aluminum containing 0.0021 pct Cu, 0.001 pct Fe and 0.003 pct Si (Alloy A) was alloyed with copper to yield aluminum containing 0.014 pct Cu (B) and 0.06 pct Cu (C). Later it was found necessary to include the lower copper Alloy K which contained 0.0001 pct Cu, 0.0004 pct Fe and 0.0004 pct Si. The upper limit for any other element can be confidently estimated as 0.0005 pct. No element other than copper appears to be present in quantities sufficient to have an effect on general corrosion as great as the observed effect of the copper in A, B, and C. The only other heavy metal detected by spectrographic examination was silver (< 0.0001 pct). The acid was made up from a selected lot of 37 1/2 pct CP hydrochloric acid containing 0.1 ppm heavy metals (mainly Pb), 0.05 ppm Fe, and < 0.008 ppm As and from water distilled from 1 megohm-cm demineralized water and believed to have contained negligible quantities of heavy metals influencing corrosion. Acid strength was adjusted to within 0.05 pct HCl of the stated value by using precision specific gravity measurements. Test blanks 10 by 41 mm were sheared from 1.65-mm cold-rolled sheet. Edges were finished by filing. The blanks were annealed in air at 645°C for 24 hr in alundum boats and rapidly water quenched. The anneal is thought to have produced a substantially homogeneous solid solution—for iron, copper, or silicon, for example, the annealing temperature was 200°C or more above the solvus-and the quench is considered to have preserved the high-temperature structure except for the condensation of lattice vacancies into dislocation loops.&apos; The 0.06 pct Cu alloy did not appear unstable in respect to slow precipitation reactions at room temperature since two pairs of tests failed to show significant differences between specimens heat treated 3 1/2 years earlier and specimens heat treated 1 or 2 days before.

Jan 1, 1962

By G. M. Gay, R. L. Perrine

This paper describes a method of numerical computation for three-dimensional, unstable, miscible displacement behavior useful for heterogeneous systems, as well as for more ideal conditions. In the method, flow equations are first linearized by a perturbation approach. The basic flow process is separated and a solution for its behavior readily obtained. The remaining problem of deviations from the basic flow caused by non-ideal conditions is then subjected to numerical analysis. Results obtained from use of the method are also presented. Although conditions assumed in the test cal-culations were severe, results show the type of dispersing flow expected and appear quite satisfactory. The method has eliminated or reduced in importance problems of oscillating values near steep fronts, excessive computer smoothing, etc. A unique advantage of the method is that the source of variations from ideal behavior can be observed. ,The one serious drawback results from the algebraic complexity of the perturbation approach, and the need for second-order terms to be retained in calculations of interest. Fewer array points are available and more computer time is required than would be desired. However, these difficulties are also experienced with other approaches to the solution of three-dimensional dispIacement problems. INTRODUCTION Prediction of behavior of the miscible displacement process within a porous medium for any system of engineering importance is plagued by a number of difficulties. With typical fluid properties the displacing fluid has the lower viscosity, and there is a natural tendency toward flow instability. The problem of predicting instability is compounded by the fact that every real system is heterogeneous. Permeability will vary from point to point — not entirely systematically, and yet not in a random fashion leading to a readily defined average. Furthermore, permeability is unlikely to be isotropic. These permeability properties, which are characteristic of any real porous medium, accentuate the effects of flow instability. Another factor to be considered is that flow dispersion accompanies the displacement process in a porous medium. Mechanistically, dispersion is due to the fact that flow between any two points in the medium follows multiple tortuous paths, each characterized by slightly different flow properties. The fact that the coefficient characterizing dispersive properties of the real medium JS a tensor with variable elements complicates matters. Further difficulties arise because the parabolic influence, while not negligible, is small. Thus, there is a tendency toward steep and uneven displacement fronts along which dispersion smoothing must be represented accurately. Yet frequently used numerical analysis schemes may tend toward either instability, oscillation or excessive smoothing, none of which gives the desired accurate picture of flow behavior. Thus, while many experimental and analytical studies of the process have been made, predictions of actual performance are still subject to considerable uncertainty and possible improvement. This paper reports on part of a study which attempts to develop improved methods for solution of the flow equations describing the miscible displacement process. Of necessity, the calculations were performed on a large digital computer. A three-dimensional system is represented, and the coefficients defining dispersion and permeability can be varied in a manner representative of a real system. Because the differential equations are nonlinear, a first step is to linearize them using a perturbation approach. Objective of this approach has been to separate the basic flow which would occur in an ideal system and for which predictions are readily made. Then the numerical analysis, which is necessarily complex, is directed toward solution of the perturbation equations. These equations state how the basic flow is modified

Jan 1, 1967

By Van Der Muelen P. A.

SOME time ago, the writer received from W. 0. Borcherdt, Manager of the, works of the Bertha Mineral Co. at Austinville, Va., several specimens of a dense yellowish slag-like material, containing cavities showing clear needle-like crystals, and representing a product formed by the fumes of zinc oxide and basic lead sulphate from the oxide furnace, attacking the firebrick lining of the flues. The matrix, which is of a greenish-yellow color, is made up chiefly of an amorphous glass, carrying numerous crystals, the latter being well-defined, clear, and transparent, but colored a very light yellow by small amounts of ferric oxide. A number of the crystals were carefully separated from the matrix and tested with the blowpipe. They have a fusibility of about 2, yield a globule of lead and a coating of lead oxide before the blowpipe on charcoal, and contain some zinc. Their hardness is 5 to 6. 1 Fig. 1. Under the microscope the crystals appear to be made up of a prism, with narrow faces, and a fairly large pinacoid, but unfortunately show no terminal faces (Fig. 1). Between crossed nicols they gave parallel extinction in. the three planes containing the crystallographic axes, and are therefore orthorhombic. The index of refraction is high. A number of the best crystals were measured on the reflection goniometer, the angles given below representing the average of measurements made on different individuals. mm = 46° 52&apos; mb = 66° 34&apos; a:b:c = 0.4334:1 :?

Jan 5, 1917

By W. A. Oates, J. A. Lambert, P. T. GaIIagher

Monte Carlo methods have been used to compute the arrangements of interstitial atoms dissolved in tetrahedral sites in bcc lattices. It is assumed that the presence of an interstitial atom "blocks " a certain number of neighboring sites and prevents their occupancy. Sites "blocked" by more than one filled site are allowed for. The computed values of. the mean occupation number (defined as the ratio of the total number of sites blocked to the number of solute atoms are used to calculate the configurational entropies of the solutions. These entropies are compared with those resulting from previous theoretical studies of this problem and also with available experin~ental data for the p Zr-H, Nb-H, V-H, and Ta-H systems. Evidence is also given that the "blocking" explanation of low limiting compositions in these systems, rather than this being due to initial limitations on the number of sites available, is probably correct. THE ideal partial configurational entropy of mixing of an interstitial solute in a metal is given by: where p is the number of interstitial sites per metal atom and Xi is the atomic fraction of the interstitial. For the bcc lattice. which we shall be concerned with in this paper, the interstitial positions are shown in Fig. 1. It can be seen that for the tetrahedral sites, p=6. whereas for the octahedral sites, p = 3. Different emphasis has been placed on the relative importance of energy and entropy effects in determining deviations from ideality in interstitial solid solutions. In some cases the same system, e.g., Fe-C, has been described by the contradictory regular and athermal solution models indicating that the enthalpy and entropy functions, derived from equilibrium data, are frequently not accu.rate enough to differentiate between these treatments. However, for certain metal-hydrogen solutions the equilibrium data is available over sufficiently wide ranges of temperature and composition to permit a reasonably accurate determination of the compositional variation of the heats and entropies. Hoch&apos; has attempted to interpret the results of interstitial solid solutions in terms of a regular solution model. In the case of the Ta-H system where 13 = 6, this model entails fitting the experimental relative partial entropies of solution, asH, to the equation: where ASgs is the relative partial excess entropy of solution of hydrogen. Hoch found that the results of Mallett and Koeh1 could be fitted to this equation with an approximately constant value of AF up to XH = 0.25. However, it is apparent from the solubility isotherms in this system which become asymptotic to the composition TaH that, since (Xh /6 - ~Xh ) becomes infinite only at TaH6, it is necessary that AS<&apos; tends to infinity at TaH. In other words, the low saturation composition of TaH, instead of the anticipated TaH,, eliminates the possibility of applying regular solution theory to such systems. Rather large negative excess configurational entropies must exist at higher hydrogen concentrations in order to explain the lower saturation values. To account for these low limiting compositions and excess entropies two distinctly different approaches have been followed. Rees and many others1-l2 have assumed that not all interstitial sites are crystallographically equivalent with respect to the interstitial addition; that is, in Eq. [I] p is less than the value anticipated from geometrical considerations. To describe, say, a bcc metal-hydrogen system with a limiting composition of MH by this approach one would consider that p = 1 in the first instance instead of p = 6.&apos;j3 In some cases, nonintegral values of B have been taken in order to improve the fit with the experimental data over limited ranges of composition. The other approach which has been used to explain the low saturation compositions is to assume that, although all sites are available for occupancy, strong repulsive interactions exist between the neighboring interstitial atoms, and hence occupancy of any site excludes or blocks a certain number of neighboring sites from being occupied. Earliest treatments of this concept considered the exclusion of an integral number, of nearest-neighbor sites from being occupied at all concentrations. In this case, the partial configurational entropy is given by: These early treatments failed to allow for the overlap of the blocked sites which will arise at all but the very lowest concentrations. More recently attempts have been made to calculate the effect of this decrease in the number of blocked sites on the configurational entropy. Using the quasichemical treatment of interstitial solid solutions as given by Lacher and assuming that an infinite repulsive interaction energy existed between the solute atoms. atom obtained an approximate configurational entropy applicable to the blocking with overlap case:

Jan 1, 1970

By A. T. Chatas

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

Jan 1, 1967

By K. Tangri, M. Chaturvedi

On water-quenching from within the (a + ß) phase region Zr-2.5 Nb and Zr-2.5 Nb-0.5 Cu alloys can undergo w transfirmation. This transformation has been attributed to the enrichment of ß Zr phase, at the solution-treatment temperature, by the solute atoms. The ß—w transformation is accompanied by a hardening 01 the alloys and in the optimum condition w-bearing specimens of the ternary alloy are about 14 pct harder than those bearing only a&apos; phase. On aging crl reverls thermally to more stable phases, which results in a considerable softening of the specirnens, as follows: w quenched -waged + ß Zr (enriched) waged - a + ß Nb The high-temperature bcc ß Zr in certain zirconium alloys can transform to a metastable w phase which has a primitive hexagonal structure. In the Zr-Nb alloy system w phase has been observed in the composition range of 7 to 15 wt pct Nb by a number of investigators.&apos;-&apos; In this composition range w can form by water-quenching ß Zr or by aging. known as athermal and thermal w, respectively. Hatt and Roberts&apos; have found that thermal w has a constant c/a ratio of 0.622 +0.002 regardless of the composition and tha: the lattice parameters are of the order of a = 5.02A and c = 3.0 ?. The actual lattice parameters depend upon composition and the temperature of aging. The lattice constants of athermal w. however, depend upon alloy composition.3 The presence of w phase in zirconium alloys can cause substantial strengthening which has been attributed as probably due to the coherency between w and parent ß Zr phase. Robinson et al.&apos; found that the presence of w can raise the strength of Zr-5.0 Nb-2.0 Sn and Zr-5.0 Mo-2.0 Sn alloys to about 170.000 psi: however. the ductility is considerably reduced. With the growing use of Zr-2.5 Nb and Zr-2.5 Nb-0.5 Cu alloys in the nuclear reactors, and the significant effect of w transformation on the mechanical properties of zirconium alloys studies were undertaken to study the w transformation in Zr-2.5 Nb and Zr-2.5 Nb-0.5 Cu alloys. This paper is concerned with the formation and reversion of w phase in these alloys. 1) MATERIALS AND EXPERIMENTAL TECHNIQUES The following two alloys were studied: 1) Zr-2.5 Nb-0.5 Cu (referred to as the ternary alloy): 2) Zr-2.5 Nb (referred to as the binary alloy). The chemical analysis of the alloys is given in Table I. The alloys, in the form of 3/8-in.-diam rods and 1/8-in.-thick sheets, were supplied by the Chalk River Nuclear Laboratories of the A.E.C.L. Initial fabrication for preparation of various specimens was carried out by cold rolling and swaging with intermediate anneals at 1000°C. All the heat treatments were carried out after the specimens were wrapped in zirconium foils and encapsulated in silica tubes under a vacuum of 5 x 10-% m of Hg. For optical metallography and hardness measurements disc specimens + in. diam by + in. thick were used. After the required heat treatments the specimens were mechanically and then chemically polished in a 45 pct HNO3, 45 pct H20. and 10 pct HF solution. The hardness measurements on the chemically polished specimens were carried out on a Vickers hardness tester using a 10-kg load. For each specimen at least fifteen indentations were made in order to obtain a representative hardness value. The phase identification and structural analysis were carried out using X-rays and electron diffraction techniques. Wires of 1.5 mm diam reduced to 0.12 mm diam by chemical etching were used for making Debye-Scherrer powder patterns using copper K-a radiations in 114.6-mm-diam camera. Thin films for transmission electron microscopy were prepared by electropolishing heat-treated a by 7 by 0.005-in.-thick strips using a modified Bollmann-Window technique. The 10 pct perchloric acid-90 pct methyl alcohol electropolishing bath was kepl at -5OCC and polishing was done at 5 to 10 v. The thinned specimens were washed in ethyl alcohol at -30" to -40°C and dried between filter papers. The thin films

Jan 1, 1970

By J. B. Newkirk

ABOUT twenty-seven years ago W. bergl discovered that interesting detail could be seen in an X-ray diffraction spot made with a rock-salt crystal if the recording photographic film were held very close to the reflecting crystal. He also showed that much of this detail could be explained in terms of the deformation structure of the specimen. In the same year, 1931, Fox and parr2 reported that X-ray reflections from a quartz crystal were much more intense while the crystal was undergoing piezoelectric oscillation than when it was at rest. This observation was soon confirmed (and explained in terms of X-ray extinction) by Barrett3 in the United States and by Nishikawa et a1.,4 in Japan. During the next few years Barrett refined Berg&apos;s method for making X-ray diffraction micrographs and applied it to the study of metallic crystals. He reported this work in the 1945 annual lecture of the AIME.5 Since 1945 relatively little use has been made of the Berg-Barrett (B-B) method for making diffraction micrographs, though considerable interest in investigating subcrystal structures by other methods has developed. Some of these methods depend upon evaporation6 or chemical etching7 and, therefore, suffer from the occasional intrusion of subtle artifacts such as those recently reported by Hooker for leada and Welsh for aluminum.9 Other methods require the introduction of foreign elements which "decorate" the substructure in some way, thereby rendering it detectable by direct" or indirect1&apos;. visual examination. Still others depend upon complex and often highly questionable analyses of X-ray diffraction line profiles. There have also been a number of related X-ray techniques, such as the Lambot,12 Schulz13 and Merlini and Guinier,14 which have proved to be quite useful for specific applications. Some of these methods give a picture of the substructure which is more or less directly related to the shape and the topography of the specimen. The Berg-Barrett method carries this desirable feature to the limit. Closely related to the B-B method are three photographic X-ray techniques recently described by Lang.15 These methods give beautiful pictures of the internal structure of crystals but are limited by the facts that the required apparatus is rather complex, exposure times are long, and the specimen crystals must be thin. An experimental appraisal which I have recently made of the original B-B method has shown that this simple technique is not only inexpensive and convenient, but can give more information than was apparently realized by its inventors. For example, first-order subgrains in iron-silicon crystals were clearly delineated and their relative orientation subject to description. Also defect structures could be seen within the subgrains and in some materials slip bands could be clearly seen with respect to the topography of the rest of the specimen. This much had already been described by Barrett. However, in the present study it was found that the actual sites of individual dislocations which emerge upon the crystal surface could be located. Furthermore, by a simple procedure, the direction (but not sign) of the dislocation Burgers&apos; vector could be defined as well.

Jan 1, 1960

By J. B. Cohen, O. F. Kimball

Aging of specimens of Au-20 at. pct Ni and Au-40 at. pct Ni has been followed with electron microscopy, resistance, magnetic measurements, and X-ray diffraction. The periodic structures observed by electron microscopy in Au-Ni alloys by Fisher and Embury are shown to be due to nickel-rich prepre-cipitates. Assuming this alloy system undergoes spinodal decomposition, our observations and available data on thermodynamic properties and elastic constants indicate that the gradient energy coefficient in this system is of the order of 10-5 ergs cm-1 . Shifts of X-ray peaks in this system have been found to be due primarily to residual stresses, not faulting. A large effect of applied stress during aging was detected by X-rays and electron microscopy for the Au-20 at. pct Ni alloy. The presence of a preprecipitation stage during aging of gold-rich Au-Ni alloys for short times below 300°C has been indicated by several previous investigations. Changes in electrical resistivity,1"3 elastic modulus,4 superparamagnetisrn,5 positions of X-ray diffraction peaks,6 and hardness7 have been studied. Electron diffraction and microscopic studies by Fukano8 and Fisher and Embury9 respectively showed that the precipitation process resulted in the formation of periodic composition fluctuations—perhaps spinodal decomposition. MOSS10 detected a satellite with X-rays in quenched specimens, but the wavelength was much smallerthan that found by Fukano. Recently Woodilla and Averbac11 detected a modulation similar to that found by Moss, but with electron diffraction. They suggested that the microstructures observed by Fisher and Embury may be simply due to a Moir&apos;e effect from a back deposit of gold during electropolishing. Satellites in electron diffraction can be enhanced by small bends or tilts of a foil. Even if this were not the case, asymmetric satellites, as are found in Au-Ni are difficult to analyze. They are principally influenced by the large differences in atomic size in this system, and obtaining information on the associated composition fluctuations depends on assumptions as to whether or not there are isolated clusters or composition "waves", and in the latter case whether or not the entire specimen has decomposed. Assuming that decomposition was complete in alloys aged at low temperatures, Woodilla and Averbach suggested that the satellites indicated a composition which varied about the average value by only± 5 pct, much less than that predicted from the phase limits, Fig. 1; they attributed this small variation to a cessation of the reaction caused by the loss of mobile vacancies. In Ref. 12 it was shown that vacancies are trapped by the preprecipitate in this system. Furthermore, after reversion, aging proceeds at a much slower rate than immediately after a quench from the solution temperature,5&apos;11,13 indicating the importance of vacancies in the kinetics. In some senses Au-Ni is a favorable system for studying preprecipitation and reversion mechanisms because of the simplicity of the equilibrium diagram14 and the availability of good thermodynamic15 and diffusion data.16,l7 The calculations necessary to relate the decomposition process to classical nucleation and growth or to spinodal decomposition, as well as those necessary for understanding reversion mechanisms, can, therefore, be made. As an example, the position of the chemical spinodal, calculated from the thermo-dynamic data of Sellars and Maak15 with the method described by Rundman and Hilliard,18 is included in Fig. 1. It is compared to the chemical spinodal calculated from the equilibrium miscibility gap by Cook and Hilliard&apos;s method.l9 Included is the coherent spinodal suggested by Golding and MOSS," using measured elastic constants for their calculations and the same ther-modynamic data. The spinodal estimated by Woodilla and Averbach,11 by determining the temperature limits for the observation of the satellites discussed above, is also given. Actually this last curve is open to question on the following grounds: As a phase boundary is approached the spacing of particles or the wavelength

Jan 1, 1970

By E. J. Nielsen, P. R. Strutt

A study has been made of the progressive changes in the distribution of etch-pit structures occurring during high-temperature creep in copper + 7 wt pct Al single crystals oriented with a [113] tensile axis. The two equally stressed glide systems with the highest Schmid factor would be expected to form subboundaries of the type predicted by Kear.2 The alignments of etch-pits on sections parallel to different (111} planes consistent with these types of boundaries were not observed. However, they were consistent with planar subboundaries (on a macroscopic scale). From an analysis of Amelinckx1 it may be shown that stable cross-grid dislocation boundaries may form in the primary slip planes. These boundaries form when dislocations with a Burgers vector not in the slip plane move into the plane by combination of climb and glide. THE geometry of subboundaries formed by the interaction of dislocations of two glide systems has been analyzed by Amelinckx,1 and the particular types produced by deforming fee crystals are predicted by ear.&apos; In this paper types of boundaries which may be formed when climb as well as glide occur are discussed as this is relevant in high-temperature creep. It is assumed in the present investigation that the etch-pits observed in Cu + 7 wt pct A1 on surfaces parallel to {111} planes delineate the sites of dislocations. Although there is no direct evidence for this previous work on a-Cu/Al single crystals by Mitchell, Chevrier, Hockey, and Mon-aghan,3 would show this assumption to be reasonable. The alignments of etch-pits which form during creep are studied on sections parallel to each {111) plane. It is then deduced that these alignments are consistent with a specific type of planar subboundary. The Cu + 7 wt pct A1 single crystals had a [113] tensile axis and Fig. 1(a) shows schematically the relation of the slip planes and slip directions (as represented by tetrahedron ABCD) with reference to the tensile axis. The two equally stressed glide systems with the maximum Schmid factor namely ß-AD and (a-BD, from the analysis of Kear,2 would be expected to form the boundaries shown in Fig. l(a) and (b), also Fig. 5(a) and (b). EXPERIMENTAL PROCEDURE The a-Cu/Al single crystals were grown and annealed in a "gettered" argon atmosphere. Chemical analysis showed the aluminum content to be uniform in each crystal and the difference between crystals was maintained to an accuracy of ± 0.25 wt pct. The initial dislocation density and mean subgrain diameter after annealing was -106 cm-2 and 250 µ, respectively. Surfaces parallel to (111) planes were produced by specially developed electrolytic machining processes. The {111} faces were next electropolished for 5 min in a solution consisting of 25 g chromium trioxide, 113 ml glacial acetic acid and 40 ml water; the applied potential was 8 v. Dislocation etch-pits were revealed using l an etchant described by 1 ml bromine, 45 ml HCl, and - 250 ml water. RESULTS In crystals strained into secondary creep at higher stresses (443 and 750 g - mm-2 at 650° C aligned rows of etch-pits parallel to slip plane traces were evident in sections parallel to the (1111, (ill), and (111) planes, see Fig. 3. As well as the longitudinal alignments in Fig. 3, well formed randomly oriented arrays indicative of an equiaxed subgrain structure are evident. At the lower stresses (100 to 230 g . mm-2) only an equiaxed structure formed during creep. The sections in Fig. 3 are from a crystal crept for 70 hr at 650°C with a CRSS of 443 g.mm-2. Two identically oriented crystals were also deformed at the same temperature and stress for 5 min and 4 hr. In the crystal crept for 5 min, the etch-pits were randomly distributed with no tendency for directional alignment, see Fig. 2(a). As shown in Fig. 2(b) aligned arrays were evident after 4 hr creep but they were not nearly so well defined as in Fig. 3. The alignments (parallel to the arrows) in Fig. 3 are consistent with the existence of boundaries in the two main slip planes a and ß. The way in which this is deduced is seen by reference to Fig. l(c), where the existence of boundaries in the a and ß planes is verified by sectioning parallel to a,ß, and d. The (111) and ß(111) planes intersect the d(111) plane along BC [101 ] and AT [011] and alignments parallel to [101] and [011] are clearly evident in Fig. 3(c) in a section parallel to the d(111) plane. Similarly the a, and ß planes in Fig. l(a) intersect each other along DC [110] and hence there will be an alignment parallel to [110 ] in sections parallel to the a-plane and the ß-plane; this is evident in Fig. 3(a) and Fig. 3(b). It is interesting to note that alignments of etch-pits consistent with the boundaries predicted by Kear2 were not observed; see Figs. l(a) and l(b). The geometry of boundaries in {111} planes as shown in Fig. l(c) is discussed later. In Fig. 4(a) the individual etch-pits are resolved and the alignments are exactly parallel to the slip trace direction [101]. However, in some areas alignments deviate away from the slip trace direction by as much as 10 to 15 deg, this is evident in Fig. 4(b), and in Fig.

Jan 1, 1969

By N. Jaffe, J. E. Dorn

This investigation is concerned with the possibility that the creep resistance of metals might be reduced as a result of recovery at the creep temperature when the applied stress is reduced. For this purpose the effect of stress on the recovery of the creep behavior of pure polycrys talline aluminum was investigated in the cross-slip temperature range where the activation energy for creep is 28.000 cal per mole. Before recovery all specimens were given the same substructure by creeping them under a stress of 7800psi to a strain of 0.218. Recovery was then accomplished by holding the specimens at the creep temperatures of 327" and 337°K under a series of reduced stresses for various periods of time. When the original stress was reapplied the specimens crept more rapidly than when the original stress was first removed, illustrating that recovery did take place. The extent of recovery could be indexed in terms of the UNDER certain conditions of test, the creep rate decreases over the primary stage of creep often by as much as a factor of 10~B. Since the external variables of stress and temperature are held constant, such pronounced and significant changes in the creep rate must be attributed to various substructural modifications that are introduced as a result of creep itself. At the secondary stage of creep, over which the creep rate remains constant and the substructure no longer changes, a balance must be achieved between the rate of increase in creep resistance due to creep straining and the rate of recovery of the substructural modifications. The general validity of these concepts on substructural changes has already been shown,&apos; including the corollary that a more creep-resistant substructure is introduced under a more rapid creep rate arising from a higher applied stress.&apos; In spite of these well-documented facts, very little progress has been additional strain undertaken by the recovered specimen to achieve the same strain rate as the unrecovered specimen. It is significant that the higher recovery stresses gave the most rapid rates of recovery. Over the high range of reduced stresses, where some creq took place during recovery, the activation energy for recovery equaled that for creep, namely 28,000 cal per mole. But for stresses during recovery below 0.80 of the original creep stress, no creep was detected during recovery and an activation energy for recovery of about 72,000 cal per mole was obtained. Whereas the activation energy of 28,000 cal per mole suggests that cross-slip of dislocations is responsible for this recovery, the mechanism for the 72,000 cal per mole process is not known. This activation energy, however, is about twice that for self-difhsion. made in determining, in satisfactory detail, precisely what substructural changes so effectively modify the creep rate.&apos; Up to the present, electron microscopy has failed to uncover, in a satisfactory quantitative manner, the mechanistic processes involved. Since, however, creep, at least at moderate and high temperatures, arises from one or more types of dynamic recovery processes, it should be possible to shed additional light on the total problem by investigating the recovery of creep resistance. The importance of such an investigation becomes apparent when the fact is recognized that there has yet been no direct experimental proof that the rate of softening due to substructural changes as a result of recovery equals the rate of hardening due to creep-induced substructural strengthening over the secondary stages of creep. Furthermore the problem of what happens to the creep-strain induced strengthening when the material recovers either at zero or some reduced operating stress has significant engineering consequences. It was for these reasons that the present investigation on the recovery of creep behavior was initiated. Inasmuch as a substantial background has now developed on the creep of high-purity aluminum, this material was again selected for the current investigation. Furthermore extensive evidence is now

Jan 1, 1963

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 Lyman Johnson

Single crystals qf copper in "single slip" orientatiorzs have been deformed in compression. During defortnation all of the independent deformation parameters have been measured. These parameters consist of thefive strain components and three components descrihing the lattice rotation. By a finite strain analysis these pararmeters , forrming a deformation gradient martrix, are related to the amounts of slip on each of the twelve slip systems. The results show that the amount of secondary slip is about equal to the amount of primary slip. This is an order of magnitude larger than has been believed previoutsly. ACCORDING to early theory and experiments, when a single crystal of a fcc metal is deformed in tension or compression it should deform by slip on only one slip system until the stress axis reaches a symmetrical orientation.&apos; However. the observation of a large increase in the secondary dislocation density during ..single slip" makes it clear that some slip does occur on secondary systems. Knowledge of the amount and distribution of this secondary slip is essential to a complete understanding of the mechanisms of single-crystal deformation. Ahlers and Haasen 2 and Mitchell and Thornton1 have tried to detect the amount of secondary slip in single crystals of silver and copper, respectively. Each simultaneously measured the angle A, between the tensile axis and the primary slip direction and the length 1 of a gage section of the specimen after incremental amounts of deformation in tension. The measured A, was then compared with the theoretical single slip angle hp. given by sin Ap = j sin . hO where ?o was the initial angle between the tensile axis and the primary slip direction and lo was the initial gage length. In both sets of experiments a small but systematic difference between ?e and ?p was found. This difference must be due to the occurrence of secondary slip. However, as Mitchell and Thornton1 pointed out. nothing quantitative can be said about the amount and distribution of this secondary slip from the measurements that they made. The reason that no quantitative conclusions could be made is because no unique solution for the distribution of slip on the twelve fcc slip systems can be determined from only two measured deformation parameters such as A and 1. There are, in fact, eight independent macroscopic deformation parameters that can be measured when a single crystal undergoes a homogeneous deformation. Physically these can be thought of as the five finite strain components and the three angles describing the crystal lattice rotation. All eight of these parameters were measured by Taylor4,5 for aluminum deformed in tension and compression. At that time the concern was to show that slip occurs on {111 (110) systems in fcc metals, and the mathematics were not available to determine what slip distributions were compatible with the measurements. In this paper the mathematics6,7 are developed that allow the slip distribution to be determined from these measurable macroscopic deformation parameters. The analysis is applied to the measurements of the strain and lattice rotation of copper single crystals deformed in compression. The results show that the amount of secondary slip is an order of magnitude larger than had previously been thought. CRYSTALLOGRAPHIC DESCRIPTION OF A HOMOGENEOUS DEFORMATION The deformation of a solid body can be represented by a transformation matrix F that transforms the un-deformed state into the deformed state. Consider a vector X connecting two material points in the unde-formed material and the vector x connecting the same two material points after deformation, where both vectors are referred to the same set of Cartesian axes. The final vector x is related to the initial vector X by the equation: X = FS. [2] Eq. [2] can be considered as the equation defining F, which is called the deformation gradient matrix. Its components are: If the deformation is homogeneous, the transformation is linear and the components of F are constants. Using subscript notation, if P is the unit vector in the initial direction of a material line, the components of the unit vector p in the direction of the same material line after deformation are given by:

Jan 1, 1970

By Mayumi Someno, Kazuhiro Goto, Masahiro Kawakami

The apparent rates of decarburization of liquid alloys of Fe-C, Fe-C-S, Ni-C, and Co-C systems and the rate of oxidation of solid graphite with pure carbon dioxide gas and with gas mixtures of carbon monoxide and carbon dioxide have been measured in the temperature range of 1000° to 1600°C. The cotnposition of carbon dioxide and carbon monoxide gas at the reaction surface has been measured by oxygen concentration cells with the ZrO,-CaO solid electrolyte. 1) The apparent rates of the carbon removal are essentially the same for all the cases of solid graphite, Fe-C, Fe-C-S, Ni-C, and Co-C systems under the same experimental conditions. 2) The apparent rates are independent of the carbon content in the high carbon concentration range but very much affected by the flow rate and the gas composition of the CO-CO2 reactant gas mixture. The ratio of the gas consumed by the reaction to the total quantity of the supplied gas is very large under the present experi~nental conditions. 3) There is a concentration gradient of&apos; carbon dioxide in the vicinity of the reaction surface and the content of CO, becomes extremely small at the reaction surface. 4) A large time fluc-tuation of the gas composition was observed. This jluctuation implies the presence of unstable flow in the gas phase in the vicinity of the reaction surface. THE decarburization of molten steel by an oxidizing gas or by slag may be one of the most important chemical reactions in steelmaking processes. Nevertheless, the kinetics of this heterogeneous chemical reaction do not seem to be well-solved even with the previous studies. Although the conditions for the reaction in steelmaking processes are quite different from those in the laboratory scale, some critical experiments may give information on the mechanism of the decarburization. From the previous work,&apos;-&apos; it is known that the rate of the decarburization is independent of the carbon content in liquid iron with more than about 0.2 wt pct C when the oxidizing gases are supplied to the surface of liquid Fe-C alloys on a laboratory scale. Two rate-controlling steps have been proposed for the decarburization of liquid iron with the high carbon content: one is the surface reaction control proposed by Swisher and Turkdogan;&apos; the other is that the rate is controlled by the gaseous diffusion through the gaseous stagnant layer. proposed by Baker. Warner, and Jenkins.7 and also by .Ito and Sano.2 In the present study, some experiments have been carried out for the evaluation of these rate-controlling steps in the decarburization of liquid iron with high carbon content. The apparent rate of decarburization of liquid iron has been compared with the rates of carbon removal of liquid Ni-C, Co-C, and solid graphite under the same experimental conditions. The composition of carbon dioxide and carbon monoxide gas at the reaction surface has been measured by oxygen concentration cells. I) EXPERIMENTAL PROCEDURE Fig. 1 shows the schematic diagram of the reaction chamber. Solid graphite and liquid metals were contained in an alumina or magnesia crucible of 32 mm ID and 35 mm in height. The samples were heated by high-frequency induction and the temperature was measured by the calibrated optical pyrometer. The temperature was held constant to within 10°C. The re-actant gases were supplied to the surface of the samples through the quartz tube of 8.0 mm ID. The distance from the end of the quartz tube to the surface was 20 mm. The block of high-purity graphite was cut and shaped to the inner profile of the crucible. The height of the shaped graphite was 18 mm, which corresponded to the depth of the liquid iron of 100 g. About 100 g of Fe-C alloy (4.20 to 4.40 pct C), Ni-C alloy (1.84 pct C), Co-C alloy (1.85 pct C), and Fe-C-s alloy (4.35 pct c, 0.5 or 1.0 pct S) were melted in the crucible. The reactant gases were pure CO, and gas mixtures of CO-CO,: the flow rates were controlled by capillary flowmeters with bleeders.

Jan 1, 1970

By Ben-Zion Weiss, Melvin R. Meyerson

Chromium diffusion coatings on commercial Armco iron lead to carbide precipitation at the grain boundaries in and below the coatings. High compressive stresses are introduced into the coating and, as a result, tensile stresses are introduced into the base material below the coating. In coated samples, fatigue cracks form at the grain boundaries below the coating after only a limited portion of the total lifetime (5 to 10 pct). Residual tensile stresses and a stress concentration caused by precipitated carbides seems to be responsible for early crack nucleation. Stage I of Propagation may be divided into two substages: I-a in which the crack Propagates along grain boundaries, and I-b in which the crack propagates along slip boundaries according to a shear mode. In uncoated samples, the cracks form at slip bands after 40 to 50 pct of the total lifetime. In coated samples the Propagation process takes longer than in uncoated samples because of the moderate rate of crack extension until the crack breaks through the coating. The chromium-diffusion coating causes little if any increase in fatigue life. CHROMIUM diffusion is one of the most popular processes used to apply coatings to iron alloys. This popularity stems from the fact that the diffusion treatment is comparatively easy to carry out, and it improves certain surface properties such as corrosion and wear resistance as well as some other properties. Very little effort has been devoted to the investigation of the effect of chromium-diffusion coatings on mechanical properties and particularly on fatigue properties. Since the chromium coating usually represents a small fraction of the total volume of the base material it is generally assumed that the macro properties are dependent principally on the base material rather than the coating. Insofar as fatigue is concerned, there are indications that the fatigue life of the chromium coating is not less than that of the basic material and the fatigue properties are generally not reduced by it.&apos; But chromium diffusion causes drastic structural and chemical changes in the region of the material surface2-4 which introduce additional residual stresses. Such changes must affect fatigue crack nucleation and may sometimes affect the initial stages of propagation. Furthermore, it is conceivable that the secondary stages of the crack propagation could be influenced by some irreversible structural changes in the core of the material, which stem from the long-time heat treatment at high temperatures required to apply the coating. This paper analyses the structural and compositional changes produced in commercial Armco iron by a chromium-diffusion coating and the effect of the coating on fatigue crack initiation and propagation. EXPERIMENTAL PROCEDURES The investigations were carried out with commercial Armco iron. The chemical composition of the iron, the number of specimens tested, the grain size produced during chromizing and heat treating, and the thickness and microhardness of the chromized layers are shown in Table I. Samples were chromized by a gas diffusion process. The temperature and time of the diffusion coating process were chosen according to a previous grain growth study. It had been found that at a temperature of 970°C and a holding time of between 3 and 9 hr there are practically no changes in grain size. To describe the grain size more accurately the ASTM method5 was supplemented by a statistical analysis of the mean volume diameter as outlined by Fullman6 and The The mean volume diameter was deter- mined from D = p/2m where D is the mean volume diameter, and m is the mean reciprocal value of grain diameters as seen on the micrograph. The shape and the size of the asymmetrical sample, see Fig. 1, were chosen so as to facilitate the study of fatigue crack initiation. Before coating, the samples were machined and the surface of the groove was polished. After chromizing, the two side surfaces were machined and polished so as to retain the coating only on the top surface of the groove. Residual stress in the chromized layer was measured on the top of the grooved section of the specimen after the diffusion process. The inclined incident X-ray beam procedure was used.9&apos;10 The computations were performed with the initial assumption of a zero surface normal stress component. An electron microprobe analyzer was used to obtain the relative amounts of chromium and iron con-

Jan 1, 1970