Search Documents

By Peter Soo

Pvestrained Ferrovac E iron has been neutron-irradiated at approximately 90°C to an integrated flux of 1020 nut (E > 0.82 mev]. The irradiation was found to produce an incveased temperature dependence of the flow stress in addition to a greatly increased athemal stress. Measurements of the flow stress and stress relaxation, from which the activation volume and activation energy for slip were deduced, show that neutron irradiation changes the rate -controlling slip process to one based on dislocation interactions with tetragonal distortions which are Produced around submicroscopic interstitial loops in the lattice. The study indicates that without prestraining prior to irradiation the chances of detecting a change in the rate -controlling slip process are greatly reduced because in the initial stages of slip a substantial fraction of the radiation defects are swept out of the slip plane by gliding dislocations. Thus, activation parameters which are subsequently measured are representative of a greatly reduced defect density and would not differ appreciably from those for unirradi-ated material. The large increase in the athermal component of the flow stress is probably connected with the presence of depleted zones in the lattice which are introduced by irradiation. ALTHOUGH fast neutron-irradiation has not been observed to markedly alter the activation parameters for slip in bcc metals,' small but significant changes do occur. Most experimenters agree that irradiation predominantly increases the athermal component of the yield stress.'-= In addition to this, Laidler and smidt7 have shown that in iron irradiated to 5 X 10" nvt and molybdenum irradiated to 10" nvt, changes occur in the activation volumes for slip. A similar conclusion has been reached by Milasin and Malkin8 for irradiated iron. Work by Ohr et a1.5 shows that for Ferrovac E iron, irradiated to 1.2 X 1016 nvt, small increases in the activation energy for slip also occur. So far these changes in the activation parameters have not been explained on a firm theoretical basis. One important factor which would minimize the chances of detecting a change in the slip mechanism upon irradiation is the presence of "channeling" which has been observed in molybdenum,9 niobium,10 and iron.11These channels are formed by gliding dislocations which sweep irradiation defects out of the active slip planes and thereby create zones in which continued dislocation motion is encouraged. The activation parameters for the dislocations gliding in the defect-free channels would, therefore, be similar to those for unirradiated iron and a change in the rate-controlling slip process would be difficult to detect. In the present work, an attempt has been made to reduce the effect of uneven deformation on the measured activation parameters for slip in neutron-irradiated Ferrovac E iron polycrystals, so that a more realistic assessment of the effects of neutron-irradiation could be made. Primarily, the experiments involve the irradiation of specimens which had been prestrained to 9 pct elongation at room temperature prior to insertion into the reactor. It was hoped that the introduction of a large number of evenly distributed dislocations would substantially decrease any channeling effect which might otherwise occur. MATERIAL AND EXPERIMENTAL PROCEDURE The starting material was vacuum-melted Ferrovac E iron, an analysis of which is given in Table I. The standard tensile specimen had a gage length of 1.125 in., a cross-sectional diameter of 0.120 in., and a re-crystallized grain size of 1.2 x 10-3 in. All tensile tests were conducted on a floor model "Instron" tensile machine at a strain rate of 3 x 10-4 per sec. The irradiation of the prestrained specimens was performed in the Brookhaven High Flux Beam Reactor to an integrated flux of 1020 nvt (E > 0.82 mev) at a temperature of about 90°C. All specimens were excap-sulated in high-purity aluminum sheaths which were lightly swaged around the samples to ensure good thermal contact. Subsequent measurements on the irradiated specimens showed that within experimental accuracy the swaging had not deformed them. EXPERIMENTAL RESULTS Fig. 1 shows the flow stresses for a series of unirradiated control samples. In order to produce a comparable dislocation substructure throughout the test sm range, all specimens were prestrained

Jan 1, 1970

By J. M. Tobin, L. H. Cadoff, L. M. Adelsberg

The reaction between liquid zirconium and graphite at temperatures above 2000 °C has been investigated. The reaction products were found to be carbon-saturated zirconium metal and ZrC which formed between the graphite and the metal. Parabolic growth behavior was observed for the ZrC Phase at all temperatures of this investigation. The parabolic growth constant at temperatures between 2000° and 2860°C was measured to be 1.83 exp 84,300/RT) sq cm per sec. The reaclion mechanism was proposed to be the rapid carbon saturation of the liquid metal and the formation of ZrC at the metal-carbide interface with diffusion of carbon through the ZrC, the "rate - determining step" of the reaction. The concentration-independent diffusion coefficient of carbon in ZrC, (DZrC), was expressed as 0.95 exp? 78,700/RT) sq cm per sec. This value mas calculated using the temperature-invariant ZrC phase fields proposed in the literature. The [ZrQiq)] — [ZYQiq) + ZrC] phase boundary over the temperature range 2000° to 2800°C was determined and the ZrC + C eutectic temperature was found to be 2890° ± 50°C. ThE Group IV and VB transition-metal refractory carbides are of interest because of their high melting points, high temperature strength properties, and relative inertness in certain corrosive environments. The present-day understanding of these materials, however, is limited by the general unavailability of accurate and reliable physical and chemical property data. This is due primarily to the difficulties associated with the preparation of suitable, high-density, high-purity carbide samples, and the achievement and control of uniform high-temperature environments. Accurate measurement of temperature is also an important factor limiting the reliability of reported data. The Zr-C reaction was selected for investigation because of the high melting point (>34000C) and favorable nuclear properties of the reaction product, ZrC. The direct reaction method afforded an opportunity to obtain kinetic data on the fully dense carbide. In this paper, layer-growth techniques were used to estimate the diffusivity of carbon in ZrC and to investigate the phase equilibria in the Zr-C system at temperatures above 2000°C. EXPERIMENTAL PROCEDURES AND DATA Crystal bar zirconium (99.9 pct Zr), purchased from the Nuclear Materials and Engineering Corp., and ZTA graphite (99.9+ pct) were used in this study. The analyses of the materials are listed in Table I. A schematic of the high-temperature carburization apparatus is shown in Fig. 1. The sample was a zirconium metal charge in a 1/2 -in.-ID graphite crucible which was capped with a tightly fitting graphite stern. The crucible and stem were designed to closely approximate black-body conditions. The graphite crucible was packed in lampblack (outgassed 1 hr at 2100°C) which provided insulation and thermal stability to the system. The inert atmosphere was maintained by a constant flow of high-purity argon or helium gas. The nozzle-diffuser section on the top flange was sufficient to prevent back-diffusion of air into the system. Chemical analysis of the carburized samples revealed oxygen and nitrogen concentrations of less than 44 and 20 ppm, respectively. The sample was heated by induction with a Westing-house 5 kw-450 kc power source. Temperature was measured with a Milletron two-color pyrometer which was sighted into the crucible by reflection from an overhead front surface mirror. Optical losses due to glass absorption and light reflection at air-glass interfaces were minimized by the employment of the dynamic-gas seal. It was found that no correction was required for the front surface reflector. The pyrometer was calibrated periodically against a U.S. Bureau of Standards tungsten ribbon secondary standard. Temperatures inside the crucible were also calibrated against the Zr-ZrC eutectic (1850°C)1 and the Nb-Nb2C eutectic (2335°C)2,3 temperatures and agreed to within ±10°C of these values. The temperature was controlled by manually adjusting the power input; variations of

Jan 1, 1967

By Y. C. Liu, G. A. Alers

The effect of rolling reduction on the textures of Cu-A1 alloys has been investigated both by pole figure and by modulus methods. In alloys which exhibit complete copper or brass types of rolling texture, the rolling reduction has little effect on the texture except to increase the degree of preferred orientation. In alloys which exhibit a transition texture, however, increased rolling reduction increases the amount of brass-type texture at the expense of the copper-type texture. The present experimental results show that there is no one-to-one correspondence between the SFE and the rolling texture of fcc metals. Additional data taken from the literature for fcc metals also support this conclusion. On the other hand, the present and previous experimental results are shown to be in good agreement with the suggestion that the texture transition occurs at a critical value for the separation distance between two partial dislocations—a consequence of the "dislocation interaction" hypothesis for texture. formation. This critical separation occurs when the parameter .r/ub is 3.75 x 10'3. From this, a value for the SFE of 39 ergs per sq cm may be deduced for a Cu-2.85 at. pct A1 alloy. ThE correlation between the rolling texture of fcc metals and the stacking fault energy, SFE, was one of the first attempts to relate atomistic properties with the type of rolling texture.' This correlation gives a qualitative explanation for the experimental observation that the addition of alloying elements, which generally lower the SFE, changes the copper-type texture to a brass-type texture. The simplicity of this correlation had led to its general acceptance and even its quantitative use.' However, it is only a correlation and is largely based on descriptive features of pole figures, and on the poorly known SFE values in dilute alloys. Quantitative verification of this phenomenologi-cal correlation is, in fact, completely lacking. One purpose of the present study is to test this correlation. Another atomistic description for the formation of rolling texture is the "dislocation interaction" hypothesis of texture formation.3 In this hypothesis, the factor controlling the type of rolling texture depends on whether or not the separation distance between two partial dislocations exceeds a critical value. Materials having a separation of less than the critical value are supposed to exhibit a copper-type texture while those with a separation above the critical value are supposed to have a brass-type texture. At the critical value, it is expected that the material should show equal amounts of copper- arid brass-type orientations in their textures, i.e., a 50 pct transition texture. The SFE appears in this hypothesis as only one of several factors which determine the separation distance between partial dislocations. It is possible to test the validity of these two concepts by studying the rolling texture as a function of rolling reduction. Since the SFE per se is an intrinsic property of the metal, it should not, by definition, be influenced by local irregularities, such as variable stress conditions. Thus, no change in texture-type is expected to occur with changes in rolling reduction. On the other hand, according to the "dislocation interaction" hypothesis, any factor that effectively influences the separation distance of partial dislocations would be expected to change the rolling texture. Since the separation distance between partial dislocations is known to depend upon local stresses,4-6 it is anticipated that there would be an effect of the degree of reduction on the texture-type. Also, since applied stresses are more likely to increase, rather than to decrease, the separation between partials,4'5 the overall effect would be to increase the amount of material in the brass-type orientations as rolling reduction is increased. Furthermore, this reduction dependence would be most prominent in alloys exhibiting the transition texture since the distance between partials in those alloys is thought to be close to the critical value. Experimental data in the literature is insufficient to distinguish between these two alternatives. Haessner studied the effect of rolling reduction on textures in a series of Ni-Co alloys by means of the X-ray intensity-ratio technique,' and found that while one texture parameter indicated no reduction dependence the other indicated a slight dependence of the rolling texture on reduction in the range of 96 to 99 pct. As has been noticed previously, the intensity-ratio technique is a convenient but controversial method7 because there is no a priori reason to suggest which intensity-ratio would describe the texture most meaningfully. A more quantitative method of describing textures is found in terms of the orientation dependence of Young's modulus. Here, the type of modulus aniso-tropy associated with the copper-type texture is sufficiently different from that observed for the brass-type texture to allow the two types to be easily distinguishable and a quantitative measure of the amount of each can be deduced from the numerical results. This ability to provide quantitative data is particularly valuable when the two textures occur simultaneously in one alloy as is the case for the transition textures. In this paper the modulus method, supplemented by pole figure data, is used to look for an effect of roll: ing reduction the texture. Also by combining the texture measurements with recent determinations of the SFE in Cu-A1 alloys'0'" it should be possible to test for a relationship between the SFE and textures.

