Search Documents

By J. S. Levine, M. Prats

A homogeneous and uniform cylindrical reservoir containing oil and gas is fractured vertically on completion and is produced at a constant bottom-hole pressure. The fracture has an infinite flow capacity, is of limited lateral extent and is bounded above and below by the impermeable strata defining the vertical extent of the reservoir. Results show that such a fractured reservoir can be represented by a reservoir of circular symmetry having very nearly the same production history. The well radius of this circular reservoir is about 1/4 the fracture length and is essentially the same as that obtained previously for a single fluid of constant compressibility. At the same value of cumulative oil production, gas-oil ratios of fractured reservoirs producing at constant terzinal pressure are larger than those of reservoirs having no fractures. This leads to more inefficient use of the reservoir energy in fractured wells and results in lower reservoir pressures for the same cumulative oil production. The reduction in operating life due to fracturing a reservoir is not as great as that for a slightly compressible fluid. This diflerence can be accounted for by the lower reservoir pressure in the fractured reservoir and its adverse effect on the average mobility and compressibility of the oil. As anticipated, the reduction in operating life increases czs the reservoir permeability decreases. The type of results presented in this report can be used to determine the economic attractiveness of fracture treatments per se, to setect the initial spacing to be used in developing a field, and to compare the relative merits of fracturing available wells and infill drilling. INTRODUCTION The effect of vertical fractures on a reservoir producing either an incompressible or a compressible liquid has already been discussed in the 1iterature.l,2 Those results indicate that the production history of such a reservoir is essentially the same as that of a circular reservoir having an effective well radius of approximately one-fourth the fracture length. The present work reports on the effect of a vertical fracture on a reservoir producing two compressible fluids —oil and gas—by solution gas drive. Because of the empirical nature of the PVT and relative permeability data used to obtain the performance of such reservoirs, results can only be obtained numerically and with the aid of high-speed computers. Since reservoirs lose their radial symmetry when fractured vertically, pressure and saturation can no longer be given only in terms of distance from the well. Two coordinates (such as x and y) must now be used to describe the pressure and saturation within the reservoir, and, since we are dealing with compressible fluids, time is also a variable. Thus the solution of a vertically fractured reservoir requires finding two unknowns (pressure and saturation) in two space variables (say x and y) and in time (t). Since no means are readily and generally available for solving such problems at the present time, we have used the results of previous work1,2 to approximate the effect of a vertical fracture on a reservoir producing both oil and gas by depletion. The purpose of the present wmk, then, is to investigate the possibility of using available numerical techniques (limited at the moment to one space variable) to study the two-space-variable flow behavior resulting from a vertical fracture. Results obtained in the course of this investigation are also reported and discussed. Input and output data of the numerical methods used are given in practical units: BOPD, feet, psi, cp, and md. Results are discussed fist in terms of specific reservoir and crude properties and geometries. Later, dimensionless parameters are introduced in order to extend results to different values of some of the reservoir and fracture properties. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The reservoir has an impermeable, circular, cylindrical outer boundary of radius r,. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below' the producing layer (reservoir). It is assumed that there is no pressure drop in the fracture due to fluid flow. 1 indicates the general three-dimensional geometry of the fractured reservoir. Gravity effects and the effects of differential depletion resulting from variations in hydrostatic head (pressure) will be neglected. Thus, the flow behavior in the fractured reservoir is described by the

By E. F. Thomas

THE late George S. Rice was active in the inves--I- tigation of bumps, particularly in the last ten years of his career as chief mining engineer of the U. S. Bureau of Mines. Since most of his investigation was carried out in Great Britain, continental Europe, and—to a lesser extent—Canada, his thinking on prevention was influenced considerably by the experience of those countries. It is not surprising, therefore, that when he was called upon a few years before his retirement to investigate bumps in the U. S. and suggest ways to prevent them, he turned to longwall mining. A longwall method had been most successful in combating the bump hazard in mining coal under deep cover, especially in Great Britain, but the prevailing method there at the time was advancing longwall mining, which he knew was uneconomical under U. S. mining conditions. For this reason he proposed a modified retreating longwall system that he believed included the best features of the advancing method. As brought out by Rice,' if the cover is 2000 ft and 50 pct of the coal is extracted, the static load on the remaining pillars will be about 4000 psi, which exceeds the ultimate crushing strength in most instances. If the pillar coal is overloaded before a pillar line is established, then the abutment zone preceding a line of extraction is no place to split pillars or extract them by any method other than an open-end system. Rice therefore advocated open-end mining, preferably by longwall, but he was willing to compromise with long-face mining if the longwall method was not acceptable. Rice's system was put into operation in a mine in Harlan County, Kentucky,3 but subsequent experience has shown that it did not take into account two important factors—avoidance of pillar-line points and maintenance of adequate development in advance of the pillar-line abutment area. For ten years after Rice's retirement the USBM did little investigation and research on bumps, chiefly because so few were occurring that there was not much cause for alarm. But in 1951 there were three occurrences involving fatal injuries, and the Bureau began a statistical survey in that year. C. T. Holland, head of the department of mines at Virginia Polytechnic Institute, was retained as a consultant. The resulting study' of 117 case histories brought out these important conclusions: 1) Almost invariably the bump occurred in a locality affected by the abutment zones of one or more pillar lines. 2) In most cases the locality of the bump was influenced by the abutment zones of more than one pillar line. The term pillar-line point has been used for many years in the Appalachian region for such a situation. Point is used in the geographical rather than the mathematical sense. 3) In pillar-line extraction the following practices are safest in preventing bumps: a. The mine layout should provide for pillars of uniform size and shape along the extraction line. b. The mine layout should be planned so that no development need be done in the abutment zone of a pillar line. c. The layout should permit open-end extraction of pillar lines from the next goaf, so that it will not be necessary to resort to pocket mining, splitting pillars, or any practice that will involve driving in the direction of the goaf within the abutment zone. d. Pillars should be large enough to support area without undue roof and floor convergence before establishment of a pillar line. These are, of course, generalities, and while they are useful in laying out areas where bumps can be expected, they are of limited help in many mines that were committed to a system of mining before it was realized that they were subject to bumps. Under such conditions it becomes necessary to choose between the following alternatives: 1) Abandon the territory, except for pillars that offer no extraction problems. 2) Through experience select the pillars that are most heavily loaded, and, by augering, induce bumps from a safe vantage point so that impinged loads are relieved. This method was first developed at the Gary, W. Va. mines of U. S. Steel Corp. and later adapted to mining thick coal beds at Kaiser Steel's Sunnyside mine in Utah. No scientific method is available to determine where to drill within a loaded pillar. Although this method of unloading has worked very successfully at Gary—with one exception—

Jan 1, 1959

By C. E. Richards, C. N. Reid

The influence of preferred orientation on the incidence of defbrtnation tuinning has been studied. High-purity iron with almost vandonz grain orientation was cotnpared uitll iron of the sa)ne grain size and composilion lza,ing a strong (110) fiber texture. As expected from published work on single crgslfls, /he ))lean stress for the onset of luitzning-, and the l,olu)nt. fraclion of twinned nzaterial obserlled in lension differed fron the 1-a1ue.s it2 co?nPression for tnolerial with a slrong texlure. The llinning stress of "rctndorrl " )zalerial did not 17ary with the sense of the aPPlied unin.via1 stress, but sirprisinglg the incidence of 1c)i)zning- was about three 1i))zes greater ill conzp?&apos;ession Illon in lension. These results (Ire attributed entirely to ovienbation and may be nderslood in ler?ns of the shear slresses acting on the allowed twinning syster)is. J. HE twins most commonly formed in bcc metals may be described as regions of the crystal in which a particular set of (112) planes is homogeneously sheared by 0.707 in the appropriate ( 111) direction. A similar twin-related crystal could be produced by a shear of 1.414 in the reverse (111) direction but twinning by this large displacement has never been reported. Thus, twinning is unidirectional and a shear stress which produces twinning does not do so when its sense is reversed. The sense of a shear Stress is reversed when the loading is changed from tension to compression, or vice versa. Consequently, for a given orientation of a crystal relative to a uniaxial stress, only a fraction of the twelve (112) twinning systems are geometrically capable of operating in tension, and the remaining systems may operate only in compression. Therefore, when twinning is involved, there are expected to be differences in behavior between crystals tested in uniaxial tension and those tested in compression. This has been verified experimentally by Reid et 01.&apos; and Sherwood el al.,&apos; although a critical stress criterion was not encountered. Furthermore, twinning stresses in colmbium," tungten, tantalum,&apos; irn,&apos; i-Fe,\ nd molybdenum7 single crystals have been shown to depend critically on orientation, although again twinning did not occur at a critical value of the macroscopic shear stress. However, when twinning occurs, it generally does so on the most highly stressed systems, 1--4&apos;6&apos;8&apos;9 implying that the stress level does have some relevance to twin formation. In view of the large orientation dependence of twinning in bee single crystals, it might be expected that such an effect would be present in poly crystalline material which possesses a recrystallisation texture. Indeed, riestner" showed that the twinning stress in tension is very orientation-sensitive it1 <&apos;grain-oriented, silicon-iron;" this material possessed a very strong t c m^ii a nnr x_____k . i-_ii__ ri_______j. _x r»i_._:__i preferred orientation obtained by secondary recrystallisation. Reid et a/.&apos; observed a marked difference in the tensile and compressive yield stresses of polycrys-talline columbium which was rationalised in terms of the effect of a preferred orientation on twinning. No other such illformation is known to the authors. Several investigations of twinning in polycrystalline bcc metals have been reported in which the possible existence of a preferred orientation was not even mentioned. It is the purpose of this paper to show that there is a strong effect of texture on twinning in polycrystalline iron, and to poilt out the difficulty in eliminating preferred orientation in recrystallised metals. 1. EXPERIMENTAL METHOD Material and Specimen Preparation. Low-carbon, high-purity iron was obtained from the National Physical Laboratory in the form of $-in. diam rod which had been cold-swaged from a diam of 1 in. The composition of the material is given in Table I. The as-received bar was cold-swaged directly to 0.185 in. diam from which cylindrical tensile and compression specimens were machined. Specimen geometry is illustrated in Fig. 1. The gage length was 0.30 in. long and 0.10 in. diam; it should be noted that, apart from the extra heads which are necessary for tensile loading, the geometry and dimensions of the two types of specimen are identical. The specimens were heat treated either by sequence A or B outlined in Table 11. The essential difference between these two treatments is that in one case the material was repeatedly cycled through the y- to a-phase change in order to produce grains of almost random orientation ("random" iron)

