Search Documents

By Ralph Hultgren

IN solid solutions atoms of differing sizes occupy the same crystalline lattice, requiring that some of them be compressed and others expanded. The energy involved has been called misfit strain energy and is an important concept of crystal chemistry. If the atomic sizes and elastic constants of interatomic bonds are known, the misfit energy may be calculated,' provided certain simplifying assumptions are allowable. Usually, isotropic crystals are assumed and interatomic distances are taken to be the statistical average determined from X-ray diffraction. Such calculations yield values of the misfit energy of the order of 1 or 2 kcal per atom in alloys such as Au-Cu at compositions of 50 atomic pct. However, evidence has accumulated in recent times that atoms change their sizes with composition of alloys, implying electronic rearrangement of the bonds. The size changes have been found particularly by application of the X-ray method developed by Warren, Averbach, and Roberts.' Thus, Averbach, Flinn, and Cohen3 determined radii in Au-Cu alloys. Oriani4 showed that these new radii led to a calculated misfit energy in disordered AuCu, which was decreased from the values calculated by the usual theory more than twenty-fold, to only 80 cal per g atom. Thermodynamic calculations from the phase diagram5 also show misfit energy to be no more than a few hundred calories per g atom in this alloy. The question of what electronic rearrangements are possible therefore becomes compelling in estimating misfit energy. In the following pages the results of certain calculations on the AuCu tetragonal superlattice are submitted. Conclusions drawn from these should be applicable in large degree to disordered solid solutions. As in all ordered states, bonding distances in the superlattice are individually known, rather than being merely average distances as found from lattice constants of disordered states. Moreover, only the Au-Au and Cu-Cu distances are strained; the elastic constants of these are known in the elementary state. In the usual calculation it is necessary to assume elastic constants for Au-Cu bonds. Misfit energy has thus been calculable without the need of many simplifying assumptions usually made. It is still assumed that equilibrium bond lengths and elastic properties of the bonds are the same in the alloy as in the pure metals. As previously discussed, this is probably not correct. Also assumed is that the bonds are not affected by strain of neighboring bonds. A calculation of Young's modulus from compressibility data shows this to be far from true; extensive electronic rearrangements take place. It would seem that misfit energy cannot be calculated from elasticity data for the elements. The usual methods may, however, give an upper limit which is often much higher than the true value. The question of electronic rearrangement is, of course, a complex one. Pauling's theory gives a simple, approximate treatment of the relation between type of bond and bond distance. This has been applied with some success to the Au-Cu system, as will be shown in a later section. Misfit Energy in Au-Cu Alloys Hume-Rothery and Raynor6 discuss the Au-CU system as a type example of strain energy. The gold atom is 12.8 pct larger in diameter than the copper atom, near the size factor limit beyond which solid solubility is severely restricted. They therefore consider the misfit energy to be large, a conclusion for which they believe they find evidence in the phase diagram. Gold and copper are completely miscible in the solid state, but the alloy has a minimum melting point at an intermediate composition. From this Hume-Rothery and Raynor conclude that the strain energy is nearly large enough to prevent miscibility; the phase diagram tends toward a eutec-tic type. In Ag-Cu, which has almost identical size relationships, solid miscibility is quite limited; whereas in Au-Ag, where atomic sizes are nearly the same, there is complete miscibility without a minimum in the melting point. From their arguments the heat of formation of Au-Cu would be expected to be endothermic or only slightly exothermic, that of Ag-Cu to be endothermic, and that of Au-Ag to be exothermic. Deviations, from Ve-gard's law of additivity of atomic radii support these conclusions, since Au-Cu and Ag-Cu both have pronounced positive deviations, and Au-Ag has a negative deviation. Nevertheless, Au-Cu alloys form exothermically; indeed, considerably more exothermically than Au-Ag, Table I. Hence, strain energy must be much less important in this case than Hume-Rothery and Raynor have supposed.

Jan 1, 1958

By F. J. Fayers, R. L. Perrine

For many years the problem of increasing ultimate recovery of oil from a reservoir has been a subject of interest to the oil industry. At present, a standard secondary recovery technique is to flood the reservoir with water once natural forces of the system are depleted. The permeability characteristics of the rock, together with interfacial forces which exist between the two fluids and the porous material. allow only part of the oil in the rock formation to be recovered by the water-flooding process. Laboratory experiments have shown that when detergents are added to the flooding water, the flow properties of oil in a porous rock arc improved. This improvement results in greater oil recovery. Even though higher recoveries are obtained using detergents, the extra cost associated with the method may make it economically unsatisfactory. To establish economic feasibility requires experimental laboratory measurement of fluid and rock properties together with a calculational procedure to predict field results. Pilot tests in the field may also be required. This paper presents a mathematical procedure for predicting detergent flood behavior. A number of reasonable simplifying assumptions are made to obtain a tractable system of equations. The resulting method may be regarded as a generalization of Buckley-Lev-crett' theory to solve simultaneously for water saturation and detergent concentration profiles. Consequently this type of calculation can be used in much the same manner as Buckley-Leverett theory is used to evaluate a water flood. The mathematical method employed is the so-called method of characteristics for solving hyperbolic partial differential equations. This application of the method is instructive because several other problems in reservoir analysis may be analyzed in a similar manner. FLUID FLOW EQUATIONS The flow behavior of oil and water without detergent through porous reservoir rock is usually expressed by Darcy's law and a continuity equation for each phase. If gravitational forces are neglected, vrg = kr (s)µr v p rg.........(1) vrg = kr (s)µr v p rg.........(1) ? pa va = F ?/?t (sd) where v is the volumetric flux, k represents permeability, S is saturation, and P the pressure. Fluid viscosity is denoted by 11, p is the density, 4 is porosity, and I the time. The subscripts w and o refer to water and oil. The capillary forces in the system lead to a difference in pressure between phases. This may be expressed by the relation: Pm -P0 =Pc(S).........(5) The capillary pressure? P,, is determined experimentally as a function of saturation. DETERGENT BALANCE EQUATION Addition of detergent to the water phase affects the flow properties of the system. Capillary pressure and relative permeability become functions of detergent concentration, c, as well as fluid saturation. That is: P--P (S,c) and k,,,.- - -- k,,,. At low oil saturation values, the presence of detergent reduces capillary pressure and reduces the relative permeability ratio. Detergents are surface-active and adsorb on rock surfaces. The quantity adsorbed, a, will depend on the concentration c of detergent solution and may also depend on the saturation; thus, a = a(S,c). In this analysis we have assumed that the detergents used are water soluble but not oil soluble, and their movement is therefore restricted to the water phase. The dispersion of detergent in water is neglected by comparison with transfer by flow within the system. Detergent movement is determined by the water flux and adsorption: ? cv m = -F ?/?t (cs + a) We will assume that the quantity of detergent in solution is too small to significantly change water density or viscosity. When the permeability and capillary pres-