Jan 1, 1970

By W. A. Backofen, R. L. Fleischer

Single crystal, bicrystal, and polycrystal tensile tests of aluminum at 4.2°K, 77°K, and 300°K have been used to examine the role of grain boundaries in the deformation process. Results indicate that a grain boundary may affect the extent and slope of easy glide. The stage II hardening rate, on the other hand, is independent of the presence or absence of grain boundaries. This conclusion allows the size of the region of multiple slip caused by an incompatible grain boundary to be determined. For the size of bicyystal sample used in this study, multiple slip occurs in about half of the cross section. PREVIOUS studies of the stress-strain characteristics of bicrystals of face-centered-cubic metals have been limited to aluminuml-5 at room temperature. Recent results, however, indicate that the stress-strain curves of single crystals of such metals may be separated into at least three stages6 in which different deformation processes are occurring7 provided testing is done at sufficiently low temperatures.' Since for aluminum a well-defined stage II develops only below room temperature, previous studies have not been able to relate effects of grain boundaries to all of the three stages of deformation. It is therefore to be expected that low-temperature deformation of aluminum single crystals and bicrystals should clarify the effects of grain boundaries on the different processes of deformation. EXPERIMENTAL PROCEDURE Single crystals and bicrystals were grown from the melt by the standard techniqueg with aluminum reported by Alcoa to be 99.993 pct pure. Ridges in the boat were used to guide the grain boundary during growth, assuring that the boundary would bisect the sample.10 The rate of furnace motion during growth was 1.0 cm per hr. During growth zone purification resulted, as evidenced by the ability of the first material to freeze to recrystallize at room temperature following severe deformation. Samples were approximately 4.4 X 6.6 mm in cross section and 103.5 mm in length between grips. Samples were annealed at 635" i 5°C for 40 hr and furnace cooled over a 7-hr period. They were then electropolished in a solution of 5 parts methanol to 1 part perchloric acid at a current density of 15 amp per sq dm for about 30 min at temperatures below 0°C. Tensile testing was performed at 295" (room temperature), 77" (sample in liquid nitrogen), and 4.2"K (sample in liquid helium) on the hard-type machine indicated schematically in Fig. 1. The machine con- sists basically of a tube surrounding a rod; one end of the sample is attached to each member, and the rod is pulled up the tube to extend the sample. The rod is rigidly mounted and is moved vertically by a system described by asinski." The pulling force is measured continuously by an electrical strain gage load cell, and the relative displacement of the tube and rod is also recorded continuously by a soft cantilever beam with electrical strain gages. Maximum stress and strain sensitivities were ±2g per sq mm and * 3-10-5. In all tests the strain rate was approximately 5.10-5 per sec. The thin wires in the tensile apparatus introduce softness, which may be corrected for, however, by measuring load vs displacement with the sample replaced by an elastic member. For loads greater than 15 kg the spring constant is 1.875.106 g per cm. The flexible wires also served to reduce substantially the large shearing forces which may arise in the case of grips having horizontal rigidity.'' As in any gripping system, however, bending moments will arise in the course of deformation by single slip. Engineering stress, s = (load)/(original cross-sectional area), and strain, E = (increase in length)/ (original length), are used for stress-strain curves unless otherwise indicated. Tables list resolved shear stress, T=mo and shear strain ? = dm, where m is the usual Schmid resolved shear stress factor for the primary slip system at the start of deformation. The first group of samples to be described forms an isoaxial set, all of the crystals making up the single crystals or bicrystals having the same tensile axis, the orientation of which is indicated by the cross in Fig. 2. For this orientation the primary slip plane and slip direction make angles of 45 deg with the tensile axis and the Schmid factor m has its maximum possible value of 0.5. Rotations about the tensile axis are indexed by means of an angle 0 between the small-area surface of the samples and the projection of the primary slip direction onto the cross section, as defined in Fig. 3. In single crystals, values of 0 were 0 and 90 deg, while in bicrystals 0 values were (0 deg, 180 deg), (90 deg, 270 deg), and (0 deg, 90 deg) as indicated in Fig. 4.

Jan 1, 1961

By Fred H. Poettmann

