Search Documents

By N. Mungan

Formation damage, i.e.. reduclion in permeability, has been generally attribuled to clay minerals which expand or disperse upon contact with water that is less saline than the connate water. Luboratory, studies show that penneahility reduction can also occur in formalions containing only nonexpandable clays such as illite or kaolinite, and can be caused also by changes in pH. Furthermore, pH changes can damage even formations that are essentially free of clays. It is suggested that permeability reduction is due to the small passages being blocked by particles, which may be dispersed clays, cemenlion material or other fine parricles. These particles are dislodged by dispersion of clays due to changes in salinity or by dissolution of calcareous cement by acids, or of silicaceous cement by alkaline solutions. In working with reservoir cores, it was found that extracted cores damaged more easily and extensively than nonextracted cores. The extent of damage depended also on tenlperatltre. INTRODUCTION Permeability is an important property of porous media and has been the subject of many studies by engineers and geologists. Many of these studies are conccrned with formation damage, i.e., reduction in permeability, resulting from exposure of oil-producing formations to water substantially less saline than the connate water. This effect causes understandable concern since during drilling, completion and production phases formations are often exposed to fresh water. The damage resulting from contact with relatively fresh water has been attributed to expansion and dispersion oF clay minerals. During laboratory investigation of the use of NaOH as a wettability reversal agent to increase oil recovery from oil-wet reservoirs, several cores used in the displacement studies suffered loss in permeability. Despite the traditional usage of NaOH for conditioning aqueous mud systems, the role of the caustic filtrate in wellbore damage seems to have been overlooked. Browning2 as recently reported on the effects of NaOH in dispersing clay minerals but he was concerned only with complications that may arise in drilling massive shale beds. The following study was made to examine the role of pH and salinity changes in core damage. Where cores from reservoirs were used, tests were performed with extracted and nonextracted cores both at room and reser- voir temperatures, since it was felt that the test environment and core condition may affect the results. Because of its limited coverage and exploratory nature, this study is not intended to provide answers to field formation damage problems. It is hoped that it will encourage research into new aspects of the permeability reduction problems, particularly those allied to new recovery and production processes. PROCEDURE In all permeability tests, fluids were pumped through the cores at a constant volumetric rate. Only deaerated fluids and reagent grade chemicals were used. The fluids were passed through two ultrafine filters before injection to remove any entrained particles. The cores, with the exception of the unconsolidated cores, were mounted in Hassler holders. Water was used to transmit pressure to the sleeve. The inlet endpiece had two entry ports which permitted scavenging one fluid with another to avoid any mixing in the small holdup volume. The cores were flushed with CO2 gas, evacuated for 5 to 6 hours and saturated with the first liquid at a pressure of 1,000 psi for 24 hours to eliminate any free gas from the cores. Pressure differences up to 20 psi were measured by transducers, calibrated in inches of water and continuously recorded. For greater pressure drops, gauges were used. All reservoir cores were cleaned with a light refined mineral oil, then with heptane, and finally dried with CO2. Compatibility tests showed that no precipitates formed when mineral oil and the crudes were mixed. Some cores were extracted in Dean-Stark-type solvent extractors using xylene and trichloroethane and dried in a vacuum oven at 450F. Each test consisted of a sequence of water, test solution, and again water flow. RESULTS AND DISCUSSION STUDIES IN BEREA CORES Salinity Contrast Berea cores 2-in. in diameter and 12-in. long were cut from sandstone quarried in Cleveland, Ohio. The clay minerals were identified by X-ray diffraction to be chlorite, kaolinite, illite and incerlayered illite. Flow of fresh water or 30,000 ppm brinc does not cause any permeability reduction (Fig. 1). However, after injection of brine the core is readily damaged by fresh water. Damage starts almost instantly as the fresh water injection is begun, and at a cumulative injection of 1.2 PV fresh water, the permeability has dropped from 190 to 0.9 md. Upon continued injection, the effluent contains clay minerals dislodged from the core. The final core per-

Jan 1, 1966

By G. R. Leverant, C. P. Sullivan

To evaluate the effect of a dispersion of nondeform-able, incoherent, second-phase particles on high-strain cyclic deformation and fracture, recrystallized TD-nickel (Ni-2ThO2) and a commercially pure nickel, Ni-200, were fatigued under strain control at total strain ranges varying from 0.009 to 0.036. Relative to the Ni-200, the slip at the surface of the TD-nickel was more wavy and discontinuous due to the presence of the thoria particles. This made crevice formation (incipient cracking) within slip bands more difficult in TD-nickel than in Ni-200. Both materials cyclically hardened to a constant (saturation) flow stress which increased with increasing plastic strain amplitude. Cellular substructures were developed in both materials during cycling. The cell size in TD-nickel was controlled by the thoria particle distribution and was independent of plastic strain amplitude over the range investigated. The cell size in Ni-ZOO was larger than that in TD-nickel at similar plastic strain amplitudes and was a function of plastic strain amplitude. These results, together with the cyclic stress-strain curves for both materials, are discussed in terms of a model for fatigue strain accommodation at saturation recently proposed by Feltner and Laird. NUMEROUS fatigue investigations have considered the interrelation of slip character, dislocation substructure, and cracking in pure metals and solid-solution alloys. However, except for the studies of the low-strain fatigue of internally oxidized copper alloys1 and cast, dispersion-strengthened lead,&apos; little is known about the effects which small, incoherent, nondeform-able, second-phase particles have on cyclic deformation and cracking processes. Effects due to the particles alone are often obscured by a dislocation substructure introduced during thermomechanical processing of dispersion-strengthened metals. In the present study, recrystallized TD-nickel and a commercially pure nickel, Ni-200, were employed to evaluate the effect of a thoria dispersion on high-strai fatigue deformation and cracking at room temperature. I) MATERIAL AND EXPERIMENTAL PROCEDURE The TD-nickel was supplied by DuPont as a 5/8-in.-thick stress-relieved plate which had been subjected to a proprietary schedule of thermomechanical treatments, and the Ni-200 as 3/4-in. bar which was subsequently annealed for 2 hr at 850°C in argon resulting in an average grain diameter of 0.05 mm. The compositions of these materials are given in Table I. The microstructure of the TD-nickel consisted of elongated grains parallel to the primary working direction with an average width of 0.16 mm, Fig. l(a). Many fine annealing twins were present indicating that the starting material was in a recrystallized condition; this supposition was confirmed by the absence of of any extensive dislocation substructure, Fig. l(b). Sheetlike stringers parallel to the rolling direction were occasionally seen both within grains and at grain boundaries. Some approximately spherical particles about 2 in diam, which may correspond to exceptionally large thoria particle aggregates, were also present. The average Young&apos;s modulus of the plate material in the rolling direction was 21.8 X 106 psi which is consistent with a {100}<001>recrystalliza-tion texture3&apos;* being prominent. In transmission microscopy, the 2.3 vol pct of thoria particles generally appeared to be uniformly distributed although some clusters, 0.1 to 0.3 in diam, of larger particles were observed as previously reported for TD-nickel sheet,5 and stringering of particles was present in some areas as welt. The average diameter of the thoria particles was 450A with a calculated mean planar center-to-center spacing of 2100A, as determined by quantitative metallographic analysis.= The 0.2 pct offset yield stress was 36,000 psi which agrees with the value predicted by the modified Orowan relation7 for edge dislocations bowing between thoria particles of the size and spacing observed in the present investigation. Fig. 2 illustrates the specimen design employed for the axial high-strain fatigue testing. Adapters were screwed onto the threaded portions of each specimen so that testing could be performed in the same manner as that reported for buttonhead specimens.8 Stressing was coincident with the working direction for both materials. The gage section of each specimen was electropolished and lightly etched prior to testing. The total strain was controlled, being varied between zero and a maximum tensile strain ranging from 0.009 to 0.036. In addition to these tests, a circum-ferentially notched TD-nickel specimen was cycled over a total strain range of 0.0075. The same strain

Jan 1, 1969

By D. G. Westlake