By M. R. Tek, M. L. Katz

The flow of fluids in stratified porous systems, in which adjacent layers of the system are in communication over part or all of their common interface, in many instances involves crossflow between layers. A theoretical study has been made of the unsteady-state flow of slightly compressible fluids during depletion of bounded stratified porous systems. Consideration has been restricted to single-phase flow, with negligible capillary and gravity effects. Analytical expressions have been derived for pressure distribution and for (luxes at the producing face in two-layer systems with complete fluid communication between layers. These expressions occur in the form of double Fourier and Fourier-Bessel series. The method of extending the solution to any number of layers is demonstrated by deriving the equation for pressure distribution in a three-layer system. On the basis of the exact results for two-layer systems, certain generalizations regarding the behavior of stratified systems in general are presented. The performance of stratified systems, in terms of cumulative flux as a function of time, is shown to lie at all times between bounds established from single-layer theory. The upper bound is given by treating the system as a single layer with arithmetic average physical properties. The corresponding lower bound is the summation of the fluxes from each layer as treated individually. Variation in performance of stratified systems between the upper and lower bounds is determined by the extent to which crossflow occurs in the system. The principal parameters affecting the amount of crossflow are the vertical permeability of the system, the ratio of system thickness to "radius of drainage" and, for radial systems, the ratio of system thickness to interior radius. At very early times, the system performance is near the lower bound. As time and the radius of drainage increase, the per- formance tends toward the upper bound. INTRODUCTION One of the principal problems confronting the reservoir engineer today is the incorporation of the effect of reservoir heterogeneity into predictions of reservoir performance. The first part of this problem is the characterization of the actual reservoir heterogeneity in a mathematically useful manner.1, 2 It is then necessary to calculate the effect of this heterogeneity on the system performance. It is frequently possible to estimate this effect by comparison with the known behavior of an idealized system. One such idealized system is the stratified or "layer-cake" model. Many studies have been made of the behavior of stratified systems, from both a theoretical and experimental standpoint. The great majority of the theoretical studies have been concerned with systems of separated layers, with no fluid crossflow. Among the more noteworthy of such theoretical papers is that recently published by Lefkovits, et al,3 in which a comprehensive analysis was presented for the behavior of bounded stratified systems which are unconnected except at the wellbore. Three recent references have dealt with stratified systems including the effect of crossflow. Jacquard4 has presented mathematical derivations pertaining to depletion of two layer systems producing at constant rate. A doctoral dissertation by Katz5 contains analytical solutions and tables of results for pressure distribution and fluid flux during constant-terminal-pressure depletion of two-layer and three-layer systems with crossflow, together with tables of eigenvalues necessary in the computations. The dissertation also presents results from experimental studies of stratified systems using a heat-transfer analog model. This paper is largely based on the material in this dissertation. A third paper has been presented by Russell and Prats6 since completion of the work resented herein, containing mathematical derivations for and discussion of the behavior

By F. E. Jaumot, R. L. Smith

Self-diffusion in zinc has been used as an instrument for comparison of the absorption and sectioning techniques as a means of studying diffusion. Single crystal as well as poly-crystal samples were used asanda the temperature range of diffusion extended from 200° to 415°C. For temperatures above 200°C, the data indicate that the results obtained from the absorption technique agree with those obtained from the sectioning technique. The effect on the values of the diffusion coefficient of electroplating vs evaporation as a means of applying the tracer was investigated and no significant difference observed. It was found that an excess or deficiency of tracer did not materially affect the results obtained from the sectioning technique, but invariably caused errors with the absorption technique. NEARLY all self-diffusion experiments and many alloy diffusion experiments have employed radioisotopes as the instrument for following the diffusion. Of the various techniques in which radioisotopes are applicable, the two most widely used are the sectioning and absorption techniques. In the former, the diffusion coefficient is obtained from a curve of log concentration vs the square of the mean penetration depth. In the latter, the diffusion coefficient is calculated from a knowledge of the decrease in activity at the surface after diffusion. In spite of the importance of knowing the reliability of the diffusion data and the widespread use of these two techniques, they had not been directly compared experimentally on the same samples. Thus, such a comparison was worthwhile. The specific system chosen for this comparison was the self-diffusion of zinc, since the metal was readily obtainable in high purity, single crystals were easily grown, and the anisotropic features of the diffusion added an additional factor of interest to the comparison. Self-diffusion in both single and polycrystals was investigated, since the possibility existed that grain boundaries and grain boundary diffusion might produce different effects for the different techniques. Also, the effects of electroplating vs evaporation as a means of applying the radioactive tracer, and the effects of applying an excess or deficiency of tracer material were studied. The general experimental procedures and the mathematical analyses of both techniques have been described previously in some detail, and only those points essential to an understanding of the present work will be discussed here. The mathematical analysis of the sectioning data requires only a knowledge of the solution of Fick's equation. where C,, is the initial concentration, C is the concentration at depth x, D is the diffusion coefficient, and p is the absorption coefficient of the solvent metal for the radiation from the solute. Having obtained F from the experiment, the value of z can be read from a plot of F vs In z, and D can then be calculated if is known. The necessity of knowing the absorption coefficient causes most of the criticism of the absorption technique. It is very difficult to measure accurately under the conditions existing in the diffusion experiment, but it is absolutely necessary, as will be indicated briefly later. Portions of this problem have been discussed in detail elsewhere." The difficulties involved in the measurement of the absorption coefficient arise from three principle factors: the physical characteristics of the absorber, the geometry, and complex radiation spectra. The first two are heavily dependent upon the third. If

Jan 1, 1957

By Z. F. Oszter

This paper discusses the materials handling problems encountered in cement plants in the areas of raw and finished product grinding and cement storage. It illustrates by examples taken from actual installations, how optimum feed rate is maintained in grinding circuits with properly applied bin dischargers (vibrating hoppers, Mexican feed boxes and aerated bin bottoms) for materials with high moisture content and/or difficult handling characteristics. Initial operating problems with the equipment will be discussed and maintenance requirements will be indicated for comparison purposes. In a modern cement plant, properly functioning bins and feeders are as essential as process control, for no process control system, however sophisticated, can perform efficiently if raw materials do not flow into the plant processing facilities at the required time and rate. Most materials used in the manufacture of portland cement - such as silica sand, iron ore fines, crushed limestone, gypsum, clinker or cement itself-have difficult handling characteristics under specific conditions. High moisture content or even dampness due to condensation may create bin hang-ups, arching and ratholing, all of which result in improper functioning of the feeders or proportioning equipment. The occurrence of these problems in existing cement plants spurred a search for the corrective devices discussed in this article. The equipment investigated included bin dischargers in the form of vibrating hoppers for moist raw materials, Mexican feed boxes for clinker and aerated bin bottoms for cement. Operating experience with different types of bin dischargers showed that the actual performance of some did not come up to the expected level under prevailing operating conditions. The aim of this paper is to record the results of this experience and there- by aid designers and/or operators in selecting the most suitable equipment for their requirements. BIN DISCHARGERS (GENERAL) A discharger is a device that causes stored solids to flow out of a hopper without metering the rate of flow as a feeder does. The investigation discussed here encompasses the following types of bin dischargers used in a cement plant: (1) vibrating hoppers (horizontal motion), (2) whirlpool hoppers (horizontal and vertical motion), (3) aerated bin bottoms (air fluidization), and (4) Mexican feed boxes. Most of the bins used in the cement industry are of the plug flow type. Raw materials with difficult handling characteristics such as high moisture content may cause flow problems when fed through small openings (12 x 18 in.), even if two walls of the bin are vertical.' The most appropriate solution to these flow problems would be the installation of separate mass flow bins for the different materials such as silica sand, iron ore fines, crushed limestone and gypsum, with each bin correctly designed from test data obtained by measuring the flow properties of representative samples. However, limited storage capacity and plant space usually preclude this approach. Quantitative design of bins for flow based on measured flow properties of representative samples of the bulk solid to be handled is no longer an "engineering art," thanks to the basic research done on the flow of bulk solids in bins during the last decade.2, Where space limitations exist, however, or where improvements are being sought for equipment already installed, bin dischargers offer a satisfactory solution to flow problems. They can be applied effectively under flat-bottom bins or mounted to existing plug-flow bins with very little modification to these facilities. In both the foregoing cases, plug flow will be maintained, but the flow channel will be considerably increased so that it is larger than the discharger's bin-side diameter, and material flow can be sustained if the discharger is correctly selected. The main purpose of a bin discharger is to draw the bulk solids across the whole area of the bin outlet, to make the opening fully effective. As long as the bin-side opening of the discharger is large enough

Jan 1, 1969

By G. R. Belton, Y. K. Rao

A galvanic cell, using liquid MgCl2 or MgC12-CaC12 mixtures as the electrolyte, has been used to determine activities in Mg-A1 liquid alloys between 700' and 880°C. The incovporation of a chlorine electrode in the cell also allowed measurements to be made of the standard free energy of formation of MgCl2(l). The results are shown to be in good agreement with thermo-chetrlical values from the literature, and this is taken as evidence that the small, known solubility of magnesium in iMgCl2 introduces no significant error in galcanic cell measurements. Within experimental error, the activity coefficients and relative partial molar enthalpies at 800°C are shown to be represented by the following "subregular" solution equations: logy~~ =-0.68(1 -xMgj3 log yAL =-1.02(1 - XMf + 0.68(1 - XMf H.M = -4400(1 - %)3 / cal Mg 7Ag' HZ =-6600(1 - xA1)' + 4400(1 - XMf cal SCHNEIDER and toll' have used a transpiration technique to measure the vapor pressures of magnesium over Mg-A1 alloys between 544" and 594°C. amsstad,' however, has since suggested that the extrapolation to zero flow rate, used by these authors in interpreting their apparent pressure vs flow rate data, gives unreliable results. Rogers, Tomlinson, and Richardson,3 in interpreting the results of solution equilibria between Mg-A1 alloys and liquid MgC12, also considered the measurements of Schneider and Stoll to be unreliable and preferred to derive activities for the alloys from the earlier boiling point determinations of ~eit~ebel~ and the partial molar heats recommended by Kubaschewski and ~atterall.~ These latter heats were based substantially on the early calorimetric work of Kawa-kami but, unfortunately, his work on other systems has sometimes been found to be inaccurate.7 Rogers et al., in the above-mentioned paper,3 tentatively concluded that the most likely species responsible for the limited solubility (0.3 mole pct at 800°C) of magnesium in MgC12 were Mg° (neutral) and Mg2++. Two more recent studies8,9 have supported Mg2++ as the soluble species. In the present study, activities in Mg-A1 alloys have been determined by means of a galvanic cell involving liquid MgCl2 or MgC12CaCl, mixtures as the electrolyte. Since the reactive nature of magnesium precluded simple Faraday yield experiments, a chlorine electrode was incorporated in the cell in order that the performance of the cell could be checked by measurements of the heat and free energy of formation of MgC12. This procedure was considered necessary since it has been suggested1' that the solubility of a metal in a molten salt might introduce electronic conductivity; also, previous determinations of the standard electrode potential for MgC12 differed by as much as 70 mv11-13 EXPERIMENTAL Materials. Analyses of the materials used in preparing the alloys and the electrolyte mixtures are presented in Table I. The alloys were prepared by induction melting weighed amounts of the metals in a graphite crucible held under an argon atmosphere. Pure anhydrous magnesium chloride was prepared by heating the mixture MgCl, . 6H20 +NH4C1(1:1) to 650°C, followed by melting under dry argon. The melting point of the dry MgC12 was found by differential thermal analysis to be 714.8oC, which compares well with the accepted value of 714C.14 This agreement was taken to be an indication of the high purity of the dried salt. Table I. Compositions of Materials, wt pn Impurity Mg Al MgCI2 CaCI, Ba - - 0.005 Ca - - 0.010 Cu 0.02 0.02 Fe - 0.10 0.001 0.010 Pb 0.01 - 0.001 0.005 Mn 0.15 0.001 Si - 0.10 Sr - - 0.005 MBSO* - 0.040 ARGON TUNGSTEN CHLORINE LEAD—411 11 II ALUMINA 'A I / ave GRAPHITE on irA~_ ( ^"l ROD MAGNESIUM S>ILICA~^ /, OR ALLOY-. I /A] ___________ -«^- _ ELECTROLYTE I y FRITTED DISCS Fig. l—Arrangement of chlorine and metal electrodes in electrolytic cell.