The vaporization characteristics of carbon dioxide in a League City natural gas - Billings crude oil system were studied at three temperatures, 38°. 120°, and 202°F and for pressures ranging from 600 to 8,500 psi. Variation of carbon dioxide concentration up to 12 mole per cent in the composite showed no effect on the equilibrium vaporization ratios (K values) of the hydrocarbon constituents or on the K value of carbon dioxide itself. It was shown that carbon dioxide is more soluble in crudes than in distillates which is contrary to the behavior of methane. A working chart of carbon dioxide K values is presented. INTRODUCTION The study of the equilibrium vaporization ratios of mixtures of paraffin hydrocarbons has been rather thorough.2,6,7,8,9 In the past few years considerable attention has been paid to the vaporization characteristics of the so-called noncondensable gases such as nitrogen, carbon dioxide, and hydrogen sulfide in mixtures of hydrocarbons. since they usually occur to some extent in most crude oils and natural gases.1,3,4,5 Knowledge of this behavior is useful to both the production and refining phases of the petroleum industry. This paper reports the equilibrium vaporization ratios (K's) of carbon dioxide in a mixture of League City natural gas and Billings crude oil, and compares them to those obtained in a natural gas-distillate system. The equilibrium vaporization ratios for the hydrocarbon components in this system had previously been studied by Roland.' In addition to the determination of the K values for carbon dioxide, the K values for methane and ethane were also determined in order to observe what effect, if any, the presence of carbon dioxide had on these K values. The concentration of carbon dioxide was also varied in order to observe the effect of this variable on the carbon dioxide K values. EXPERIMENTAL PROCEDURE The apparatus used in this study cotlsisted of a stainless steel equilibrium cell of about 2 liters capacity. The cell was mounted on trunions permitting rocking in a thermostatically controlled oil bath. Two high pressure valves fitted with steel tubing were mounted on the top of the cell. one was used for sampling the equilibrium gas phase and the other for sampling the equilibrium liquid phase by means of an induction tube within the cell. Stainless steel tubing from the bottom of the cell led to a mercury reservoir and manifold which was connected to a free-piston type pressure gauge manufac- lured by the American Instrunlent Ctr. and to a volumetric. putrip. The temperature of the oil bath was measured by means of a ralibrated mercury-in-glass thermometer. The recorded temperatures are believed to be accurate to ±0.5 °F. The pressures are correct to 22 psi. The crude oil used in this study was stock tank oil obtained from the Wilcox formation in the Billings Field, Noble County. Okla. The natural gas was obtained from the League City Field. Galveston County, Tex. The oil was treated with anhydrous calcium sulfate in order to remove the last traces of water. To insure a supply of constant composition gas at room temperature the cylinders of League City gas were cooled to about 30°F, inverted, and the condensed liquid was allowed to drain from the cylinders. The analysis of the gas and crude are tabulated in Table I. The carbon dioxide came from Pure Carbonic, Inc., and was .stated to have a purity of 99.5 per cent or better. The procedure used to obtain samples of the equilibrium liquid and vapor was similar to that employed by others making use of the rocking type equilibrium cell.6,7,8 The equilibrium cell was evacuated and calculated quantities of carbon dioxide, natural gas, and crude oil were charged to the cell to the desired pressure. In charging the equilibrium cell an attempt was made to maintain the ratio of the natural gas to crude oil as close as possible to that employed by Roland. After the cell was charged, samples of

Jan 1, 1951

By Fred H. Poettmann

The vaporization characteristics of carbon dioxide in a League City natural gas - Billings crude oil system were studied at three temperatures, 38°. 120°, and 202°F and for pressures ranging from 600 to 8,500 psi. Variation of carbon dioxide concentration up to 12 mole per cent in the composite showed no effect on the equilibrium vaporization ratios (K values) of the hydrocarbon constituents or on the K value of carbon dioxide itself. It was shown that carbon dioxide is more soluble in crudes than in distillates which is contrary to the behavior of methane. A working chart of carbon dioxide K values is presented. INTRODUCTION The study of the equilibrium vaporization ratios of mixtures of paraffin hydrocarbons has been rather thorough.2,6,7,8,9 In the past few years considerable attention has been paid to the vaporization characteristics of the so-called noncondensable gases such as nitrogen, carbon dioxide, and hydrogen sulfide in mixtures of hydrocarbons. since they usually occur to some extent in most crude oils and natural gases.1,3,4,5 Knowledge of this behavior is useful to both the production and refining phases of the petroleum industry. This paper reports the equilibrium vaporization ratios (K's) of carbon dioxide in a mixture of League City natural gas and Billings crude oil, and compares them to those obtained in a natural gas-distillate system. The equilibrium vaporization ratios for the hydrocarbon components in this system had previously been studied by Roland.' In addition to the determination of the K values for carbon dioxide, the K values for methane and ethane were also determined in order to observe what effect, if any, the presence of carbon dioxide had on these K values. The concentration of carbon dioxide was also varied in order to observe the effect of this variable on the carbon dioxide K values. EXPERIMENTAL PROCEDURE The apparatus used in this study cotlsisted of a stainless steel equilibrium cell of about 2 liters capacity. The cell was mounted on trunions permitting rocking in a thermostatically controlled oil bath. Two high pressure valves fitted with steel tubing were mounted on the top of the cell. one was used for sampling the equilibrium gas phase and the other for sampling the equilibrium liquid phase by means of an induction tube within the cell. Stainless steel tubing from the bottom of the cell led to a mercury reservoir and manifold which was connected to a free-piston type pressure gauge manufac- lured by the American Instrunlent Ctr. and to a volumetric. putrip. The temperature of the oil bath was measured by means of a ralibrated mercury-in-glass thermometer. The recorded temperatures are believed to be accurate to ±0.5 °F. The pressures are correct to 22 psi. The crude oil used in this study was stock tank oil obtained from the Wilcox formation in the Billings Field, Noble County. Okla. The natural gas was obtained from the League City Field. Galveston County, Tex. The oil was treated with anhydrous calcium sulfate in order to remove the last traces of water. To insure a supply of constant composition gas at room temperature the cylinders of League City gas were cooled to about 30°F, inverted, and the condensed liquid was allowed to drain from the cylinders. The analysis of the gas and crude are tabulated in Table I. The carbon dioxide came from Pure Carbonic, Inc., and was .stated to have a purity of 99.5 per cent or better. The procedure used to obtain samples of the equilibrium liquid and vapor was similar to that employed by others making use of the rocking type equilibrium cell.6,7,8 The equilibrium cell was evacuated and calculated quantities of carbon dioxide, natural gas, and crude oil were charged to the cell to the desired pressure. In charging the equilibrium cell an attempt was made to maintain the ratio of the natural gas to crude oil as close as possible to that employed by Roland. After the cell was charged, samples of

Jan 1, 1951

By E. W. Filer, L. S. Darker

ONE of the primary purposes of this investigation was to determine how far blast-furnace metal and slag depart from equilibrium, particularly with respect to sulphur distribution. In studying the equilibrium between blast-furnace metal and slag, there are two approaches that can be used. One method is to use synthetic slags, as was done by Hatch and Chipman;' the other is to equilibrate the metal and slag from the blast furnace by remelting in the laboratory. In the set of experiments here reported, metal and slag tapped simultaneously from the same blast furnace were used for all the runs. The experiments were divided into two groups: 1—a time series at each of three different temperatures to determine the t.ime required for metal and slag to equilibrate in various respects under the experimental conditions of remelting, and 2—an addition series to determine the effect of additions to the slag on the equilibrium between the metal and slag. An atmosphere of carbon monoxide was used to simulate blastfurnace conditions. The furnace used for this investigation was a vertically mounted tubular Globar type with two concentric porcelain tubes inside the heating element. The control couple was located between the two porcelain tubes. The carbon monoxide atmosphere was introduced through a mercury seal at the bottom of the inner tube. On top, a glass head (with ground joint) provided access for samples and a long outlet tube prevented air from sucking back into the furnace. The charge used was iron 6 g, slag 5 g for the time series, or iron 9 g, slag 7 % g for the addition series. This slag-to-metal ratio of 0.83 approximates the average for blast-furnace practice, which commonly ranges from about 0.6 to 1.1. A crucible of AUC graphite containing the above charge was suspended by a molybdenum wire in the head and, after flush, was lowered to the center of the furnace as shown in Fig. 1. The cylindrical crucible was 2 in. long x % in. OD. The furnace was held within &3"C of the desired temperature for all the runs. The temperature was checked after the end of each run by flushing the inner tube with air and placing a platinum-platinum-10 pct rhodium thermocouple in the position previously occupied by the crucible; the temperature of the majority of the runs was much closer than the deviation specified above. The couple was checked against a standard couple which had been calibrated at the gold and palladium points, and against a Bureau of Standards couple. The carbon monoxide atmosphere was prepared by passing COz over granular graphite at about 1200°C. It was purified by bubbling through a 30 pct aqueous solution of potassium hydroxide and passing through ascarite and phosphorus pentoxide. The train and connections were all glass except for a few butt joints where rubber tubing was used for flexibility. The rate of gas flow was 25 to 40 cc per min. As atmospheric pressure prevailed in the furnace, the pressure of carbon monoxide was only slightly higher than the partial pressure thereof in the bosh and hearth zones of a blast furnace—by virtue of the elevated total pressure therein. Simultaneous samples of blast-furnace metal and slag were taken for these remelting experiments. The composition of each is given in the first line of Table I. There is considerable uncertainty as to the significant temperature in a blast furnace at which to compare experimental results. This uncertainty arises not only from lack of temperature measurements in the furnace, but also from lack of knowledge of the zone where the slag-metal reactions occur. (Do they occur principally at the slag-metal interface in the crucible, or as the metal is descending through the slag, or even higher as slag and metal are splashing over the coke?) The known temperatures are those of the metal at cast, which averages about 2600°F, and of the cast or flush slag, which is usually about 100°F hotter. To bridge this uncertainty, remelting temperatures were chosen as 1400°, 1500" (2732°F), and 1600°C. For the time series the duration of remelt was 1, 2, 4, 8, 17, or 66 hr; crucible and contents were quenched in brine. The addition series were quenched by rapidly transferring the crucible and contents from the furnace to a close-fitting copper "mold." Of incidental interest here is the fact that the slag wet the crucible

Jan 1, 1953

By A. F. van Everdingen

The pressure drop in a well per unit rate of flow is conrolled by the resistance of the formation, the viscosity of the fluid. and the additional resistance concentrated around the well bore resulting from the drilling and completion technique employed and, perhaps, from the production practices used. The pressure drop caused by this additional resistance is defined in this paper as the skin effect. denoted by the symbol S. This skin effect considerably detracts from a well's capacity to produce. Methods are given to determine quantitatively (a) the value of S, (b) the final build-up pressure, and (c) the product of average permeability times the thickness of the producing formation. INTRODUCTION Equations which relate the pressure in a well producing from a homogeneous formation with pressures existing at various distances around the well are generally used within the industry. The relation ii quite simple when the fluid flowing is assumed to be incompressible. It becomes somewhat more complicated when the flowing fluid is considered compressible so that the duration of the flow can he considered. In each case the major portion of the pressure drop occurs close to the well bore. However analyses of pressure build-up curves indicate that the pressure drop in the vicinity of the well bore is greater than that computed from these equations using the known, physical characteristics of the formation and the fluids. In order to explain there excessive drops it is necessary to assume that permeability of the formation at and near the well bore is substantially reduced as a result of drilling. completion and, perhaps. production practice. This possibility has been recognized in the literature. A method to compute the pressure drop due to a reduction of the permeability of the formation near the well bore. which is designated as the skin effect. S, is given in the following paragraphs. To start, equations normally used to describe flow in the vicinity of a well are given without considering this effect. These equations then are modified to include the effect of a skin on the pressure behavior. Finally a method is given to estimate the effect of the skin on the pressure and production behavior of a well. PRESSURE EQUATIONS Incompressible Fluid Flow If p is defined as the flowing pressure in a well of radius the pressure at distance r from the well has been shown to be:" The total pressure drop between the drainage boundary, and the well bore is given by These equations are valid only if the flow towards the well occurs in a horizontal homogeneous medium and the fluids are incompressible. The assumptions imply that all fluid taken from the well enters the system at r a condition rarely encountered in practice. Compressible Fluid Flow, Steady State A more realistic equation is obtained if it is assumed that the compressibility, c, of the flowing fluids is small and has a constant value over the pressure range encountered. After the well has been producing for some time so that its rate has become constant and steady state is reached, the pressures throughout the drainage area are falling by the same amount per unit of time, and the pressure differences between a point in the drainage area and the well are constant. When these conditions are met. the rate of production, q, from a well is equal to where dp/dt is the pressure drop per unit time. The fluid flowing at a distance from the center of the well is equal to From the last equation and from Darcy's law it can be shown that The equation holds for a depletion-type reservoir of radius drained by a well located in its center, provided the compressibility of the fluid per unit pressure drop is small and constant, and no fluid moves across the boundary Compressible Fluid Flow — Nonsteady State Table 111 of reference (5) shows the relationship between the pressure at the well bore and the reduced time, The pressure-drop function, p represents the drop below the original reservoir pressure, p caused by unit rate of production for several values of R, the ratio of drainage boundary radius to well radius, r In most reservoirs the values of approach infinity. and under these conditions the values of p shown in Table I of reference (5) can be used where p then signifies the difference between the pressure in the well and the prevailing reservoir pressure per unit rate of flow. The total pressure drop below prevailing reservoir pressure amounts to where the factor converts the cumulative pressure drop per unit rate of production to cumulative pressure drop for actual rate. q. For values of T > 100 the P function may be written (equation VI-15 of reference 5) as Using the time conversion the difference in pressure between reservoir and well becomes If values for the physical constants of the formation and the fluids are inserted, it is found that T exceeds 100 after a few seconds of production (or closed-in time), so that the approximation becomes valid almost at once. A simple relation between the pressure in the well and in the reservoir can also be derived by considering the well as a point source"" '" instead of a unit circle source, that is, by using Lord Kelvin's solution instead of the unit circle source

Jan 1, 1953

By A. S. Odeh

The mechanism of two-phase flow in porous media has been a subject of wide controversy. One of the properties essential for understanding the dynamic behavior of two-phase flow is relalive permeability. Relative permeability to a certain phase is defined as the ratio of the effective permeability of that phase to its permeability when it is the only fluid present and powing. In this research, a theoretical analysis was made to determine the effect of viscosity ratio between the non-wetting and the wetting phase on relative permeability. Experimental work was conducted to test the validity of the derived equations. The experiment was conducted on four natural cores. Four oils were used as the non-wetting phases with a viscosity range of 0.42 to 71.30 cp and two wetting phases with a viscosity range of 0.86 to 0.96 cp. Oil and bring were made to flow simultaneously at various ratios, and relative permeability curves were constructed. A total of eight relative pertileability cycles representing eight viscosity ratios were run oil each sample. It was found that relative permeability to the non-~cletting phase varies with viscosity ratio. The relative effect of this variation on relative permeability values was a function of the sample's single-phase permeability, decreasing with its increase. It was concluded that, for .samples of single-phase permeability over I darcy. the effect of viscosity ratio could be disregarded, and relative permeability would be, in effect, a function of satrtration only. INTRODUCTION Two-phase as well as multiphase flow occurs in many fields of science. This type of flow is of particular interest in petroleum production. The knowledge of relative permeability, which describes the dynamic behavior of two-phasc as well as multiphase flow, is essential for solution of problems arising in that field. Thc relative permeability ot a porous medium to a given phase in multiphase flow. is generally considered to be only a function of the saturation of that phase, independent of the properties of fluids involved and ranging in value from zero to unity. Work by Leverett' and Leverett and Lewis' apparently supports this concept. In his experiments Leverett used a clean, packed unconsolidated sand of high permeability (3.2 to 6.2 darcies) with two phases (water and oil) flowing and a viscosity ratio range of 0.057 to 90.0. His results showed that the wide range of viscosity had practically no effect on relative permeability-saturation relationship. Recently accumulated evidence from work performed by several laboratories and a paper by Nowak and Krueger,2 in which relative permeability to oil of a few core samples in the presence of interstitial water was considerably greater than single-phase permeability to water, cast some doubt on the conclusions reached by Leverett' and subscribed to by a large number of individuals in the oil industry. One explanation advanced to explain this behavior states that it is caused by the variable extent of hydra-tion of clay minerals present in the sand. The greater the water saturation, the greater will be the area of contact between water and clay minerals; therefore, the greater will be the extent of swelling with corresponding reduction in permeability. Yuster4 presents another explanation for the recently accumulated evidence. Utilizing Poiseuille's law, he analyzed concentric flow in a single capillary where the non-wetting phase flows in a cylindrical portion of the capillary and concentric with it. The wetting phase flows in the annulus between the non-wetting phase and the capillary wall. The equations obtained indicate that relative permeability to the non-wetting phase is a function of saturation and viscosity ratio. Although Yuster's equations show that fractional rel-ative permeability to oil could be greater than unity, as was indicated by the data of Nowak and Krueger,1 they failed to present an explanation to the experimental data of early investigators such as Leverett.1 Due to the importance of relative permeability in understanding the flow behavior of petroleum reservoir fluids, this work—theoretical as well as experimental —was undertaken to determine whether relative permeability is a function of saturation only as was concluded by Leverett1 or a function of saturation and viscosity ratio as was theorized by Yuster.4 THEORETICAL ANALYSIS An equation will be derived for the rate of oil flow through a porous medium that is initially filled with water. Based on this equation, an analytic expression for relative permeability will be developed. The porous medium will be assumed to consist of .straight circular capillaries of different radii. It will also be assumed that there are no interconnections among the capillaries and no mass transfer across the oil-water interface. Consider a porous sample initially saturated with a wetting phase (water). As a non-wetting phase (oil) is

By Hugh H. Waesche

Of all the military services, the Signal Corps is the most concerned with piezoelectric minerals because of its function as a supply service to the strategic and tactical military forces. Consequently this paper is written from the point of view of one associated with that organization. The Signal Corps is responsible for the research, development, and supply of communications, radar, and components to the using services of the Department of the Army and to some extent the Other branches of the National Defense Department. Their work therefore includes the study of the sources* characteristics, and application of quartz and other piezoelectric materials. These materials have become a vital consideration in strategic planning and are essential for efficient tactical operation by all the Armed Forces. The Signal Corps at the beginning of world War 11 Was respon-sible for both Army Ground and Air Force electronic equipment. Since that time this Army service organization has probably done more in the development of frequency control devices using piezoelectric materials than any other group. The U.S. Department of the Interior, Bureau of Mines, Minerals yearbook of 1945, shows that during the four war years, 1942 through 1945, 9,598,-410 Ib of quartz crystal were imported for all uses and of this total, 5,168,000 lb were consumed to produce 78,320,-000 crystal units for electronic application. Other government records confirm these data which conclusively show that approximately 53 pct of the crystalline quartz imported was consumed in the production of electronically applied quartz crystal units. It may be assumed that some effort was made to maintain a stockpile over demands for all purposes. and this would mean that the actual percentage of quartz used electronically was considerably over the 53 pct figure. These data only emphasize that electronic application of crystalline quartz was the greatest requirement, and per- haps the actual value in this application to national defense is many times greater in importance than is apparent on first inspection. Current electronic research and development programs of the Armed Forces are planned around the fundamental use of piezoelectric minerals for frequency control and this at present, at least, means quartz. Definition and Early Development The word piezoelectricity is formed from combination of the Greek word "piezein". meaning "to press," and "electricity." It is that property shown by numerous crystalline substances whereby electrical charges of equal and opposite value are produced on certain surfaces when the crystal is subjected to mechanical stress. It appears to be intimately associated with the better known property, pyro-electricity and in fact, the two may be manifestations of the same phenomeuon. This property was discovered by Pierre and Jacques Curie in quartz, tourmaline, and other minerals in 1880 while studying the symmetry of crystals. The converse effect, that is, mechanical strain in the crystal when placed in an electrical field, was predicted by the French physicist, G. Lippman, in 1881, and verified by the Curies almost immediately. As has been the case with many discoveries of similar character in the basic sciences, not much attention was paid to this property for man)- years except as an entertaining curiosity. Between 1890 and 1892 a series of papers was published by W. voigt in which the theoretical physical properties were put into mathematical form. The first practical application of piezoelectricity occurred during World War I when professor P. Langevin of France used quartz mosaics to produce underwater sound waves. The same mosaics were used to pick up the sound reflections from submerged objects which were in turn, amplified by electronic means and used to determine the distances to such objects. This device was intended for use as a submarine detector but development was not completed in time for war service although it was used later for determining ocean depths. About the same time, A. M. Nicholson, of Bell Telephone Laboratories, developed microphones and phonograph pickups using Rochelle salt crystals. A major step in the application of piezoelectric quartz came in 1921, when professor W. G. Cady, of wesleyan university, showed that a radio oscillator could be controlled by a quartz crystal; from that date, this use of quartz has increased steadily, reaching its peak in world war 11 as is shown by the figures previously given. Essentially all American electronic equipment for communication, navigation, and radar, utilized quartz crystals for the exacting frequency control required. Crystalline Minerals with piezoelectric Properties QUARTZ Hundreds of piezoelectric crystalline materials are known, most of which are water soluble. Of these, quartz appears to be without a peer for electronic frequency control. Unfortunately, the quartz must be of superior quality. It must be free of mechanical flaws, essen-tially optically clear, free of both Brazil and Dauphiné twinning and must be, for average uses, over 100 g in weight. Because of these stringent requirements, raw quartz of the quality desired is of rare occurrence. In addition to quartz, several other naturally occurring crystalline materials are known to have the piezoelectric property and could perhaps be substituted for quartz in some applications. These

Jan 1, 1950

By Michel Dupuy, Daniel Calais

The study of self-diffusion of plutonium in E phase has been carried out by the welded couples method. The tracer used was puZ4O which is detected by its X-ray emission (conversion lines of uranium which are computed between 13 and 21 kev). Intensities were measured with a scintillation counter. Each layer was removed in a direction parallel to the original interface with a grinding machine and a thickness measured with a pneumatic comparator. The concentration-penetration curves obtained were corrected for the effect of heating time from room temperature to annealing temperature and for the expansion due to phase transformations of plutonium. They were analyzed by the generalized Gruzin method. Self-diffusion of plutonium in E Phase is very fast cm per sec between 500" and 620°C) and the diffusion zones are 2 to 3 mm wide for annealing times ranging from 30 min up to 10 hr. The Arrhenius law gives the temperature dependence in the form: From the point of view of self-dqfusion, PUE phase falls into the anomalous bcc metals category (Tip , Hfp, Zrp, Uy) with a low-frequency factor and an activation energy lower than those provided by standard correlations. No theory proposed hitherto to explain these anomalies (influence of dislocations, of extrinsic vacancies bonded to inlpurities, of bi-vacancies) can clearly explain the self-diffusion coeffzcients of plutonium. DIFFUSION in bcc metals is a present-day problem. A recent symposium (Gatlinburg, 1964), followed by a book,' has been devoted to it. A great many experiments seem to show that diffusion in certain bcc metals obeys unexpected laws. The activation energies measured are sometimes strangely low (B hafnium, y uranium). For certain metals (0 zirconium, p titanium) the curves of log D (D = diffusion coefficient) as a function of 1/T (T = absolute temperature) are not linear. The frequency factors Do, which are of the order of 1 sq cm sec-' in fcc metals, vary from 1 to 10~6 sq cm sec-'. Various theories have been put forward to explain these anomalies; none is yet satisfactory. We wished to introduce a new experimental result by studying the self-diffusion in c plutonium. This allotropic phase, stable from 475°C up to the melting point (640°C), is in fact bcc. Unfortunately, nothing is known of the characteristics of the point defects in this phase, which limits the scope of the hypothesis which can be made about the mechanism(s) of self-diffusion in plutonium. 1) EXPERIMENTAL METHODS 1) Principle. We used the welded couple method. The two pellets of the couple initially have different 240 isotope contents (X emitter). After diffusion, the concentration/penetration curves are drawn up by the generalized Gruzin method. 2) Gamma Spectrography. The metal used in our study is plutonium, either low in puZ4O (isotopic content 1 pct) or high in puZ4O (8 pct). The latter also contains plutonium 241 (-1 pct) and 300 ppm of ameri-cium produced by the reaction Pu2U-AmM1 + 8-. The emission spectra of these two plutoniums placed in leak-tight vinyl bags have been studied by y spectrograph~. The detector is a thin crystal of thallium-doped sodium iodide. The activity of the plutonium rich in 240 is about twice that of the plutonium low in 240 in the energy band of 17 kev (L conversion lines of uranium); this band was used in these measurements. 3) Preparation and Examination of the Diffusion Couples. Diffusion couples were made from plutonium with a high and low PU"' content. Pellets (6 6 mm. thickness 3 mm) mounted on a polishing disc with ground parallel faces were polished mechanically on both sides. In this way, pellets with two parallel faces were easily obtained. The polished pellets were joined by a 6 phase anneal (420°C, 1 hr) in a small screw press (pressure of 20 kg per sq mm cold); a centering ring enabled the two pellets to be pressed coaxially. The couples then were subjected to the diffusion treatment by annealing in the E phase in sealed silica ampules in argon at atmospheric pressure. The annealing temperatures and times are given in Table I. The couples were encased in a mild steel ring, the joint interface being thus parallel to the ground face of the ring. The diffusion couple/ring assembly underwent successive abrasions by means of a magnetic plate grinder. The thickness of the abraded layer was measured with a Solex pneumatic comparator when it was less than 0.1 mm (accuracy 0.2 p) or with a mechanical micrometer (accuracy 3 p) for passes of the order of 0.2 mm. All these operations were done in glove boxes, as plutonium is particularly toxic. After each abrasion we determined the emission spectrum of the ground face. The emissive surface is defined by means of a diaphragm 3 mm in diam. We noted more particularly the X activity in the 17-kev

Jan 1, 1969

By Harry A. Wahl

Precent fracture design procedures are bared on the total fracture area created. A method to distinguish beI,,.ecn total area and [he propped or effective fracture area has not been available. This paper presents a solution to this problem, applicable to horizontal fractures. The difierences between effective fracture area and torn] area are demonstrared in example calculations. This work is hayed on experimentally determined transport efficiencier of solids in sand-liquid slurries. Newtonian and non-Ne~vtoninn systems are considered. INTRODUCTION Fifteen years after commercial introduction, hydraulic fracturing remains the most successful stimulation technique in the oil field. This success is primarily due to ability of induced fractures to penetrate and alter permeabilities deep within formations. Many fields producing today could not have been developed without the hydraulic fracturing process. Because of wide usage, fracture-treatment design has received a great deal of engineering and research effort. This work, resulting in improved equipment and materials, has increased the benefits from fracture treatments as well as the applicability of the process. A major contribution was the development of fluid-loss additives. Necessarily, the number of parameters to be considered in treatment design has steadily increased, resulting in more complicated design techniques. Almost all present design procedures are based on the precepts set forth by Howard and Fast. Relating the fluid volume lost into the formation, the volume required in extending the fracture, and the total slurry volume injected, they developed an expression for the total fracture area created in terms of pertinent treatment parameters. Fluid loss during treatment was expressed as a function of time for three flow mechanisms. Although modifications of fluid loss equations have been made, the total fracture area concept has remained unaltered. A vast amount of field data indicate that induced fractures must be propped and held open to be effective. A notable exception is the Mesa Verde formation in the San Juan basin. However, analysis of these treatments shows that improved well productivities are obtained when propping agents are incorporated in the treating fluid. Although propped fracture area has been recognized as an important design parameter, a method to distinguish between total area and effective fracture area has not been available. The necessary information on slurry-sand transport in fractures has been lacking. Interest in the propped region of induced fractures is not restricted to areal extent alone. The distribution of sand within fractures is important from the standpoint of fracture flow capacities. Flow capacity affects the increase in well productivity after stimulation. The work of Huitt and Darin4 hows that partial monolayers of sand have large flow capacities compared to thick* dense sand packs. It has been postulated that gelled fluids have the ability to transport sand within the fractures at the deired low concentrations. An early contribution in the area of sand placement in fractures was made by Kern et al.' They studied sand movement in a transparent vertical fracture model. It was observed that the sand tended to settle out in the bottom of the model before moving very far. When the fluid velocity exceeded a certain critical value, all of the sand injected began moving through the crack even though it settled to the bottom. This critical velocity was determined under several flow conditions. Some work on sand movement in horizontal fractures has been reported in Russian publications. Sand movement was studied by Izyumova and Shan'gin' using a transparent "pie-shaped" flow model to simulate a horizontal radial flow system. However, the data were limited, especially in a quantitative sense. Dorozhkin, Zheltov and Zheltow studied the behavior of sand-liquid slurries in a horizontal linear flow model. The quantitative data were restricted primarily to the thickness of sand deposits formed at the bottom of the fracture. An earlier paper provided basic data on the flow of sand in horizontal fractures. This study was designed to yield specific quantitative information on rate of advance of sand particles and pressure behavior under various flow conditions. A comprehensive photographic study was undertaken in a 10-ft windowed flow cell to provide the necessary qualitative and quantitative data. Since the number of potential variables far exceeded the capacity of the initial study, emphasis was placed on the effects of sand concentration, oil viscosity and oil flow rate. A detailed description of these experiments and the results are described in Ref. 9. However. the implications of this work on the fracture design calculations were not discussed. An analysis of these data as well as new data is provided in the following sections. EXPERIMENTAL RESULTS The primary objective of the experimental investigation was to provide information on the rate of advance of the solids in sand-liquid sturries. A 10%-ft long transparent

Jan 1, 1966

By L. A. Wilson, P. J. Root

The relative costs of heating a reservoir by steam injection and by combustion have been examined. The comparison was based on a model similar to that proposed by Chu.' The cost of boiler feed water, the price of fuel, pressure and plant capacity were parameters in determin-ing the costs of air compression and steam generation. The analyses compare the cost of heating to the same radius by the two methods. Results suggest that the two primary factors for comparison are the price of fuel and the amount of crude burned during underground combustion. The cost of fuel has a greater effect on the cost of heat from steam than it does on its cost by combustion. As a result, analyses indicate that when the price of fuel is low, steam may be unequivocally cheaper than air. The influence of heat loss is such, however, that as the heated radius increases combustion becomes relatively more competitive depending upon the amount of crude burned. This implies that steam may be cheaper for small stimulation jobs (huff and puff) but combustion may be more economically attractive for heating large areas (flooding). INTRODUCTION Use of thermal methods of recovery is an accepted fact today. After an induction period of several years, processes are being widely used that involve reservoir heating to augment recovery. Of the several techniques, steam injection and forward combustion appear to be destined to dominate the field. Although the objectives of both are the same, the basic differences between generating heat in situ and injecting heat after surface generation influence the cost in different ways. This study compares the cost of heating a reservoir either by steam injection or by forward combustion. There has been no consideration of recovery. Presumably, recovery from the swept region would be high in either case. The sole consideration was the cost of heating to the same radial distance by either process. PROCEDURE THE MODEL The basis for comparison was a mathematical model similar to that used by Chu' for combustion. The model simulates a radial heat wave in two-dimensional cylindricaI coordinates. It includes heat generation, conduction and convection within the reservoir and conduction in the bounding formations. Thus, heat losses from the formation are considered. Three significant modifications were made. 1. Equal logarithmic increments rather than equal increments were used for the mesh spacing in the r direction. By this technique large distances were simulated with relatively few mesh spaces. 2. A backward difference approximation to the convection term was used to avoid troublesome oscillations which result from a central difference approximation when the convection term is large. 3. The radial increments of the combustion zone motion were not necessarily uniquely related to the mesh configuration. The cumbersome step function introduced by the heat of vaporization of steam was circumvented by assuming the enthalpy of the steam to be a linear function of temperature between reservoir temperature and steam temperature. This is equivalent to assuming an average heat capacity numerically equal to the difference between the enthalpy of saturated steam and the enthalpy of water at reservoir temperature divided by the difference between the two temperatures. Heat losses obtained by this model are in essential agreement with those obtained by the analytical solution of Rubenshtein.' A detailed description of the model is presented in the Appendix. Using the model, the times required to heat to particular radial distances were obtained as a function of injection rate and other physical parameters. For the steam case, injected fluid was assumed to be saturated steam at pressures of either 500, 1,000 or 1,500 psia. The corresponding temperatures are 467, 544 and 596F, respectively. Thickness ranged from 10 to 50 ft and injection rate ranged from 100,000 to 1 million Ib/D. Reservoir and overburden temperatures at the injection well were assumed to be that of saturated steam at the injection pressure. The effect of maintaining the overburden temperature at the well at a different temperature (initial reservoir temperature) was examined with no significant change in behavior. The influence of wellbore heat losses for the steam case was determined in the following manner. The rates of heat loss as a function of time were estimated using an approach similar to that suggested by Ramey." he data were based on injection through 2%-in. tubing in 7-in. casing. Integration of these data over the entire iniection period yielded the total heat loss. Total heat losses were then corrected to their equivalents in steam (this number resulted from dividing the total heat loss by the latent heat). This was considered additional steam required to accomplish the reservoir heating and the total cost was increased accordingly.

Jan 1, 1967

By G. P. Airey, R. F. Mehl, T. A. Hughes

The coarsening of cementite particles in a ferrite matrix has been studied in a series of steels with 0.15 pct C only and 0.15 pct C plus 1 pct Ni, Mn, and Cr, respectively. Two initial states were employed: quenched nartensite, and quenched and cold-rolled martensite. A series of tempering temperatures between 500' and 700" and tempering times of up to 190 hr were used. The structures were studied by replica and transmission electron microscopy. Particle size distribution curves were determined. From the average size value coarsening curves were obtained. These were plotted in accordance with the Wagner analysis assuming diffusion control. A discussion of the significance of the results is given. L HE reactions occuring upon the tempering of martensite have long engaged the attention of metallurgists. The latter stages, when cementite particles coarsen in a ferrite matrix, have been studied both qualitatively and quantitatively. Studies of such coarsening processes have recently been spurred by the publication of the Lifshitz-Wagner theory1, and the extension of this to the a Fe-Fe3C system by Oriani3 and by Li, Blakely, and einold. Following Wagner the coarsening process is often designated as "Ostwald Ripening". The only quantitative data on the rate of coarsening, except for the work of Hyam and uttin,' in the a Fe-Fe3C system are those of Bannyh, Modin, and odin' for a commercial eutectoid steel and those of Heckel and ereorio" using a pure eutectoid steel. The data of Bannyh, Modin, and Modin have been employed by 0riani3 to derive the a Fe-FeE interface energy. The reaction is one of the most important ones in steel and is worthy of detailed study. This is the purpose of the present study. Laboratory heats were prepared; these were steels with approximately 0.15 pct C, selected so that the num ber of carbide particles would be relatively small and thus so that the overlapping diffusion fluxes would be minimized, presumably a desirable circumstance.'-3 In addition to Fe-C alloys, comparable heats containing 1 pct of Ni, Mn, and Cr, respectively, were included with a view of appraising the effect of alloying elements. This report includes an account of the micro-structures observed, primarily with the electron microscope, and of kinetic data and their interpretation. MATERIALS AND TREATMENT The alloys were prepared from electrolytic iron ("Plastiron") and high-purity graphite; these were melted in a zirconia crucible using a vacuum furnace. The alloy steels were made by adding electrolytic nickel, electrolytic manganese, and "vacuum grade" chromium, respectively, under a partial pressure of argon. Each melt was poured into a mold within the vacuum furnace and cooled in the mold. The ingots were 2 in. in diam. and 8 to 10 in. long. The analysis of the alloys is given in Table I. These ingots were hot-rolled to strip 0.1 in. thick, then cold-rolled to 0.05 in. and each alloy split into two batches. One batch was austenitized at 1200 for 1 min, quenched in cold brine, then cold-rolled to 0.02 in.; samples given this treatment are hereinafter designated as "worked". The other batch was cold-rolled to 0.025 in., austenitized at 1200" for 1 min, and quenched in cold brine; such samples are hereinafter designated as "quenched". These two batches were then tempered, as below. The purpose of the treatment given the first batch was to provide an initial structure of cold-worked martensite, with the expectation that the additional defect structure created by cold work would encourage a higher rate of nucleation of cementite on tempering and hence a more uniform distribution of cementite particles. Individual specimens were sealed in evacuated quartz or Pyrex tubes, then tempered in a muffle furnace. The temperature control was better than 3'C at 700. Tempering treatments wer: performed at 400°, 500°, 550°, 600°, 65o°, and 700C for time periods between 15 min and 190 hr. PREPARATION OF SPECIMENS Specimens for optical and replicalelectron microscopy were mounted, polished conventionally, and etched with 2 pct nital. For electron microscopy, single-stage "formvar" replicas were made, dry-stripped and rotary-shadowed with chromium at an angle of 30 deg. Carbide extraction replicas were prepared from electropolished specimens usirig the method described by Smith and uttin.' Thin foils for electron transmission microscopy were prepared by chemical thinning in an H202-HF bath prior to electropolishing in a chromium trioxide-acetic acid solution. The most

Jan 1, 1969

By P. A. Witherspoon, R. W. Donovan, T. D. Mueller

The use of petroleum-barren aquifers for underground storage has become extremely important to the natural-gas industry. A critical problem in assessing the feasibility of a specific aquifer for such use is the permeability determination of the caprock over the proposed storage project. The approach used here is to conduct both static and dynamic field tests on the aquifer being analyzed. Valuable information on the possibility of communication between the storage aquifer and any other aquifers above can be obtained by measuring hydrostatic water levels and water analyses. Significant differences in such data give evidence of the lack of communication between the intended storage reservoir and other horizons. The dynamic approach requires that one well be pumped in the storage aquifer, and changes in fluid levels recorded in both the aquifer and its caprock. The interpretation of the data from such pumping tests involves the solution of nonsteady radial flow in an infinite aquifer and the influence on such flow of a leaky caprock. A finite-difference model has been used to investigate this problem, and the transient behavior has been solved numerically with a digital computer. It has been found that the pressure transients in the storage aquifer are not affected significantly by moderate caprock leakage. The pressure behavior of the caprock is a much better indicator of the degree of leakage, and generalized solutions for this behavior are included. Field data are presented to demonstrate both the static and dynamic approach. If is concluded that appropriate investigation of the groundwater hydrology in an aquifer-type gas-storage project can provide much valuable information for determining the effectiveness of the caprock to hold gas. INTRODUCTION Underground storage of natural gas in the United States has been developing at a rapid rate over the past few years. In 1955, the total gas-storage capacity was about 1.6 trillion cu ft; by 1961, this figure was almost 3.2 trillion cu ft, an increase of 100 per cent in six years.' This trend un- doubtedly will continue because the economics favor the development of gas storage, as opposed to the construction of new pipelines, to meet the inherent cyclic demand for fuel in the metropolitan areas of this country.' About 15 per cent of the current underground gas storage has been developed in petroleum-barren aquifers, i.e., geological domes or anticlines in which no commercial quantities of oil or gas had been produced prior to the storage operations. The necessity for using barren aquifers outside many metropolitan areas of this country has been due to the lack of depleted oil or gas fields that were near enough and large enough to meet the demands of such consuming areas. Pipeline companies have developed aquifer storage along their transmission lines to meet the fluctuating needs of their complex systems. Considerable thought has also been given to the problem of storing gas in a structureless aquifer, both in this country' and in the Soviet Union outside the city of Leningrad.'," Conditions such as these have led to the development of aquifer gas-storage projects in many parts of the U. S. Most of these developments have centered in the Mid-Continent area, and the greatest amount of activity has been concentrated in Illinois.6 Thus, the use of petroleum-barren aquifers for gas-storage purposes has become extremely important to the natural-gas industry. There are three basic problems in developing aquifer-type storage: (1) finding an adequate geologic structure, (2) finding a suitable storage reservoir within the structure and (3) determining the tightness of the caprock over the intended storage zone. The first two problems can be solved by applying conventional methods of exploration geology, but once these problems are solved, the question arises as to why no oil or gas is present in an otherwise favorable setting. Two situations are possible: (1) an adequate source bed was never present, or (2) a source bed was present but the petroleum seeped away because of a leaky caprock. Determining the tightness of the caprock is one of the most critical problems in assessing the feasibility of a specific aquifer for storage purposes. In attacking this problem, one usually takes cores of the caprock and subjects them to a rigorous investigation. Such core data are desirable, but they only detail the matrix properties and cannot be expected to reveal the gross characteristics of the caprock. Several gas-storage projects in the U. S. have had considerable leakage where

By J. G. Richardson, J. W. Graham

The purpose of this paper is to present the results of theoretical and experimental studies of water imbibition. The imbibition processes are involved in recovery of oil from stratified and fractured-matrix formations in natural water drives and water flooding. An understanding of the role of inhibition in implementing the recovery of oil from such formations is deemed essential to proper control of these reservoirs to achieve maximum recovery. The theoretical studies involved development of the differential equations which describe the spontaneous imbibition of water by an oil-saturated rock. The dependence of the rate of water intake by the rock on the permeability, interfacial tension, contact angles, fluid viscosities and fluid saturatiorls is discussed. A few experiments were performed using core samples to determine the effects of core length and presence of a free gas suturation. The role of water imbibition in recovery of oil from a fractured-matrix reservoir by water flooding was investigated by use of a laboratory model. This model was scaled to represent one element of a frac-tured-matrix formation. Water floods were made at various rates with several fracture widths. Interpretations were made of the behavior expected in a system containing many matrix blocks. The presence of a free gas sntu.ration was found to reduce the rate of water imbibition. In the reservoir prototype of the fractured-matrix model, water imbibition rather than direct displacement by water was the dominant mechanism in the recovery of oil at low rates. INTRODUCTION Imbibition may be defined as the spontaneous taking up of a liquid by a porous solid. The spontaneous process of imbibition occurs when the fuid-filled solid is immersed or brought in contact with another fluid which preferentially wets the solid. In the process of wetting and flowing into the solid, the imbibing fluid displaces the non-wetting resident fluid. Common examples of this phenomenon are dry bricks soaking up water and expelling air, a blotter soaking up ink and expelling air and reservoir rock soaking up water and expelling oil. As increasingly better lithological descriptions have been made of the characteristics of petroleum-bearing formations, it has become obvious that imbibition phenomena which were once considered laboratory curiosities are of practical importance. For instance, in reservoirs composed of water-wet sand strata of different permeability in intimate contact, the tendency of water to channel through the more permeable stratum is offset by the tendency for water to imbibe into the tight sand and expel oil into the coarse sand. Also, in fractured-matrix formations the tendency of water to channel through the fractures is offset by water-wet matrix blocks. As some imbibition of the water into the of the largest fields in the world are fractured-matrix reservoirs, it has become increasingly important to understand all the factors involved in the imbibition process. Examples of fractured-matrix reservoirs are the Spraberry field in West Texas which produces from a fractured sandstone', the giant Kirkuk field in Iran', the Dukhan field in Qatar, Persian Gulf2, and the Masjid-I-Sula-main and the Haft-Kel fields in Southwestern Iran, which produce from fissured limestone3. Research into recovery of oil from fractured-matrix formations was stimulated by the rapid decline of oil productivity of wells in the Spraberry formation. One result of this research was the water imbibition process developed by the Atlantic Refining Co.4 Another idea was that much of the Spraberry oil could be recovered by conventional water-flooding procedures5. Subsequently, pilot floods were conducted in this field to test the feasibility of these ideas. It was felt that an understanding of the role played by imbibition processes in displacement of oil from a fractured-matrix reservoir could not be obtained from field data alone because of the many complicating factors and uncertainties involved. Therefore, theoretical and laboratory studies were undertaken to provide this understanding. Study of the equations which describe the linear, countercurrent imbibition process provided an insight into the role of various factors in the process, such as the permeability of rock and inter-facial tension. In addition to the theoretical studies, imbibition experiments were conducted with core samples to determine the effect on the rate of imbibition of such variables as core length and free gas saturation. The principal experimental studies were conducted by water flooding a scaled model of an clement of a frac-tu red-matrix reservoir to evaluate

By H. R. Thresh

An oscillational vis cometer has been constructed to measure the viscosity of liquid metals and alloys to 800°C. An enclosed cylindrical interface surrounds the molten sample avoiding the free surface condition found in many previous measurements. Standardization of the apparatus with mercury has verified the use of Roscoe's formula in the calculation of the viscosity. Operation of the apparatus at higher temperatures was also checked using molten lead. Extensive measurements on five different samples of zinc, of not less than 99.99 pct purity, indicate i) impurities at this level do not influence the viscosity and ii) the apparatus is capable of giving reproducible data. The variation of the viscosity ? with absolute temperature T is adequately expressed by Andrade's exponential relationship ?V1/3 = AeC/VT , where A and C are constants and V is the specific volume of the liquid. The values of A and C are given as 2.485 x 10-3 and 20.78, 2.444 x 10-3 and 88.79, and 2.169 x 10-3 and 239.8, respectively, for mercury, lead, and zinc. The error of measurement is assessed to be about 1 pct. Prefreezing phenomena in the vicinity of the freezing point of the zinc samples were found to be absent. AS part of an over-all program of research on various phases of melting and casting nonferrous alloys, a systematic study of some physical properties of liquid metals and their alloys was undertaken in the laboratories of the Physical Metallurgy Division.1,2,3 The most recent phase of this work, on zinc and some zinc-base alloys, was carried out in cooperation with the Canadian Zinc and Lead Research Committee and the International Lead-Zinc Research Organization. One of the properties investigated was viscosity and the present paper gives results on pure zinc; the second part, on the viscosity of some zinc alloys, will be reported separately. Experimental interest in the viscosity of liquid metals has virtually been confined to the past 40 years. The capillary technique was already established as the primary method for the viscosity of fluids in the vicinity of room temperature; all relevant experimental corrections were known and an absolute accuracy of 1 to 2 pct was possible. Ap- plication of the capillary method to liquid metals creates a number of exacting requirements to manipulate a smooth flow of highly reactive liquid through a fine-bore tube. Consequently, the degree of precision usually achieved in the high-temperature field rarely compares with measurements on aqueous fluids near room temperature. However, the full potential of the capillary method has yet to be explored using modern experimental techniques. As an alternative, many investigators in this field have preferred to select the oscillational method. Unfortunately, the practical advantages are somewhat offset by the inability of the hydrodynamic theory to realize a rational working formula for the calculation of the viscosity. In attempting to overcome this restriction many investigators have employed calibrational procedures, even to the extent of selecting an arbitrary formula for use with a given shaped interface. However, where calibration cannot be founded on well-established techniques, the contribution of such experiments to the general field of viscometry is questionable. A critical appraisal of the viscosity data existing for pure liquid metals reveals a somewhat discordant situation where considerable effort is still required to establish reproducible and reliable values for the low-melting point metals. The means of rectifying this situation have gradually evolved in recent years. Here, the theory of the oscillational method has undergone major advances for both the spherical and cylindrical interfaces. The basic concepts of verschaffelt4 governing the oscillation of a solid sphere in an infinite liquid have been adequately expressed by Andrade and his coworkers.5,6 Employing a hollow spherical container and a formula, which had been extensively verified by experiments on water, absolute measurements on the liquid alkali metals were obtained. The extension of this approach to the more common liquid metals has been demonstrated by culpin7 and Rothwel18 where much ingenuity was used to surmount the problem of loading the sample into the delicate sphere. Because of the elegant technique required to construct a hollow sphere, the cylindrical interface holds recognition as virtually the ideal shape. On the other hand, loss of symmetry in one plane increases the complexity of deriving a calculation of the viscosity. The contributions of Hopkins and Toye9 and Roscoe10 have markedly improved the potential use of the cylindrical interface in liquid-metal viscometry. The relatively simple experi-

Jan 1, 1965

By C. Laird, C. E. Feltner

May different kinds of dislocation structures have been observed in strain-cycled metals and alloys. In order to understand their pattern and causes, an experimental program has been carried out to determine the influence on the dislocation structures of the three variables: 1) slip character of the material, 2) test temperature, and 3) strain amplitude. The results show that at high strain amplitudes cell structures me formed when the slip character is wavy, and that these are progressively replaced by uniform distributions of dislocations as the stacking fault energy is decreased. At lower strains, dislocation debris is formed which consists primarily of dipoles in wavy slip mode materials and multipoles in planar slip mode materials. Temperature merely acts to change the scale of the structure, smaller cells, and clumps of dislocation debris being associated with lower temperatures. It is shown that the results for many metals fit this pattern, which Parallels that occurring in unidirectional deformation. DISLOCATION structures produced by cyclic strain (fatigue) have been examined in a number of metals by transmission electron microscopy. These studies have produced a variety of interesting and often seemingly conflicting results. For example, different investigators have reported such structural features as cells.le4 bands of tangled dislocations,4'5 dense patches or clusters of prismatic dislocation loops, planar arrays,4'10 and various combinations or mixtures of these different structures. Most of these observations have been made on materials which were initially annealed and cyclically strained at low amplitudes resulting in long lives. Recently we have reported observations of the dislocation structures produced in copper and Cu-7.5 pct Al cycled at large amplitudes, resulting in lives of less than 104 cycles.4 These results, examined in combination with those in the literature, have suggested that a common or consistent structural pattern exists. Variations in this pattern appear to be determined chiefly by the three variables, namely, the slip character of the material,4,11 test temperature. and the strain amplitude. To verify this interpretation, we have studied [he influence of the above three variables (in different combinations) on the resultant structures in cyclically strained metals. Copper, fatigued at room temperature, was chosen as a reference state to which all other observations can be compared. The effect of slip character has been investigated by employing fcc metals of different stacking fault energy. Thus aluminum which has a more wavy slip character than copper, and Cu-2.5 pct A1 having a more planar slip char- acter, have been examined. The aluminum samples were fatigued at 210°K thus making their homologous temperature equal to that of copper at room temperature. The influence of temperature has been evaluated by examining the structures in copper at room temperature and 78°K. Finally the effect of strain amplitude was studied by looking at the structures at amplitudes giving lives ranging from 104 to 107 cycles. All of the specimens were examined at the 50 pct life level at which stage the structures have reached a stable configuration.12 I) EXPERIMENTAL PROCEDURE Strip specimens, 0.006 in. in thickness, were prepared from base elements of 99.99 pct purity or greater. Specimens were fatigued by cementing the strips to a lucite substrate which was subjected to reverse plane bending. This method of testing has been described e1sewhere.7 After fatiguing, specimens were thinned and examined in a Philips EM 200 which was equipped with a goniometer stage capable of ±30-deg tilt and 330-deg rotation of the specimen. On the basis of separate calibrations,13 allowances were made for the relative rotation and inversions between the bright-field images and the diffraction patterns. II) RESULTS AND DISCUSSION The life behavior of the materials under different test conditions is shown in Fig. 1 in the form of plots of total strain range vs cycles to failure. Comparisons of structures produced in the different materials were made at amplitudes which produced equal numbers of cycles to failure. The influence of strain amplitude on the structures produced in the reference state material (copper tested at room temperature) is shown in Fig. 2. At the 104 life level the structure produced comprises cells similar to those previously observed.3,4 They are approximately 0.5 p in diam and the cell walls are generally more regular or sharper than those produced by unidirectional deformation.14 At the 10' life level the