The electrical resistance of a series of V-H alloys (0 to 3.5 at. pct H) has been measured over the temperature range G° to 360°. Interstitial impurities made contributions to the residual resistivity, but not the ideal resistivity. The contribution of hydrogen in solid solution is expressed by Ap = 1.12 microhm-cm per at. pct H; but the contribution of precipitated hydride was negligible. A portion of the so1vu.s for the V-H phase diagram is presented. The solubility limit is given by In N (at. pct H) = (5.828 i 0.009) - (2933 i 44)/RT. Comparison of critical temperatures joy hydride precipitation and published critical temperatures for hydrogen embrittlement suggests the two are related. ThiS study was initiated as part of an investigation of the mechanism by which small concentrations of hydrogen embrittle the hydride-forming metals at low temperatures. It has already been shown that, in the case of hcp zirconium, a reduction in ductility accompanies the strengthening resulting from precipitation of a finely dispersed hydride phase.&apos;&apos;&apos; Our attempts to detect a similar precipitation of a second phase at low temperatures in V-H alloys by transmission electron microscopy have been thwarted because we have been unable to prepare thin foils that are representative of the bulk material with respect to hydrogen concentrati~n.~&apos;~ The present investigation establishes the solvus of the V-H system at subambient temperatures. Subsequently, we hope to be able to determine whether the embrittlement temperature is related to the critical temperature for precipitation of the hydride in a given V-H alloy. veleckis5 has proposed a partial phase diagram for the V-H system based on extrapolations of the pressure-composition relations he measured at higher temperatures. Kofstad and wallace&apos; conducted a similar study of single-phase alloys but did not attempt to establish the phase diagram. Zanowick and wallace&apos; and ~aeland&apos; have studied a portion of the phase diagram by X-ray diffraction, but they investigated no alloys in the hydrogen concentration range 0 to 3 at. pct, the range of interest to us. EXPERIMENTAL PROCEDURE The vanadium was obtained from the Bureau of Mines, Boulder City, Nev., in the form of electrolytic crystals. The analyses supplied with them listed 230 ppm by weight metallic impurities, 20 ppm C, 100 ppm N, and 290 ppm 0. The crystals were electron-beam-melted into an ingot that was rolled to 0.64 mm. Strips, 60 mm long and 4.2 mm wide, were cut from the sheet, and both rolled surfaces were ground on wet 600-grit Sic paper to produce specimens 0.4 mm thick. They were wrapped in molybdenum foil, vacuum-encapsulated in quartz, and annealed 4 hr at 1273°K. The specimens were annealed in a dynamic vacuum of 2X lo-&apos; Torr for 30 min at 1073°K for dehydrogenation, and charged with the desired quantity of hydrogen by allowing reaction with hydrogen gas at 1073°K for 2 hr and cooling at 100°K per hr. Purified hydrogen was obtained by thermal decomposition of UH3. Sixteen specimens were studied: two contained no hydrogen and the others had hydrogen concentrations between 0.5 and 3.5 at. pct (hydrogen analyses were done by vacuum extraction at 1073°K). Electrical resistances were measured by the four-terminal-resistor method on an apparatus similar to the one described by Horak.~ The specimen holder was designed so that both current and potential leads made spring-loaded mechanical contact with the specimen. The potential leads were 30 mm apart, and the current leads were 55 mm apart. The current was 0.10000 amp. We used the following baths for the indicated temperature ranges: liquid nitrogen, 77°K; Freon 12, 120" to 230°K; Freon 11, 230" to 290°K; and ethanol, 290" to 340°K. Temperatures lower than 77°K were achieved by allowing the specimen to warm up after removal from liquid helium. Temperatures above 77°K were measured by a calibrated copper-constantan thermocouple (soldered to the specimen holder) and below 77°K by a calibrated carbon resistor. The temperature of the bath changed less than 0.l0K between duplicate measurements of the resistance. RESULTS AND DISCUSSION Typical plots of resistivity p vs temperature T are shown in Fig. 1. In the interest of clarity, only five curves are presented and the data points have been

Jan 1, 1968

By A. McLean, D. W. Gomersall, R. G. Ward

An experimental study has been made of the solubility of nitrogen in liquid iron and liquid Fe-C alloys using levitation melting and a rapid quenching device. Iron alloy droplets were equilibrated with nitrogen gas at 1 atm pressure, quenched, and analyzed. Previous techniques for studying the Fe-C-N system have produced data which me in marked disagreement. This disagreement is due largely to errors caused by reaction between the molten alloy and the crucible material. With the levitation procedure, errors from this source have been eliminated and precise solubility data obtained for temperatures between 1450° and 1750°C. C-N interactions in molten iron have been expressed in terms of first- and second-order free energy, enthalpy, and entropy parameters. ALTHOUGH the solubility of nitrogen in iron base alloys is in general small, the effects of nitrogen on the properties of steel may be quite profound. For most purposes nitrogen in finished steels is undesirable, particularly in the low-carbon grades, since on cooling to room temperature the solubility limit of nitrogen in the steel may be exceeded and this can lead to embrittlement and loss of ductility on aging. On the other hand, nitrogen can improve the work-hardening properties and machinability of steels while in certain stainless grades nitrogen is important in order to stabilize the austenite phase. It is, therefore, desirable that one should be able to predict the solubility of nitrogen in liquid iron alloys. To do this, information is required concerning the interactions between nitrogen and the various alloying elements which may be present in liquid iron. There have been several investigations of these effects in recent years1"7 and the interactions between nitrogen and many elements dissolved in liquid iron are now known to a high degree of precision at steel-making temperatures. Unfortunately, a number of iron alloy systems which are of interest in steelmaking have been difficult to deal with by the experimental techniques generally used for this type of investigation. Among the most important of these are the Fe-C alloys. In the past, two methods have been widely used for determining nitrogen solubilities: the Sieverts&apos; technique, in which the amount of nitrogen required to saturate a given mass of liquid metal at a particular temperature and pressure is measured volumetrically, and the sampled-bath technique in which liquid metal held in a crucible is equilibrated with a gas phase containing a known partial pressure of nitrogen, and samples drawn from the melt are quenched and analyzed. These two methods have been discussed in detail elsewhere.5,8 With the Sieverts&apos; technique, errors may be introduced from the following sources: i) Gas adsorption on metal films which have condensed on the cooler parts of the reaction chamber. ii) Uncertainty in the determination of "hot volume" calibrations. iii) Crucible-melt interaction, particularly if a gaseous reaction product is formed or if the melt becomes contaminated with material from the crucible walls. The sampled-bath method may also suffer from errors due to reaction between the melt and the crucible material. In addition, there is the possibility that gas may be lost from the sample during solidification and cooling. In the present investigation, the solubility of nitrogen in liquid iron alloys has been studied by means of a new technique based on the use of levitation melting equipment and a rapid quenching device. In addition to the fact that problems of the type outlined above are avoided, this particular approach has the following advantages: i) The high-frequency current induces vigorous stirring within the levitated droplet so that gas-metal equilibration is rapidly attained. ii) The gas phase surrounding the melt can be changed very quickly and is easily controlled. For example, a droplet may be levitated in helium, deoxidized in hydrogen, equilibrated with nitrogen, and quenched, within a period of 15 min. ii) Melts can be readily under cooled or superheated, thus extending the effective temperature range of an investigation and allowing temperature-dependent data to be determined with a high degree of precision. Excellent reviews of levitation melting techniques and their application to physical-chemistry studies at high temperature have been published recently by Jenkins et al.,9 Peifer,10 and Rostron.11 In the present investigation a levitation melting technique has been used to obtain data for the solubility of nitrogen in pure liquid iron and liquid Fe-C alloys at temperatures between 1450° and 1750°C. The solution of nitrogen in liquid iron can be described by the reaction:

Jan 1, 1969

By D. Mills, G. B. Craig