Jan 1, 1970

By H. C. Fiedler

The development of cube-on-edge secondary re crystallization texture in Si-Fe strip depends upon the ability of inclusions, such as manganese sulfide, to restrain nomal grain gvowth. The ability of inclusions to perform this function is shown to depend upon their size and distribution. The desired uniform dispersion of many small inclusions is obtained by controlling the cooling rate from above the solution tempmature of the inclusions or, alternatively, by drastic 'quenching from above the solution temperature followed by heat treating at a lower temperature. The development of the cube-on-edge secondary recrystallization texture in Si-Fe strip depends upon both establishing the proper primary recrystallization texture and having inclusions to restrain normal grain growth. It has been suggested that the first requirement is met with a texture which, though weak, is near (110)[001]. A mechanism for secondary recrystallization originally proposed by ielsen' is consistent with the presence of this primary texture. It is visualized that a secondary grain appears where two adjacent grains happen to be sufficiently similar in orientation to be capable of behaving as a single large grain. This pair of grains will grow at the expense of its smaller neighbors as a consequence of its having lower grain boundary free energy. Other explanations of the origin of the secondary recrystallization texture have been offered,' and it can be concluded that the mechanism is not fully understood. Nevertheless, it is apparent that operationally the necessary primary recrystallization texture is obtained simply by interposing an anneal during the cold reduction to final gage. The role of inclusions in the development of the secondary recrystallization texture has been demonstrated by May and Turnbull4 in a comparison of a high-purity heat and a heat containing manganese sulfide. They showed that the driving force for secondary recrystallization is inversely proportional to the diameter of the matrix grains and that the function of the inclusions is to maintain a small matrix grain size. Other inclusions, such as various sul- fides,')' silicon nitride,8 and vanadium nitride,' titanium carbide,'' and vanadium carbide," have been shown to be effective in promoting secondary recrystallization to the cube-on-edge texture. There is also evidence that silica inclusions in Si-Fe strip made by powder-metallurgy techniques perform this function.' This report is concerned with the second of the two requirements for secondary recrystallization, namely the importance of the metallurgical control of the size and therefore number of inclusions. While the concepts involved are simple, they are fundamental to the making of this material. EXPERIMENTAL PROCEDURE A heat containing 3.27 pct Si, 0.026 pct S, 0.06 pct Mn, and 0.004 pct C was made in an air induction furnace from electrolytic iron, 98 pct ferro-silicon, electrolytic manganese, and iron sulfide. The heat was divided between graphite and baked-sand molds with a cavity 2-5/8 by 5 in. and with walls 2-3/4 in. thick. Slabs, 1 in. thick and 2-5/8 by 5 in., were cut from the ingots and hot-rolled either in the as-cast condition or after one of the heat treatments described later. After pickling, the 80- to 100-mil-thick hot-rolled band was, unless otherwise noted, heat-treated for 5 min at 900°C before being cold-rolled to 25 mils. After a heat treatment of 1 to 3 min at 860°C, the strip was cold-rolled to 12 mils, the final thickness. The percent cube-on-edge texture was determined by dividing the maximum torque as averaged from four samples by the maximum torque of a 3-1/4 pct Si-Fe single crystal having the (110) plane in the plane of the sample. The torque measurements were made in a field of 1000 oe on 1-in.-diam disks. RESULTS AND DISCUSSION The importance of the thermal history of the ingot is evident in Fig. 1 from the appearance of strips heat-treated in a temperature gradient. The sample

