Search Documents

By J. P. Hirth, G. M. Pound

The flux of evaporation J from {loo) , {111), and (110) planes of silver single crystals and from poly crystalline specimens was measured at 744°C under a high vacuum using a time average weight loss method. The evaporation coefficient a was computed from the measured flux and vapor pressure data through the relationship J = aPe/21rmkT. The evaporation coefficients for the close packed planes and the poly crystalline specimens were found to be about 0.4 providing tentative confirmation of the theoretical limiting law for evaporation of metal crystals, a = 73(P/Pe) + % . The (110) planes faceted to form (111) planes, resulting in a greater evaporating area and larger a . The evaporation morphology of the plane surfaces and the surfaces of spherical single crystals is described in some detail. RECENTLY a detailed theoretical analysis of the mechanism and kinetics of metal crystal evaporation was made.1'2 A mechanism for vaporization involving the recession of monatomic ledges across a metal surface was deduced, and evaporation rate equations were derived by considering the dynamics of these ledges. In most metals, the process of evaporation is one of dissociation of atoms from ledges, surface diffusion, and desorption. The gross vaporization rate is controlled by the surface diffusion and desorption rates. It was shown that the monatomic ledges, under the action of surface diffusion gradients, tend to assume a terminal velocity and spacing. Further, it was shown that crystal edges and grain boundaries are a ready source for monatomic ledges. In the case of evaporation of large, perfect crystals with clean, low-index surfaces, the kinetic treatment1 leads to a limiting law for crystal evaporation: where J is the gross flux of vaporization, a is the evaporation coefficient, Pe the equilibrium vapor pressure, m the atomic mass, k Boltzmann's constant, T the absolute temperature, and P the actual pressure of vapor over the crystal surface. Imperfections such as grain boundaries, pores, cracks, and macroscopic steps should have the effect of causing positive departures from this limiting law. Such an imperfection is a source for monatomic ledges, which, near their origin, are closely spaced and thus give rise to higher evaporation rates there. However, under many experimental conditions,' this effect should be negligible due to the relatively low concentration of imperfections. In general, a more important cause of positive departure should be screw dislocations, and the effect of these has been considered in some detai1.2,3 Monatomic ledges emanate from points where the surface is intersected by a dislocation with a screw component normal to the surface. The dynamics of ledges arising from a screw dislocation are similar to those of ledges arising from a crystal edge. There is a perturbed region near the origin where the monatomic ledges are, under certain experimental conditions, closely spaced and the

Jan 1, 1960

By M. J. Marcinkowski, A. S. Sastri

A detailed compressive stress-strain analysis and transmission electron microscopy investigation has been made of the deformation behavior occurring in a 50 at. pct Ni-Ti (hypoeutectoid) alloy and a 54.5 at. pct Ni-Ti (hypereutectoid) alloy. In the case of the hypoeutectoid alloy, three stages of work hardening are observed. Stage I occurs at a very low stress and is associated with plastic deformation via martensite formation. Stage 11 is characterized by very rapid work hardening and is due to difficulties in causing further deformation in the fine martensite aggregate produced in Stage I. Stage III which occurs at very high stress levels is characterized by smaller work hardening rates and is due to the plastic deformation arising from alternate reconversions of the original martensites to martensites of varying orientation. Rapid quenching of the hypereutectoid alloy leads to very high yield strengths and is related to a fine precipitate dispersion that such treatment brings about. The present investigation represents the final phase of a three-part study directed toward an understanding of the solid-state transformations in near equi-atomic Ni-Ti alloys as well as the deformation mechanisms associated with these alloys. In the first part,"2 to be henceforth referred to as I, it was found that alternate simple shears on {112} planes and in (111) directions convert the parent B2 structure in the equiatomic NiTi alloy into two distinct close-packed monoclinic martensites. All of the marten-sites were of this type, whether they were formed by cooling or by plastic deformation, whether induced to form in bulk samples or in thin foils, or whether examined in the electron microscope at room temperature or below. On the other hand, in the second part of this investigation,3 to be reffered to as 11, it was shown that upon slow cooling to about 640°C. alloys in the neighborhood of NiTi which possess the B2 structure transform eutectoidally into their equilibrium phases Ti2Ni and TiNi3. However, preceding the formation of these equilibrium phases a series of metastable intermediate phases are formed. This paper will set as its goal the elucidation of the remarkable deformation behavior exhibited by NiTi. In particular, Buehler and Wiley4 have found equiatomic NiTi to be surprisingly soft, while Buehler et al.5 have shown this alloy to possess a memory effect: i.e., upon bending at room temperature it will revert to its original shape when heated to above about 50°C. In I it was shown that NiTi was soft in the sense that the yield stress was low; nevertheless, the alloy work-hardened at an extremely rapid rate to very high stress levels. On the other hand, the hypereutectoid alloys with somewhat higher nickel, say 54.5 at. pct (60 wt pct) have enormously increased yield strengths compared to those of the equiatomic alloys. In order to determine the atomistic processes giving rise to the above behavior, it was decided to examine samples that were wafered from bulk specimens deformed in compression to various strains using the techniques of transmission electron microscopy. EXPERIMENTAL TECHNIQUE All of the alloys used in the present investigation contained either 50 at. pct Ni (55.06 wt pct) or 54.5 at. pct Ni (60 wt pct) and were arc-melted in the form of a finger using the same techniques described in I and II. The finger was capsulated in a stainless-steel jacket and swaged at 850°C into rods. Compression specimens 0.300 in, long and 0.200 in. in diam were machined from these rods. In order to completely re-crystallize the samples and remove residual stresses, all of them were capsulated in evacuated quartz, annealed for 1/2 at 1050°C. and then furnace-cooled. Compression tests were carried out in an Instron tensile testing machine covering a range of temperatures from —196° to 200°C using procedures described previously.6'7 In all cases crosshead speed was 0.02 in. per min. Wafers 0.015 in. thick were spark-cut from the cylindrical samples at 45 deg to the compression axes after they had been deformed to the desired strain. These specimens were then spark-planed to about 0.005 in. and then electrochemically thinned for examination by transmission electron microscopy as described in I.

Jan 1, 1969

By J. S. Levine, M. Prats