ZIRCONIUM has been float-zone-refined in an electron-beam furnace and the redistribution of oxygen, iron, and tungsten has been measured. The iodide zirconium used in the present experiments initially contained both oxygen and iron in the range 100 to 200 ppm by weight, and tungsten in amounts less than 10 ppm. Float-zone refining of zirconium, using induction heating, has previously been attempted by Kneip and Betterton1 and Langeron.2 Kneip and Betterton were primarily interested in the removal of iron and nickel. They achieved some purification, with respect to both these elements, in material given up to six passes at 150 mm per hr under an argon atmosphere. Langeron reported purification for a large number of elements after four completed passes in a static vacuum system operating at pressures less than 10-6 mm of Hg and with rates of zone travel between 30 and 40 mm per hr. He did not report oxygen concentrations, but stated that there was inverse segregation of this element. westlake5 has also used a minimum of four zone passes to remove iron from zirconium. No rates or conditions were given. Belk6 has shown that tungsten contamination can occur during electron-beam melting. He reported an increase from 0.01 to 0.05 pct by weight tungsten in molybdenum after four complete zone passes. The vacuum system for the electron-beam unit used in the present investigation consisted of a single-stage rotary pump backing a liquid-air-trapped oil-diffusion pump. A pressure of less than 10-8 mm of Hg was obtained with Dow-Corning 705 fluid in the diffusion pump. In order to avoid contamination of the zirconium by evaporation from the tungsten filament, a special focusing cage4 was employed. Three rates of zone travel, 114, 38, and 4 mm per hr, were investigated, with the liquid zone moving from the bottom to the top of the bar. The bars used were 3 mm diam and the total melted length was 130 mm. Oxygen was analyzed by neutron-activation analysis using a neutron flux of 10&apos; neutrons per sq cm per sec to form the isotope N16 by the 016 (n,p) N16 reaction. The standard deviation is quoted for each oxygen determination. The iron and tungsten analyses were performed spectrographically, and the precision is estimated to be ±6 pct. Analyses for tungsten were all below the detectable limit of 10 ppm, confirming the protection given by the focusing cage. No significant redistribution of oxygen was found at the two higher rates of zone travel. The redistribution of oxygen and iron obtained after five passes at a rate of zone travel of 4 mm per hr (1.1 x 10-4 cm per sec) is recorded in Tables I and 11. Burris, Stockman, and Dillon3 have estimated the distribution of solutes during multipass zone refining. Using their curves for an effective distribution coefficient of Keff = 1.5 and 2 for oxygen and Keff = 0.3 for iron, the expected concentrations were estimated for the present material. These are shown in Tables I and 11, along with the experimentally measured values. The calculated concentrations are based on a molten-zone length of 10 pct of the total melted length, whereas in the present experiments the molten zone was approximately 5 pct of the total melted length. The effect of zone length on solute redistribution is most pronounced after a large number of zone passes. Comparison with Pfann&apos;s published data8 for solute redistribution for various Keff&apos;s and zone lengths indicates only small differences at low numbers of completed zone passes. It is evident that the expected distribution has not been realized in the case of iron. Scrapings of material deposited on the inner surfaces of the electron-focusing cage were found to be magnetic. It is, therefore, concluded that redistribution of iron is masked by evaporation. Deposition of iron inside the focusing cage was observed at all rates of zone travel. The results of the investigation may be summarized as follows: 1) Segregation of oxygen is typical for a solute with a distribution coefficient, Keff > 1. 2) No redistribution of oxygen takes place at high rates of zone travel. 3) The distribution coefficient for oxygen lies between 1.5 and 2.0. 4) Purification with respect to iron occurs mainly by evaporation. 5) The focusing cage is effective in preventing tungsten contamination.

Jan 1, 1967

By H. Conrad, V. V. Damiano, G. J. London

The slip mode activated during the c axis compression of single crystals of commercial-purity ingot SR beryllium, high-purity (twelve-zone-pass) beryllium, and Be-4.4 wt pct Cu and Be-5.2 wt pct Ni alloys in the temperature range of 25° to 364°C was determined using two-surface slip trace analysis, slip-step height analysis, and electron transmission microscopy. All three techniques indicated the occurrence of copious pyramidal {1 122) (1123) slip in the alloys over the entire temperature range, the amount increasing with temperature. Pyramidal slip was also indicated in the high-purity beryllium by slip trace analysis and electron transmission microscopy, but the amount was somewhat less than in the alloys. For the commercial-purity ingot crystals, only a very small number of pyramidal slip lines were observed, and these were in the immediate vicinity of the fracture surface. No pyramidal dislocations could be detected by electron transmission microscopy in this material. Dislocatransmissiontions with Burgers vectors [0001] and +(ll20) were identified by electron transmission microscopy inthe (1122) slip bands, as well as those with the j (1123) vector. This was interpreted to indicate that the edge components of the 3(1123) vector dislocations activated during c axis compression dissociate upon unloading according to the reaction i (1123) — [0001] + 3(1120) THE microstrain c axis compression of single crystals of commercial-purity ingot SR beryllium (99.6 pct), high-purity twelve-zone-pass beryllium (99.98 pct), Be-5.24 pct Ni and Be-4.37 pct Cu alloys was described in a previous paper.1 This paper covers in detail the analysis of slip traces observed on two mutually perpendicular lateral surfaces of these specimens, and a detailed description of transmission electron microscopy studies performed on foils cut from the bulk crystals after they had been deformed to fracture in the c axis compression. Observation of slip traces on single surfaces of deformed single crystals are generally insufficient to positively identify slip or twinning modes. The use of two carefully cut and oriented perpendicular surfaces can greatly aid in the positive identification and index- ing of slip traces, although even this technique may be quite inadequate if more than one type of slip system operates and if an insufficient number of traces are observed on the surfaces. The problem is greatly simplified for symmetric cases like that for c axis compression of an hep crystal such as beryllium, in which the operating slip systems are all equally inclined to the direction of the applied stress, and each slip system of a given slip mode has an equal chance of operating. For such cases, the traces of any given slip mode observed on the surfaces cut parallel to the c axis are symmetrically tilted about the c axis. It is therefore possible to quickly determine whether one or more slip modes are operating. Confirmatory evidence in support of the observations made on the external surfaces can be obtained from foils cut from the deformed crystals and examined by transmission electron microscopy. This latter technique serves to identify not only the operating slip plane but also the Burgers vector of the dislocations which participate in the slip. For this purpose, a simplified technique based upon a double tetrahedron notation is used in the present paper. The planes and directions in the hep lattice are all designated by letters rather than indices and extinction conditions are easily determined if the Burgers vector lies in the plane contributing to the diffraction. RESULTS 1) Slip Trace Analysis. The standard (0001) stereo-graphic projection of beryllium is shown in Fig. 1. The two mutually perpendicular, lateral surfaces of the compression specimen are represented by the diametrical planes AA&apos; and BB&apos;, also referred to as surface A and surface B. For the specific case represented (a Be-5.24 pct Ni specimen deformed by c axis compression at room temperature), the A surface is tilted 5 deg to the (10i0&apos;) plane and the B surface is tilted 5 deg to the (1120) plane. Two surface trace analyses may be facilitated by examining in turn the intersection of various great circle traces of specific pyramidal planes with two surfaces and comparing the angles made with the (0001) plane with those actually observed on the two surfaces. One then identifies the slip traces by trial and error on a best-fit basis. The (1122) type planes (it was found that slip occurred on these planes) are shown plotted on the stereographic projection in Fig. 1. One obtains directly the angles between the (0001) plane and the {1122) traces by measuring the angle from the periphery to the point of intersection along the lines

Jan 1, 1969

By John L. Miles