Jan 1, 1964

By AIME AIME

IN a recent technical paper of the Central Committee of the French Coal Mines (Note technique No. 198, by L. Lahoussay) the author points out that continuous in¬crease in depth of mine shafts makes it more and more necessary to look for appropriate means to reduce the weights and dimensions of the hoisting ropes. Increasing depth demands a larger section and weight for the cable if the useful load is not to be cut down at the same time. This naturally leads to increased cost per ton, and at the same time creates serious mechanical difficulties. The dimensions of the hoisting engine must be increased, acceleration at the start must be decreased so as to avoid periods of excessive peak power loads, and a longer time be allowed to slow down the speed of the hoist when applying brakes. The hoisting rope, upon which depend the dimensions and hence the mass of the hoisting engine, presents the greatest obstacle to the development of hoisting from greater depths. Evidently, for a given depth, the total weight of the rope can be notably reduced by utilizing a wire rope of decreasing section, approaching the theoretical logarithmic form. This is a solution frequently adopted for great depths. Of course, the weight of the rope could be easily' diminished by reducing the useful load, but it is certain that with increasing depth of mining, more effort must be made to confine the hoisting to a small number of shafts, and increase the load hoisted per trip.

Jan 1, 1933

By E. T. Turkdogan

This is a compilation and a critical review of the data on the Fe-Ca-0 ternary system. Using the results on the reductiorz-equilibria, an oxygen potential diagram is drawn for the greater part of the system. A number of univariant systems and invariant points, not determined experimentally, are evaluated by making use of the theorem on univariant curves intersecting at an invariant point. The formation of solid solutions and ternary compounds and the crystal structures of some of the phases in the composition-diagram are given for a feu! temperatures. DURING the past thirty years there has been a continued interest in the study of the Fe-Ca-0 system. A detailed literature survey made by white1 in 1943 indicated that further work was required to clarify the phase relationships; in fact, the work carried out in the subsequent years, reviewed by Burdese2 in 1954, made a valuable contribution to a better understanding of this subject. Most of the data were obtained by the step-wise reduction of mixtures of iron oxides and calcium ferrites. However, there is some uncertainty about the phase relationships in this system, arising mainly from insufficient use of thermodynamic concepts in the construction of phase diagrams. In this paper an attempt is made to establish the invariant points from the available data. OXYGEN POTENTIAL DIAGRAM Schenck and coworkers3,4 were the early investigators who studied part of this system by the step- wise reduction method and suggested the formation of two ternary compounds, CaO.FeO.Fe2O3 and 2Ca0.5Fe0.2Fez03, the latter being referred to as point X. They also deduced from their results that the phases Y (8.74Ca0.32.5Fe0.8.74Fe2O3 and Z (2.02CaO-97.98FeO) were formed near the w?stite corner of the system. Martin and voge15 postulated the existence of the compounds CaO.9FeO and 4Ca0.3Fe3O4; the former is very close to the composition of point Z and the latter to that of CaO.FeO-Fe2O3. The most valuable contribution to the presen knowledge on this system was made by Cirilli and Burdese6-8 who showed by X-ray and chemical anal yses that there were two ternary compounds, CaO-Fe0.Fe2O3 and CaO.3FeO-FezO3; the former is the same as that suggested by Schenck et al. and the latter is close to their point X. The results of Cirilli and Burdese show that the phases Y and Z postulated by Schenk et al. are solid solutions of calcium oxide in w?stite. In most of the above-mentioned reduction-equilibria measurements carbon monoxide-carbon dioxide mixtures were used and in a few cases hydrogen-water vapor mixtures were employed. Using the

Jan 1, 1962

By L. S. Darken

It has been known for sometime that in an ionic lattice, such as that of Ag,S or FeO, the migration velocity of the anion may differ markedly from that of the cation, the cation being usually the more mobile. It has been shown, for example, that during the formation of silver sulphide tarnish on silver, the tarnish layer grows principally by the motion of silver cations through an essentially stationary silver sulphide lattice. Similarly, the growth of wüstite scale on iron proceeds principally by motion of iron ions, this being accomplished with the aid of vacant iron lattice sites into which the iron may readily move. Likewise investigations of solid solutions of various salts show that there is usually considerable difference in the mobility of the different cations. Measurements of the mobility of various metals in liquid mercury indicate that the various metals move at different rates relative to the mercury. In view of the foregoing it seems reasonable to suppose that when any single phase solution is acted upon by a field, the atoms of different elements respond in different ways; in particular, that the force arising from a particular concentration gradient in a binary alloy would cause the atoms of one component to move with a drift velocity which is in general different from that of the atoms of the other component. The recent experiments of Smigelskas and Kirkendall' support this viewpoint, since they found, in essence, that during diffusion in alpha brass zinc moves appreciably faster than copper. This was evidenced by the relative motion of inert molybdenum wires. Such a phenomenon is incompatible with the place interchange theory of diffdsion proposed by several authors, most recently by Bir-chenall and Mehl.2 This theory has, however, also been questioned3v4 for other reasons, principally because of the high energy required for direct atomic interchange. The rather large difference between self-diffusivity, as determined by tracer technique, and mutual diffusivity, as determined with the aid of another element at low concentration, might also be interpreted as favoring the above idea; however. Birchenall and Mehl have pointed out that there is another possible interpretation of this difference; this is discussed later. A rather interesting and pertinent observation was made by Schwartz, who noticed that tungsten lamp filaments, heated by direct current, commonly fail by necking down in the region at the temperature gradient near the positive end, with a concomitant swelling at the corresponding position near the negative end. He introduces considerable auxiliary evidence to show that this phenomenon is occasioned by the motion of tungsten cations which migrate in the high temperature region from the positive end toward the negative end. The mechanism of this migration is not clearly understood. If, however, this experimental evidence be accepted, we are forced to conclude that alloying elements in the tungsten may migrate by the same or a Similar mechanism, whatever it may be, at a greater Or lesser rate, thus giving rise

Jan 1, 1949

By L. S. Darken

It has been known for sometime that in an ionic lattice, such as that of Ag,S or FeO, the migration velocity of the anion may differ markedly from that of the cation, the cation being usually the more mobile. It has been shown, for example, that during the formation of silver sulphide tarnish on silver, the tarnish layer grows principally by the motion of silver cations through an essentially stationary silver sulphide lattice. Similarly, the growth of wüstite scale on iron proceeds principally by motion of iron ions, this being accomplished with the aid of vacant iron lattice sites into which the iron may readily move. Likewise investigations of solid solutions of various salts show that there is usually considerable difference in the mobility of the different cations. Measurements of the mobility of various metals in liquid mercury indicate that the various metals move at different rates relative to the mercury. In view of the foregoing it seems reasonable to suppose that when any single phase solution is acted upon by a field, the atoms of different elements respond in different ways; in particular, that the force arising from a particular concentration gradient in a binary alloy would cause the atoms of one component to move with a drift velocity which is in general different from that of the atoms of the other component. The recent experiments of Smigelskas and Kirkendall' support this viewpoint, since they found, in essence, that during diffusion in alpha brass zinc moves appreciably faster than copper. This was evidenced by the relative motion of inert molybdenum wires. Such a phenomenon is incompatible with the place interchange theory of diffdsion proposed by several authors, most recently by Bir-chenall and Mehl.2 This theory has, however, also been questioned3v4 for other reasons, principally because of the high energy required for direct atomic interchange. The rather large difference between self-diffusivity, as determined by tracer technique, and mutual diffusivity, as determined with the aid of another element at low concentration, might also be interpreted as favoring the above idea; however. Birchenall and Mehl have pointed out that there is another possible interpretation of this difference; this is discussed later. A rather interesting and pertinent observation was made by Schwartz, who noticed that tungsten lamp filaments, heated by direct current, commonly fail by necking down in the region at the temperature gradient near the positive end, with a concomitant swelling at the corresponding position near the negative end. He introduces considerable auxiliary evidence to show that this phenomenon is occasioned by the motion of tungsten cations which migrate in the high temperature region from the positive end toward the negative end. The mechanism of this migration is not clearly understood. If, however, this experimental evidence be accepted, we are forced to conclude that alloying elements in the tungsten may migrate by the same or a Similar mechanism, whatever it may be, at a greater Or lesser rate, thus giving rise

Jan 1, 1949

By T. M. Allen, W. H. White

THE Greenwood-Grand Forks area of southern central British Columbia, known as the Boundary District, has a long history of mining exploration and production. At the turn of the century this was the premier copper mining camp in the British Empire, its total production amounting to some 20 million tons. Most of this ore came from the great Granby mines at Phoenix, but the Motherlode mine at Deadwood camp, 6 miles to the west, and several mines in Summit camp, 5 miles north of Phoenix, made important contributions. The large deposits were exhausted in 1918 and the district since has seen only desultory exploration and salvage operations. The orebodies are mineralized skarn zones in limestone members of a thick series of Upper Paleozoic sedimentary and volcanic strata. Chalcopyrite is the primary ore-mineral. Copper carbonates and silicates occur sparingly in outcrops, but the oxidized zone generally is very shallow. Much of the surface is mantled by glacial drift which in most places ranges in thickness from 2 to 15 ft. In some of the hanging valleys, however, the glacial drift may be as much as 100 ft thick and may assume drumlin-like forms. In 1951 an ambitious program aimed at the discovery of new orebodies and important extensions of abandoned deposits was launched by Attwood Copper Mines, Ltd. In this district so thoroughly searched by an earlier generation of prospectors, any orebody which had remained undiscovered must have little or no surface indication. Consequently, in addition to the basic detailed geological work, the program of exploration included magnetometer and self-potential surveys. Geological bets and geophysical anomalies were tested further, prior to diamond drilling, by a study of copper distribution in tree twigs and/or in the soil. The soil sampling and analytical methods used and some of the results seem of sufficient importance to warrant this paper. The authors had done some plant sampling in this and other districts, using the dithizone neutral-color-end-point method (Warren and Delavault, 1948, 1949; White, 1950),1-3 but they were unfamiliar with its soil application. Finally, after much experimenting in the field, they adopted the methods described here. These methods are not entirely original or defensible on theoretical grounds, but under field conditions of rapid sampling and analysis the results are reliable enough to be of use. Fig. 1, which shows the results of duplicate analyses of duplicate soil samples taken at 50-ft intervals across an anomalous zone, indicates the relative dependability both of the sampling and analytical methods. Sampling and Analytical Equipment A 2-ft piece of 1-in. solid drill steel, one end sharpened to a broad, conical point. The steel is marked at 1 ft from the point. A 2-ft piece of ½-in. black iron pipe, one end filed to a bevelled cutting edge. The pipe is marked at 1 ft 3 in. from the cutting end. A 3-lb hammer. A plastic or rubberized sheet about 18 in. square. Moisture-proof assay pulp envelopes. A 10-mesh seive made from window screen with the paint burnt off. A small assay spatula. A pan balance sensitive to 10 mg. Two ignition trays about 4 in. square, made of sheet iron turned up along the edges. A Coleman two-burner gasoline stove. An asbestos board about 5x8 in., used as a hot plate on the gasoline stove. A circular aluminum rack to hold 8 test tubes while refluxing (design of Almond and Morris). Pyrex Glassware Large refluxing test tubes, 25x200 mm, marked at 40 ml volume. Breakers, 20 ml. Pipettes, 1, 5, and 10-ml capacity. Graduate, 50 ml. Shaking cylinders, 100 ml, glass stoppers. Burette, 25 or 50-ml capacity, with holder. Chemical Supplies 1 N sulphuric acid. Hydroxylamine hydrochloride, solid crystals. Fisher Alkacid test paper. Copper standard solution. Dithizone standard solution 60 mg per liter. Water reasonably free of metals. Soil Sampling Method: The problem of how to take a soil sample is extremely crucial. The method outlined below, adopted after a number of tests, has the advantages of uniform pattern, uniform depth, and uniform size of sample. The area to be tested was marked off by chain and compass lines 100 ft apart, normal to the strike of possible ore deposits. Numbered stakes were set at 50-ft intervals along these lines and a soil sample was taken at each stake in the following manner. The drill steel was driven into the ground normal to the slope of the surface to the marked depth of 1 ft, moved slightly from side to side, then carefully withdrawn. The iron pipe was inserted to the bottom of this hole, tapped down to the marked depth of 1 ft 3 in. and withdrawn; the 3-in. soil plug in the

Jan 1, 1955

By R. M. Rose, D. P. Ferriss, J. Wulff

Single crystals of tantalum, tungsten, and the binary alloys thereof were grown by electron-beam zone melting and tested in tension between 77° and 373°K. The ductile-to-brittle transition temperature increased sharply with solute concentration, for both tungsten-rich and tantalum-rich alloys; apparently, the maximum transition temperature exists near the 50 pet concentration. Twinning was observed to be a consequence of cleavage-crack propagation. Cleavage occurred on(100)planes, and extensive tearing and "river markings" were present. The onset of brittleness could not be associated with any change in deformation mode, the yield point phenomenon, or the large rise in yield stress with decreasing temperature. Cleavage was observed to initiate at either the specimen surface, or subgrain boundaries. The observation was made that the grain boundaries of Polycrystalline materials are probably even better as nucleation sites for fracture, and are responsible for brittleness at high temperatures and lower solute concentrations than those of the corresponding single crystals. TUNGSTEN and tantalum, the two high-melting neighbors in the third transition period of the periodic table, present striking contrasts. Although both are bee, tungsten has far higher elastic moduli, and the highest ductile-to-brittle transition (DBT) of all the pure metals.' Tantalum, however, is ductile at liquid helium temperatures. Large strain aging and yield point phenomena are observed in tantalum,2,3 but only to a lesser extent in tungsten.4 Tungsten work-hardens at a high rate,5 whereas this effect is negligible in tantalum.' Since the mutual solid solubility of these two metals is complete,7 and no ordering has yet been reported in the binary system, the possibility of depressing tungsten's DBT by tantalum additions was considered. Also, tungsten additions to tantalum could produce hard, high-melting alloys which might be ductile at room temperature. Such alloys offered excellent bases for complex alloy development, particularly with regard to the effect of rhenium on Group VI metals.8"10 In addition, there was the opportunity to examine flow and fracture modes in the vicinity of the DBT, with particular emphasis on the correlation between brittleness and fundamental deformation mechanisms. Myers11 has reported that recrystallized alloys with up to 10 pet W in tantalum are ductile at room temperature, but alloys with 20 pct W could not be swaged at temperatures as high as 500°C. The whole spectrum of Ta-W alloys has been surveyed by Braun, Sedlatschek, and Kiefferl2 using sintered powder compacts. In addition to some physical properties they reported a maximum hardness at about 50 at. pet W, dropping to much lower hardness levels for the pure components. They were able to cold work tantalum-base alloys with as much as 30 pet W and to warm work tungsten-base alloys with up to 20 pet Ta by procedures similar to those used for pure tungsten. EXPERIMENTAL In order to study deformation modes and to keep the impurity concentrations at a minimum, electron-beam zone melting was used to prepare single crystals of high purity. The electron-beam zone melting apparatus is similar to that described by Calverley, Davis, and ever" and incorporates a grid control rectifier to maintain control of emission current. Pure metals were refined by making three zone passes at 4 mm per min up the rod. Alloys were prepared by melting together wire of one constituent on a previously melted rod of the major constituent. Subsequently, two zone passes were made in opposite directions to level the composition after the manner suggested by pfann.14 Purification of the starting materials was accomplished by evaporation of the several impurities, and the degree to which this was accomplished is indicated in the chemical analyses of Table I. The uniformity of composition along a zone levelled rod was within the i 1 pet accuracy of the fluorescent X-ray technique used to determine composition. Crystallographic orientation of the

Jan 1, 1962

By Pol Duwez

The effect of the rate of cooling on titanium, zirconium, and thallium has been measured. For titanium-molybdenum alloys, it has been shown that for molybdenum concentration up to 8 pct the rate of cooling has no effect on the transformation temperature, but this temperature decreases progressively with increasing molybdenum content. HIGH strength to weight ratio and remarkable corrosion resistance have often been advanced as the main reasons for the present widespread interest in titanium. For the physical metallurgist, the existence of an allotropic transformation in titanium is another important factor. As in the case of iron, an allotropic transformation may be the fundamental property upon which the response of Ti-base alloys to heat treatment will depend. It is therefore important to study the mechanism of the allotropic transformation in pure titanium and the effect of alloying elements on this transformation. The existence of an allotropic transformation in titanium was reported in 1931 by Schulze.' Later, De Boer et al.' determined the transformation temperature from a discontinuity in the electrical resistivity at 882"42O°C. According to recent work," the transformation temperature is 882.5"&1°C. The crystal structure of the room temperature a phase is hexagonal close-packed and the high temperature fi phase is body-centered cubic.' It is well known that the high temperature /3 phase of pure titanium cannot be retained by quenching to room temperature. Several investigators'-' have shown that in alloys of titanium with metals which are soluble in the @ phase (for example, nickel, molybdenum, chromium, or manganese), the /3 structure can be retained, providing the concentration of the alloying element is greater than a certain critical value. Below this critical value, the alloys transform during quenching, at least partially, into a supersaturated solid solution having the a. structure. This behavior is very similar to that found in some iron binary alloys (Fe-Ni and Fe-Mh, for example) and is an indication that the transformation might occur by means of a shear mechanism analogous to that observed in the austenite-martensite reaction. The present investigation is mostly concerned with the effect of the rate of cooling on the /3 to a transformation. In the first portion of this paper, the results obtained on pure titanium are described. Two other pure metals, zirconium and thallium, were also investigated, because they both have an allotropic transformation like titanium, from an hexagonal close-packed to a body-centered cubic structure. For purposes of comparison, the effect of rate of cooling of the y to a transformation of iron was included in the experimental program. The second part of this paper deals with Ti-Mo alloys. The equilibrium phase boundaries involving the solid state reactions on the Ti-rich side of the system were determined and a study was made of the effect of the rate of cooling on the transformation of the high temperature /3 solid solution. Experimental Procedure Since the main object of the present investigation was the determination of the effect of the rate of cooling on a transformation temperature, the best technique available was that used by Greninger" in his study of the martensite arrest in carbon steels. This method consists of heating a very small sample of metal to which fine thermocouple wires are .attached, cooling the sample rapidly by means of a blast of helium, and recording the temperature vs. time curve on an oscillograph. The apparatus used in the present experiments was very similar to that described by Greninger, except that the heating coil was made of molybdenum wire (0.040 in. in diam) and the assembly was placed in a bell jar under vacuum of 10" mm Hg or less. The specimen con-

Jan 1, 1952

By B. M. Larsen

Attempting to review the fundamental aspects of segregation in steel ingots of all types in a paper of reasonable length, we encounter two difficulties: (I) the fact that a large number of different physical and chemical effects are involved, and (2) rather incomplete and confusing evidence on certain aspects of detail and mechanism in the freezing process in which the segregation effect is involved. The best solution seemed to be to link together most of the underlying factors into a simplified theoretical picture that appears to the writer to explain most of the observed data, in the hope that this will serve as a reasonably good target for discussion and criticism. Segregation effects, especially in steel ingots or castings, are often divided into " microscopic" and "macroscopic," depending on whether the variations in concentration of solutes or dissolved elements in the iron occur as: (I) a periodic variation over the small distances between the nuclei and the surfaces of the primary grains, or (2) as a larger scale effect involving concentration gradients over distances of the order of several inches, more or less. The reason for this somewhat artificial division is, of course, that it is only the macroscopic scale effect that has much practical importance. About the only effect of practical importance that results from microscopic segregation is the banded structure often left in rolled sections, and even this is relatively harmless in most cases. However, the causes of microscopic segregation also untlerly the larger scale effect, so we must first consider matters such as the formation and growth of crystal nuclei, selective freezing and diffusion, and others that are fundamental to all sepregation effects. Undercooling and Nuclei Formation It is probable that although the degree of undercooling is small in iron, as in metals in general, a certain amount of it takes place, and also that the number of crystal nuclei increases very rapidly with increased undercooling below the liquidus temperature. In practical terms this means that the nuclei number illcreases very rapidly with increased cooling or freezing rate, thus being dependent on the rate of heat flow out from the zone of crystallization. Temperature Gradients and Heat Dissipation during Freezing Each pound of steel evolves about 100 B.t.u. when freezing. In addition to this, the temperature gradient in the solidified portion (required to dissipate the heat as freezing progresses into the axis) necessitates a loss of heat required to cool the solid metal below its freezing range, which is equivalent to an amount of not so much less in magnitude. In a 5-ton ingot, the total would include about one million B.t.u. for freezing alone and about 0.4 to 0.5 million B.t.u. for the gradient, up to the end of the freezing period. Thus the rate of freezing depends primarily on the rate of heat dissipation. In the ordinary ingot shapes, the axial rate of growth of the solid

Jan 1, 1945

By R. G. Hawthorne

This report is concerned with fluid displacement in porous media, in those cases where viscous and gravitational forces control the displacement. Such a system would usually be found in a sand body of large physical dimensions such as an oil reservoir, although it is possible to create such a system in the laboratory. It is shown that the position of the fluid interface can be predicted by numerical calculations using a basic idea presented by Dietz. Fluid flow is considered in a vertical plane in a homogeneous, porous medium of sufficient thickness that the capillary transition zone is small in comparison with the total reservoir. A theory developed by Dietz' is used to make numerical calculations of the position of the fluid interface. The results for several conditions are compared with scaled model experiments. The results show that, for gas drive in a reservoir of steep dip, a relatively low flow rate can displace large volumes of oil before gas breakthrough. On the other hand, water injection at favorable mobility ratio and low dip may show best performance at high rates. Water tends to underride the oil and, given sufficient rime, will break through without much oil displacement. For certain conditions, which include relatively low flow rate, the interface is a straight line and its behavior is simple to calculate. At higher flow rates, the interface is unstable, and a numerical solution was programed for an automatic computer. In general, good agreement is shown between the fluid model and the computed results so long as gravitational forces have control. For a water drive at very unfavorable mobility ratio, many small water fingers appear. These viscous fingers are not controlled by the relatively small gravitational forces. When viscous fingering becomes the controlling factor, the mathematical model is oversimplified, and results do not check the fluid flow model. INTRODUCTION Present methods of reservoir analysis depend upon certain simplifying assumptions to obtain mathematical descriptions of practical use. Material-balance methods (Muskat2 or Tarner2) assume uniform fluid saturations in the entire reservoir, or in a few subdivisions of the reservoir. An unsteady-state flow calculation by West, er at considered pressure and saturation changes in flow to a well during solution gas drive, and neglected gravity effects. Results showed only a 4 per cent difference for ultimate oil recovery by the Muskat method, even though the case chosen for study was one in which unsteady-state effects should be high. The Buckley-Leverett5 method commonly assumes a one-dimensional flow system. It is applicable at high flow rates where viscous forces predominate over gravity forces. Simultaneous, parallel flow of the two fluids is assumed, and the concept of a fluid interface is not introduced. Permeabilities to each fluid for a given saturation must be known. The method is not applicable for a two-dimensional system where cross flow becomes possible. Less well known is the displacement equation derived by Dietz. This method is designed for two-dimensional flow systems and assumes a definable fluid interface within the porous medium. Dietz showed that, for a range of low flow rates, the interface would be stable. straight and at an angle of inclination which could be simply calculated. At a certain critical flow rate, the calculated interface tilt would equal the formation dip. For higher flow rates, a finger of displacing fluid would invade the displaced fluid. Dietz indicated that his method applied only to macroscopic reservoir behavior, while the Buckley-Leverett method applied to the small transition zone at the fluid interface. The examples worked out in this report are based on the fluid-displacement theory of Dietz. It is shown that the Dietz theory may be used to derive equations analogous to the Buckley-Leverett equations. In contrast to the Buckley-Leverett method, flow is considered in a plane rather than being limited to a line. Rather than a frontal advance, the movement of a fluid interface is followed. For flow rates substantially exceeding the critical rate and for high viscosity ratio, many fingers of invading fluid occur-—rather than the single finger assumed by Dietz. On the other hand, so long as some gravitational influence remains, the flow is not entirely parallel to the bedding planes as assumed by Buckley and Leverett; therefore, both methods fail to give an adequate descrip-

By William Hurst, A. F. van Everdingen

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst' shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.' The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy's equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By A. F. van Everdingen, William Hurst

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst' shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.' The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy's equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By R. F. Steigerwald

The effects of heat transfer on the corrosion behavior of type 304 stainless steel in boiling water have been studied. Heat transfer conditions increase the tendencies of the stainless steel toward stress-corrosion mucking when the water is contaminated with Cl-. Surface preparation is the most important variable in determining the severity of the stress-corrosion problem in water with a given Cl- content. Heat treatment, chemical cleaning, and degree of wet film boiling also affect the corrosion of stainless steels used as heat transfer surfaces in boiling water. An up-paratus for corrosion testing under controlled heat transfer conditions is described. THE principal objective of this investigation was to study the corrosion behavior of AISI Type 304 stainless steel in boiling high-purity water under conditions of heat transfer. The principal variables in the study were: a) boiling conditions: limited boil, agitated boil; b) metallurgical structure: annealed or sensitized; c) surface condition: rolled, pickled, or ground. The results of the study should elucidate the role of heat transfer in the corrosion of stainless steel in boiling high-purity water and provide a base line from which to assess the severity of this type of corrosion. Although there has been an investigation of the effect of heat transfer on the corrosion behavior of stainless steels in alkaline waters,' work on the corrosion of stainless steels under heat transfer has been largely confined to various acid solutions.'-4 Furthermore, the study of the corrosion of stainless steels in high-purity water has been generally restricted to auto-clave and instream testing in superheated (e.g., 300°C) water and steam.5'6 Cooling water problems led to this study of the corrosion behavior of type 304 stainless steel in high-purity water at its atmospheric boiling point. In addition to the effects of heat transfer, this study also considered the influences of the type of boiling,7 metallurgical structure, and surface condition on the corrosion behavior of stainless steel in high-purity water. Two boiling conditions were investigated in this study: a) limiting boiling, i.e., when the heat-transfer surface had just reached the boiling point and boiling nucleated randomly over the specimen surface and b) agitated boiling, i.e., when a mild degree of superheat had been achieved on the specimen surface and boiling occurred generally over the specimen surface. Both boiling conditions are broadly classified as wet-film boiling. Originally, this program included the study of a third type of boiling, dry-film conditions. At dry-film boiling the temperature of the heat transfer surface is high enough so that it is covered with a continuous film of steam. Experimental difficulties made it impossible to study this kind of boiling with the available apparatus. Since surface condition affects heat transfer rates, three typical surfaces were chosen for study: a) as-rolled, the 2B finish received from the supplier, b) as-pickled, 30 min in 8 pct HNO3-1 pct HF, and c) as-ground on 120 grit paper. In order to evaluate the effect of improper heat treatment or welding on hot wall corrosion, both annealed and sensitized mi-crostructure were included in this study. EXPERIMENTAL Test Apparatus. The particular problems of interest in this investigation required a knowledge of the thermal gradient through the test specimen. From this gradient, the heat flux to the surface and the surface temperature could be determined. In order to satisfy this need, an apparatus similar to the boiling-disk apparatus of Fisher and whitney2 was designed. The solution is contained in a thermally-insulated, steel pipe lined with TEFLON resin. A stainless steel _______;________________________________________