Several methods are available for calculating the performance of solution-gas-drive reservoirs from the PVT properties of the oil and from the relative permeability and other properties of the formation. These methods require a number of simplifying assumptions. The present method of computation has made use of a high-speed computer to solve simultaneously the nonlinear partial differential equations that describe two-phase flow by solution-gas drive in order to calculate the performance of a reservoir. Some of the results obtained by the nonlinear partial differential equation solution are compared with those obtained with an approximate method, which has been called the semisteady-state solution. The pressure and saturation profiles from the wellbore to outer boundary calculated by the two methods are compared for one constant-terminal-rate case and two constant-teminal-pressure cases. The agreement in these profiles, as well as in the values of average reservoir pressure and cumulative recovery, leads to the conclusion that, for most engineering calculations, the semisteady-state method will give a reasonable approximation to the numerical solution of the differential equations describing solution-gas drive. An unfavorable (as regards ultimate oil production) set of relative permeability curves was used in the calculations in the belief that the effect of the parameters which were studied would be emphasized to a greater degree. Furthermore, the reservoir was assumed to be completely homogeneous, and these results should not be considered applicable to any other type of reservoir. Gravity effects are not considered. The absolute permeability was varied from 25 to 0.5 md. At an economic limit of 2 B/D, the recovery for a 25-md reservoir is about 1.8 times as great as that for a 0.5-md reservoir. The effect of permeability on the producing gas-oil ratio is minor. Once PVT properties of the oil and the relative penneability properties of the reservoir are fixed, the pro- during gas-oil ratio is found to be a function of the fraction of oil recovered. Well spacings of 10, 40 and 80 acres were considered. For the assumed homogeneous-reservoir properties, the effect of spacing on recovery at an economic limit of 2 B/D was very slight. Certain dimensionless groups can be used to extend the results to other fields having different penneabilities, spacings, reservoir thicknesses, well radii and porosities, so long as the PVT and relative permeability properties are similar to those used in this paper. INTRODUCTION An important aspect of reservoir engineering is the prediction of the performance of the reservoir based on the limited information normally available early in the life of a field. Usually, soon after a field has been discovered, it is necessary to decide upon the spacing to be employed and the production program to be used for most efficient utilization of reservoir energy. Where these decisions have not been dictated by legal or political considerations, they have frequently been based on experience factors which indicate that two fields with similar characteristics will probably have similar performances. To the extent that the industry has been successful in economically exploiting oil reservoirs, this rule-of-thumb method has been found to have merit. A large amount of theory has been developed with which it is possible to predict the performance of a reservoir by using certain known properties of the oil and the formation. Because of the difficulties of the mathematics involved for a drainage area having square or rectangular boundaries, the problem usually has been simplified and idealized by assuming each drainage area to have radial symmetry and the field to be represented by a number of these drainage areas. Muskatl was one of the first to formulate the theory for two-phase flow and to solve the equations for a few cases by numerical integration, using the relatively slow means of computing available at that time. His equations describing

By W. C. Hauber

Techniques previously published to predict the production performance of a water flood when the mobility ratio is not unity have been primarily restricted to a five-spot well pattern. When other types of well patterns are considered, the prediction of performance has required either a model study or a complex numerical solution of a multiple-fluid system with a high-speed digital computer. Simple methods are developed in this report for predicting the areal sweep efficiency and injectivity history of a water flood for arbitrary well patterns and mobility ratios. The effect of an initial gas saturation can be incorporated into this method of prediction. With the layer superposition concepts published by Prats, et al.,1 these methods can be used to predict the production performance of a water flood in a reservoir having a wide range in permeability. A reasonably good agreement is shown by comparisons of the results calculated by this method with those obtained either by model studies or by complex numerical solution of multiple-fluid systems. Therefore, it has been concluded that a reasonably accurate method has been developed for predicting the production performance of a water flood employing an arbitrary type of well pattern for any given mobility ratio. A medium-speed digital computer can be used in applying this method. INTRODUCTION Production performance of multi-fluid five-spot water floods can be predicted by the techniques published by Prats, et al.1 This technique takes into account the effects of initial gas saturation, mobility ratio and vertical variations in horizontal permeability. A basic assumption in applying this technique is that unit displacement efficiency does not vary with time after the passage of the flood front. Prats, et al. demonstrated that this assumption is valid for normal waterflood applications with small mobility ratios. To extend Prats' technique to other types of well patterns has required either an experimental model study or a numerical solution of a multiple-fluid system to describe the performance of an individual homogeneous layer. Sheldon and Dougherty2 presented an exact numerical solution employing any arbitrary well pattern. This technique is complex to apply, requiring (1) numerical solution of the single-fluid system, (2) conformal mapping of the streamlines and iso-potential lines to form a new grid network, and (3) a numerical solution of the multi-fluid system on the new grid network. A method was developed in this paper to approximate areal sweep efficiency and injectivity as functions of time for an individual homogeneous layer for arbitrary well patterns and mobility ratios. Analytic expressions were derived to directly apply this method to five-spot and direct line drive well patterns. For other types of well patterns, this method is similar to the approach of Higgids and Leighton,6 in that it utilizes the stream tubes and iso-potential lines obtained by the solution of differential equations when the mobility ratio is unity. However, this method is easier to apply than the Higgins-Leighton method since it does not require knowledge of the variation between oil and water relative permeabilities and saturation. Using the relationship between areal sweep efficiency, injectivity and the volume of water injected into an individual homogeneous layer as obtained by the method presented in this paper, and Prats' technique for the superposition of individual layers, one can predict the production performance of any water flood. ASSUMPTIONS GENERAL Basic assumptions made in this report are that fluids are incompressible and the layers are horizontal, homogeneous and uniformly thick. It is assumed that as water is injected into the layer, three distinct regions are formed —a water bank, an oil bank, and an unflooded region. Each region is assumed to displace other regions in a pis-tonlike manner, with only one fluid flowing in each region, so that there is no change in unit displacement efficiency with time after the flood front has passed in the reservoir. Therefore, the residual oil saturation in the water bank and the oil saturations in the unflooded region are constant throughout the layer. The ratio (F) of the bulk volume of a layer occupied by the oil and water banks at any time before oil breakthrough, to the corresponding bulk volume of a layer occupied by the water bank, is given' by the equation BREAKTHROUGH AREAL SWEEP EFFICIENCY To predict the behavior of any type of waterflood pattern, one must know the breakthrough areal sweep efficiency Eb(Mc,o, F) for the pattern for any water:oil mo-

Jan 1, 1965

By D. B. Hunter, H. W. Rosenberg