John L. Miles (Arthur D. Little, 1nc.)—Pollack and orris" have reported measurements on electron tunneling through A1-A12O3-A1 sandwiches in which the oxide was formed by gaseous oxidation in a glow discharge. From these measurements they deduced the asymmetry of the barrier and, since this is small, conclude that the mechanism suggested by Mott19 for the growth of oxide in thin A12O3 films is inapplicable. In earlier papers20 Pollack and Morris report similar work for oxide films grown thermally. In this case they find a greater asymmetry and conclude that the Mott mechanism is valid. I would like to point out that both these conclusions are quite unjustified. Mott suggests that the growth of the oxide film on aluminum results from the passage of ions through the already present film of oxide under the action of an electric field. This field results from a constant voltage which is in effect a contact potential between metal on one side of the barrier and adsorbed oxygen ions on the other side of the barrier. The theory does not require that the oxide grown is nonuniform either in stoichiometry or structure. It does however specifically assume that the partial layer of ionized oxygen on the surface remains adsorbed on the surface of the growing oxide. In other words, the so-called "built-in field" remains in the oxide only as long as the ionized oxygen is present. When a counter electrode of aluminum is deposited on the oxide, it will react with the adsorbed oxygen on the surface of the oxide, thus forming a small additional amount of oxide. It is clear, then, that there is no requirement in the Mott theory of oxide growth which would necessitate tunneling currents through an Al-A1203-A1 sample to be different when the polarity is reversed. Neither does the theory eliminate the possibility that some additional mechanism could cause the tunneling barrier to be asymmetric and hence tunneling currents to be a function of polarity in such a sandwich. Thus these tunneling-currents measurements are not germane to the question of whether the Mott mechanism is the true method of growth of aluminum oxide films. In fact, it is not surprising that there should be a difference between the oxide properties at the two interfaces (with resulting asymmetry in the tunneling barrier) since the growth conditions and growth rates must have been quite different at these two positions. S. R. Pollack and C. E. Morris (authors&apos; reply)— The point raised by Miles above is one has caused some confusion in the past. The following is an attempt to clarify this point. The built-in field which is responsible for the growth of the thermal oxide at low temperatures arises, according to Mott, because of the passage of electrons from the Fermi surface of the oxidizing metal to surface states introduced by the adsorbed oxygen. It is assumed that the energy of these surface states lies below the Fermi energy of the metal. Electrons therefore continue to flow from the metal to the surface until the built-in electric field raises the potential energy of the surface states to the value of the Fermi energy in the metal, at which time equilibrium is obtained between the surface states and the metal. That is in equilibrium as many excess electrons pass from the metal to the surface per unit time as vice versa. The surface of the oxide prior to deposition of a metallic counterelectrode can then be pictured as follows. The Fermi energy lies in the energy gap of the oxide and is essentially pinned at the energy of the oxygen surface states. The vacuum work function of the oxide is then given by the sum of the electron affinity of the oxide (i.e., the difference in energy between the vacuum and the conduction-band minimum) plus the energy difference between the conduction-band minimum and the Fermi energy. The deposition of a metal onto the surface of the oxide can result in a transfer of electrons across the extremely thin oxide only if there is a contact potential difference between the deposited metal and the parent metal or oxide. That is if the vacuum work function of the deposited metal differs from that of the parent metal, then charge can be redistributed across the oxide in order to equilibriate the Fermi energy across the structure. (It should be

Jan 1, 1965

By R. M. Waterstrat

The phases occurring in the V-Mn system were studied by means of X-yay diffraction and metallo-paphic techniques, using are-melted alloy specimens annealed in the temperature range 800° to 1150°C and quenched. The bcc solid solution extends at 1250°C all the way from vanadium to 6-manganese. Below 1050°C the a-phase is formed, and the terminal a-manganese phase is stabilized up to about 900°C by vanadium in solid solution. IN the only previous general survey of the V-Mn system Cornelius, Bungardt and Schiedtl reported the existence of three intermediate phases corresponding to the approximate compositions VMn,, VMn, and V5Mn. The phase VMn8 has recently been identified as a o phase2 but the alloy VMn was found to have a bcc structure2 corresponding apparently to the vanadium solid solution rather than to the large cubic unit cell reported by Cornelius et al. 1 Subsequent work by Rostoker and Yamamoto3 has shown that the vanadium-base bcc solid solution extends to at least 15 pct Mn at 900°C. An alloy corresponding to the composition VMn, was examined by Elliott,4 who reported that the as-cast sample as well as samples annealed at 1200o and 1300°C had bcc structures, but that annealing at 1000°, 800") and 600°C produced two phases. One of these phases was apparently the bcc solid solution and the other resembled the o phase structure. Hellawell and Hume-Rothery5 established the phase relationships in manganese-rich alloys above 1000°C, and showed that the o phase in this system is replaced by the 6 Mn (bcc) solid solution at temperatures above 1050°C. These results suggest that a continuous bcc solid solution may exist above 1050°C between vanadium and 6 Mn. The present investigation was undertaken in order to develop more complete information in regard to this system. EXPERIMENTAL METHODS The alloys used in the present work were prepared by arc-melting electrolytic manganese having a minimum purity of 99.9 pct and vanadium lumps with a purity of 99.7 pct. The major impurities present in these metals were carbon, nitrogen, and oxygen and this would account for the small percentage of nonmetallic inclusions observed metal-lographically. The arc-melting was at first performed under a helium atmosphere and it was necessary to keep the melting times as short as possible in order to minimize the loss of manganese by vaporization. It was later found that the evaporation of manganese was considerably reduced when the melting was done under argon atmosphere. The final composition of each alloy was calculated by assuming that the total weight loss during melting was due to evaporation of manganese. Compositions which were calculated in this manner agreed reasonably well with the results of chemical analysis, as shown in Table I. Spectrographic analysis revealed the presence of contamination by tungsten, but in no case was the percentage of tungsten greater then 0.4 at. pct. The specimens were in each case broken in half and the fractured section was examined visually and microscopically for evidence of inhomogeneity. Each specimen was homogenized at temperatures near l100°C, as shown in Table I. After this treatment most specimens consisted of large columnar grains of the bcc vanadium solid solution. The etchant used in most of the metallographic work consisted of 20 pct nitric acid, 20 pct hydro-flouric acid, and 60 pct glycerine. It was found that this etchant would clearly delineate the phases present in these alloys although it does not produce any striking contrast between the phases. For certain manganese-rich alloys, a 1 pct aqueous solution of nitric acid was used. This etchant gave a brown color to the a-manganese phase, whereas the o phase was virtually unattacked and appeared very light as shown in Fig. 1. The etchants used by Cornelius et a1.l were found to produce spurious effects in some of these alloys. In particular, the vanadium-rich alloys etched in hot sulfuric acid often appeared to consist of two phases when both X-ray diffraction and etching with the glycerine-acid mixture indicated the presence of single phase bcc solid solution. A few percent of what appears to be an oxide or nitride phase was found at the grain boundaries and in the interior of the grains, especially in the vanadium-rich alloys. All alloys were annealed in sealed silica tubes containing 1 atm of pure argon and these tubes were then quenched in cold water. Although some manganese loss occurred during annealing, the loss seemed to be confined to the surface of the speci-

Jan 1, 1962

By H. T. Green, R. M. Brick