Jan 1, 1969

By H. R. Cooper, T. S. Mika

T. S. Mika (Department of Mineral Technology, University of California, Berkeley, Calif.) - Dr. Cooper&apos;s attempt to establish a correlation between process behavior and operational variables on the basis of a statistical analysis after imposing a reasonable process model is a very commendable improvement on the use of standard regression techniques. However, it must be recognized that the imposition of a model has the potential of yielding a poorer representation if its basic assumptions or mathematical formulation are invalid. It appears that at least two aspects of his treatment require some comment. First, the limitations on the kinetic law where xta represents a hypothetical terminal floatable solids concentration (cf. Bushell1), should be mentioned. Most current investigations2-9 appear to utilize the concept of a distribution of rate constants rather than a single unique value, k, to describe flotation kinetics. A distributed rate constant is certainly a more physically meaningful concept than that of a terminal concentration. The study of Jowett and safvi10 strongly indicates that xta is merely an empirical parameter, whose actual behavior does not correspond to that expected from a true terminal concentration. Rather than being a strictly mineralogical variable, as Dr. Cooper&apos;s treatment implies, it apparently represents the hydromechanical nature of the test cell as well as the flotation chemistry. The extension of batch cell kinetic results to full-scale continuous cell operation is a suspect procedure if the effect of such nonmineralogical influences on x,, remain unevaluated. There is evidence that introduction of a terminal concentration is necessitated by the inherent errors which arise in batch testing and are eliminated by continuous testing methods.&apos; Possible lack of validity of the author&apos;s use of Eq. 1 is indicated by two unexpected results of the statistical analysis of his batch data. The first is the apparent corroboration of the assumption that the rate constant, k, is independent of particle size, i.e., of changes in the size distribution of floatable material. This assumption directly contradicts numerous results 2,4,11-l8 for cases where first order kinetics prevailed and ignores the phenomenological basis for the analysis of flotation in terms of a distribution of k&apos;s. It must be recognized that, if the rate constant is size dependent, the lumped over-all k would be time dependent; Eq. 1 would then no longer be valid. Cooper&apos;s x,, is determined by batch flotation of a distribution of sizes for an arbitrary period of time. If the size dependence of k is artificially suppressed, x,, will become a function of the experimental flotation time used in its determination. Upon reviewing the rather extensive literature concerning batch flotation kinetics, there appear to be few instances where constant k and x,, adequately adsorb variations in floatability due to particle size. The second surprising result is the low values of the distribution modulus, n, determined. Contrary to Cooper&apos;s assertion, most batch grinding (ball or rod mill) products yield values of n > 0.6, which increase as the material becomes harder.&apos;&apos; It is likely that the values of n = 0.25 and n = 0.42 for Trials 1 and 2, respectively, are completely unreasonable, and even the value n = 0.54 obtained for Trial 3 is unexpectedly low. Possibly, this indicates inherent flaws in the three trial models considered, in particular the assumed particle size independence of the rate constant, k. The above does not necessitate that Eq. 1 (and the terminal concentration concept) is invalid; it could constitute a good first approximation. However, the qualitative arguments used by Dr. Cooper in its justification are somewhat frail and require verification, particularly since much of the flotation kinetics literature is in opposition. Apparently, no effort was made to test these hypotheses on the actual data; in fact, since they pertain to a single batch test time, his data cannot be utilized to evaluate the kinetics of flotation. To evolve a control algorithm on the basis of this infirm foundation seems a questionable procedure. Another difficulty in his analysis arises in consideration of the froth concentrating process. As Bushel1 &apos; notes, for Eq. 1 to be valid it is necessary that the rate of recycle from the froth be directly proportional (independent of particle size) to the rate of flotation transport from the pulp to the froth, a restrictive condition." Harris suggests that it is more realistic to assume that depletion occurs in proportion to the amount of floatable material in the pertinent froth phase volume (treating that volume as perfectly mixed).12,21,22 The physical implications of

Jan 1, 1968

By Larry Kaufman, Martin Schatz

The effect of pressure on the extension of the ? loop in the FeSi system has been determined by means of metallogvaphic studies and hardness measurements performed on a series of high-purity Fe-Si alloys containing 7.5, 11.0, and 13.9 at. pct Si, respectively. These mensurements, performed at 42 kbar and temperatures up to 1200oC, indicate that the ? loop is expanded to about 10 at. pct Si at 42 kbar as opposed to a maximum extension of 4 at. pct Si at 1 atm. Comparison of the experimental results with thermodynamic predictions of the pressure shifts yields satisfnctory results. DURING the past few years, several studies have been performed in our laboratory1-&apos; in order to determine the effect of high pressure on phase equilibrium in pure iron and iron-base alloys. The purpose of these studies has been to elucidate the effects of high pressure experimentally and to compare the observed results with predicted pressure effects derived on the basis of known thermody-namic and volumetric data at 1 atm. These studies have included work on pure iron2,5,7 as well as Fe-Ni,1,5 Fe-cr,l,5 and Fe-c4-6 alloys. In addition, Tanner and Kulin3 have reported results of pressure studies on two Fe-Si alloys containing 2.0 and 6.25 at. pct Si. At the time of this latter study, no detailed information was available concerning the difference in volume between the a (bcc) and ? (fcc) phases in the Fe-Si system as a function of silicon content. In order to compare their observations with calculated pressure shifts, Tanner and Kulin were forced to assume that silicon had no effect on the difference in volume between a and ? iron. The resulting discrepancy between their calculation of the a/? phase boundary at 42 kbar and the observed results led them to the conclusion that silicon additions probably decrease the difference in volume between a and ? iron. Recently: Cockett and Davis8,9 have reported de- tailed studies of the lattice parameters of a series of Fe-Si alloys at temperatures ranging from 20" to 1150°C. These measurements, performed on alloys in the bcc and fcc range, show that silicon does indeed decrease the difference in volume between a and ? iron. By correcting the calculations of Tanner and Kulin in line with the observed effect of silicon they were able to show improved agreement between computed and observed pressure shifts.&apos; The present measurements were undertaken to provide additional corroboration of this effect, by extending the range of composition, in addition to exploring a situation where large extensions of a ? loop could result in impingement of the ? field with an ordered bcc phase (based on Feo.75Sio.25). I) EXPERIMENTAL PROCEDURES AND RESULTS The alloys investigated were obtained from Dr. F. Kayser of M.I.T. They were prepared at the Ford Scientific Laboratory by vacuum melting electrolytic iron and high-purity silicon. The melts were poured under an argon atmosphere into hot-topped steel molds. Subsequently the ingots were hot-worked down to 1/2-in.-diam rods. Three alloys containing 7.5, 11.0, and 13.9 pct Si were studied. Carbon, regarded as the principal impurity, analyzed at, or below, 0.001 wt pct for all of the alloys. Prior to pressure-temperature treatment, the rod was annealed for 24 hr in vacuum at 1000°C, water-quenched, and subsequently machined into 0.100-in.-diam by 0.100-in.-long specimens. Subsequent to machining, the specimens were again annealed and then examined metallographically. They were found to exhibit a clear coarse-grained ferrite similar to Figs. 10 and 110 of Ref. 1 and Fig. 2 of Ref. 3. Subsequently, specimens of each alloy were equilibrated at 42 kbar at various temperatures in supported piston apparatus.1,3,4,6 Three specimens, one of each alloy, were wrapped in platinum and exposed simultaneously. The pressure-temperature cycle consisted of increasing the pressure from ambient to 42 kbar at 25oC, heating rapidly to the desired temperature, holding for 15 min, and quenching to 100°C, followed by slower cooling to 25°C and pressure release. The temperature was measured with a Pt/Pt-13 pct Rh thermocouple which was not corrected for pressure effects. Subsequently, specimens were examined metallographically and by

Jan 1, 1964

By R. J. Reynik