The titanium-rich portion of the Ti-Pd system was investigated from 0 to 75 wt pct Pd by metallo-graphic and X-ray techniques. A 0 eutectoid occurs at 24 wt pct Pd and 1190°F. Two compoutzds are indicated in the region below 75 wt pct Pd, Ti,Pd and Ti2Pd3. The solubility of palladium its a titanium is low, probably less than 1 pct. In 1960 Rudnitskii and Birunl published a complete version of the Ti-Pd phase diagram. However, their work was in disagreement with the earlier literature in that they reported only one inter metallic compound, whereas three had been reported earlier. In view of these discrepancies, it was therefore necessary to redetermine those portions of the diagram of immediate interest. The following account describes our work on the system over the range of 0 to 75 wt pct Pd. MATERIALS AND METHODS Distilled titanium sponge and elemental palladium were used in the formulation of the alloys; the chemistry of these materials is detailed in Table I. The alloys were prepared as 10 to 50 g blended compacts that were melted into buttons by arc melting under gettered argon on a water-cooled copper hearth. Weighing of the ingredients before and after melting showed that negligible weight changes occurred. Therefore, no analyses were undertaken and the compositions of all alloys are nominal. All alloys were fabricated by hot rolling at 1700°F to 0.070-in.-thick sheet. Scale was removed by sandblasting and pickling in a 5 pct HF-35 pct HNO,, balance H20 solution. For metallographic examination, specimens were mounted after heat treatment transverse to the rolling direction, ground on silicon carbide papers of increasing fineness to 600 grit, and then electro-polished using a solution containing 600 ml me-thanol, 60 ml perchloric acid, 360 ml butyl cello-solve, and 2 ml "Solvent X". Unless otherwise specified, etching of alloys containing up to 42 pct Pd was carried out by swabbing with a 12 pct HN03-1 pct HF aqueous etch where a bright etch was required, or by a 1 pct hydrofluoric in saturated oxalic acid solution where contrast between phases was required. The Ti-52.8 Pd alloy was etched with a solution of 25 ml HF, 40 ml glycerine, 35 ml methanol, and 18 g benzalkonium chloride. For X-ray examination, 1/2-in.-square speci- mens of sheet were mounted flat in a standard 1-in. metallographic mount and ground and polished as above. X-ray diffraction was performed using a Norelco type 12045 Diffractometer, employing CuKa radiation with a nickel filter at 40 kv and 20 ma. Specimens were rotated about the sheet normal during exposure. Although this procedure did not remove the effects of sheet texture from the relative intensities, it had the advantage that oxidation or contaminants entering during preparation of powder samples could not confuse the patterns obtained. RESULTS AND DISCUSSION Fig. 1 illustrates the Ti-Pd phase diagram according to Rudnitskii and Birunl with the work of the present authors superimposed. Both interpretations agree that the system is of the 0 -eutectoid type with an extensive 0-phase field, and that the eutectoid temperature is just below 1200°F. There is also agreement that the solubility of palladium in a titanium is restricted. Our work would indicate that the a solubility of palladium is low, probably less than 1 pct. However, whereas Rudnitskii and Birunl place the eutectoid composition at 41 wt pct Pd, this investigation shows it to be at about 24 wt pct Pd. Moreover, this investigation confirms the existence of compounds at Ti2Pd and Ti2Pd3, whereas Rudnitskii and Birun report only a single Berthol-lide phase covering the TiPd to TiPd, range. Laves et a1.' and Wallbaum, whose work was summarized by Maykuth , reported the existence of Ti2Pd3 and TiPd, in addition to Ti2Pd. More recently, Nevitt and Downe~' have reported the structure of

Jan 1, 1965

By A. R. Brown, D. Silverman

EquipmENT and methods for locating the points of entry of salt water into oil wells are described. These techniques make it possible to delineate accurately the top and bottom boundaries of water zones whether bottom water or intermediate zones of production, Data are given on tests by this survey method on a number of producing wells, Complete results were obtained on wells shelving water percentages ranging from high values to as Little as 10 pct salt water. A number of the wells tested were plugged back and the results are reported in detail. While not all the work-overs were successful, evidence, including the results of reruns by this survey method, points to failure of the plugback operations. The data provided in the tests discussed show, contrary to general opinion, that a large percentage approaching 50 pct of the wells showed intermediate water instead of bottom water. Furthermore, approximately a third of the wells tested showed multiple zones of water production. The equipment comprises a long pipelike assembly carrying 10 pairs of electrodes spaced apart a distance of 4 ft. A surface-controlled solenoid-operated switch is provided to connect each of the individual pairs of electrodes in turn to a single-conductor steel-armored cable by means of which the conductivity of the fluid in the vicinity of each of the electrode pairs can be continuously recorded at the surface. Any number of electrode pairs and any spacing can be used. This assembly is lowered on tubing and placed in position opposite the section of the well to be studied. A conventional rod pump is used to produce the well for the survey. The pump inlet is placed below the lowest electrode pair. The well is conditioned by pumping Out the salt water normally standing in the well and introducing fresh water into the annulus at the surface. When the resistivity of the water standing in the well has been increased to a sufficiently high value, the inflowing fresh water is stopped and the drawdown operation of the survey is started. As the fluid head is lowered by pumping, fluids flow into the well and the entry Of salt water is indicated by the change in resistivity of the fresh water opposite the 'Ones Of water entry, The process can be repeated as Often as necessary with little effort and without moving any well equipment to be sure that the fluids produced from each zone are representative Of the true reservoir content. Introduction The oil industry has been concerned from its inception with the control of fluids entering its wells. In order to effect the greatest economy of operation and to get the greatest total production of petroleum liquids out of the ground, the flow of gas and water must be controlled and often extensive workovers are planned and carried Out for this purpose. This has become even more important in recent Years because of the need for the conservation of natural resources. Whereas producing operations concerned with the original completion of a well are assisted to a great extent by a wide variety of instruments developed and perfected for that purpose, workover operations on oil wells have been greatly hindered by a lack of suitable instrumentation. The industry has searched for many years for a satisfactory instrument and process that

Jan 1, 1948

By M. R. Tek

A new method for correlating the data on multiphase flow through vertical pipe is presented. The correlation is based on a "two-phase f factor" concept which was developed and successfully applied to horizontal multiphase flow by Bertuzzi, Tek and Poettmann. Field data previously published on several flowing and gas-lift wells have been used as the basis for the developed correlation. The application of the method to actual two-phase flow problems indicate that this method is capable of predicting the pressure distribution in vertical multiphase flow strings well within the accuracy range usually desired for common engineering and design calculations. A new working chart developed for calculation of two-phase pressure gradients and a graphical step-by-step procedure for the computation of pressure distribution are presented along with an example problem. INTRODUCTION Multiphase flow through vertical pipe is encountered in many engineering installations. In petroleum, chemical process, nuclear engineering and many other industries, problems associated with simultaneous flow of two or more phases through vertical pipe have been of interest for a long time. This interest has increased considerably during recent years due to applications to new processes in petroleum production and refining and to problems of steam generation and heat removal from nuclear reactors. One prominent example of vertical two-phase flow is provided by the gas-lift process where oil, water and gas flow simultaneously. If the pressure profile in a gas-lift well can be predicted within reasonable accuracy, it would be possible to get good estimates of the power required to lift the oil, the optimum depth, the pressure and the rate at which to inject gas. Furthermore, the effect of production rate and tubing sizes on these quantities can be evaluated before any design decision is made on the installation and operation of the flow string. The majority of experimental work available in the literature deals with two-component systems where individual phase flow rates in and out of the pipe remain constant. The general problem of prediction of pressure drops in multiphase flow systems is very complicated. The co-existence of numerous flow patterns of widely different geometry and mechanism, conditions of surface instability and the nature of force fields acting upon the system are among the major difficulties commonly encountered. The classical approach of fluid dynamics which would be based upon the formulation and solution of Navier-Stokes equations has been found by many investigators completely devoid of any hope of success-—not only because of inherent nonlinearities but also because of insurmountable analytical difficulties standirng in the way of setting up the boundary conditions. The presence and effect of interfacial forces on mu1tiphase flow systems further complicate the theoretical approach. For these reasons, many investigators choose to adopt semi-empirical if not purely empirical approaches in order to obtain solutions of engineering utility. A prominent practical solution has been given by Poettmann and Carpenter' in the form of an empirical correlation. In their paper, total flowing densities of fluids and solubility effects of gas in oil have been taken into account for the correlation of field data covering a wide range of operating conditions for oil wells. They treated the gas. oil and water as a single phase of combined propcrties and correlated the multiphase friction factor as a function of the product (pvd), i.e., density X velocity X diameter of the flow string. However, because the product (pvd) is dimensional while the frictior factor is dimensionless, the generality of their result's is somewhat restricted. It seems that the omission 01 viscosity effect may be one of the reasons for the scatter of data as shown in their correlation chart. The concepts of the "two-phase f factor" and the "two-phase Reynolds' number function" were recently developed and successfully applied to correlate horizontal multphase flow by Bertuzzi, Tek and Poettmann.' Recently, two new methods of correlation by Tek and Chan' have been presented on simultaneous flow of liquid and gas through vertical pipe. These two methods of correlation and the working equations in field units necessary for their application are included in this paper. The extension of the conccpt of two-phasc Reynolds' number function successfully developed for horizontal flow into vertical multiphase flow systems, development and evaluation of working charts permitting calculation of two-phase pressure gradients and