True stress-strain data on alloys of pure iron with up to 2.4 pct Al were obtained in the temperature range +100° to —185°C. Alumi-num was found to reduce yield and flow stresses of iron at low temperatures but to have little or no effect on ductility. The effects of temperature and composition on strain hardening are discussed. SEVERAL independent studies of the behavior of high purity iron binary alloys at low temperatures are now in progress in attempts to evaluate systematically the variables affecting the low temperature brittleness of ferritic steels. This paper reports the results of one such investigation in which the tensile properties of aluminum and aluminum plus silicon ferrites were measured from 100" to —192°C. True stress-natural strain data have been obtained in order to evaluate as many as possible of the parameters which describe the behavior of the materials involved. In comparable studies at the National Physical Laboratory in England, iron and iron alloys of high purity have been produced&apos; and tested at subat-mospheric temperatures.&apos; True stress-natural strain curves were obtained there also. The purest iron contained 0.0025 pct C and 0.001 pct O and N. Even this, as normalized at 950°C following hot rolling, showed little ductility at -196°C. The grain size was ASTM No. 3, and the room-temperature yield strength was 17,800 psi (which seems too high for pure iron). Some of the NPL irons contained considerably more oxygen and demonstrated intergran-ular fracture at —196°C. The authors2 carefully differentiated between intergranular fractures associated with excessive oxygen content and transcrys-talline cleavage with little ductility encountered at —196°C in the purer material. The cleavage stress was half again as great as that associated with inter-granular fracture. Test Material, Preparation, and Procedures Of a number of Fe-A1 alloys produced, eight were considered to be sufficiently pure for testing. Partial chemical analyses (Table I), low observed yield points, and high ductilities indicate these alloys to be comparatively pure for vacuum-melted irons of sizable ingots, 5 Ib or more. To produce the binary Fe-A1 alloys, electrolytic iron was melted in air, cast into slabs, and rolled to strips 0.010 in. thick. These strips, joined into a continuous ribbon and wound into 2 1/2 in. diameter spools, were subjected for four weeks to a moving atmosphere of purified dry hydrogen in a stainless-steel tube at 1050" to 1150°C. Charges of these spools were melted in beryllia crucibles under good vacuums (1 micron), and aluminum (99.97 pct Al) was added to the melts. Compositions of these alloys are recorded in Table I. The ingots were hot forged and then cold rolled at least 65 pct to 3/8 in. rods which were vacuum annealed to the desired grain size, approximately ASTM No. 4, prior to machining into tensile test bars. All tensile specimens had gage sections 1 in. long, with a fillet of 1.5 in. radius to the shoulder. Gage diameters were 0.250 in, except for a few rods where additional cold work required use of a 0.200 in. gage section. After machining, 0.002 in. was removed from the gage diameter using 240, 400, and 600-grit metallo-graphic papers. The final polish with 600 grit left the fine scratches running in the longitudinal direction. By this means, surface metal strained during machining was removed. A few specimens heat treated after machining were similarly reduced 0.004 in. to remove any material affected chemically by the atmosphere during heat treatments, as is discussed in a later section. Tensile tests of the eight alloys at constant temperatures from +100° to —185°C were performed in apparatus which has been described." The essentials include a double-walled insulated metal vessel which contained the liquid heat-transfer medium surrounding the test specimen. A constant temperature was maintained by means of a pyrometer which regulated the pressure of dry air driving liquid air through a copper coil. Temperature variation was less than ±2°C during a specific test. For axial straining, two lengths of case-hardened chain, terminating in simple shackles, loaded the specimen through threaded grips. The lower grip bar passed through a hole in the bottom of the test vessel to which it was joined by a thin-walled

Jan 1, 1955

By Herman Dykstra, R. L. Parsons

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.&apos; THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy&apos;s Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. &apos;Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If&apos; .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951

By R. L. Parsons, Herman Dykstra

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.&apos; THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy&apos;s Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. &apos;Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If&apos; .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951

By J. H. Swisher

The solubility of sulfur and rate of solution of sulfur in pure Lron were measured in H2S + H2 and H2S + H2 H2O gas mixtures. The solubility and diffusivity of sulfur at 1410°Care 0.13 pet S and 1.0 x 10-5 sq cm per sec, respectively. The solubility iS the same, but the rate of sulfurization is slower in the presence of H2O in the reacting gas. Under these conditions, the over-all rate is controlled jointly by a slow surface reaction and by solid-state diffusion; the mechanism for the surface reaction has not been identified. KNOWLEDGE of the behavior of sulfur in solid iron is desirable for the metallurgy of such products as free machining steel, where a high sulfur level is required, and inclusion-free high-strength steels, where the sulfur specifications are very low. The present investigation was undertaken to check previously reported values for sulfur solubility and diffusivity in 6 iron, and to study the poisoning effect of chemisorbed oxygen on sulfurization kinetics in H2-H2S-H2O gas mixtures. All of the experiments were performed at 1410°C. The thermodynamic behavior of sulfur in 6 iron was the subject of a paper by Rosenqvist and Dunicz.&apos; The sulfur solubility at 1400" and 1500°C was determined by equilibrating pure iron specimens with H2-H2S gas mixtures. The maximum solubility of sulfur in 6 iron was alsc determined by Barloga, Bock, and parlee2 by reacting iron wires with sulfur in sealed capsules. In another investigation, the diffusion coefficient of sulfur in 6 iron at temperatures up to 1450°C was measured by Seibel.3 The method used was to measure sulfur concentration profiles in diffusion couples containing radioactive sulfur EXPERIMENTAL Apparatus. A vertical resistance furnace wound with molybdenum wire and containing a recrystallized alumina reaction rube was used for the experiments. The hot zone in the furnace was approximately 2 in. long with a temperature variation of ±3oC. The hot zone temperature was automatically controlled to within ±2°C, and the test temperature was measured with a pt/Pt-10 pet Rh thermocouple before and after each experiment. Flow rates of the reacting gases were obtained using capillary flow meters. Materials. The source of H2S in the gas train was a premixed cylinder containing 5 pet H2S in H2. This mixture then was diluted with additional hydrogen and argon. In some experiments, water vapor was introduced by passing hydrogen and argon through a column containing 10 pet anhydrous oxalic acid and 90 pet oxalic acid dihydrate. The vapor pressure of water above this mixture is well-known.4 Argon was used as a diluent to minimize thermal segregation of H2S in the furnace5 and to reach higher H2O:H2 ratios than could be obtained in mixtures of H2 and H2S alone. Argon was purified by passage over copper chips at 350°C and subsequently over anhydrone. Hydrogen was purified by passage over platinized asbestos at 450°C and then over anhydrone. The H2-H2S mixture was purified by passage over platinized asbestos and then over P2O5. The specimen stock was made by melting and vacuum-carbon deoxidizing electrolytic "Plastiron" in a zirconia crucible. The principal impurities are listed in Table I. In some of the equilibrium experiments, six-pass zone-refined iron was used to minimize impurity side effects. This zone-refined iron had a total impurity level of about 25 ppm. Procedure. Specimens were annealed in hydrogen for a period of at least 2 hr at the beginning of each experiment. The specimens were held in the reacting gas for times varying between 10 min and 17 hr, and cooled to room temperature in a water-cooled stainless-steel block at the bottom of the furnace. The pH2S/pH2 ratios reported are those for gas equilibrium at 1410°C. Calculations based on available thermodynamic data8 showed that the only other gaseous8 species that formed in significant amounts in the furnace were S2 and S. Even when water vapor was introduced into the gas mixture, the concentrations of SO2, SO, and so forth, were negligible. The initial partial pressure of H2S was therefore corrected for its partial dissociation to S2 and S in determining the equi-

Jan 1, 1968

By G. R. Belton, R. J. Fruehan