A semiempirial small flunctation theory of diff- sion in liquids is presented, which employs a fluctuation energy assumed quadratic for a small atomic or molecular displacement and Einstein&apos;s random-iralh model. The resulting diffusion equation is given by In these equations. D is the diffusivity, is the average liquid shite coordination number (at interatomic distance d. cm. T is the absolute temperature, xu. em, is (the diffusive displacement. K, is the quadratic fluctuation energy force constant, and rg, cm, are the radii oj diffusing atoms A and B, respectively. The quantities Xn and K are calculated from the computer-filled values of the slope and intercept. respectively. The radius of self-diffusing atom or radii and of diffusing atoms A and B are eta United and compared with values reported in the literature.. The predicted linear variation of diffusivity with. It tempera lure htm been observed in approximately thirty-iire metallic liquid systems, and in over seventy-fiee other liquid systems, including the organic .alcohols, liquified inert gases, and the molten salts, ALTHOUGH the average density within a macroscopic volume element of liquid is constant for fixed total number of atoms. pressure. and temperature, there exist microscopic: density fluctuations within the respective volume element. As such the microscopic volume available to an atom and its Z first nearest neighbors at any instant of time fluctuates above and below the average volume available to these atoms. If one assumes that liquid state atoms vibrate as in a solid. and further postulates that the mean position of any atom in the liquid state is not stationary. but shifts during every .vibration a distance 0 5 j 5 xo. then every atom in the liquid state continuously undergoes diffusive displacements which vary in the range 0 5 j 5 ro. Mathematically. for a binary liquid system consisting of atcrms A and B. the maximum diffusive displacement. .YO, is defined by the equation: where d is the average liquid state interatomic distance at specified liquid state coordination number Z. and v~ \ and vg are the effective radii of diffusing atoms A and B: respectively. For self-diffusion. r^ equals rg , and Eq. [I.] reduces to: It is interesting to note that Eq. [l] or [2] can be used to compute the radii of the diffusing atoms, provided one had an experimental evaluation of xo. As such. the computed radii could be compared with metallic or crystallographic ionic radii to ascerlain the electronic character of the diffusing atoms. Thus it is proposed that in the liquid state the n~otion of an atom relative to its original equilibrium position of oscillation represents the thermal vibration of any atom and its Z first nearest neighbors. while the small and variable displacements. 0 5 1 5 xc,. of the centers of oscillation represent the complex diffusive motions of the atoms at constant temperature and pressure. This is consistent with data obtained from slow neutron scattering by liquids1 &apos; and resembles an itinerant oscillator model of the liquid state.&apos;" It is further postulated that the atomic displacements characterizing the liquid state diffusion process are essentially a random-walk process. As such. it nlay be described by Einstein&apos;s equation:&apos; where D is the diffusivity. sq cm sec-&apos;. j2 is the mean square value of the diffusive displacement. and i> is the frequency of density fluctuations giving rise to diffusion. FORMULATION OF DIFFUSION EQUATION The effective spherical volume occupied by an atom, as a consequence of a microscopic density fluctuation which enlarges the volume available to any atom, exceeds its average liquid state atomic volume by an amount: where AV is the enlarged spherical volume, v is the radius of the diffusing atom. and j is the elementary displacement distance from the original center of oscillation of the vibrating atom to a new center of oscillation position. For small atomic displacements. where c is a constant whose value depends upon the assumed geometry of the enlarged volume. For a spherical increase in volume, c equals 4nr2. Following the treatment of Furthl&apos; and ~walin." assuming the enlarged volu~nes AL7 for the diffusing atoms are distributed in a continuunl. the probability of finding a fluctuation in the size range 0 5 j 5 xo defined by Where c includes the geometric constant cl and Eij) is the fluctuation energy causing the volume change. But the proposed model assumes all the Z first nearest-neighbor atoms are centers of oscillation. and hence the probability that any of these atoms is adjacent to a fluctuation of magnitude 05j5xo is unity. Thus:

Jan 1, 1970

By J. H. Moran, E. E. Finklea

The pressure build-up technique is a recognized method of determining permeability from conventional drillstem tests. In this paper an effort is made to extend such techniques to the interpretation of data obtained from the wireline formation tester. Such a study is necessary because of the differences, for this case, in the magnitude of the flow parameters (rate of flow, amount of recovered fluids) and in the flow geometry (flow through a perforation vs flow across the face of the wellbore, etc.) involved in the solution of the equations of flow for compressible fluids. The perforation is replaced by a spherical hole, and the effect of the borehole is neglected, so that the flow can be considered to be radial in a spherical co-ordinate system. Arguments are presented to justify this idealization. Assuming single-phase flow, general relations between pressure and flow rate are developed for a homogeneous medium. The study is then extended to permeable beds of finite thickness. It is shown that the early stages of pressure build-up tend towards spherical flow, while the later stages tend towards cylindrical flow. The thinner the bed, the more quickly flow approaches the cylindrical model. The prevalence of thin beds in practical work makes this analysis quite important. Cases involving permeability anisotropy are treated. INTRODUCTION From wireline formation tester operation, two types of data are obtained: (1) the nature and amount of recovered fluids, and (2) the pressure history recorded during the test. A number of papers have been written dealing with the interpretation of formation production on the basis of the recovered fluids.&apos;.&apos; In general, the methods described have been quite accurate for both high- and low-permeability formations. The present paper will deal with an analysis of the pressures observed. An analysis of the pressure build-up curves obtained in hard-rock country has already been attempted on the basis of the formula proposed by Hor-ner. Although this approach has met with success in many instances, some questions have been raised as to its validity. It is the aim of the present study to place the analysis of pressure build-up in the formation tester on a firmer basis, from which more detailed methods of interpretation can evolve. Because of the great differences between the operation of the wireline formation tester and the conventional drillstem test, modifications are necessary in the interpretation. The major difference relates to the flow geometry. Once the flow geometry has been established other features such as multiphase flow, skin effect, afterflow, etc., well described in the literature, can be introduced. It will be assumed that the mechanical operation of the formation tester is already known to the reader.6 t will suffice here merely to state that the tester provides the means for taking a relatively small sample of the fluid immediately adjacent to the borehole, and for recording the subsequent pressure response. In comparison with conventional drillstem tests, the time required for a satisfactory pressure build-up response is much shorter, because of the relatively small quantity of fluid withdrawn by the wireline tester. This feature is highly desirable in the case of low-permeability formations. For an analysis of the pressure response within the formation, three simple flow geometries are considered— linear, cylindrical and spherical. The spherical and cylindrical flow geometries are most pertinent to the formation tester; therefore, they will receive the major emphasis. Since the configuration of the borehole and the perforation made by the tester complicate the flow geometry, it is necessary to allow for them in the drawdown response. However, because of the volume of formations contributing to the pressure-response, the details of the perforation shape are unimportant in the build-up period. Since relatively small amounts of fluid are withdrawn from the formation, in contrast to a conventional drill-stem test, a study of the "depth of investigation" and the significance of drawdown as well as build-up data will be included. Because the "depth of investigation" will be shown to be rather large, the effect on the build-up curves of the finite thickness of the permeable bed is considered. It is this consideration that leads to the importance of cylindrical flow geometry. Also included is a discussion of permeability anisotropy and its effect on the interpretation of the tester results. The pressure curves recorded by the formation tester will follow two general patterns, depending upon whether the formation is of high or low permeability. Fig. I (a and b) schematically illustrates these two responses. In Fig. 1(a), the high pressure recorded during fill-up of the tool is essentially the pressure differential across the choke in the system. In Fig. l(b), the flow rate is

By W. R. Hibbard

THE classical work on the electrical conductivity of alloys was carried out by Matthiessen and his coworkers1 in the early 1860&apos;s. He attempted to correlate the electrical conductivity of alloys with their constitution diagrams, but the information regarding the latter was too meager for success. Guertler2 reworked Matthiessen&apos;s and other conductivity data in 1906 on the basis of volume composition (an application of Le Chatelier&apos;s principle with implications as to temperature and pressure effects), and obtained the following relationships between specific conductivity and phase diagrams (plotted as volume compositions) : 1—For two-phase regions, electrical conductivity can be considered as a linear function of volume composition, following the law of mixtures. 2—For solid solutions, except intermetallic compounds, the electrical conductivity is lowered by solute additions first very extensively and later more gradually, such that a minimum occurs in systems with complete solid solubility. This minimum forms from a catenary type of curve. Intermetallic compound formation with variable compound composition results in a maximum conductivity at the stoi-chiometric composition. Landauer" has recently considered the resistivity of binary metallic two-phase mixtures on the basis of randomly distributed spherical-shaped regions of two phases having different conductivities. His derivation predicts deviations from the law of mixtures which fit measurements on alloys of 6 systems out of 13 considered. Volency (Ionic Charge) Perhaps the first comprehensive discussion of the electrical resistivity of dilute solid-solution alloys was presented by Norbury&apos; in 1921. He collected sufficient data to show that the change in resistance caused by 1 atomic pct binary solute additions is periodic* in character. The difference between the period and/or the group of the solvent and solute elements could be correlated with the increase in resistance. Linde5-7 determined the electrical resistivity (p) of solid solutions containing up to about 4 atomic pct of various solutes in copper, silver, and gold at several temperatures. He reported that the extrapolated"" increase in resistance per atomic percent addition is a function of the square of the difference in group number of the solute and solvent as follows: ?p= a + K(N-Ng)2 where a and K are empirical constants and N and Ng are group numbers of the constituents. This empirical relation was subsequently rationalized theoretically by Mott,8 who showed that the scattering of conduction electrons is proportional to the square of the scattering charge at lattice sites. Thus, the change in resistance of dilute alloys is propor-t,ional to the square of the difference between the ionic charge (or valence) of the solvent and solute when other factors are neglected. Mott&apos;s difficulty in evaluating the volume of the lattice near each atom site where the valency electrons tend to segre-gate: limited his calculations to proportionality relations. Recently, Robinson and Dorn" reconfirmed this relationship for dilute aluminum solid-solution alloys at 20°C, using an effective charge of 2.5 for aluminum. In terms of valence, Linde&apos;s equation becomes ?P= {K2 + K1 (Z8 -Za)2} A where K1 and K2 are coefficients, A is atomic percent solute, Z, is valence of solvent, and Zß, is valence of solute. Plots of these data for copper, silver, gold, and aluminum alloys are shown in Fig. 1. The values of K1 and K2 are constant for a given chemical period (P), but vary from period to period. The value of K, increases irregularly with increasing difference between the period of the solvent and solute element (AP), being zero when AP is zero. The value of K, appears to have no obvious periodic relationship. All factors other than valence that affect resistivity are gathered in these coefficients. Because of the nature of the coefficients, Eq. 1 is of limited use in estimating the effects of solute additions on resistivity unless a large amount of experimental data are already available on the systems involved. It is the purpose of the first part of this report to investigate the factors that may be included in the coefficients of Linde&apos;s equation. On this basis, it is hoped that the relative effects of solute additions on resistivity can be better estimated from basic data, leading to a more convenient alloy design procedure. It is well 10,11 that phenomena that decrease the perfection of the periodic field in an atomic lattice, such as the introduction of a solute atom or strain due to deformation, will also increase the electrical resistivity. Thus, in an effort to relate changes in electrical resistivity to alloy composition, it appears appropriate to consider the atomic characteristics related to solution and strain hardening

Jan 1, 1955

By Thomas R. Mager

An in.t!estigation has been made of the equilibrium conditions at 1200°C in the reaction between hydrogen sulfide gas and sulfur dissolved in Fe-Si alloys From this the equilibrium constant, activity coefficient, and activity of sulfur in solution were calculated. A number of studies of the equilibrium of sulfur with iron and iron alloys have given closely agreeing results from which the activity and free energy of the dissolved sulfur may be found. Sherman, El-vander, and chipman1 discussed the significant researches of dilute solutions of sulfur in liquid iron prior to 1950, and the results of this study indicated that the relationship between the ratio of PH2S/PH2 in the environment and the percentage of sulfur in solution is not a linear one. Morris and williams2 studied the equilibrium conditions in the reaction between hydrogen sulfide gas and sulfur dissolved in liquid iron and Fe-Si alloys, and reported that silicon dissolved in iron has a pronounced effect on the equilibrium conditions. They found that the activity of sulfur in iron is increased by the addition of silicon. At a silicon content of 4 pet the activity coefficient of sulfur was about twice that for sulfur dissolved in pure iron. Sherman and chipman3 investigated the chemical behavior of sulfur in liquid iron at 1600°C through the study of the equilibrium: H2 + S = H2S; K = PH2S/PH2 . 1/as [1] From the known equilibrium constant of the reaction between H2, H2S, and S and the experimental data, the activity of sulfur in the melt was determined. They found that the activity coefficient of sulfur defined as fs = as/%s is increased by silicon and decreased by manganese. Morris4 and Turkdogan5 also reported that manganese decreases the activity coefficient of sulfur in liquid iron and iron-base alloys. A recent technique of sulfur analysis developed by Kriege and wolfe6 of the Westinghouse Research Laboratories permits an accurate sulfur analysis of 0.5 * 0.2 ppm in the range of 0.1 to 3 ppm, whereas in the range of 3 to 50 ppm the accuracy is ±1 ppm. This technique of sulfur analysis was utilized in this experiment. Previous unpublished data reported that sulfur analysis by the combustion technique was not accurate below 20 ppm. EXPERIMENTAL PROCEDURE Five 5-lb ingots of high-purity Fe-Si were prepared. Three of these ingots were prepared without the addition of manganese but with a variation of silicon contents from 2 to 4 pet. The remaining two ingots contained 3 pet Si with the addition of manganese. Ingots were made at each of three silicon levels: 2, 3, and 4 pet. No alloys were made with less than 2 pet Si since below approximately 1.8 pet Si the binary alloy exhibits a to ? transformation. The two additional ingots of 3 pet Si-Fe were made at each of two manganese levels: 0.20 and 0.50 pet. To minimize the effects, if any, of impurities on the activity of sulfur on Si-Fe, the best metals available were used for melting. All ingots were vacuum-melted in magnesium oxide crucibles. After obtaining samples for chemical analyses, the ingots were processed. This consisted of hot rolling and subsequently cold rolling the alloys. Each ingot was hot-rolled at 1000°C, reheating between every pass to minimize grain growth. All heating was done in a protective argon atmosphere. The slabs were hot-rolled to strips 50 mils thick. After hot rolling, all the material was pickled to remove the scale formed on the surface of the strip during hot rolling. The material was then cold-rolled to 12-mil strips. Single strips of the material used in this experiment were hydrogen-annealed at 1200°C for 16 hr in an alumina tube. Chemical analyses of strips M-1, M-3, M-4, M-7, and M-8 are given in Table I. Sulfur, silicon, and manganese analyses were made from the millings from the cold-rolled 12-mil strips. The oxygen analyses were made from slugs of the as-cast material. The hydrogen sulfide used in these experiments was supplied from cylinders containing a mixture of argon and 1 pet hydrogen sulfide. The parts per million of hydrogen sulfide were determined from the analysis of the exit gas of the annealing furnace during each anneal. The flow rate of hydrogen was approximately 1 liter per min in all anneals. The

Jan 1, 1964

By Rex O. McLellan

The formation of a monovacancy in a metal is simulated in an elastic model by the displacement of the surface of a small spherical cavity in a large elastic continuum. The application of linear elasticity to this distortion results in a well- known formula for the energy and an expression for the concomitant entropy change due both to the shear strain in the continuum and also to the dilation of the solid resulting from the boundary conditions at the surface of the solid. Elastic data (the sliear modulus and its temperature coelficient) are used to calculate the entropy and energy of formation for many metals. Despite the simplicity of the assumptions involved, the agreement between the calculated entropies and energies and experimental values is remarkably good. In recent years there has been a large increase in measurements of the absolute concentration of mono-vacancies in metals as a function of temperature. Hence new data for both the energy and the noncon-figurational entropy of formation of monovacancies has become available. Recent measurements&apos; of the anomalous (non-Arrhenius) self-diffusion in many bcc metals has also focused interest on the prediction of the thermodynamic parameters of mono- and multi-vacancies in those metals for which no data are available. Damask and Dienes&apos; have discussed the various theoretical calculations of the energy of formation EL, of a monovacancy. These include simple models involving the breaking of atomic bonds on moving atoms from the interior of a crystal to the surface, models combining elastic calculations with surface-energy terms and detailed quantum mechanical calculations. The simler models give the correct order of magnitude of &, but tend to overestimate it by a factor of about two. The quantum mechanical calculations4"7 have been carried out for the noble and alkali metals with generally reasonably good agreement with the available Ef data. The calculation of entropy of formation Sfv14 lnvolves a fundamental calculation of the perturbation of the phonon spectrum caused by the creation of a vacancy. Huntington, Shirn. and wajda8 have given an approximate evaluation of sJV by considering an Einstein model for the localized vibrations in the immediate neighborhood of the defect and then using elastic theory to calculate the entropy associated with the shear stress field in the distorted crystal (as originally proposed by Zenerg). They also included a term due to the dilation of the crystal. They obtained a value of 1.47k for copper, in good agreement with the experimental value (1.50k). However, Nardelli and Tetta- manzi1° have recently shown that neglecting the coupling between atoms (Einstein Model) may lead to a serious error so the agreement may be somewhat fortuitous. In this work simple linear elastic theory is used to calculate the entropy and energy of formation of mono-vacancies. Despite the simplicity of some of the assumptions involved, the agreement with the available experimental data is remarkable. However. the reasonable degree of success in the application of linear elastic calculations to the excess entropy of a solute atom in a dilute solid solution1&apos; indicates that the application of elastic theory to vacancies. where the interaction of different atomic species is not involved, may not be inappropriate. THE ELASTIC MODEL The metal is assumed to be a spherical elastic continuum. A small spherical cavity of volume V = 4i;v:&apos;/3 is cut from the center. removed. and dissolved rever-sibly in the bulk of the material. TO a good approximation no net atomic bonds are broken and the material does not undergo a volume change although the externally measured volume of the body would increase by V. The radius of the sphere of metal is much larger than r Next a negative pressure is applied to the cavity causing its surface to be displaced inward by an amount simulating the relaxation of the lattice around a monovacancy. In this model the energy and entropy accompanying the distortion are taken as 4, and <. As a first approximation the equation of state for the solid is taken as: r = ro(i + *~D LiJ where K is the bulk modulus. P the hydrostatic pressure. Vo the volume of the material at 0°K and zero pressure. and d+/dT = 30. where 0 is the linear thermal expansion coefficient. The variation of entropy with hydrostatic pressure is given by the Maxwell equation: These equations give the entropy change resulting from increasing the hydrostatic pressure from 0 to P as: and since • we have: This is the entropy arising from the dilation resulting

Jan 1, 1970

By H. C. Gatos, A. D. Kurtz

WITH the growing interest in the mechanism of self-diffusion of metals, the study of accurate and convenient methods for determining self-diffu-sion coefficients appears highly desirable. It was with this objective in mind that the present investigation was undertaken. Gatos and Azzam1 employed an autoradiographic technique for measuring self-diffusion coefficients of gold. This method involved sectioning of the specimen through the diffusion zone and recording the radioactivity directly on a photographic film. Because of the very short range of the emitted ß rays in gold, the activity recorded on the film was essentially the true surface activity. With proper choice of the sectioning angle, sufficient resolution could be obtained and the entire concentration-distance curve recorded in one measurement. For the boundary conditions of the experiment, where an infinitesimally thin layer of radioactive material diffuses in positive and negative directions into the end faces of a rod of infinite length, the solution of the diffusion equation is C/Cn = 1/v4pDt exp (-x2/4Dt) where C is the concentration of diffusing element (photographic density in this case), C,, is the constant (depending upon amount of radioactive material), x is the diffusion distance, D is the diffusion coefficient, and t is the time. Thus, by plotting the logarithm of the concentration vs the square of the diffusion distance, a straight line results and the slope contains the diffusion coefficient. In this manner, the self-diffusion coefficient of gold can be obtained as a function of temperature. In the present investigation the results reported by Gatos and Azzam1 have been verified, and the autoradiographic technique has been further developed and applied for the determination of the self-diffusion coefficient of gold at a number of temperatures. Furthermore, the energy of activation for the self-diffusion of gold has been conveniently determined. . Experimental Techniques Preparation of Specimens: The inert gold of high purity was received in the form of a rod from which cylinders were cut and machined to a diameter of 0.500 in. The specimens were annealed to a suitably large grain size and the faces were surface ground prior to the deposition of the radioactive layer. The radioactive isotope Au198 was chosen. It was produced in the Brookhaven pile by means of the reaction Au197 + n ? Au108. It decays by ß emission (0.96 mev) to Hg108 with the subsequent emission of a y ray (0.41 mev). 70Au 108 ? 80Hg 108 + -1e°. The half life of the Au108 is 2.7 days so that a strict time schedule had to be maintained in order to secure sufficient activity until the end of the experiments. For this reason, initial activities as high as 10,000 millicuries per gram were used. The gold arrived in the form of foil and was evaporated onto one face of each gold specimen cylinder to a thickness of about 100A. A sandwich-type specimen was formed by welding two such cylinders together. Evaporation of Gold: The gold was evaporated under vacuum from heated tantalum strips which were bent in such a way as to limit the solid angle through which the gold was allowed to vaporize, thus insuring a more efficient utilization of the gold. The specimens rested on flat brass rings which had an inner diameter of 0.475 in. The entire specimen-holding assembly could be manipulated from outside the vacuum system by means of a magnet which attracted a slug of soft iron attached to the assembly. By evaporating inert gold on glass slides under conditions identical to those employed for the radioactive gold, it was found that the thickness of the films was about 100A. Welding: The welding was performed by hot pressing in a stainless steel cylinder. The inside of the cylinder was threaded and fitted for two plugs. The specimens to be welded were placed in the middle of the cylinder and two pressing disks, one at each end, were inserted to avoid shearing stresses as the plugs were tightened. Mica disks were placed between the pressing disks and the specimens to prevent them from welding. The plugs were then tightened with a hand wrench and the entire unit was placed in an argon stream for about an hour to remove the oxygen. The unit was then inserted in the center of an argon atmosphere furnace maintained at about 700°C and left there for about an hour. Because of the difference in the temperature coefficient of expansion of the two metals, as the temperature rose. the pressure on the specimen-rollple increased and a weld resulted Welding was generally satisfactory under the conditions described.

Jan 1, 1955

By Vittal S. Bhandary, B. D. Cullity

Swaged iron wire has a cylindrical {001} <110> texture. The texture is also cylindrical after re-crystallization in the absence of a magnetic field, but <111> and <112> components are added to this texture when recrystallization occurs in a field. The mecizanical properties in tension and in torsion are not greatly altered by these changes in texture. AS shown in a previous paper,1 cold-worked wires of the two fcc metals copper and aluminum can be made relatively strong in torsion and weak in tension, or vice versa, by proper control of preferred orientation (texture). The deformation texture can be controlled by selection of the starting texture (texture before deformation), because certain initial orientations are stable during deformation. The present paper reports on similar work performed on bcc iron. In this case it was clear at the outset that there was no hope of controlling the deformation texture, which is one in which <110> directions are aligned parallel to the wire axis. (1t has usually been regarded as a fiber texture, but Leber2 has recently shown that it is a cylindrical texture of the type {001} <110>. In either case, <110> directions are parallel to the wire axis.) There is general agreement on this texture among a large number of investigators, which in itself suggests that the starting texture has no influence on the deformation texture. More direct evidence was produced by Barrett and Levenson,3 who reported that iron single crystals of widely varying initial orientations all had a single <110> texture when cold-worked into wire. Thus <110> is a truly stable end orientation for iron and probably for other bcc metals as well. Under these circumstances attention was directed to the possibility of controlling the recrystallization texture. This texture is normally <110> in iron,4 just like the deformation texture. However, it is conceivable that this texture could be modified by a proper choice of the time, the temperature, and what might loosely be called the "environment" of the recrystallization heat treatment. In the present work the environmental factor studied was a magnetic field. The effect of heating in a magnetic field ("magnetic annealing") on recrystallization texture has been investigated by Smoluchowski and Turner.5 They found that a magnetic field produced certain changes in the recrystallization texture of a cold-rolled Fe-Co alloy. The texture of this material is normally a mixture of three components, and the effect of the field was to increase the amount of one component at the expense of the other two. Smoluchowski and Turner concluded that the effect was due to magnetostriction. With the applied field parallel to the rolling direction, the observed effect was an increase in the amount of the texture component which had <110> parallel to the rolling direction. In the Fe-Co alloy they studied, the magnetostriction is low in the <110> direction and high in the <100> direction. Thus nuclei oriented with <110> parallel to the rolling direction will have less strain energy than those with <100> orientations and will therefore be more likely to grow. In a later paper on the same subject, Sawyer and Smoluchowski6 ascribed the effect to magneto-crystalline anisotropy and made no mention of magnetostriction. In the papers of Smoluchowski et al. the intensity of the magnetic field was not reported but it was presumably large, inasmuch as it was produced by an electromagnet. In the second paper6 it is specifically mentioned that the specimens were magnetically saturated. But if magnetostriction has a selective action on the genesis of stable nuclei during recrystallization, that selectivity must depend only on differences in magneto-strictive strains between different crystal orientations and not on the absolute values of those strains. Thus the saturated state does not necessarily produce the greatest selectivity, because the relative difference in magnetostrictive strains between different crystal directions may be larger for partially magnetized crystals than for fully saturated ones.7 In the present work the specimens were subjected to relatively weak fields (0 to 100 oe) produced by solenoids. MATERIALS AND METHODS Armco ingot iron rod (containing 0.02 pct C and 0.19 pct other impurities) was swaged from 0.25 in. in diam. to 0.05 in., a reduction in area of 96 pct. The mechanical properties in tension and torsion were measured as described previously.&apos; Textures were measured quantitatively with chromium or iron radiation and an X-ray diffractometer,8,1 and

Jan 1, 1962

By A. U. Seybolt

In part 111&apos; of this series it was shown that during the selective oxidation of a 3 1/4 pct Si-Fe alloy in damp hydrogen, only silica, (observed at room temperature) as low cristobalite or low tridy-mite or both, was formed as an oxidation product. In some in- „ stances where the film was fairly thin (probably well under 100A) there was some suggestion of an amorphous form of SiO2. The present investigation of oxidation rate showed that the selective oxidation of silicon-iron can be rather complicated, and apparently impossible to rationalize in an unequivocal manner. In some temperature regions, notably near 800" and 1000°C, the data seem to obey the familiar parabolic rate law. However, at intermediate temperatures complications were noted, some of which are possibly due to the order-disorder reaction in the silicon-iron solid solution. IN an earlier report&apos; it was shown that during the oxidation of 3 1/4 pct Si-Fe alloys in H2O-H2 atmospheres only silica films were formed in the temperature range from 400° to 1000°C in hydrogen nearly saturated with water at room temperatures, or at dew points as low as -45°C. In the work to be reported here, some observations are made on the rate of oxide film formation. As in the earlier investigation, electron diffraction patterns generally showed either low tridymite or low cristobalite or both, except for some very thin films. These sometimes showed diffuse rings, presumably due to a very small crystallite size, or in a few cases, diffuse bands probably caused by an amorphous film. EXPERIMENTAL PROCEDURE Vacuum-melted silicon iron made of high-purity materials was rolled into strips 0.014 in. thick, and cut into samples 1/2 in. wide by 1 in. long. Chemical analysis showed 3.2 pct Si and 0.002 pct C. All samples were surface abraded with 600-grit paper, were solvent cleaned, and then placed in an paper,apparatus containing a "Gulbransen type"2 micro-balance. Here the gain in weight of the samples of about 5 sq cm area could be followed as a function of time during the oxidation caused by the water in atmospheres of various controlled water-hydrogen ratios. The water-hydrogen ratios can most easily be described as varying from a dew point of 0°C (PH2O-p^2 = 6.2 x 10-3 , to K (P j -40°C (PH2O/PH^= 1.3 X 10-* Most of the experiments were conducted at the 0°C dew-point atmosphere because drier atmospheres caused so little gain in weight that the accuracy of measurement was poor. Because of this, only the data obtained at PH2O,/P,,,= 6.2 x X3 will be reported. The temperature range extended from 800" to 1000°C; and most of the oxidation runs lasted for about 24 hr. The reproducibility of any reading was about ± 1 ?, but the sensitivity of the balance was about 0.2 ?. The atmosphere, flowing at 200 cm per-min, was preheated to the furnace temperature before contacting the specimen. While the gas flow caused a measurable lift on the sample, it was ordinarily sufficiently constant so that it was not an appreciable source of error. X-ray and electron diffraction checks of the samples before and after oxidation showed no evidence of preferred orientation, either on the metal samples or on the silica films formed. EXPERIMENTAL RESULTS The data obtained are summarized in Table I, and some are given in detail in Figs. 1 to 4. In the fourth column of Table I, kp refers to the parabolic rate constant in the expression (?/cm2)2 = kpt + c [1] where ? = micrograms gain in weight kp = parabolic rate constant in units r2 /cm4 t = time in minutes c = constant It will be noted that in many cases no value for kp is given; this is because in these instances the data did not obey the parabolic rate law. The silica film thicknesses given in the last columns are values calculated from the weight gain, an average tridy-mite-crystobalite density, and by assuming a perfectly plane surface. Fig. 1 shows the data plotted in the form of Eq. [I], hence a linear plot indicates parabolic behavior. It has been frequently observed in the literature that

Jan 1, 1960

By E. Teghtsoonian, E. M. Schulson

The tensile behavior of bcc ordered P&apos; AuZn single crystals (CsCl structure) has been investigated under varying conditions of temperature, composition, and orientation. Between -0.2 and 0.4 T, multi-stage hardening occurs fm stoichiometric and nonstoichio-metric crystals oriented near the middle of the primary stereographic triangle. At higher and lower temperatures, parabolic type hardening occurs, followed by work - softening at the higher temperatwes. Deviations from stoichiometry give rise to increased flow stresses. Multi-stage hardening was observed for most orientations, except along the [loll-[lll] boundary and near the [001] corner of the stereo -graphic triangle, where parabolic type hardening occurs. Along two slip systems, (hk0)[001] and (, operate simultaneously while in the [001] comer, slip occurs mainly on the system. Electron microscopy of deformed crystals revealed bundles of edge dislocations forming walls approximately Perpendicular to the glide plane. In general the plasticity of 4&apos; AuZn closely resembles the plasticity of bcc crystals. In recent years, considerable interest has arisen concerning the mechanical properties of the CsCl type intermetallic compounds Ag Mg,&apos;- Fe co,&apos; and Ni Al.&apos;-&apos; The compound P&apos;AuZn is structurally similar. It has a low and congruent melting point of 725"~,&apos;" remains ordered up to the melting point,16 and pos-esses a range of solid solubility from 47.5 to 52.0 at. pct Au at room temperature.15 The present paper reports the results of an investigation on the general tensile behavior of material in single crystal form. Some dislocation configurations characteristic of the deformed state are also reported. The results of a detailed study of the slip geometry in AuZn are presented in a separate paper.17 PROCEDURE Alloy preparation, crystal growing techniques, and the procedure followed in selecting specimens of minimum composition variation are reported elsewhere.17 Dumb-bell shaped tensile specimens were prepared by carefully machining single crystals in a jewellers&apos; lathe to a gage length of 0.80 in. and diam of 0.090 in. Back-reflection Laue X-ray patterns and room temperature tensile tests revealed that machining damage could be eliminated by electrochemically polishing 0.005 in. from the machined surface followed by annealing at 300°C for 1 hr. Specimens were polished in fresh 5 pct KCN solution (40°C, 12 v). Experiments were performed by gripping specimens in a self-aligning pin-chuck and threaded collet system, then straining in a floor model Instron tensile machine. All tests were performed in duplicate. Experimental variables included temperature, composition, and orientation. Unless otherwise stated the strain rate was 2.5 x 10"3 per sec. Liquid testing environments included nitrogen (WOK), nitrogen cooled petroleum ether (133" to 293"K), and silicone oil (293" to 488°K). Resolved shear stress-shear strain curves were electronically computed from autographically recorded load-elongation curves. Stress and strain were resolved on the macroscopic noncrystallographic (hkO) [001] system operative under the specific test conditions of temperature, strain rate, and orientation reported earlier.17 RESULTS The temperature dependence of the work-hardening curves is shown in Fig. 1 for gold-rich crystals of 51.0 at. pct Au oriented near the center of the stereo-graphic triangle. Over the range of intermediate temperatures from -200" to 400°K, they are very similar to those classically observed for fcc metals (reviewed by Nabarro et al.).&apos;&apos; The beginning of deformation is characterized by a region of decreasing hardening rate, stage 0, which is followed by a region of low linear hardening, stage I, and then a region of higher linear hardening, stage 11. At the higher temperatures, stage 111 is observed, a region of decreasing hardening rate. Over the intermediate temperature range, the extent of stage 0 and of the slow transition between stages I and I1 decreases with increasing temperature. Total ductility is large, often greater than 300 pct shear. As the temperature is either increased or decreased, the extent of stage I is decreased, giving rise to parabolic type flow and reduced ductility. Similar temperature effects have been reported for bcc ~r~stals.~~-~~ Below -14O°K, hardening is terminated in brittle fracture while above -400°K. initial hardening is followed first by work-softening and then by chisel-edge type ductile fracture. Stoichiometric (50.0 at. pct Au) and Zn-rich (51.0 at. pct Zn) crystals were also tested from 77" to -500°K. The effect of composition on the flow behavior is illus-

Jan 1, 1970

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

The mechanical properties of unidirectionally solidified A1(rich)-Mg and A1(rich)-Cu castings containing up to 15 wt pct solute have been determined with re -spect to the volume fraction of interdendritic eutectic. Pioperties were determined in the directions pumllel and Perpendicular to that of solidification; the volume fraction of eutectic was varied between the "as-cast" and equilibrizcm amounts by approperiate heat treatment following solidification. The principles of fiber strengthened composites and dispersion strengthened materials are adapted to explain the mechanical properties of these castings. It is generally accepted that castings often have inferior mechanical properties when con~pared to wrought products. However, there is little quantitative data available concerning the factors which make apparently sound castings weak and/or brittle. The relative ease and inexpensiveness of the casting process have always been attractive and, therefore, an understanding of the factors which contribute to the mechanical properties of castings would seem desirable. Such an understanding may lead to an improvement in the mechanical properties to an extent where castings would become competitive in applications where presently only wrought products are considered to have the requisite properties. Such an understanding could also improve the reliability of present cast products. Much of the recent research on castings has centered about determining the extent of segregation in cast alloys. Macrosegregation, particularly inverse segregation, has been studied in some detail 1-8 and a considerable understanding of microsegregation has been obtained.9&apos;10 The effect of solidification rate on dendrite spacing and on the amount of interdendritic eutectic in binary alloys has been established, particularly for Al(rich)-Cu alloys.""0 However, the extension of these ideas to relate the amount of interdendritic eutectic, concentration gradients, micro-segregation, dendrite spacings, and so forth, to the rnechanical properties has been limited. Dean and spear" have related the mechanical properties of an Al-Si-Mg alloy, A356-T62, to the dendrite spacing and have shown that the mechanical properties improve with decreasing dendrite spacing. Passmore et al.12 have shown that annealing at high temperature improves the mechanical properties of Al(rich)-Cu al- loys and Archer and Kempf 13 have shown that an Al-1 pct Mg-1.75 pct Si alloy behaves in a similar manner. Ahearn and Quigley 14 have shown that high temperature homogenization also enhances the mechanical properties of an SAE 4330 steel. However, in the above investigations, no underlying reasons were suggested for the improvement in mechanical properties. The purpose of the present investigation is to relate the mechanical properties of castings to some of the solichfication variables and to derive some equations by which calculations of the mechanical properties may be attempted. In particular, the effect of the amount of interdendritic eutectic and the effect of stress direction with respect to that of solidification on the mechanical properties will be considered. The Al(rich)-Mg and Al(rich)-Cu binary alloy systems were chosen for study. The A1-Mg system was chosen because its constitutional relationships are such that large volunles of eutectic (up to 24 vol pct) may be obtained in the as-cast condition and then be completely dissolved by subsequent heat treatment at about 440°C. This allows a comprehensive study relating the mechanical properties of castings to the amount of interdendritic eutectic. Also the Al(rich)-Mg eutectic is almost a single phase 15 which should make the experimental results more amenable to theoretical interpretation and calculation. The A1-Cu system was chosen for study because of the large amount of related information available concerning segregation, dendrite spacing, and so forth. Unidirectionally solidified castings were used throughout the investigation so that the effect of solidification direction with respect to the direction of applied stress could be determined. THEORETICAL It is well known that upon solidification of binary alloy castings, the nonequilibrium amount of eutectic which forms is given by 10 where fe o is the weight fraction of eutectic, Cs is the solid solubility of solute at the eutectic temperature, k is the equilibrium partition coefficient, and C, is the average composition of the alloy. In the development of Eq. [I], it is assumed that the effects of inverse segregation and diffusion in the solid are negligible, and that no porosity is present. If the casting is homogenized at a high temperature for a long period of time, some (or all) of the eutectic is dissolved and the amount of eutectic for this "equilibrium" condition may be calculated directly from the constitutional diagram. By appropriate intermediate annealing, the

Jan 1, 1970

By S. G. Epstein

Using the capillary-reservoir technique, electromi-gvation rates of cadmium and indium in liquid bismuth were measured at several temperatures. The electric mobility of cadmium Jrom 305° to 535°C and indium from 310° to 595°C can be expressed as a function of temperature by the equations UIn = 1.52 x 10-3 exp sq caz per v-sec Migraion of both solutes was cathode-divected at a rate rnore than four tiMes tHAt previously found for siluer in liquid bisnmth. The electric mobilities of cadmium and indiulrz in liquid bismuth at 500° C are nearly identical with their respective mobilities in mercury at room temperature. AS part of a systematic study of the variables which are considered to control electromigration in liquid metals, the electromigration rates of cadmium and indium in liquid bismuth have been measured. Mass transport properties of silver in liquid bismuth have been reported previously,&apos; and measurements of tin and antimony in liquid bismuth are forthcoming. Comparisons will be made with literature values for these same solutes in mercury.2&apos;3 This series of solutes was selected to determine the effect of the solute valence on its electromigration. Silver, cadmium, indium, tin, and antimony have nearly equal atomic masses but have chemical valences ranging from +1 to +5. They are all fairly soluble in bismuth above 300°C and all have radioactive isotopes, which are an aid in making analyses. EXPERIMENTAL TECHNIQUE Electromigration of cadmium and indium in liquid bismuth was measured by the modified capillary-reservoir technique previously described.&apos; In this method irradiated cadmium or indium is added to bismuth to form alloys containing about 1 wt pct solute (<2 at. pct solute). Several quartz or Pyrex capillaries: 1 mm ID and 5 cm long, vertically oriented, are simultaneously filled in the reservoir of the liquid alloy. A direct current is passed through two of the capillaries, which contain tungsten electrodes sealed in the upper end. The other capillaries sample the reservoir during the experiment. After a measured time interval the capillaries are removed from the reservoir and rapidly cooled. The glass is then broken away from the solidi- fied alloy, which is then weighed, dissolved in acid, and analyzed for solute content by chemical and radiochem-ical techniques. An electric mobility (velocity per unit field) can be calculated from the amount of solute entering or leaving each capillary by the simplified expression1 in which Ui is the electric mobility of the solute, ?mi the solute weight change, Ci the solute concentration of the reservoir, I the current, p the alloy resistivity, and l the duration of the experiment. This expression is valid as long as the experiment is terminated before a concentration gradient develops across the capillary orifice. Earlier experiments showed that the concentration gradient formed initially at the electrode changes with time and eventually reaches the orifice, due to back-diffusion. This condition produces a solute exchange between capillary and reservoir by diffusion or convection, opposing the electromigration, which results in a lower measured value for the electric mobility. To determine if the concentration gradient had reached the orifice, the capillaries used in some of the experiments were sectioned at 1-cm intervals and the solute content of the alloy from each section was radiochemically determined. A typical concentration profile for an experiment with indium in bismuth is shown in Fig. 1; cadmium in bismuth showed similar behavior. As illustrated in the graph, very little back-diffusion has occurred in the capillary containing the cathode, since the concentration gradient is confined to the upper 1 cm of the capillary. In the capillary containing the anode, however, the concentration gradient is much broader, extending nearly to the orifice, even though the net change in solute concentration is nearlv the same in both capillaries. Since cadmium and indium probably lower the density of bismuth when alloyed, depletion of the solute from the alloy adjacent to the anode would increase the density of the liquid in the uppermost region of the capillary. This would give rise to convective mixing within the capillary, causing the broadened concentration gradient. Conversely, the alloy adjacent to the cathode should have a reduced density as the solute concentration is increased by migration, explaining the "normal" concentration profiles found in these capillaries. This disparity was not found for electromigration of silver in bismuth. Both metals have similar densities at the operating temperatures, and nearly symmetrical concentration profiles were found in the two capillaries of each exueriment. This density effect was also apparently encountered when an attempt was made to measure diffusion coefficients for indium in liquid bismuth by the same technique which was successfully used to measure diffusion of silver in bismuth.&apos; Capillaries 1 mm ID and 2 cm

Jan 1, 1968

By Tuncel M. Yegulalp, Malcolm T. Wane

In general, many problems relating to the exploitation of mineral deposits are probabilistic in nature. This derives from the fact that the geologic universe is inherently random. Probability theory and statistics have been found useful for forecasting the behavior of natural events that occur in the geologic universe. The objective of this paper is to illustrate the application of the theory of extremes to this fore-casting problem. For example, it is customary for design purposes to determine the rupture strength of geologic materials. The theory of extremes is exceedingly useful in describing that portion of the frequency distribution of rupture strength which contains the least strengths. Parameters describing the distribution of the least strengths are more important to the designer of mining excavations than parameters describing the total distribution. The basic principles of the theory of extremes will be detailed and illustrated. Any person required to work in the laboratory of nature is aware that uncertainty is a salient feature of all mining enterprises. A mining engineer required to plan the most efficient, practicable, profitable, and safe mine finds himself face to face with numerous ill-understood and often unquantifiable states of nature. Basic information necessary for adequate planning is often lacking or derived from incomplete tests on samples or experience of doubtful validity. The planning procedure usually takes the form of determining a feasible layout with the intent of determining an optimal layout when and if the necessary details and information become available. The crux of the entire procedure is the choosing of numbers to put into the operational and structural models which encompass the plan. Many times these numbers must be assigned qualitatively from past experiences and are called the "most probable ones." At other times, load records, performance records and material tests provide a basis for extrapolation. In any event, the numbers are chosen from a distribution or set of all numbers. Since each number in the distribution represents a possible state, the choice of any particular value is based upon a decision rule. To illustrate, consider the design of an underground structure or the design of a rock slope. The initial step is the formulation of the various possible structural actions which result from the geometry of the layout. For a given structural model various intensities of behavior are possible depending upon the load, deformation, and material characteristic spec-trums, respectively. Of particular interest to mining people is the failure behavior or condition, i.e., when there is a complete collapse of structural resistance by either structural instability or fracture. A necessary feature of the analysis is the "rupture strength" of the material. Information on the rupture strength is derived from testing either in situ or in the laboratory and the usual outcome is a variation in the test results. The methodology used to overcome this variation is to construct a frequency distribution of rupture strengths, and then determine a measure of central tendency and variability. The main idea involved is that the central tendency number will be used in the failure calculations and the measure of dispersion will be used to estimate the probability of failure. In particular if the distribution of rupture strength is normal, the mean rupture strength is the central tendency number and the standard deviation of the rupture strength is the measure of variability. Suppose the mean value of rupture strength is 1000 psi and the standard deviation is 200 psi. Insertion of 1000 psi into the failure calculation produces results that are unsafe, hence a common decision rule is to reduce the mean value by a "factor of ignorance" so that the failure calculation will produce a "safe result." If two is chosen as a factor of ignorance, this means the value inserted in the calculation is 500 psi or 2.5 times the standard deviation. The next step is to determine the percentage chance that failure will occur from a design created on this basis. Tables on the normal distribution function show that this percentage chance is 0.621% or approximately 7 times out of 1000. In practice, however, the situation is more complicated than represented by the foregoing illustration. The laboratory or field testing program usually constitutes a pathetically small sample of the geologic universe of interest and not enough testing is carried out to determine the exact form of the distribution of the test results. The normal, Cauchy and Student&apos;s T distributions are strikingly similar, and it becomes a matter of mathematical convenience to assume the normal law for phenomena which follow other laws.

Jan 1, 1969

By W. M. Baldwin, J. T. Brown

The effect of hydrogen on the ductility, c, of SAE 1020 steel at strain rates, i, from 0.05 in. per in. per rnin to 19,000 in. per in. per rnin and at temperature, T, from +150° to —320°F was determined. The ductility surface of the embrittled steel reveals two domains: one in which and the other in which The usual "explanations" of hydrogen embrittlement are in accord with the first of these domains only. THE purpose of this investigation was a fuller A characterization of this of the investigation effects of varying temperature and strain rate on the fracture strain of hydrogen-charged steel. To be sure, it is known that low and high temperatures remove the embrittlement that hydrogen confers upon steels at room temperature,1 * see Fig. la and b, and that high strain rates have a similar effect,&apos;-&apos; see Fig. 2a, b, and c. However, the general effect of these two testing conditions on the fracture ductility of hydrogen-charged steels is not known, i.e., the three-dimensional graphical representation of fracture ductility as a function of temperature and strain rate is not known—only two traverses of the graph are available. The need for such a graph is not pedantic. To demonstrate this point, Fig. 3a, b, and c shows three of many three-dimensional graphs, all possible on the basis of the two traverses at hand. The important point (as will be developed in the Discussion) is that each of them would indicate a different basic mechanism for hydrogen embrittlement. It will be noted that the four types of ductility surfaces in Fig. 3a, b, and c may be characterized as follows: Material and Procedure Tensile tests were made at various temperatures and strain rates on a commercial grade of % in. round SAE 1020 steel in both a virgin state and as charged with hydrogen. The steel was spheroidized at 1250°F for 168 hr to give the unembrittled steel the lowest possible transition temperature. The steel was charged cathodically with hydrogen as follows: The specimen was attached to a 6 in. steel wire, degreased for 5 min in trichlorethylene, rinsed with water, and fixed in a plastic top in the center of a cylindrical platinum mesh anode. The assembly was placed in a 1000 milliliter beaker containing an electrolyte of 900 milliliters of 4 pct sulphuric acid and 10 milliliters of poison (2 grams of yellow phosphorous dissolved in 40 milliliters of carbon disulphide). A current density of 1 amp per sq in. was used which developed a 4 v drop across the two electrodes. All electrolysis was carried on at room temperature. Temperatures for tensile tests were obtained by immersing the specimens in baths of water (+70° to + 150°F), mixtures of liquid nitrogen and isopen-tane (+70° to —24O°F), and boiling nitrogen (-240" to-320°F). Specimens were tested in tension at strain rates of 0.05, 10, 100, 5000, and 19,000 in. per in. per min. The 0.05 and 10 in. per in. per rnin strain rates were obtained on a 10,000 lb Riehle tensile testing machine, the 100 in. per in. per rnin rate on a hydraulic-type draw bench with a special fixture, and the 500 and 19,000 in. per in. per rnin rates on a drop hammer. The fracture ductility of hydrogen-charged steel at room temperature and normal testing strain rates (-0.05 in. per in. per min) is a function of electro-lyzing time, dropping to a value that remains constant after a critical time.&apos;* Under the conditions of • The hydrogen content of the steel continues to increase with charging time even after the ductility has leveled off to its saturated value.&apos; this research the saturated loss in ductility occurred at approximately 30 min, see Fig. 4, and a 60 min charging time was taken as standard for all subsequent tests. After charging the steel with hydrogen, the surface was covered with blisters. These have been described by Seabrook, Grant, and Carney.&apos; The original diameter of the specimen was not reduced by acid attack, even after 91 hr. Results The ductility of both uncharged and charged specimens is given as a function of strain rate in Fig. 5, and as a function of temperature at four different strain rates in Fig. 6. These results are assembled into a three-dimensional graph in Fig. 7. It is seen that the locus of the minima in the ductility curves of the charged steels divides the ductility surface into two domains. At temperatures below the minima,

Jan 1, 1955

By W. M. McKewan

Magnetite pellets were reduced in flowing hydrogen at pressures up to 40 atm over a temperature range of 350° to 500°C. The rate of weight loss of oxygen per unit area of the reaction surface was found to be constant with time at each temperature and pressure. The reaction rate was found to be directly proportional to hydrogen pressure up to 1 atm and to approach a maximum rate at high pressures. The results can be explained by considering the reaction surface to be sparsely occupied by adsorbed hydrogen at low pressures and saturated at high pressures. PREVIOUS investigation1,2 have shown that the reduction of iron oxides in hydrogen is controlled at the reaction interface. Under fixed conditions of temperature, hydrogen pressure, and gas composition, the reduction rate is constant with time, per unit surface area of residual oxide, and is directly proportional to the hydrogen pressure up to one atmosphere. The reduction rate of a sphere of iron oxide can be described3 by the following equation which takes into account the changing reaction surface area: where ro and do are the initial radius and density of the sphere; t is time; R is the fractional reduction; and R, is the reduction rate constant with units mass per area per time. The quantityis actually the fractional thickness of the reduced layer in terms of fractional reduction R. It was found in a previous investigation2 of the reduction of magnetite pellets in H2-H,O-N, mixtures, that the reaction rate was directly proportional to the hydrogen partial pressure up to 1 atm at a constant ratio of water vapor to hydrogen. Water vapor poisoned the oxide surface by an oxidizing reaction and markedly slowed the reduction. The enthalpy of activation was found to be + 13,600 cal per mole. It was also found that the magnetite reduced to meta-stable wüstite before proceeding to iron metal. The following equation was derived from absolute reaction-rate theory4,8 to expfain the experimental data: where Ro is the reduction rate in mg cm-2 min-&apos;; KO contains the conversion units; Ph2 and PH2O are the hydrogen and water vapor partial pressures in atmospheres; Ke is the equilibrium constant for the Fe,O,/FeO equilibrium; Kp is the equilibrium constant for the poisoning reaction of water vapor; L is the total number of active sites; k and h are Boltzmann&apos;s and Planck&apos;s constants; and AF is the free energy of activation. Tenenbaum zind Joseph5 studied the reduction of iron ore by hydrogen at pressures over 1 atm. They showed that increasing the hydrogen pressure materially increased the rate of reduction. This is in accordance with the work of Diepschlag,6 who found that the rate of reduction of iron ores by either carbon monoxide or hydrogen was much greater at higher pressures. He used pressures as high as 7 atm. In order to further understand the mechanism of the reduction of iron oxide by hydrogen it was decided to study the effect of increasing the hydrogen pressure on rebduction rates of magnetite pellets. EXPERIMENTAL PROCEDURE The dense magnetite pellets used in these experiments were made in the following manner. Reagent-grade ferric oxide was moistened with water and hand-rolled into spherical pellets. The pellets were heated slowly to 550°C in an atmosphere of 10 pct H2-90 pct CO, and held for 1 hr. They were then heated slowly to 1370°C in an atmosphere of 2 pct H2-98 pct CO, then cooled slowly in the same atmosphere. The sintered pellets were crystalline magnetite with an apparent density of about 4.9 gm per cm3. They were about 0.9 cm in diam. The porosity of the pellets, which was discontinuous in nature, was akrout 6 pct. The pellets were suspended from a quartz spring balance in a vertical tube furnace. The equipment is shown in Fig. 1. Essentially the furnace consists of a 12-in. OD stainless steel outer shell and a 3-in. ID inconel inner shell. The kanthal wound 22 in. long, 1 1/2, in. ID alumina reaction tube is inside the inconel inner shell. Prepurified hydrogen sweeps the reaction tube to remove the water vapor formed during the reaction. The hydrogen is static in the rest of the furnace. The sample is placed at the bottom of the furnace in a nickel wire mesh basket suspended by nickel wire from the quartz spring. The furnace is then sealed, evacuated, and refilled with argon several times to remove all traces of oxygen. It is then evacuated, filled with

Jan 1, 1962

By Robert E. Green, Wolfgang Sachse

Simullaneous ultrasonic attenuation measurements of both quasishear waves propagating in single cryslals of aluminum indicate that, in the undeformed annealed state, the dislocation density is generally not uniform on all slip systems. Change oof attenuation measurements made during plastic defortnation of crystals , which possessed specific orientations ideal for studying the orientation dependence of dislocation damping, indicate that, for low strain levels, dislocation motion occurs on additional slip systems besides the primary one, even for crystals oriented for plastic deformation by single slip. THE sensitivity of internal friction measurements permits such measurements to be used successfully in studying the deformation characteristics of metal crystals. On the basis of experimental observations, T. A. Read1 was the first to associate internal friction losses with various dislocation mechanisms. Since that time further work2-&apos; has been performed and a dislocation damping theory has been formulated by Granato and Lucke.6 In the amplitude independent region, this theory predicts the attenuation a to be dependent on an orientation factor O, a dislocation density A, and an average loop length L. if is a constant, independent of crystallographic orientation. For a given crystallographic orientation, changes in dislocation density and loop length give rise to the observed attenuation changes accompanying plastic deformation. The Granato-Liicke theory suggests the investigation of the orientation dependence of attenuation measurements in hopes of obtaining information to separate dislocation motion losses from other losses.7 An experimental study of the orientation dependence of attenuation in undeformed annealed single crystals should yield an insight into the uniformity of dislocation distribution throughout the entire specimen. A similar study on crystals plastically deformed in a prescribed fashion should give information about the alterations in the dislocation distribution on the slip systems activated during plastic deformation. The possible modes of elastic waves which can be propagated in aluminum,8 copper,9 zinc,10 and other hexagonal metals" have been calculated. Associated with each mode of wave propagation are dislocation damping orientation factors, which are based on the resolution of the stress field of that particular elastic wave onto the various operative slip systems in the material. These orientation factors have also been calculated as a function of crystallographic orientation in the papers cited above. Einspruch12 obtained agreement between predicted and observed attenuation values of longitudinal and shear waves in (100) and (110) directions of two undeformed aluminum crystal cubes. He ascribed the slight deviations between predicted and observed values to a nonuniform dislocation distribution, or to other loss mechanisms. In shear deformation of zinc crystals, Alers2 found that the attenuation of shear waves having their particle displacements in the slip plane was very sensitive to the deformation, while the longitudinal wave attenuation was affected only when the wave propagation direction was not normal to the slip plane. Using aluminum single crystals oriented for single slip, Hikata3 et al. found that during tensile deformation the change of attenuation of the shear wave (actually quasishear) having particle displacements nearly perpendicular to the primary slip direction exhibited the easy-glide phenomena, while longitudinal waves did not. Similar results were reported by Swanson and Green5 during compressive deformation of aluminum crystals. These results are in qualitative agreement with the calculated orientation factors for specimens of this orientation. In well-annealed, undeformed aluminum crystals, the damping is expected to be due to dislocations vibrating on all twelve slip systems. The orientation factors associated with this initial damping will be designated by O2 and O3, where a, represents the average orientation factor for the slow shear (or quasishear) wave and O3 represents the average orientation factor for the fast shear (or quasishear) wave. The calculation of these values for aluminum crystals by Hinton and Green8 shows that they vary very little as a function of crystallographic orientation—at most, by a factor of 2.47. If the dislocation density and loop length are uniform, then in the initial undeformed state, Here the subscript zero refers to the initial value of the attenuation. Also for aluminum, the calculations8 show that the orientation factors for primary slip only, associated with each shear wave, exhibit a sharp minimum for particular crystallographic orientations. A composite plot of the two shear wave orientation factors for primary slip only is shown in Fig. 1. Since these orientation factors are associated with dislocation motion occurring on the primary slip system only, the proper condition to check these factors might be attained by slightly deforming a single crystal oriented for primary slip. For dislocation motion on the primary slip system only,

Jan 1, 1969