Jan 1, 1970

By R. G. Davies, R. H. Richman

Several techniques, X-ray diffraction, transmission electron microscopy, conventional, metallography, and mechanical property measurements are combined to investigate decomposition in the suPersaturated solid solution Fe-27 at. pct Be. The following sequence of structures for aging temperatures below 400°C is established: massive short-range order, pseudo-FeBe (L2o) order, Fe3Be (DU3) order, sinusoidal modulation of composition, general or continuous precipitation, and discontinuous or cellular precipitation. In addition to this sequence itself, several features of the individual stages are noteworthy; for example, when coherent modulations occur, their development is parallel to (100) everywhere in the specimen apparently without incubation, and both the wavelength and the composition amplitude of the modulations increase with aging time. Two types of discontinuous reaction are observed, depending upon the aging temperature. At low temperatures (<550°C) the boundary-nucleated variety of discontinuous precipitate predominates; but at higher temperatures (>600°C) the crystallographic variety occurs instead. Hardness and flow-stress measurements clearly reflect the quoted structure sequence. Some qualitative aspects of the structure changes and of the relationship between mechanical properties and structure are discussed. BERYLLIUM in solid solution markedly strengthens iron.&apos;" A recent study demonstrated that not only were supersaturated single crystals of Fe-27 at. pct Be very strong, but they deformed primarily by mechanical twinning.3 Such behavior, coupled with the decrease in solid solubility of beryllium in a Fe from about 33 at. pct at 1165°C to about 10 at. pct at 500°C4 as shown in Fig. 1, leads naturally to questions about property changes with aging in this system. In 1929 roll,1 and later Seljesater and Rogers,5 reported that useful increases in strength could be obtained by aging quenched Fe-Be alloys, but it was not until quite recently that any of the details of precipitation in these alloys were studied. Richman6 observed that both an Fe3Be superlattice (with a DO3 crystal structure the same as Fe3A1, illustrated in Fig. 2) and a modulated structure preceded the appearance of a discrete second phase in the 20 at. pct Be composition. Here, a modulated structure refers to a periodic composition fluctuation within one crystal, while its lattice remains coherent. Subsequent work on a Fe-Be alloys by Heubner7 and von Heimendahl and Heubner8 confirmed the existence of the modulated structure, although only an FeBe (L20) ordered structure was detected (also shown in Fig. 2). Both Rich-man and Heubner corroborated the large increases in hardness during aging reported earlier, but other mechanical properties were not examined. Since modulated systems are encountered relatively infrequently, the influence of coherent modulation upon mechanical properties is not well-known. With the accumulation of information about age-hardening systems, it has become clear that many precipitation reactions cannot always be described by the concept of heterogeneous nucleation at some lattice discontinuity followed by diffusional growth. Rather, complex sequences involving prior "preparation" of the matrix are quite common; that is, a transition structure occurs in the parent phase, instead of in some "nucleus" form. Delineation of such subtle sequences appears to be a worthwhile step toward the understanding of solid-state precipitation. Since not many structural progressions of comparable nature have been reported, the objectives of the present survey are to identify and elucidate the various

Jan 1, 1967