The Knudsen cell-mass spectrometer combimtion has been used to study the Fe-Al and Ag-Al liquid alloys. By application of the recently developed integration technique to the measured ion-current ratios, activities have been derived for the Fe-A1 system at 1600° C and for the Ag-Al system at 1340"C. The results are partially represented by the following equations: Internal consistency between the data on silver-rich and iron-rich alloys is demonstrated by application of the literature measurements on the distribution of aluminum between the nearly immiscible liquids iron and silver. The usual restrictions on the ratio of the mean free path of the escaping atoms to the orifice diameter of the Knudsen cell are shown not to be limiting in this technique. DESPITE the importance of a knowledge of the activity of aluminum in understanding deoxidation equilibria in molten steel, no direct studies have been made of activities in liquid Fe-A1 alloys at steel-making temperatures. Lower-temperature direct studies have, however, been carried out on aluminum-rich liquid alloys by Gross, Levi, Dewing, and Eilson&apos; at 1300°C and by Coskun and Elliott&apos; at 1315°C. Apart from phase diagram calculations by Pehlke, other determinations have been indirect and were made by measurement of the distribution of aluminum between iron and silver475 and combination of these data with extrapolated activities in the Ag-A1 system.~-% ecently, however, Woolley and Elliott have made a significant contribution by directly measuring heats of solution in the Fe-A1 system at 1600°C. The present authorslo have recently employed a Knudsen cell-mass spectrometer technique in a study of activities in iron-based liquid alloys. In this technique activities and heats of solution are determined from a series of measurements of the ratio of ion currents of the components; and since ion-current ratios are used, problems caused by changes in instrument sensitivity or cell geometry are overcome. Results obtained for the Fe-Ni system were found to be in excellent agreement with previous work, thus demonstrating the reliability of the method. The present paper describes a similar study of activities in the liquid Fe-A1 and Ag-A1 systems, this latter system being included in order that a meaningful comparison can be made with the above-mentioned indirect studies. INTEGRATION EQUATIONS A detailed derivation of the equations used to determine the thermodynamic properties from the measured ion current ratios has been given elsewhere;&apos;&apos; however it is useful to summarize them here. By the combination of the Gibbs-Duhem equation with the direct proportionality between ion-current ratios and partial pressure ratios, it was shown that for a binary system at constant temperature and pressure: where al is the activity of component 1 with pure substance as the standard state, N, is the atom fraction of component 2 in the solution, and I; and t&apos;2 are ion currents of given isotopes of the components. The activity coefficient is given by: this latter equation being more suitable for graphical integration. Combination of Eq. [l] with the Gibbs-Helmholtz equation gives an expression for the partial molar heat of mixing: EXPERIMENTAL A Bendix Time-of-Flight mass spectrometer model 12! fitted with a 107 ion source and a M-105-G-6 electron multiplier, was used to analyze the vapor effusing from the Knudsen cell. The arrangement of the Knudsen cell assembly was essentially that of the commercial instrument (Bendix model 1030) but with several modifications. Instead of heating with a single tungsten filament, a cylindrical tantalum-mesh heater was employed. Up to 1400°C simple resistance heating was used but above this temperature electron bombardment between the tantalum mesh and the tantalum cell susceptor was necessary. The temperature was measured by means of a Leeds and Northrup disappear ing-filament type optical pyrometer sighted on an essentially black-body hole in the side of the cell. Details of the temperature control, temperature measurement, and in situ calibration of the optical pyrometer can be found elsewhere.I0 In the investigation of the Fe-A1 system the Knudsen cells were constructed of thoria crucibles with fitted thoria lids (Zircoa). The cells employed in investigating the Ag-A1 alloys were made up of high-purity alumina crucibles (Morganite) with lids of recrystal-lized alumina (Lucalox). The cells were 0.370 in.

Jan 1, 1970

By D. W. Rudd, D. W. Vose, S. Johnson

The permeability of Mo-0.5 pel Ti to hydrogen was investigated over a limited range of temperature and pressuire (709° to 1100°C, 1.i and 2.0 atm). The resulting permeability, p, is found to obey the The experimental data justifies the permeation mechanism as a diffusion contl-olled pnssage of Ilvdrogen atoms through the metal barrier. 1 HE permeability of metals to hydrogen has been investigated by a number of workers and their published results have been tabulated by Barrer&apos; up to 1951. Since most of the work on the permeability has been accomplished prior to this date, the compilation is fairly complete. Mathematical discussion of the permeability process has been reported by Barrer, smithells, and more recently by zener. From these efforts several facts are observed. First, the permeability of metals to diatomic gases involves the passage through the metal of individual atoms of the permeating gas. This is evidenced by the fact that the rate of permeation is directly proportional to the square root of the gas pressure. Second, the gas permeates the lattice of the metal and not along grain boundaries. It was shown by Smithells and Ransley that the rate of permeation through single-crystal iron was the same after the iron had been recrystallized into several smaller crystals. Third, it has been observed that the rate of permeation is inversely proportional to the thickness of the metal membrane. Johnson and Larose5 verified these phenomena by measurirlg the permeation of oxygen through silver foils of various thicknesses. Similar findings were noted by Lombard6 for the system H-Ni and by Lewkonja and Baukloh7 for H-Fe. Finally, it has been determined that for a gas to permeate a metal, activated adsorption of the gas on the metal must take place. Rare gases are not adsorbed by metals, and attempts to measure permeabilities of these gases have proved futile. ~~der&apos; found negative results on the permeability of iron to argon. Also, Baukloh and Kayser found nickel impervious to helium, neon, argon, and krypton. From what was stated above concerning the dependence of the rate on the reciprocal thickness of the metal barrier, it is seen that although adsorption is a very important process, at least in determining whether permeation will or will not ensue, it is not the rate determining process for the common metals. A case in which adsorption is of sufficient inlportance to cause abnormal behavior has been noted in the case of Inconel-hydrogen and various stainless steels.&apos;&apos; APPARATUS The apparatus used in this study is shown in Fig. 1. The membrane is a thin disc (A), but is an integral part of an entire membrane assembly. The entire unit is one piece, being machined from a solid ingot of metal stock. When finished, the membrane assembly is about 5 in. long. Two membrane assemblies were made; the dimensions of the membranes are given in Table I. The wall thickness is large compared to the thickness of the membrane, being on the average in the ratio of 13 to 1. There exists in this design the possibility that some gas may diffuse around the corner section of the membrane where it joins the walls of the membrane assembly, If such an effect is present, it is of a small order of magnitude, as evidenced by the agreement of the values of permeability between the two membranes under the same temperature and pressure. A thermocouple well (B) is drilled to the vicinity of the membrane. The entire membrane assembly is then encased in an Inconel jacket and mounted in a resistance furnace. The interior of the jacket is connected to an auxiliary vacuum pump and is always kept evacuated so that the membrane assembly will suffer no oxidation at the temperatures at which measurements are taken. The advantages of this configuration are: 1) there are no welds about the membrane itself, so that the chance of welding material diffusing into the membrane at elevated temperatures is remote. 2) It is possible to maintain the membrane at a constant temperature. Since the resulting permeation rate is very dependent upon temperature, it is advisable to be as free as possible from all temperature gradients. 3) It is possible to obtain reproducible results using different specimens. The only disadvantage to this configuration is the welds (at C) in the hot zone. The welding of molybdenum to the degree of per-

Jan 1, 1962

By R. G. Treuting, W. C. Ellis

The twinned state is characterized by a lattice of coincidence sites. Imperfections are required at stable lateral twin interfaces. Twinned regions can occur with relative ease in the diamond cubic IN recent contributions1,2 on the origin and growth of cubic annealing twins, attention has been directed to the orientation relations between such twinned components and their parent matrix. There are some aspects of twinning which may be illuminated by a more detailed consideration of the twinned state" alone. As an extreme example, the dense twinning in cast ingots of germanium,&apos; as contrasted with the rarity of twins in cast face-centered cubic metals, is yet to be accounted for. It has been this that has led us to the present work, which, it will be noted, uses methods and constructions in many respects similar to those of Kronberg and Wilson.&apos; In the cubic systems, a 70" 32&apos; rotation about a <110> axis is angularly equivalent, as to twinning, to the more usually considered 180" rotation about a <111> axis. Figs. 1 and 2 show a (110) projection of a twinned face-centered cubic lattice and a twinned diamond cubic lattice. In both figures, the two adjacent planes A and B, shown by the larger and smaller circles, are sufficient to represent the entire array. In each case a section of lattice, the original atom sites of which are shown by open circles, has been rotated as indicated through 70" 32&apos; to bring an original. [112] direction into coincidence with the [112] diiection. The latter is the intercept on the (110) projection of the (ill) plane normal thereto, the twinning plane. In the face-centered cubic case the rotation can be performed about an axis passing through an atom-site; the mirror plane then is also a composition plane containing atoms common to both twinned and untwinned lattices. The diamond cubic lattice may be construed as two interpenetrated face-centered lattices. Its (111) planes recur in a sequence of alternately short and long interspacings. Consequently a mirror plane for twinning cannot be a composition plane, but must be the bisector of one of the spacings. When the longer spacing is selected, the closest distance of approach across the mirror plane in the [ill] direc- tion is identical with that in the untwinned structure. In each case periodically recurring (ill) planes (parallel with the twinning plane) are found, on which there is coincidence of atom sites of the pre-twinned and twinned orientations; these are indicated by the cross-hatched circles. In the face-centered lattice there is such coincidence every third (ill) plane; in the diamond cubic lattice, on two adjacent planes in every six. At the twinning interface in the latter, there is on each side of the mirror plane a (ill) plane of atoms common to both twin components. Conceivably, there is little influence on a plane of atoms about to be adhered to such a pair of coincidence planes, whether it be laid down in a normal or in a twinned position with respect to the previously formed structure. Slawson% as attributed the high incidence of twinning in diamond to this boundary state. Further examination shows that the motion of intermediate planes can consist of various pairs of equal and opposite translations, for example of (ill) planes in the [l&apos;;i2) direction, the familiar twinning shear, indicated in the small schematics in the figures. Since the translations form a system of shears of alternating sign between coincidence planes, twinning could take place by such a mechanism over an extended region without extensive shear; in fact, in this case any atom moves but the distance in the [1i2] direction. One alternative construction for the face-centered cubic lattice leading to the same end result is illustrated in Fig. 3. The plane (711) with respect to the pretwinning orientation (the twinning plane of Fig. 1) is given, the twinned region arbitrarily bounded by <110> and <112> directions. The coupled shear is identical to that of Fig. 1. The "rotational" movement about coincidence sites generating the same twinned position could consist as shown of the translation a,/d% for each atom of a group of three in the B layer in a different one of the three <112> directions, and a similar translation of the underlying three atoms in the C layer in either the same or the opposite sense. This is not dissimilar to Kronberg and Wilson&apos;s construction for their 22" rotation of three adjacent (111) planes.