Jan 1, 1969

By D. W. Lillie

It has been demonstrated that considerable bend ductility exists in bulk specimens of polycrystalline high-purity silicon. The possibility of hot-forming at 1200°C is suggested. EXCELLENT corrosion resistance in many media and low cross section for absorption of thermal neutrons (0.13 barn) would make silicon of interest to nuclear engineers were it not for extreme brittle-ness and the difficulty of fabrication by any reasonable means. The use of silicon for structural purposes also has been considered in view of its light weight and oxidation resistance. Johnson and Han-sen' have investigated the properties of silicon-base alloys and concluded that there was no way of making pure silicon or silicon-rich alloys ductile at room temperature. In view of reports of appreciable ductility in germanium single crystals above 550°C'." and some plastic deformation in single-crystal silicon above 900oC,' the present investigation was undertaken to define more precisely the limits of high-temperature ductility in pure silicon. After this investigation was begun torsion ductility in both germanium and silicon was reported by Greiner." Through the courtesy of F. H. Horn, a small bar of cast extra high-purity silicon was obtained and small bend specimens were made from it by careful machining and grinding. All of the reported tests results were obtained from samples from this bar (bar No. 1) and one other of similar source (bar No. 2). No complete analysis was obtained but, based on analysis of similar semi-conductor grade material, metallic impurities were under 0.01 pct total. Vacuum-fusion analysis for oxygen showed a value of 0.0018 2 0.0003 pct for the first bar tested and metallographic analysis showed no evidence of a second phase. Bend tests were carried out on an Instron tensile machine using a bend fixture with a 1 -in. span loaded at the center. Supporting and loading bars were 0.250 in. round and the load was applied by downward motion of the pulling crosshead of the machine. Specimen thickness and width were approximately 0.10 in. and % in. respectively. Loading rate was controlled by holding crosshead motion constant at 0.02 ipm. In some cases a smaller specimen was used on a 5/8-in. span with a 0.129-in.-diam loading bar. The entire bend fixture was surrounded by a hinged furnace and all heating was done in air atmosphere. Temperature measurement was made with thermocouples fastened directly to the bend fixture within less than 1 in. from the specimen. Autographic stress-strain curves were recorded during each test, and breaking load, total deflection, and plastic strain could be obtained from these curves. Stress was calculated from the beam formula S = 3PL/2bh2, where P is the load in pounds, L the span in inches, b the specimen width in inches, and h the specimen thickness in inches. This formula is strictly correct only in the elastic range but has been used to calculate a nominal stress for convenience in the plastic range. The stress given is the maximum stress in the specimen. Results The results of the complete series of tests are shown in Table I. The first group of tests (specimens Nos. 1-6) showed the beginning of plastic flow at a test temperature of 900°C, so two additional tests (Nos. 8 and 9) were made at 950°C on small-size specimens from bar No. 2. Specimen No. 8 was tested in the as-machined condition, and No. 9 was heat-treated in hydrogen at 1300°C for 2 hr, cooled to 1200°C and held 1 hr, cooled to 1000°C and held 1 hr, cooled to 900°C and held 1 hr, and finally cooled to a low temperature before removal from the hydrogen. It is apparent that the heat-treatment had a significant effect on yield strength and ductility. In addition, the magnitude of the yield point was conslderably reduced in the heat-treated specimen as is shown m Fig. 1 by tracings of the stress-strain curves. After obtaining a furnace capable of reaching higher temperatures specimens Nos. 10 to 13 were tested at 1100 and 1200°C. Strain rate was increased by up to a factor of 10 to see whether the ductility observed was excessively strain sensitive. Specimen NO. 10, strained at 0.02 ipm and 1100oC, was still bending at a deflection of 0.322 in. when the load rate was increased to 0.2 ipm, resulting in immediate