By J. P. Hirth, A. H. Clauer

A simple graphical method is presented for the determination of climb forces on dislocations exerted by uniaxial stresses. In conjunction with standard stereographic projections, the technique is applicable for any crystal system. Climb forces on glide dislocations in bcc crystals are determined as an application. In connection with a study1 of elevated temperature creep of molybdenum single crystals, a detailed analysis of the climb force on a dislocation was necessary. In general, such forces can be calculated for an arbitrary dislocation segment from the Peach-Koehler equation2 (or one of its variants3"6) where Fk are the components of the force per unit length of dislocation line, ?i of the unit sense vector tangent to the line, ojl of the stress tensor, bl of the dislocation Burgers vector, and is the Einstein permutation operator with components ?123 = ?231 = ?321 = —?132 = — ?321 = —?—213 = unity, and all other components zero. It is understood in the suffix notation of Eq. [I.] that repeated indices are to be summed over 1, 2, and 3. For the most general coordinate system, the matrix manipulations required to solve Eq. [I] are fairly complex. However, a simplification can be achieved by specifying a coordinate system based on the dislocation in question. The pertinent system involves cartesian coordinates xi, x;, xi, with attendant unit vectors e'i; ei being parallel to the Burgers vector b, and e'3; being parallel to the climb directions (b x ?), i.e., normal to the slip plane. In this coordnate system b = (bi, 0, 0) and ? = (?'1,?'2, 0). Hence, Eq. [I] yields for the climb force per unit length F3 = ?'2s'11b'1—?'1s'12b'1 [2] The first term on the right side of Eq. [2] is the climb force on the edge component of the dislocation, while the second is the (cross-slip) force normal to the glide plane acting on the screw component. Thus, the problem of determining the climb force reduces to that of resolving the stress components —il and s'12 in the xi coordinate system. Several authors7- l0 have shown, for the glide case, that such stress resolution can be performed most readily by graphical methods. Therefore, we present a similar method for the climb case. Of the suggested methods, we follow that of Hartley and Hirth,9 which can be performed directly with a standard stereographic projection for the crystal system in question. We treat the case of simple tension, or compression, this being the most common creep test in which climb is important. A set of cartesian coordinates (x1, x2, x3) is fixed with x3 parallel to the specimen axis. This set is connected to the x'i set by the factors defined in Table I. A perspective view of the coordinate systems relating them to the glide ellipse, is given in Fig. 1, while a stereographic projection of the coordinate systems is presented in Fig. 2. The transformation matrix connecting the xi and xi coordinates is x'i =aijXj [3] where an a11 a13 aij = a21 a22 a23 a31 a23 [4] s The corresponding transformation for the stresses is i [5] The only nonzero component of okl for the simple tension case is s33,. Hence, using Eqs. [4] and 151, we find that Eq. [2] becomes F3 = b'1s33 sin2 cos ?[cos ? sin a - sin ? cos a] [6] = m1b'1s033 sin a- rn2b's33 cos a where sin a = 5: and cos a = ?i, with a the angle measured from b to ?, i.e., a polar angle in the x1, x2 plane. The same equation with a sign change applies for compression. A graphical plot of the "Schmid" factor m1, relating to the climb force on the edge component, is presented in the stereographic projection of Fig. 3. In this projection the glide plane normal points toward the top of the page and the glide direction points normal to the plane of the page. To use the graph, one

Jan 1, 1970

By D. E. Peacock, B. Harris

The mechanical behavior of a niobium (columbium)TUNG alloy containing 20 wt pet Ta. 15 wt pet W, and 5 wt pct Mo has been studied in the temperature range 77° to 423°K. All specitrzens tested, apart from those into which oxygen was introduced deliberately, were fully recrystallized at 1500°C, had an average grain diameter of approximately 0.03 mm, and a combined carbon, ox-vgen. nitrogen, and hydrogen content of not not more than 250 ppm The meckatnical behavior of this alloy was found to be similar in seine respects to that of unalloyed bcc transilion metals. Both yield and flow stresses were sensililae to changes in deformation rate and both increased rapidly as the temperature was lowered below 150°C. A detailed examination revealed, however, that most of the temperature dependencc of the yield stress could be accounted for by the temperature dependence of the elaslic, modulus, whish is not true of the unualloyed metals. In lension at moderate strain rates, the trunsition fvom ducfile to hvittle behavior occurred just belour room teunperature bur incveasing the oxygen content from 190 to 200 pptn raised the transition temperature to about 4000K. TUNGSTEN and molybdenum have been found to be very effective high-temperature solid-solution strengtheners of niobium, but the amounts of these elements that can be added are limited if the alloy is to retain useful room-temperature ductility. It has been shown, however, that partial replacement of niobium with tantalum lowers the ductile-to-brittle transition temperature of such alloys and allows larger additions of molybdenum and tungsten to be made. Some of the high-temperature properties of niobium alloys containing these elements have been reported by Bartlett et al.,' and the work presented here describes the low-temperature behavior of one of these, a solid-solution alloy containing, by weight, 20 pct Ta, 15 pct W, and 5 pct Mo. The elastic modulus, the temperature and strain-rate dependence of the yield and ultimate tensile stresses, and strain-hardening effects are discussed and compared with those of unalloyed niobium. EXPERIMENTAL MATERIAL AND PROCEDURE The material used in this investigation was produced from a 6-in.-diameter, double arc-melted ingot. This was can-extruded at about 1150°C to a 3-in.-diameter billet, and after annealing for 1 hr at 1345°C the billet was heated to 1315OC and part rolled and part forged to 0.5-in.-thick plate. The resulting material exhibited a cold-worked structure

Jan 1, 1965

By O. D. Sherby

In 1938, Kanter' revealed that the steady-state creep rate of low-carbon iron alloys could be correlated by an activation energy expression, where the activation energy for creep, Qc, was found equal to about 90.0 kcal per mole for most of the alloys investigated. He suggested that creep at high temperatures might be controlled by some sort of solid self-diffusion process since the value of Qc was about equal to the known value for self-diffusion, Qd. Correlations were not further attempted along these lines until many years later, when new evidence2-8 was presented suggesting that Qc = Qd for many metals. This equality now appears to be established for ten metals (Al, Pb, Ni, Mg, Cd, In, a Fe, y Fe, Cu, and Au). This experimental fact prompted several investigators to postulate various mechanisms of dislocation motion which would yield an activation energy for creep equal to that for self-diffusion.7-9 Possibly the most accepted theory is that based on dislocation climb7 where a dislocation is envisioned to surmount a barrier by an atom-vacancy exchange process, thereby "climbing" from one slip plane to another. Such a model would suggest that the steady-state creep rate, i, should be directly proportional to the rate of self-diffusion, D, for a given creep stress. In order to more clearly reveal the relationship between creep and diffusion, and to eventually determine the correct model, it would be desirable to obtain data showing correlations between i and D, rather than with a comparison of the temperature coefficients, Q. It has been known for some time that the self-diffusivity of iron atoms is greatly decreased when body-centered-cubic (a) iron transforms to face-centered-cubic (y) iron. This has been ascribed to the more close-packed structure of y iron. The diffusivity ratio at the transition temperature is Da/Dy =350. Correlations of creep data for a and y iron5 reveal that the creep rates under a given stress at the transition temperature differ by about the same ratio (ea/ey 2 200). It would be difficult to reconcile these results on any other basis than that the diffusivity of atoms probably controls creep of iron in some very direct way. No other published data on other allotropic metals appear currently available for further checking such a possible correlation between i and D. Thallium transforms from a hexagonal close-packed (a) structure to a body-centered-cubic (ß) structure at 230°C (503°K) and the self-diffusivities of the two structures have been determined by shirn.10 His results are plotted in Fig. 1 where it will be noted that Dß/Da = 25 at the transition temperature. It would appear desirable, therefore, to determine the creep rates of polycrystalline a and ß thallium at the transition temperature, as well as determine

Jan 1, 1959

By P. S. Rudman

Crystal-structure occurrences in binary phase diagrams of fifth- and sixth-period elements, zirconium through palladium and hafniun through platinum, are surveyed with the conclusion that average .group number GN is the dominant correlative. This is interpreted in terms of curves of enthalpy differences between competitiue structures that are continuous func -tions of GN. It is further reasoned that the absolute formation energies may also be continuous functions of GN and if so can be simply obtained from the regatives of the heats of sublimation of the elements. It is shown that from such a heat of formation-GN plot the heats of mixing of alloy systems can be quantitatively predicted. For example, for sixth-period alloys very large exothermic heats of mixing are predicted for systems that span GN = 6 (i.e., Hf-Pt, Ta-Os, ...). Thermodynamic data with which to test these predictions do not appear to be available. If in the latter transition (A subgroup) elements we eliminate from consideration the fourth period where magnetic terms contribute to crystal-structure stability, as shown in Table I, there is a very regular sequence of structures. Using the group number (GN) assignments of Table I, which differ from the usual assignments of 8 to the elements listed as 8 through 10, it is perhaps trivial to say that the phase occurrence is a unique function of GN. However, the interesting aspect of alloys from among these elements is that the alloy structures also appear to correlate with GN. The GN for an alloy is given as GN = SixiGNi, where xi and GNi = atomic fraction and GN of component i, respectively. It is difficult to credit anyone in particular with having discovered this correlation, so many people have implicitly or explicitly employed it in the last 10 years. However, Bloom and Grant' in their correlation of s-phase occurrence with total number of d + s electrons were certainly one of the first to appreciate the correlation and Hume -Rothery2 appears to be responsible for the use of average GN as defined above. Only in the last few years have the phase diagrams of enough alloy systems been determined for this correlation to be surveyed with any degree of completeness. We shall here, in an abbreviated format, summarize the available phase-diagram literature and discuss some consequences of the GN correlation. BINARY DIAGRAM SURVEY In Figs. 1 to 3 the phase occurrences are summarized. In Fig. 4 the occurrence limits for the more common structures are summarized. That the GN of an alloy falls within the occurrence limits of a particular structure is a necessary but not a sufficient condition for the appearance of that structure. The preceding figures certainly illustrate that GN is a dominant correlative for structures of these elements. Many of the diagrams have complications, with ordering and other phase changes occurring in intermediate phases, that cannot be clearly accommodated in such a schematic survey and have hence been left out. Such omissions cannot, however, negate the GN correlations discussed here. Relative atomic size is also an important crystal-structure determinant and since atomic size varies quite regularly with GN for these elements it is possible that the apparent GN correlation is basically a size correlation. In Fig. 5 relative size occurrence limits are given for intermediate phases and extended terminal solutions. The relative volume, ?v/(v) - 2(vA - vB)/(vA + uB), has been taken as the measure

Jan 1, 1968

By J. N. Anderson

Developments in reverberatory furnace practice at Noranda over the period 1928 to 1953 are described. Features of interest are increasing furnace tonnage from 700 to 2000 tons per furnace day, the use of the suspended basic roof, and improvements in furnace fuel ratio. Points of comparison and contrast between Canadian and American furnace practice are given. OPERATIONS at the Noranda smelter have been described in a series of papers published between 1930 and 1945.1-4 In the present paper the development of reverberatory furnace practice over the 25 year period from the beginning of operations up to the present time will be reviewed and current furnace practice will be described in some detail. Particular attention will be given to those features which are characteristic of Canadian practice: first, the use of the suspended basic roof; second, the deliberate driving of furnaces to give tonnages in the range of 1400 to 2000 tons of solid charge per furnace day; and third, the continuous operation of furnaces at high tonnage for many years without any extended shutdown for repairs or rebuilding. Finally, points of comparison and contrast between reverberatory furnace practice in Canada and the United States will be discussed. The Noranda smelter, which started in December 1927, was built for the reduction of ore from the Horne mine where early development had shown a comparatively small tonnage of rich cop-per-gold ore of a grade and character which would require direct smelting. The smelter was designed for a capacity of 1000 tons of ore a day and originally consisted of eight roasters, two reverberatory furnaces, and two converters, together with necessary auxiliaries. The discovery of a large new body of ore in 1928 made it necessary to increase both the size and the scope of the reduction plant. A mill which was ultimately brought to a capacity of 3000 tons a day was built for concentrating low grade ore and in addition the capacity of the smelting plant had to be greatly expanded, so that for the next fifteen years every effort in the smelter was directed toward achieving maximum tonnage. Because the layout and design of the smelter made it impracticable to add a third furnace, it was necessary to increase the tonnage on the two original furnaces by improvements in furnace design and operation, roasting and converting equipment being added as required to meet the demands of increased furnace capacity. The smelter tonnage, which reached a peak in 1941, was curtailed in 1944 when it became necessary to operate on a one-furnace basis because of a reduction in the output of the Horne mine due to a shortage of labor. Two-furnace operation at a moderate tonnage was resumed in 1949. In the early years the smelter tonnage was made up almost entirely of Noranda ores and concentrates. However, since 1937 an increasing tonnage of ores and concentrates from other mines has been treated in the smelter. In recent years, with decreased production from the Horne mine, copper concentrates from several mines in the district, treated on a custom basis, have formed an important part of the smelter tonnage. A review of developments in reverberatory furnace practice at Noranda falls naturally into the following four periods: 1—1928 to 1935, the early years; 2—1935 to 1944, the period of greatest change and development; 3—1944 to 1949, one-furnace op-

Jan 1, 1955

By Robert E. Hoffman, John B. Hudson

The self-diffusion coefficient of pure lead has been measured at five pressures between atmospheric and 40 kb. over a temperature range of about 150°C near the melting point at each pressure. Measurements were made using a radiotracer technique, samples being sectioned electrochemically after the diffusion anneal to determine the diffusion penetration curve. Pressure was applied in a "belt" type pressure device similar to that used by Hall et al., in the diamond synthesis. The self-diffusion coefficient, D, can be represented by an equation of the form at constant pressure, and an equation of the format constant temperature. An apparent correlation between the effects of pressure on self-diffusion and on the melting point of the material observed by Nachtrieb et al., at pressures below 10 kb. was found to be only approximate as measurements are carried to higher pressures. OVER the course of the past 30 years a large number of studies1-' of the rate of atom movements in solid metals have been carried out. These studies have in general been aimed at increasing the understanding of the detailed atomic process by which atom movements in metals can take place. The largest portion of these studies consisted of measurements of the self-diffusion coefficient, D, as a function of temperature. From such studies con- siderable information regarding the energetics of the motive process has been deduced. Recently such measurements have been extended to the pressure dimension. However, the existing work in this area is scanty and has in the past been limited to the pressures attainable in oil or gas bombs—i.e., up to about 10 kilobars (kb). In the present work, the self-diffusion behavior of pure lead has been investigated using the "belt" type high-pressure apparatus designed by Hall7 over a pressure range from atmospheric to 40 kb at temperatures extending over about a 150°C range close to the melting point at each pressure.

Jan 1, 1962

By P. Somasundaran, D. W. Fuerstenau

The length of the collector is found to influence the flotation of the mineral even at incipient conditions, which are below the concentration at which interaction at the solid-liquid interface begins to take place to form hemi-micelles. To study this dependence, concentration for incipient flotation of quartz was determined as a function of pH with collectors of various chain lengths. The observed effect of chain length on flotation is ascribed to that of collector adsorbed on the bubble surface. In previous studies, it was shown that at low concentrations the alkyl collector ions adsorb at the solid-liquid interface as individuals.''2 At higher concentrations, the collector ions adsorbed at the solid-liquid interface associate with each other to form two-dimensional aggregates called hemi-micelles. Above the hemi-micelle concentration, the length of the hydrocarbon chain is extremely important since the hydrocarbons are in effect removed from water during the association, making the energetic conditions more favorable for adsorption at the interface. Because of this enhanced adsorption, one observes a very rapid increase in flotation associated with the hemi-micelle formation at the solid-liquid interface. However, a dependence of flotation on the chain length at concentrations below that required for hemi-micelle association was also observed,' and this cannot be explained by the above mechanism which postulated hydrocarbon chain interactions only at the solid-liquid interface. This prompted an investigation into other possible reactions of the hydrocarbon chains and an examination of the conditions at the bubble surface involved in the flotation system and how these observations might explain the reactions at the solid-gas interface which cause the particle-bubble attachment required for flotation. To obtain more information on chain length effects, flotation, under incipient conditions, was tested by vacuum flotation techniques. The collector-concen-tration-pH relationships for flotation of quartz with alkyl ammonium acetate collectors was delineated by observing the pH at which quartz particles begin to float to the liquid surface. By investigating flotation as a function of pH, it was also possible to study the effect of neutral molecules on incipient flotation conditions, since the aminium ions hydrolyze to amine molecules at higher pH values. EXPERIMENTAL WORK Brazilian quartz specimens were crushed and sized, and the 270 x 400 mesh fraction was used for flotation studies. The samples were leached with concentrated hydrochloric acid until no coloration of the acid occurred. The leached material was washed free of chloride ions and stored in distilled water. The vacuum flotation technique developed by Schuhmann and prakash3 was used to determine the critical pH-concentration curves. This method, which can be used to delineate conditions for incipient flotation, is fairly simple and rapid. About 0.5 gm of 270 x 400 mesh quartz was placed in a 100 ml graduated cylinder which was then filled to the 100 ml mark with the collector solution made from high-purity alkyl ammonium acetate salts. The water used for the test was conductivity water saturated with air that had been passed through a cleansing train consisting of Drierite, Ascarite, a water wash bottle, and a trap. After the pH was adjusted, the cylinder was then conditioned for thirty minutes. In the tests where an acid pH was desired, sufficient acid was added before the collector solution to avoid any effect due to slow desorption of collector from the quartz surface. After conditioning, vacuum was applied to the system and the flotation or nonflotation of the quartz was noted. The pH at which the quartz particles began to float to the liquid-gas interface was taken as the critical PH. Critical pH curves were thus determined for different concentrations of the various collectors. Hallimond tube flotation data were taken from the authors' previous publication1 for correlation with that from the vacuum flotation. RESULTS AND DISCUSSION The results of vacuum flotation studies for determining critical pH-concentration curves, i.e., curves which delineate conditions for incipient flotation, are shown in Fig. I. In generaI, all the curves exhibit an upper and lower pH limit between which flotation will

Jan 1, 1969

By E. D. McMurry

INTRODUCTION Before discussing the use of the free piston; it may be well to explain briefly the type of gas lift valves with which the loo] is designed to function. The gas lift valve ii nothing more or less than a small back pressure regulator. It; construction consists of an outer hous ing covering a seamless metallic bellows in such a way as in form a sealed pressure vessel or pressure dome. To the bellow* in attached a stem and to this stem is attached a .rating member. When the valve is run into the hole, a pre-delermined pressure charge is placed in the pressure dome. As shown in Fig. 1. the valve is furnished in four different major assem-blies. to accomplish good results in different types of wells. The valve is mounted on a pup joint (Fig. 2a) : the eoin-plele assembly is entered into the tubing siring at the desired location. Each valve is furnished with a revere flow check valve (Fig. 2-b) to eliminate any back wash through the mechanism. Valve action is as follows: At any time the casing pressure is allowed to exceed the pressure forces within the pressure dome of the valve, the valve simply opens and passes gas from the rasing annulus into the tubing. When the casing pressure forces arc allowed to bleed to a force less than the pressure force within the dome, of the valve. the valve closes. With this type of construction. it is possible to open and close the valve or valves at will from the surface, the medium of the gas pressure in the rasing annnlus. The cycle of operation is governed by a tine cycle motor valve, commonly referred to as a gas lift interutitter or gas life surfacr controller (Fig. 3). In actual practice. after the well is unloaded only one valve is in operation . . . all other valve.: remain in the closed position. On pre-determined time cycles, it is possible to open and close this operating valve at will, and to pass the desired amount of gas through this valve wheneler needed. This method of gas lift is widelv used and has been found In be very efficient. However, it has long been known that even greater efficiency would be obtained if it were possible to place Piston between the injected gas and the fluid head. The free piston which will be discussed here was designed to accomplish this end and has been found to do so with very satisfactory results. DESIGN AND OPERATION OF THE FREE PISTON Basically. the free piston consist- of a telescoping assembly in which are incorporated the various control and operating parts (Fig. A set of sleeve valves which are opened and closed by movement of the telescoping section is contained in the. tool. A cross section of the free piston is shown in Fig. 4-11, 'rile main body cotain- a seamless monel bellows of the same type as that generally used in pressure charged gas lift valves (Fig. 4-c), and a valve core to permit placing a pressure charge in the bellows before the free piston is placed in service. The bellows is the main control unit of the free piston and performs two function: 1. It governs the opening and closing of the sleeve valves. a.: dictated by the producing activities of the well in which the piston is operating. 2. In its normal position, it holds a rubber packer in an elongated position to permit passage of fluid between the packer and the tubing walls. The lower section of the free piston contains a hydraulic shock absorber which cushions the return of the tool to its

Jan 1, 1953

By S. K. Tarby, W. O. Philbrook

Limited experimental data and a critical review of the literature are given to indicate that the true equilibrium distribution of manganese between carbon-saturated iron and blast-furnace type slags has not yet been satisfactorily established as a function of the composition, temperature, and pressure of the system. THE distribution of manganese between carbon-saturated iron and blast-furnace type slags has been studied by a number of investigators.1-7 Generalization of the results of these works indicates that the weight percent distribution ratio, (%Mno)/[%Mn], decreases with increasing temperature and slag basicity, in accord with thermodynamic data and practical observation. However, where comparisons among these studies can be made for similar slag compositions and temperatures, the reported values of the manganese distribution ratio are discordant. The variations are attributed to failure of the experimenters to bring their systems completely to equilibrium with graphite and CO at 1 atm, as judged by departures from silicon equilibrium. Even in the most recent work,5-7 this problem was recognized but the authors made the questionable assumption that equilibrium could prevail for manganese distribution even though other reactions were still in progress. Manganese equilibrium values for three particular slag compositions were required incidental to a kinetic study8 and were determined experimentally. Special pains were taken to approach concurrent equilibrium with respect to silicon. Although the results were not completely definitive, they indicate that the true manganese equilibrium distribution ratios are lower than have heretofore been reported. The data are presented together with a brief commentary on the literature to suggest the need for further work in this area. EXPERIMENTAL METHOD Equilibrium runs were carried out in a stationary, induction-heated graphite crucible closed by a graphite lid and flushed continuously by a slow flow of carbon monoxide at 1 atm pressure. Temperature was controlled manually within ±3 deg C with the aid of immersed thermocouples in mullite tubes protected against slag attack by graphite sleeves. The desired temperature was established by a calibrated Pt vs Pt-10 pct Rh thermocouple and control was then maintained using a more stable W-Mo couple. Runs were made at 1500º and 1575ºC for durations of 7.5 to 10.5 hr. The charges were calculated to yield 400 g of metal and 100 g of slag. The primary slag compositions and designation numbers were as follows:

Jan 1, 1963

By W. N. Hiatt

The mathematical theory of the flow of two-phase, incompressible fluid through porous media is clarified, and the development of a general fluid displacement equation for two-phase, incompressible vertical flow through porous media is outlined. The relationship between the three components of the equation — the "fluid drive effect", the "gravity effect", and the "capillary effect'' — is discussed, as well as implications of the analysis in connection with the study of certain types of petroleum reservoirs. INTRODUCTION The purpose of this paper is to clarify and summarize the mathematical basis of the theory of two-phase incompressible fluid flow through porous media, particularly as it pertains to vertical flow, and to present a qualitative* exposition of the controlling factors involved in circumstances where this type of flow may exist. Although the generally accepted theory of heterogeneous fluid flow through porous media was developed a number of years ago by Muskat,1 and Buckley and Leverett2, 3 using the basic premise of Darcy, (and since then expounded upon by many others) there is a need for an over-all summary and clarification of this theory, especially regarding certain types of applications in connection with so-called "gravity-drainage-type reservoirs". The source and derivation is outlined for what is referred to here as the "General Fluid Displacement Equation for Incompressible Vertical Two-Phase Flow through Porous Media''. This equation would apply under the condition where the pressure drop along the path of flow is negligible in comparison with the pressure level on the system, or where compressibility and solubility effects are negligible. These conditions exist ordinarily in gravity-drainage-type reservoirs. Specifically, it will be shown that this general equation may be developed directly from the basic equations, which, according to present theory, govern the general case of viscous heterogeneous flow through porous media. This general equation has appeared in the literature previously, but in most cases its development from basic theory is not indicated clearly, or the equation is presented in terms of groups of constants and variables involving reduced or dimensionless quantities. This presentation is expedient in cases of experimental or computer applications; but it may make it difficult for the practicing petroleum engineer to identify all the factors involved or to understand the relationships between them. Therefore, in the following development, only ordinary mathematical parlance will be used, and no terms involving groups of other equations or reduced or dimensionless quantities will be used. DERIVATION The following equations for two-phase oil and gas flow result from analogy with Muskat's general equations,** if solubility and compressibility effects are neglected, and if all volumes are considered in terms of reservoir conditions, and the water phase, if present, is considered to be immobile (or effectively a part of the porous medium.) In this case, the sum of the oil and gas saturations may be taken as unity, and the saturation of each phase is in terms of the total hydrocarbon saturation. Note that in Eqs. 1 and 2, and in the derivation following, z is taken as positive in the vertically downward direction,

Jan 1, 1969

By David Cornell, D. L. Katz

Procedures for computing turbulent flow of gas in steady state near the well bore and a graphical method for predicting unsteady state laminar flow at distances from the well have been combined to compute pressure gradients in gas reservoirs. Methods are discussed for predicting single and multiple transients at constant flow rate in an infinite reservoir, predicting constant How rate. in a finite reservoir, reproducing and interpreting back pressure test data, and prediction of the behavior of a closed-in gas well. An example of the graphical method is given for a single transient and the results are compared to a published analytical solution. A typical calculation of the pressure gradient in a reservoir and the production of gas is made .starting with data from a hack pressure test. INTRODUCTION The calculation of pressure gradients throughouse a ga-reservoir at any point in its production history is a complex problem involving turbulent flow near the well bore and unsteady state flow of the gas from the reservoir. The problem is complicated further by the variety of boundary conditions that may be imposed upon the flow. Such boundary conditions include finite or infinite reservoirs. constant or variable rates of production, constant or varying bottom hole pressures, complex initial pressure distributions throughout the reservoir, and pressure maintenance through cycling. In addition to these problems there are the practical difficulties associated with natural gas reservoir analysis These might include: (1) variation. in the permeability and porosity throughout the formation. (2) variations in the thickness of the formation, (3) communication of the producing formation with other producing zones. (4) radial variations in permeability due to influx of drilling fluid. acidizing. or the buildup of liquid in the formation around the well, (5) partial penetration of the producing zone. (6) water flood, and many others. These problems will be considered eventually. In the mean- time, it is necessary to present methods of handling the case of radial flow through a homogeneous, regular, producing stratum to a single well uncomplicated by other factors. The idealized case of steady state, radial flow from a natural gas reservoir was studied by Elenbaas and Katz5 and by Mac-Koberts.' An analysis of the unsteady flow of gases through porous media has been made by Aronofsky and Jenkins' for the one dimensional laminar case with the boundary conditions of constant downstream pressure and uniform initial pressure. Solutions for the partial differential equation for unsteady state radial flow of liquids through porous media have been given by Van Everdingen and Hurst" for the constant bottom hole pressure and constant production rate cases for laminar flow in a reservoir initially at a constant pressure. The analogous heat transfer equation has been treated by Perry and Berggren9 for the case of quenching a cylindrical hole in an infinite medium to a constant temperature. The back pressure curve, which involves the variation of the bottom hole pressure as it is determined by the behavior of the reservoir as a whole. has been studied by Rawlins and Schellhardt,10 Binck-lev,' Baumel and Breitung,' and others. No complete. adequate treatment has been given for the unsteady state, radial flow of gases into a cylindrical hole from a porous medium with or without the presence of turbulent flow and for any set of boundary conditions. If one relies only on analytical solutions of the basic equations, a new solution must be obtained for each new boundary condition. Furthermore, the complexity of the mathematical expressions for boundary conditions other than the very simple cases limits the use of analytical solutions. Graphical methods exist, however. that are general in nature, accurate, and readily employed. The analysis of steady and unsteady, laminar and turbulent, radial flow for various boundary conditions with application to specific natural gas well problems by means of graphical procedures forms the content of this paper. STEADY STATE RADIAL FLOW EQUATION FOR LAMINAR AND TURBULENT FLOW A steady state radial flow equation for laminar and turbulent flow through unconsolidated sands has been given by Muskat as discussion of the work by Elenbaas and Katz.5 This equation is based on the properties of unconsolidated

Jan 1, 1953

Jan 1, 1929

By E. J. Dulis, R. M. Fisher, K. G. Carroll

IT is well known that ferritic steels containing more than 15 pct Cr when subjected to temperatures in the range of 700" to 1000°F exhibit increasing hardness and decreasing ductility. The phenomenon has been widely termed the "885 °F embrittlement," after the temperature of most marked effect.'.' In view of the excellent review articles available in the literature3-6 only a brief account of experimentally established facts need be given here. The extent of changes in physical characteristics during embrittlement depends on chromium concentration and time at temperature, higher alloy content and longer time both promoting more rapid and extensive changes. In a 27 pct Cr steel, changes in impact strength and in angle of fracture in bending can be detected after only a 1 hr exposure at 885°F; after 50 hr this steel becomes quite brittle. Hardness increases slowly with time during thousands of hours exposure and may attain a maximum hardness number twice as large as that of the unexposed steel. Microstructural changes accompanying embrittlement have been described as an initial widening of grain boundaries followed by eventual darkening of ferrite grains. Embrittled steels etch more readily, e.g., the weight loss of a 27 pct Cr steel in acid solution may occur at a rate one hundred fold greater following exposure at 885 OF. Marked changes which accompany embrittlement have been observed in electrical resistivity, specific gravity, and magnetic coercive force. Changes in physical properties may be readily removed by heating at temperatures above the embrittling range, such as a treatment at 1100°F for 1 hr. It has frequently been noted that the 885 °F embrittlement suggests precipitation on a submicro-scopic scale of a chromium-rich constituent, the nature of which has not been revealed by X-ray diffraction. Progressive broadening of the body-centered cubic diffraction lines during embrittlement has been observed," and recent observations by Lena and Hawkes' upon single crystals have shown early asterism in X-ray photographs, disappearing within an hour at 900 °F. Many workers have ascribed8-13 the phenomenon to a precipitation of a phase (FeCr), which is known to cause embrittling effects at temperatures much higher than 885°F. Two general observations, however. suggest that a precipitation cannot be responsible for the 885°F phenomenon: 1—prior cold work greatly enhances a formation, whereas it scarcely affects the 885 °F embrittlement, and 2—the presence of an alloying element such as nickel or manganese may have an effect on the 885 °F embrittlement which is opposed to its effect upon a formation. The slight enhancement of a formation and 885°F embrittlement observed in the presence of elements with strong carbide and nitride forming tendencies' is probably a consequence of lessened chromium depletion of the matrix. The bar graph in Fig. 1 shows a typical example, taken from two 27 pct Cr steels used in this work, of the hardness after exposure for 10,000 hr at 900°, 1050°, and 1200°F. Steel A (0.03 pct C, 3.13 pct Mn) showed marked hardening at 900" and 1200°F, whereas steel B (0.12 pct C, 0.63 pct Mn) exhibited only the 900°F hardening. The cr phase was found in steel A at the higher temperatures but not in steel B. Presumably a formation is enhanced by the low-carbon and high-manganese concentrations in A," Thus there are two distinctly different hardening phenomena present which cannot both be ascribed to v precipitation without invoking a transition phase possessing remarkable properties. Materials A number of chromium steels exposed for long periods (5000 to 34,000 hr) at 900°F, as well as unexposed samples of one of the steels, were available for this investigation. Table I gives the chemical compositions and aging treatments of these steels. In addition to these steels exposed in the elevated temperature test furnaces of the National Tube Division, a number of high-chromium steels were heated for short periods in small laboratory air furnaces and lead baths. Supplementing these commercial steels, a sample of high-purity (0.018 pct C, 0.002 pct N) 28 pct Cr iron, exposed 1000 hr at 887°F. was furnished by the Union Carbide and Carbon Corp. In addition, an alloy of iron and chromium of high purity containing 46 pct Cr was used. This

Jan 1, 1954