Jan 1, 1952

By George B. Clark

Many basic engineering principles of all four phases of mining operations, namely, prospecting, exploration, development, and exploitation, can be analyzed better in terms of quantitative geology. Geological data from both field and laboratory will also complement scientific methods now being developed. THE geological data emphasized so successfully in prospecting for new deposits, that is, structural controls, strength of solutions, and type of mineralization, are basically those required for successful exploitation. In the mining of newly discovered deposits the most economical methods should be employed as early as possible to keep the overall cost per unit produced at a minimum and to permit maximum extraction of valuable minerals. A crucial question is: How can geological data be translated into useful quantitative results which will aid in achieving this end? H. E. McKinistry&apos; has suggested that a solution may be reached in one of two ways: 1—the usual approach, use of judgment based on experience; or 2—mathematical calculations and tests on models, both subject to certain limitations. He also suggests that in addition to better use of geology more case data and theoretical data are needed on which to base sound judgment. Further research, therefore, is necessary. Perhaps in this field the emphasis should be on more specialization in mining methods and ground movement by men with thorough training in physics, engineering, geology, and underground mining. These specialists would be equipped to point out the most economical and scientific methods of exploitation. Selection of a stoping method is governed by the amount and type of support a deposit will require in the process of being mined, or by the possibility of employing the structure of the deposit to advantage in mining the ore by a caving method. In addition to these factors there are others which almost invariably influence the choice of an economical method of mining:&apos; 1—strength of ore and wall rocks; 2—shape, horizontal area, volume, and regularity of the boundaries of the orebody, and thickness, dip and/or pitch of the deposit and individual ore shoots; 3—grade, distribution of minerals, and continuity of the ore within the boundaries of the deposit; 4—depth below surface and nature of the capping or overburden: and 5—position of the de- posit relative to surface improvements, drainage, and other mine openings. In the final analysis it is usually necessary to disregard the less important of these factors to satisfy the requirements of the more important. Because of the variation of geological conditions throughout and surrounding the deposit, no mining method will be everywhere ideally applicable to the conditions encountered in one ore deposit. The immediate problem is to interpret the above physical characteristics of deposits in terms of geological characteristics. Very few quantitative geological data are available on the factors related to a choice of mining methods. However, there are many descriptive data in mining and geological literature which collectively show how important an effect details of geology have upon all phases of mining operations. The following categories of basic mining methods were investigated to establish the geological factors that have affected their successful application: 1— open stopes with pillars; 2—sublevel stoping; 3— shrinkage stoping; 4—cut-and-fill stoping; 5— square-set mining; 6—top slicing and sublevel caving; and 7—block caving. It should be noted that the first five of these methods are listed in the order of increasing support requirements. Mines were selected as examples only where geological descriptions were complete enough to warrant their use. A study of the geological factors involved in mining operations led to a choice of the following classifications, employed in Table I: 1—structural type of orebody; 2—dimensions (geometry); 3— country rock (type); 4—faulting, folding, and fracturing; 5—alteration of ore and rock; 6—type of mineralization; and 7—geological factors determining mining method (summary). Of these factors only one yielded results that can be defined from available data in a quantitative manner, i.e., dimensions of the deposit. These are the most reliable guides that can be used in selection of suitable mining methods. They are, in general, the properties of geologic structure most difficult to evaluate by studies of models, pho-toelastic studies, and other laboratory methods, all of which are at present more limited in their applications than the geologic method. Application of geology has proved a reliable guide in other phases

Jan 1, 1955

By G. W. Thomas

A mathematical tnodel of forward combustion in an oil reservoir is treated in this paper. The model describes a radial system having a vertical section of essentially infinite thickness, all of which is permeable to gas flow. Combustion, however, is presumed initiated over a limited thickness of the total vertical section. In the interval supporting cotnbustion, the mechanisms of radial conduction, convection and heat generation are taken into account. Above and below the burning interval, heat transport in the radial direction is by cottduction and convection. Vertical heat losses from the ignited interval are accounted for by conduction alone. A general solution is presented for the temperature distribution caused by radial movement of the combustion front. The results show that no feedback of heat occurs into the ignited interval when convection and conduction are acting in the bounding media. Peak temperatures are also 5 to 10 per cent higher than in the case where heat transport in the bounding media is by conduction alone. We arbitrarily define vertical coverage to be that fraction of the total ignited interval which is at 600F above atnbient, or greater, at any given time. The radial distance at which the vertical coverage becomes zero is the propagation range of the combustion front. It was found that an increase in vertical coverage results when the oxygen concentration, fuel concentration or gas-injection rate is increased. Moreover, the combustion front can be propagated 10 to 15 per cent further than in the case where only conduction is acting above and below the ignited interval. INTRODUCTION In the theoretical treatment of forward combustion in a radial system, one of the problems encountered is the determination of the transient temperature distributions caused by an expanding cylindrical heat source. Bailey and Larkin&apos; and Ramey&apos; simultaneously presented analytical solutions to the problem assuming heat transport by conduction alone. In a subsequent publication, Bailey and Larkin3 included the effects of both conduction and convection while treating linear and radial models. In this latter work, however, vertical heat losses were largely neglected. Selig and Couch&apos; dealt with a radial model in which both conduction and convection were acting. Only a limiting case involving vertical heat losses was considered, however. Namely, temperatures on the boundary of the bed of interest were set equal to zero. Solutions thus obtained were representative of a system having a maximum vertical heat flux. Chu5 recently treated a more general case in which a permeable bed was considered bounded by impermeable media. Conduction and convection took place within the bed, and only conduction outside of the bed. The effects of vertical heat losses were included in his study. Solutions were obtained by numerical techniques. This paper is an extension of the theoretical work of other authors pertaining to forward combustion in a radial system. In particular, a mathematical model of the process is treated in which heat generation occurs over a small vertical interval of a larger permeable section. In the interval supporting heat generation, and above and below this interval, the mechanisms of radial conduction and convection are also presumed acting. Heat losses from the ignited interval are accounted for by vertical conduction. An analytical solution for the temperature distribution caused by radial movement of the burning front is presented. The effects of certain process variables are indicated and comparisons with Chu&apos;s results are made. THEORY To render the mechanism of forward combustion tractable to mathematical treatment, we idealize the problem to the extent of assuming continuous reservoir media possessing homogeneous and isotropic properties. The following additional assumptions are implicit in this analysis. 1. The thermal parameters, i.e., heat capacities, thermal conductivities and thermal diffusivities are invariant with temperature and pressure. Moreover, the bounding media possess the same thermal properties as the bed of interest. 2. The temperatures of the porous media and its contained fluids at any point and at any time are equal. 3. The reaction rate between the oxidant gas and the fuel is infinite. This assumption implies that the incoming oxygen concentration instantaneously goes to zero within an infinitesimal distance, i.e., the width of the combustion zone is negligible. 4. The rate of gas injection is constant and corresponds to the average rate throughout the lifetime of the project. 5. The fuel concentration is constant throughout the volume of rock swept out by the burning zone. 6. There is complete burnoff of fuel. This assumption demands that the rate of propagation of the burning front equals the rate of fuel burnoff. In a radial system, with a

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 J. W. Spretnak, G. W. Powell, J. H. Bucher