Jan 1, 1959

By David W. James, Howard Jones, George M. Leak

Conditions are described for producing, by primary recrystallization, a matrix suitable for the growth of large grains near (110)[001] orientation in silicon iron strip by secondary recrystallizaliun in a steep temperature gradient. The orientation distribution of these large grains is expressed in terms of rotational deviations about the cross-rolling direction, the rolling direction, and the normal to the sheet, the deviational spread increasing in that order. With the aid of cowplenientary published data on the orientation dependence of growth rate, it is shown that this observation is consistent with the oriented-growth theory of recrystallization lextures. It is conclutled that growth-rate and orientation-distribution data obtained in a steep thermal gradient should be used with caution to account for isothermally Produced recrystallization textures. SEVERAL authors have reported methods of growing large grains by re crystallization of a small-grained matrix in silicon iron 1- B and pure a cr The present study was a preliminary in the growth of single crystals and bicrystals for surface relaxation," grain boundary mobility, and grain boundary diffusion studies. The method was to control the growth of a seed crystal into a suitable primary re crystallized matrix by feeding through a steep temperature gradient. The driving energy for growth derived from the grain boundary energy released as the seed crystals grew into the matrix. Thus, stability of the matrix against normal grain growth was considered to be essential for success. It was known that the manganese sulfide dispersion present in commercial silicon iron performs this function during secondary recrystallization to the (110)[001.] texture.12 Hence commercial, rather than high-purity, material was used throughout. The paper describes the growth conditions for grains large enough to be used as seed crystals for further growth into single crystals. The orientation distribution of the seed crystals is analyzed and its significance for the theory of recrystallization textures is discussed. EXPERIMENTAL PROCEDURE Strip material was supplied by the Steel Co. of Wales, Ltd. The chemical analysis in weight percent was Si, 2.90; C, 0.015; Mn, 0.059; P, 0.011; S, 0.027; Ni, 0.032; 0, 0.009; Fe, balance. A gradient furnace of similar design to one described previously4 was loaned from B.I.S.R.A. It consisted essentially of a vertical water-cooled copper slot projecting downwards into the hot zone of a molybdenum furnace. Hydrogen was passed through the furnace to protect both heating element and specimen from oxidation. Strip specimens up to 8 cm wide and 0.2 cm thick were sealed into the furnace at the mouth of the copper slot. A coating of light oil on the strip surface maintained the seal during translation of a specimen. The maximum temperature gradient in the region just below the copper slot was 500°C per cm over 1 cm, with the hottest point controlled at 1175°C. Several large grains would usually grow by secondary recrystallization from the primary matrix when a specimen was immersed in the hot zone for about 30 min. A back-reflection X-ray camera was constructed to facilitate rapid and accurate orientation determinations of the large grains produced. It was possible to reproduce a standard geometry, with regard to strip and camera, without the tedium of careful alignment on each occasion. Specimens, typically 4 cm wide and 75 cm long, were cut with the longitudinal axis parallel to the rolling direction of the original strip. The surfaces were cleaned by immersion alternately in a hot aqueous solution containing 2 pct hydrofluoric acid plus 10 pct sulfuric acid and in cold 10 pct nitric acid. The nitric acid etch was just sufficient to reveal the grain structure. Rolling and annealing treatments to prepare the matrix (discussed below) were followed by growth of seed crystals in the gradient furnace. The matrix was transformed to a single crystal by growth of a selected seed crystal connected to the matrix by a thin neck. 4,5 Growth was promoted by controlled feeding into the gradient furnace. Several single crystals of controlled orientation were grown successfully from seed crystals by twisting the interconnecting neck in a reorien-tation jig.4 EXPERIMENTAL RESULTS AND DISCUSSION Growth Conditions. A suitable matrix for growth of large grains was prepared starting from primary re-crystallized strip 1.9 mm thick. This was cold-rolled in two stages each being followed by a recrystallization anneal at 800°C for a few minutes. Such treatment gave the required growth matrix only if the two cold-reduction stages were each performed in several passes and in the following ranges: the first, 30 to 70 pct; the second, 10 to 50 pct. Immersion in the temperature gradient otherwise resulted in an equiaxed

Jan 1, 1967