Tlze room-temperature tensile fracture oj smooth, round specitnens of three ultrnhigh- strength steels tempered to a wide range of strength levels was studied by means by light and electron-microscopic examination of the fracture surfaces. The fracture of AISI 4340 and 300 M at all the strength levels studied, and H-11, except after tempering at 1200° and 1300°F, occurs in three stages. The initiation of fracture is internal (except in some lightly tcmpeved specimers in which fracture is initiated at surface flaws), and is nucleated largely by separation at metal-second phase intevjaces. TIze voids grow and, coalesce to form a crack. When the crack has reached a sufficienl size, rapid propngutio~z ensues. Failure in this stage of fracture usually occurs by dimpled rupture of inicroshear stefis. In the case of H-11 tempered in the 1125° to 1300°F range, fracture in the shear steps is predominantly by concentrated deformation without void formation. The termination of fracture is usually occomplished by the formation of a shear lib in which fracture occurs by shear dimpled rupture. In the case of H-11 tempered at 1200° and 1300°F, no shear lip was obserued, and the radial elelments extend to the surface—a true termination slage does not exist. ThE tensile fracture of several metals and alloys has been investigated.2-4 In the case of polycrystal-line materials, cup-cone fracture usually results. The mechanism of cup-cone fracture may be summarized as follows.5 Cavities are formed in the necked region of the specimen. They usually are initiated by inclusions or second-phase particles. The cavities extend outwards by means of internal necking, and a crack lying about perpendicular to the length of the specimen is formed in the necked region. Subsequent crack growth occurs by the spread of bands of concentrated plastic deformation inclined at an angle of 30 to 40 deg to the tensile axis. Cavities are formed in the bands of concentrated deformation. The deformation bands zigzag across the bar with the net result that mac-roscopically the crack extends about perpendicular to the specimen axis. The final separation, or cone formation, appears to occur by continued crack propagation along one of the deformation bands out to the surface of the specimen. The micromechanics of the tensile fracture of ultrahigh-strength steels have not been thoroughly investigated. Larson and carr6,7 studied the tensile-fracture surfaces of AISI 4340 with a low-power microscope and reported that three stages of fracture could be observed in general. A centrally located region characterized by circumferential ridges, an annular region characterized by radial surface striations, and a peripheral shear lip were found. It was first pointed out by 1rwin8 that the central region is very probably one of fracture initiation and slow growth, and that the annular, radially striated region is one of rapid crack growth. Presumably the crack grows slowly, assuming roughly a lenticular shape, until it is large enough for the initiation of rapid propagation. In this investigation, it was attempted to determine the fine-scale aspects of the room-temperature tensile fracture of some ultrahigh-strength steels, and to relate the variation in fracture mode with microstructure. The steels studied were AISI 4340, 300M, and H-11 tempered to a wide range of strength levels. I) EXPERIMENTAL PROCEDURE The compositions of the steels studied are given in Table I. The steel was received in the form of hot-rolled bar stock 5/8 to 1 in. in diameter from which oversized specimens were machined and heat-treated. The heat treatments employed are given in Table 11. Subsequent to heat treatment, the specimens were ground to the final dimensions and stress-relieved by heating for 1 hr at 350°F (with the exception of the as-quenched steel). Standard smooth round specimens of 0.252-in. diameter and 1-in. gage length were tested in a Tinius Olsen Universal Testing Machine using a cross-head speed of 0.025 in. per min. The relatively coarse aspects of the fracture topography were determined by light-microscopic examination of sections through the fracture surface of nickel-plated specimens. A direct carbon-replication technique9 was used in the electron-microscopic study of the fracture surfaces. The replicas were examined in the electron microscope, and stereo pairs of electron micrographs were taken. The stereo pairs were then examined using a Wild ST4 Mirror Stereoscope. Carbide and inclusion particles extracted in the replicas were analyzed by selected-area electron diffraction. II) EXPERIMENTAL RESULTS The mechanical testing data are summarized in Table 111. The values reported are the average of

Jan 1, 1965

By W. R. Purcell

An apparatus is described whereby capillary pressure curves for porous media may be determined by a technique that involves forcing mercury under pressure into the evacuated pores of solids. The data so obtained are compared with capillary pressure curves determined by the porous diaphragm method, and the advantages of the mercury injection method are stated. Based upon a simplified working hypothesis, an equation is derived to show the relationship of the permeability of a porous medium to its porosity and capillary pressure curve, and experimental data are presented to support its validity. A procedure is outlined whereby an estimate of the permeability of drill cuttings may be made with sufficient acuracy to meet most engineering requirements. INTRODUCTION The nature of capillary pressures and the role they play in reservoir behavior have been lucidly discussed by Lev-rett&apos;, Hassler, Brunner, and Deah12, and others. As a result of these publications the value of determining capillary pressure curves for cores has come to be generally recognized within the oil industry. While considerable attention has been directed toward the subject in an effort to provide a reliable method of estimating percentages of connate water, it has been recognized that capillary pressure data may prove of value in other equally important applications. This paper describes a method and procedure for determining capillary pressure curves for porous media wherein mercury is forced under pressure into the evacuated pores of the solids. The pressure-volume relationships ob- tained are reasonably similar to capillary pressure curves determined by the generally accepted porous diaphragm method. The advantages of the method lie in the rapidity with which the experimental data can be obtained and in the fact that small, irregularly shaped samples, e.g., drill cuttings, can be handled in the same manner as larger pieces of regular shape such as cores or permeability plugs. Based upon a simplified working hypothesis, a theoretical equation will be derived which relates the capillary pressure curve to the porosity and permeability of a porous solid, and experimental data will be presented to support its validity. This relationship aplied to capillary pressure data obtained for drill cuttings by the procedure described provides a means for predicting the permeability of drill cuttings. METHODS FOR DETERMINING CAPILLARY PRESSURES Several techniques have so far been employed in determining capillary pressure curves and these fall into two principal categories: (1) Liquid is removed from, or imbibed by, the core through the medium of a high displacement pressure porous diaphragm (2) Liquid is removed from the core which is subjected to high centrifugal forces in a centrifuge4,&apos;. There are? however, certain limitations inherent in both methods. The greatest capillary pressure which can be observed by method (I), above, is determined by the maximum displacement pressure procurable in a permeable diaphragm which at the present time appears to be less than 100 psi. An even more serious limitation of the diaphragm method is imposed hy the fact that several days may be required to reach saturation equilibrium at a given pressure; hence, the time re- quired to obtain a well-defined curve may be measured in terms of weeks. Furthermore, to date, no suitable technique for handling relatively small, irregularly shaped pieces of rock, such as drill cuttings, has been reported and, therefore, measurements must be made, in general, on cores, or portions thereof. The centrifuge method offers the distinct advantage over the porous diaphragm method of arriving at saturation equilibrium in a relatively short time by virtue of the elimination of the transfer medium for the liquid. The calculation of capillary pressures from centrifuge speeds is somewhat tediousa, however, and the equipment required is fairly elaborate. While there exists the possibility that this method might be adaptable to the determination of the capillary pressures of cuttings, this particular ramification has not been investigated, as far as is known. In view of the limitations of the two principal methods for determining capillary pressures, the apparatus described in the following sections has been devised in order that difficulties previously encountered might be circumvented. MERCURY INJECTION METHOD FOR DETERMINING CAPILLARY PRESSURES Theory The methods described above for determining capillary pressures are characterized by the fact that one of the fluids present within the pore spaces of the solid is a liquid which "wets" the solid, i.e., the contact angle which the liquid forms against the solid is less than 90" as measured through that phase. For these "wetting" liquids the action of surface forces is such that the fluid spontaneously fills the voids within the solid. These forces likewise oppose the withdrawal of the fluid from the pores of the solid.

Jan 1, 1949