Search Documents

By J. Szekely, V. Stanek

Steady state heat transfer and fluid flow rates are predicted for the counter-current contacting of gas and solid streams. For a .fixed pressure drop across the bed the convective heat transfer rate is affected by the gas flow rate, which in turn is affected by the temperature distribution in the system. Variations in bed voidage and in particle size were found to change the predicted temperature profiles quite markedly, in particular when the thermal capacities of the gas and solid streams were equal. IN Part I we examined the effect of step and pulse changes in the inlet gas temperature on transients in the flow of gases through fixed packed beds. Here we shall discuss heat flow and fluid flow interactions in the counter-current heat exchange between packed solids and a gas stream. Counter-current heat exchange in packed beds has been extensively studied. In his now classical paper, Furnas' assumed constant property values and developed a simple analytical. expression for the temperature profiles in the gas and in the solid stream. A highly sophisticated analysis of counterflow systems, allowing for transients and interparticle conduction was developed by Amundsen. In preparing a very readable review of this topic Elliott' also presented a number of elegant solutions, including the case involving heat generation. A great deal of work has been done on the application of these concepts to the modelling of the blast furnace. However, in all this previous work the mass flow of the gas was regarded as a quantity independent of the gas-solid heat exchange process. In the following we shall consider counter-current heat exchange between a gas stream and a packed bed of solids, under conditions where the pressure drop across the bed is specified, rather than the mass flow of the gas. FORMULATION Consider a packed bed, extending from y' = 0 to y' = L, in which moving solids and a moving gas are contacted in a counterflow arrangement. Let cold solids, with an initial temperature Tsi and mass flow GA, enter at y' = L, and hot gas with an initial temperature of Tgi and mass flow Gg enter at y' = 0. In contrast to the problems considered in Part I this system has a nontrivial steady state and we shall restrict our attention to the study of the pressure drop flow relationships and the resultant temperature distributions for steady state conditions. Under these conditions the formulation of the problem requires the statement of the following: 1) heat balance on the gas 2) heat balance on the solids 3) the mechanical energy balance equation and in addition the knowledge of the equation of state, and a functional relationship between viscosity and temperature and an expression for the heat transfer coefficient in terms of the other parameters. On neglecting temperature gradients within the solids and for a constant mass flow of the gas, the solution to this problem is well known.'" The expressions for the temperature profiles of the gas and solid streams may be put in the following form: where Tg and T, denote the dimensionless gas and solid temperatures respectively, y is a dimensionless distance, and R is the ratio of the thermal masses. These quantities are defined as follows: For the particular case of equal thermal mass capacities by the gas and the solid, that is R = 1, the Eqs. [1]and[2] simplify to Eqs. (71 and[8], respectively: i.e., the resultant profiles are linear. As for most practical cases , solutions given by Eqs. [7] and [8] will contain only the inlet temperatures as parameters. For a fixed, predetermined value of C;b the above

Jan 1, 1970

By A. A. Mathews

During the past year (1967-68), there have been many instances of either partial or complete failure of a tunnel-support system. It looks like the art is in a pretty sad state, but as long as we permit the field of tunnel design to remain basically an art, that's about what we can expect. The application of scientific principles to tunnel design involves much more than assigning values to a series of factors in an equation. I am reminded of the case where the engineer calculated that steel ribs with struts were required to support a horseshoe-shaped tunnel in firm ground. The contractor had driven several hundred feet without installing the struts when the engineer visited the job. When he was informed that the struts had been left out, he fled the tunnel in alarm. As one authority has so aptly put it, "There is far too much calculating and too little observing." Initially, one might think of a tunnel support as a structural member, such as a strut or arch which is installed in order to support some definite load. Similarly, a tunnel lining might be considered as a treatment of the interior surface of the tunnel in order to produce desired aesthetic, hydraulic, or architectural results, or to serve a protective function. Actually, many linings, such as a concrete lining, serve a dual purpose and support a load as well as satisfy other requirements. Not too many years ago a tunnel, during the driving or excavating stage, was considered to be either supported or unsupported. In recent years, however, several techniques have been developed, which by helping the tunnel to support itself, permit the elimination of separate structural support. Since these techniques in many cases actually become a substitute for the conventional types of supports and linings, a state-of-the-art review would hardly be complete without including them for consideration. Thus, for the purpose of this discussion, we shall define tunnel supports as any system or technique which effectively maintains the integrity of the tunnel opening. In their chapters, both Jacobs (No. 6) and Lippold (No. 7) pointed out how support systems have advanced from timber supports to steel supports

Jan 1, 1970

By S. P. Ray, B. D. Sharma

THE activation energy values for self-diffusion in pure metals have been correlated with a number of physical properties such as melting points,1 valences,' Debye temperatures,3 and cohesive energy densities.4 Extending the melting point correlation, the self-diffusion activation energy values for the alloys can be similarly correlated with the solidus temperatures. For metals and alloys which exhibit phase transformations, such a correlation is not feasible because of the uncertainty of various physical properties of the corresponding phases. Using thermodynamical approach, Ardell5 calculated the hypothetical melting points of the low melting phases of the polymorphic metals. However, these calculated melting points do not show reasonable correspondence with the activation energy of self-diffusion. Recently, a new approach treating metals and alloys as electron concentration phases has been developed by Engel6 and Brewer.7 These workers have attempted to correlate the crystal structures with the electron configurations of the atoms existing either in ground state or in an excited state of low promotion energy. The basic postulate of this theory is that whereas the contribution to the bonding of the metal comes from the interaction of both d and (s-p) electrons existing either independently or in hybridized states, the existing crystal structure is controlled by the (s-p) electrons alone. Specifically, a metal is supposed to crystallize in bcc, hcp, or fcc phase provided the metal has one, two, or three valence electrons per atom in the (s-p) configurations. The diamond, or ZnO and ZnS structures, occur with four (s-p) electrons per atom. Secondly, the contribution of d electrons to bonding is supposed to increase with atomic number

Jan 1, 1970

By L. W. Sink, B. H. Kear, S. M. Copley, M. F. Hornbecker

IT is well known that under appropriate conditions of chill casting, many metals and alloys develop a columnar grain structure having a high degree of preferred orientation.' In cubic materials, this principle has been employed to produce shaped castings with a unidirectional cube oriented columnar grain structure as well as shaped single crystals with a cube axis parallel to the growth ais.' "he purpose of this note is to describe how this same principle can be enlployed to produce wires, rods, or more complex crystal Shapes with a variety of conlpositions in cube or random orientations by a single casting operation. The procedure for directional solidification of in Fig. 2. Examination under crossed polarizers indicated that the striations were a second phase (presumably the B19 phase) which frequently forms initially at primary grain boundaries and then grows toward the interior as thin ribbons. X-ray studies on filings of this alloy reannealed at 1300°C for 15 min indicated a strong pattern of the fcc phase together with a substantial pattern of the B19 phase. The un-annealed filings, however, showed only the pattern of a cold-worked, fcc structure. In these two-phase alloys it appears, therefore, that the B19 phase is transformed into a disordered fcc structure by cold-working. The mechanical properties of these alloys would be of some interest since a fcc structure might be expected to possess greater ductility than the or-thorhombic B19 structure from which it forms. The author wishes to thank Mr. Harvey Yakowitz and Mr. Richard Manuszewski for obtaining the photomicrograph in Fig. 2. This research was partially supported by Research Grant DE-02455 from the National Institute of Dental Research. columnar grained ingots involves pouring a melt into a preheated mold having a water-cooled copper chill at its base? and controlling the temperature gradient in the mold so as to maintain columnar grain growth from end to end of the ingot.4 This is illustrated schematically in Fig. l(a). It should be noted that the longitudinal temperature gradient promoting unidirectional columnar grain growth is maintained by heat loss primarily through the solidifying ingot to the chill, so that both the temperature gradient and the growth rate diminish with increasing distance from the chill. As the distance from the chill increases? the columnar structure coarsens, and the perfection in preferred orientation tends to increase progressively. Using this same procedure of gradient solidification, Fig. l(b) shows a simple scheme for producing a batch of crystals having the same cube orientation rather than an assemblage of columnar grains. Wires and rods are produced by cementing ceramic tubes to the wall of the mold. The steep temperature gradient produced in the central chamber is controlled by the mold heater to produce favorable conditions for crystal growth. Clearly the success of the operation lies in the selection of individual growing columnar grains by some kind of orifice at the base of each channel in the mold. Using simple constrictions, experience with nickel-base alloys has shown that 4 -in.-diam holes at a distance of -1 in. from the chill yield crystals within 7 deg of (100); closer control of orientation is possible by placing the crystal selector farther away from the chill.

Jan 1, 1970

By William G. Gates

INTRODUCTION The Heath Steele mine is located in the Bathurst-Newcastle area of northern New Brunswick. It is jointly owned by American Metal Climax, Inc. and the International Nickel Company of Canada who hold a 75% and 25% interest, respectively. The ore deposits of lead, zinc and copper sulphide occur in very fine grained pyrite and/or pyrrhotite masses. A consistent spatial relationship with a quartz feldspar porphyry is characteristic of these deposits and the many other similar sulphide deposits in the district. The two principal mineral assemblages of pyrite-sphalerite-galena, and pyrrhotite-chalcopyrite-pyrite, occupy characteristic positions relative to the enclosing rocks. This feature and the nature of the sulphide mineralization are described and the variations in metal content are noted .with special emphasis on the lead-zinc ore zones. The mine product ion rate is 3,000 tons per day of a mixed lead- zinc and copper ore which is extracted from the B-1 zone. The complex nature of the ores presents many metallurgical problems. Optimum metal recoveries are achieved when milll feed grades are relatively consistent and proportions of the contained metals are with in acceptable limits. Ore outlines and stoping blocks are delimited by close spaced diamond drilling. From this information it is possible to forecast the grade of any portion of a stoping block. Production from the various stopes may then be scheduled so as to achieve an acceptable mill feed grade.

Jan 1, 1970

By Martin E. Chenevert

The early work of Terzaghi1 showed that the true stress state of a rock can be expressed by the principle of effective stress. This principle states that effective stress is equal to total applied stress minus the pore pressure. Since the introduction of this principle, a great effort has been made to measure the pore pressure of rocks so that the true stress state can be determined. For clean, permeable sandstones, it is usually found that the effective stress principle can easily be applied. Direct measurement of the pore pressure is achieved with a hydraulic pressure device and the effective stress state thereby computed. However, for fine-grained rocks which contain large amounts of clay (referred to herein as shales), pore pressure measurements are usually quite difficult because of the adsorptive forces of the clay platelets, which produce negative internal pressures. No suitable device is known for measuring internal pressures of dry shales; the effective stress approach is therefore not feasible and data are usually presented in another form. As an example, data from compressive strength tests are presented as a function of rock wetness, rather than effective confining stress. The work of Belopoloskii2 and Colback and Wiid3 are examples of such data. It was pointed out by Lambe and Whitman4 that a better way must be found to express the effective stress of argillaceous rocks. The information presented herein is an attempt to fill that need. Various authors have been able to measure moderate negative pore pressures (0-240 psi) within rather wet, argillaceous rocks by means of various pore pressure devices. In 1961 Lambe5 resented information which shows that pore pressures in argillaceous rocks are indeed negative ; he presented values of suction pressure as low as - 14 psi. Likewise,

Jan 1, 1970

By Edward J. Felten

Binary Cr-A1 alloys containing from 2.5 to 30 wt pct Al and 0.7 wt pct Y were heated in oxygen, air, and nitrogen between 1000" and 1200°C. The reacLivity of the alloys was found to be dependent both on the alloy composition nnd the nature of t he atmosphere. In oxygen, nllojs containing up to 15 to 20 wt pct A1 reacted to produce an external scale of Crz03 and a subscale consisting Predominently of Al203. Alloys contazning 20 to 30 wt pct A1 react in oxygen to produce an A1203 external scale and little m no subscale. The latter alloys were markedly more oxidation resistant than those of low alurninum content. In air, the alloys on which an external Crz03 scale was formed were found to be permeable to nitrogen ns evidenced by the copious amomts of chromium and aluminum nilrides observed ns part of the subscale. The reactizities in nir (or nitrogen) of these alloys increase <m their aluminurn contents increase. However, alloys on which Al,O, us an external scale is formed were nol culnerable to nccelerated attack in air, and no eltldence of nitvide subscnles were observed. For all alloys, yttrium serwed pYimarily to improve oxide adhrence. THE role of chromium in the oxidation resistance of Fe-Cr alloys &apos;-&apos; and that of aluminum in Fe-Cr-A1 al10s&apos; has received considerable attention in recent years. This is understandable since many of these alloys have excellent oxidation resistance due to the formation of either a Cr203 or a-Ala03 film between the metal and the oxidizing atmosphere. Small additions of yttrium or other rare earth metals are effective in preventing spalling of the protective oxide from the metal substrate."" In contrast, little is known regarding the oxidation resistance of Cr-A1 alloys, although some work has been done by Tumarov et a1.&apos; The poor niechanical properties exhibited by Cr-A1 alloys make them undesirable for use as structural components, but their use as coatings cannot be disregarded. The use of chromium-rich aluminide coatings for refractory metal alloys is an example of the potential use of this type of sytem. The purpose of this work is to examine the oxidation behavior of Cr-A1 alloys containing 2.5 to 30 wt pct A1 and 0.7 wt pct Y. The effects of temperature, atmosphere, and thermal cycling have been determined. EXPERIMENTAL PROCEDURE The alloys used in this investigation can be divided into two groups. Those containing 2.5, 5, 7.5, and 10 wt pct A1 and 0.7 wt pct Y were extensively evaluated in the temperature range from 1000" to 1200°C. Alloys containing 15, 20, 25, and 30 wt pct A1 and 0.7 wt pct Y were tested only at 1200°C. All of the alloys were prepared by standard arc-melting techniques in the form of cylinders approximately 4 in. long and 19 in. in diam. Wafers were cut from the cylinders and subsequently subdivided into rectangular coupons. The alloys were brittle and therefore some cracks were found in almost all specimens. The coupons were prepared for oxidation by mechanically polishing through 600 grit Sic paper, and were thoroughly degreased just prior to testing. Two types of oxidation experiments were conducted, namely; cyclic tests in which the specimens were examined and weighed after each 2 hr exposure, and continuous thermal balance tests run in a controlled atmosphere (oxygen, air, or nitrogen) for 20 hr. In the former test the spalled oxide was not included when the specimens were weighed. The physical condition of a specimen was noted visually after each cycle and testing was continued either to failure or until the performance of the specimen was well characterized. Both Micro and Semi-Micro Thermal Balances (Ains-worth) were used in the continuous tests. The oxidized specimens were sectioned and prepared for metal log raphic examination. The specimens were polished through 600 grit Sic paper. After polishing through 6 and l p diamond, a final mechanical polish with Linde B-Alz03 was used. Specimens containing 2.5 pct A1 were etched electrolytically using a 10 pct oxalic acid solution at 4 v for about 2 sec. Selected specimens were examined in the electron microprobe analyzer. Oxide specimens were examined by standard X-ray diffraction techniques. EXPERIMENTAL RESULTS For convenience, the test results have been broken down according to the exposure temperature, and further subdivided according to the type of test and atmosphere employed. Because of the poor quality of the specimens a larger than normal amount of scatter was observed in the measured rate constants. Also, the evaluation of the weight gain data was done on a somewhat arbitrary basis and may not be truly representative. However, the results obtained do show a significant trend in behavior regarding both alloy composition and the nature of the oxidizing atmosphere. I) Oxidation Behavior at 1000°C. A) Continuous Oxidation estsin Oxygen. This series of experiments was run in the Ainsworth Micro-Thermal Balance using pure oxygen at a pressure of 76 mm Hg. Under these conditions all specimens oxidized in accordance with the parabolic rate law over a major portion of the exposure time; the rate constants appear in Table I. The oxide formed externally on all specimens was predominantly Cr,O,, which was generally adherent. In some cases a slight amount of spalling in the form of a fine powder was noted. a-A1203 was observed as a subscale, along with Yz03 in all alloys. Alloys containing up to 7.5 wt pct A1 oxidize more rapidly than the Cr-0.7Y alloy.

Jan 1, 1970

By R. O. Williams

Analytical representation of the thermodynamics of solutions is highly desirable from the standpoint of accuracy, compactness, and numerical manipulations. In particular, computer calculations are greatly implemented. Mathematical considerations show that previous expressions have one or more serious defects. This investigation shows a Fourier series to be satisfactory but that it is also possible to derive a new series which fits certain additional conditions. Included examples show the value of analytical expressions in giving a simple characterization of each system using some two to five parameters, the elimination of the Gibbs-Duhem integration, and the es timation of the error for the experimental function as well as derived functions. It is further shown that the present characterization provides easy comparison between systems. IN the past, thermodynamic calculations have depended to a considerable extent on tabular and graphical methods. As the volume and precision of such data increase such methods become less satisfactory. Specifically, the selection of the optimum representation and the estimation of errors require statistical methods which in turn require analytical representation. The utilization of such data require further manipulations which are best done analytically for maximum precision. For example, phase equilibria are determined by common tangents to free-energy curves: a graphical determination is normally of low accuracy. As computers are increasingly used analytical representations become almost mandatory. Insufficient mathematical consideration has been given previously to the selection of empirical expressions. Those expressions having some theoretical justification are generally too inflexible and mathematically unattractive. We consider the problem in some detail and show that a Fourier series can be effectively used. Also a new series is defined which has certain advantages. ANALYSIS We wish to consider the analytical representation of the heat of mixing, AH, the excess free energy, ?Gxs, and the excess entropy, ?sXS, as a function of composition, X, for binary solutions relative to the pure components in the same state. When a distinction is not required, we use W to denote any one of the above functions. One may use a Taylor expansion around X = 0 to generate a power series. As the derivatives are un- known we represent the series as W = A + BX + CX2 + DX3 + EX4 + ... [l] where the constants A , B, C , ..- are to be selected to provide some optimum fit. For the extremes of composition W is necessarily zero so it follows that A = 0 [2a] B +C + D + E +••• = 0 [2b] Nonelectrolytes, which we are considering, appear to satisfy the condition that d3W/dx3 = 0 [3] in the terminal regions. This is the basis of the a, ß, and Q functions used by Hultgren et al.&apos; and others. While this condition does not have a strong theoretical basis it does appear desirable that any analytical relation should satisfy this condition. Darken2 and Turk-dogan and Darken3 have shown that many systems exhibit this behavior over an extended range from each terminal region, departure being restricted to a limited intermediate region. Since we have no a priori knowledge as to where this transition occurs we can require that this condition be satisfied only as a limit at the extreme compositions as a general condition. We will show later how more restricted conditions can be included in specific solutions. Darken2 has called this behavior the quadratic formalism; we call our application the limiting quadratic formalism, LQF. This condition applied to the above power series requires that D = 0 [4a] 4-3-2E +5-4-3_F + 6 • 5 . 4G + ••• =0 [4b] The form of the power series normally used, due to Margules,4 is W=X(1-X)(A + BX + CX2 + DX3 + EX4 + •••) [5] where A, B, C, --. are a new set of constants. (Guggenheim5 has given a variation of this expression in a more desirable form. Since, however, it is contained in the above expression it does not require separate consideration.) This form is precisely what results by incorporating the conditions in Eq. [2] into the power series and regrouping the constants. The LQF requires that B =C [6a] and 4.3.2(D-C) +5-4-3(E-D) + ••• =0 [6b] Thus, the correct form of the Margules expression with two adjustable parameters is w =X(1-X)[A + B +X2-2/3x3)] 171 and the EX4 term must be included before three adjustable parameters are permitted.

Jan 1, 1970

By A. Timur

Fluid flow properties of porous media have long been of interest in such varied disciplines as geology, geophysics, soil mechanics, and chemical, civil, and mechanical engineering. This interest has resulted in numerous models of porous media that have been proposed, tested, and found to be useful under, at best, only special circumstances. A critical review of most of these models is given by Scheidegger. 1 The application of each new tool to the study of porous media has helped to test some of the existing theories and to form new ideas based on the new parameters being measured. In the application of nuclear magnetic resonance (NMR) techniques to the study of the properties of fluids in porous media, the theoretical studies of Korringa et al.2 resulted in a model for the relaxation of spin polarization of protons in a hydrogenous fluid in a pore of a solid. Seevers3,4 verified this model, established a method for measuring the surface-to-volume-ratio distribution in a porous medium, and proposed a technique for determining specific permeability in sandstones. In the present investigation, the uses of NMR methods are extended and a simple model of porous media is developed from an analysis of pulsed NMR measurements. In this model, the pore spaces of a porous medium are divided into three groups on the basis of their surface-to-volume ratios. To evaluate predictions of the NMR model, laboratory measurements of spin-lattice relaxation time, porosity, permeability, and residual (irreducible) water saturation were conducted on more than 150 sandstone samples obtained from four oil fields in North America. Applications are discussed for estimating the volume of movable fluid, ?p, and the specific permeability, k, of sandstone samples. The former parameter, which can be considered as producible porosity, is expressed by ?p = [1 -- (Swr/100)]? , .... (1) where + is porosity in percent of bulk volume and Swr is the residual (irreducible) water saturation in percent of pore volume at a capillary pressure of 50 Psi. Fast, nondestructive laboratory methods are described for determining porosity, movable fluid, and permeability of sandstone samples from measurements with a pulsed NMR (spin-echo) apparatus. Procedures are given for adapting these methods to samples of irregular shape and small size. Theory General In the Korringa, Seevers, and Torrey2 (KST) model the observed decrease in the spin-lattice relaxation time, T1, of protons of a hydrogenous liquid contained in the pore spaces of a porous solid is attributed to an increase in the correlation time for the random motion of the water molecules and to the presence of paramagnetic centers at the liquid-solid interface. This KST model predicts the observed relaxation time of a hydrogenous liquid in a pore of a solid through

Jan 1, 1970

By D. S. Pye, J. P. Gallus

Successful hydraulic pressure treatment of oil and gas wells that have long producing intervals, multiple pay zones, fluid thief zones, or natural or induced fractures has long been a serious problem plaguing the petroleum industry. High fluid-injection rates, straddle packers, ball sealers, and either liquid or solid additives to the treating fluids are being used to cope with this difficult and costly problem.l-6 The use of ball sealers and straddle packers is limited, in practice, to wells that have been completed with perforated cemented casings or liners. Well stimulation with ball sealers and straddle packers is usually successful as long as the cement remains competent. The other techniques of diverting are used with all types of completions. However, there are two significant disadvantages of the high rate method: excessive horsepower requirements, and the need of large treating-fluid volume. Consequently, most well stimulation techniques that involve diverting depend upon ball sealers or granular materials. Liquid diverting agents, and some of the solid ones in use today, frequently will cause formation damage. On the other hand, the most popular particulate diverting agents often are so soluble in the treating fluid that they become largely ineffective; their particle size is severely reduced as they dissolve while being transported through thousands of feet of pipe to the formation face down hole. In our attempts to deal with the shortcomings of these diverting agents, we experimented with a large variety of potentially suitable deformable materials the solubility of which could be controlled. Following are the criteria for the diverting agent sought: 1. The material must be effective in reducing fluid loss to thief zones, including existing fractures and very high permeability layers in the formation to be pressure treated. 2. The material must have enough strength to bridge under the treating pressure. 3. The material must be insoluble or slowly soluble in the treating fluids so that the particles retain sufficient size to be effective after they reach the formation face. 4. The material must be completely soluble in the produced oil or condensate liquid so that it can be removed readily; it must not cause formation damage. A number of products consisting of various polymers blended with waxes and other solid hydrocarbons were examined. All these materials satisfied the criteria for a diverting agent and also exhibited some degree of plastic deformability. On the basis of initial tests of strength and hardness, several blends of wax and polymer were selected for extensive testing in the laboratory and the field. These compounds have congealing temperatures ranging from 140 to 190F; they are insoluble in water and soluble in oil. They conformed, therefore, to the criteria specified for a theoretically superior diverting agent. Laboratory and field tests were conducted using these wax-polymer materials to demonstrate that the use of a deformable diverting agent can significantly improve fluid diverting performance. Data in this

Jan 1, 1970

By C. F. Weinaug, P. M. Blair

Many difference equations used to approximate reservoir flow problems treat the phase pressures implicitly but not the mobility-density coefficients. Such difference equations are neither wholly explicit nor implicit, but might be described as mixed. Mixed equations are relatively easy to apply. But the associated time truncation error is relatively large, and when used to solve problems characterized by high flow rates, these equations may be unstable for practical size time steps. This paper outlines the development of a completely implicit difference analogue for reservoir simulation, along with a Newtonian iterative method for solving the resulting nonlinear set of algebraic equations that arise at each time step. While this implicit equation requires two to three times more work than does the mixed equation, it is shown to markedly decrease the time truncation error and to yield a stable solution for much larger time steps than does the mixed equation. INTRODUCTION Calculation of multiphase, multidimensional flow in porous media is generally accomplished by numerical methods that involve approximating systems of partial differential equations by systems of partial difference equations. The early development of difference equation techniques was directed toward the solution of linear differential systems. However, the equations of multiphase fluid flow through porous media are highly nonlinear in that mobility and density often are strong functions of pressure and saturation. Thus, solution of the flow problems by difference equations involves solving sets of algebraic equations whose coefficients change from step to step of the calculations. Current practice is to evaluate most coefficients at the beginning of a time step and then to apply difference methods well suited for solving linear problems. At least two major difficulties arise from the use of this practice. First, evaluation of coefficients at the old time level and evaluation of pressures at the new time level results in larger time truncation errors near displacement fronts than if all quantities in the distance derivative were evaluated at the same time level. Second, evaluation of mobility coefficients at old time levels results in an unstable difference equation in regions of high flow rate. The first of these difficulties often results in optimistic recovery behind the displacement front, and the second requires the use of small time steps or special techniques in solving coning problems1 or problems of gas percolation.2 Evaluation of all quantities in the distance differences at the new time level results in a completely implicit difference system. Such a difference system has a lower time truncation error than equations in which mobilities are evaluated at the old level and pressures evaluated at the new level. The algebraic equations resulting from the fully implicit difference equation are, in this case, nonlinear and require some iterative method for their solution. This paper develops fully implicit difference equations for two-phase flow in porous media, describes their solution by Newtonian iteration, and gives examples of problems more easily solved by the new method than by previous methods. EVOLUTION OF DIFFERENCE EQUATIONS FOR NONLINEAR PROCESSES The completely implicit difference equation used in these calculations results from evolution beginning with explicit difference equations, proceeding through what we will call "mixed" equations, and arriving now at a completely implicit equation. Let us illustrate the evolution by considering different analogs of the nonlinear equation

Jan 1, 1970

By L. B. Davidson

The petroleum literature contains many reports on the relative permeability properties of porous media. However, only recently have studies of relative permeability at elevated temperatures been made.1 -5 Each of these studies presents relative permeability data on water-oil systems in linear sandpacks or sandstone cores under isothermal conditions. The results of these investigations are dissimilar. Edmondsonl found that in Berea sandstone a temperature-dependent effect was present, while Shilobod,2 using the same material, found no temperature dependence after normalizing for initial water saturation. Workers at Texas A&M U.3-5 found that the ratio kw/ko was temperature dependent in an unconsolidated system. The differences in results of these studies are probably a consequence of differences in the material employed, in the experimental procedures followed, and in the properties of the porous systems used. This paper describes an experimental investigation of permeability ratio temperature dependence, conducted to clarify some of the results obtained by previous workers. Each experiment described in this paper involved the isothermal displacement of Chevron No. 15 white oil by either nitrogen, steam, or distilled water and displacement of water by nitrogen. All permeability ratios were calculated using the Welge procedure. Values of the water-oil permeability ratio were obtained over the temperature range 75OF to 540°F. Results indicate the ratio is temperature dependent at low water saturations, probably due to interfacial effects. However, as the saturation is increased, temperature dependence decreases until, above some minimum saturation, the ratio becomes insensitive to temperature. Previously published work indicates that the dependence reappears at very high water saturations, a region of the ratio curve not completely examined in the present study. It appears that water-oil permeability ratios are temperature dependent at the extremes of the ratio curve where changes in interfacial tension with temperature induce variations in residual saturations. Between the extremes, there is a region in which the permeability ratio is insensitive to temperature. The size and location of this region is probably a function of the characteristics of the porous system. Gas-oil permeability ratios calculated from nitrogen-white oil displacement experiments indicate a definite temperature dependence over the range 75OF to 500°F and over the entire gas saturation range studied. Evidence discussed in this report indicates that this dependence is a consequence of molecular slippage in the gas phase. This is an extension of the Klinkenberg effect6 to two-phase systems at high temperatures (the existence of slippage in two-phase systems at room temperature has been demonstrated7). Evidence for a similar source of temperature dependence in superheated steam-white oil experiments

Jan 1, 1970

By Burt C. Mariacher

The Minerals Beneficiation Division meets each year at the Annual Meeting of AIME. All members are invited and encouraged to attend. For the information of those who did not attend the meeting in Denver last February, this column presents a summary of the reports given and topics discussed. Technical Committee Coordination Consolidation of technical programs having common interest to MBD and the Extractive Metallurgy Division is proving a considerable success, both on the technical committee level-through the efforts of the MBD-EMD Coordinating Committee-and in special symposia. MBD Chairman, R. S. Shoemaker, noted that J. D. Prater will be serving as Chairman of both the MBD Chemical Processing Committee and the EMD Hydrometallurgy Committee for 1970-71. R. P. Ehrlich commented on the cooperative effort of both divisions in organizing the Copper Metallurgy Symposium held at the Denver Meeting. This effort, he said, was important to the success of the Symposium. H. R. Spedden discussed inter-Division coordination activities on the SME level. The forthcoming MBD-M&E-EMD Joint Lead-Zinc Symposium is yet another example of progress in this direction.

Jan 1, 1970

By R. Stefanko, Y. J. Wang, H. R. Hardy, R. Y. Kim

Since 1964, the Dept. of Mining at The Pennsylvania State University has been carrying out extensive studies associated with the phenomena of microseismic activity and inelastic behavior in geologic materials. Briefly, the program to date has involved the following stages 1) Study of the microseismic noise and energy rates for selected geologic materials under incremental, uniaxial tension.1,2 2) Investigation of the frequency spectra of microseismic activity of selected geologic materials under incremental, uniaxial tension.3,4 3) Investigation of the time-dependent inelastic deformation of geologic materials and the analysis of such deformation in terms of mechanical models.5, 6 4) Simultaneous study of microseismic activity and time-dependent inelastic deformation in geologic materials under incremental, uniaxial compression.7,8 Details of the research involved in stage 4 of the program constitute the subject of this chapter. It should be noted that the majority of the experimental data quoted in this chapter forms part of a thesis to be presented by one of the authors (Kim) in partial fulfillment of the requirements for the Ph.D. degree in mining engineering at The Pennsylvania State University. During the last two or three years, a number of important papers have been published which deal primarily with the fundamental structure and mechanical behavior of geologic materials. Recent papers 9,10 have provided research workers in the rock mechanics field with a meaningful overall picture of rock as a structural material. The concept that many

Jan 1, 1970

By H. C. Gatos, J. T. A. Pollock, V. Sadagopan

PREVIOUS investigations on the effect of neutron irradiation on the superconducting properties of bulk transition metals indicated that little or no change occurred in poly crystalline disordered materials, particularly if cold-worked prior to irradiation. Only the superconducting properties of ordered intermetal-lic A15 compounds were found to be sensitive to irradiation, and the observed reduction in Tc and increased magnetic hysteresis or critical current densities were attributed to the introduction of flux pinning vacancy clusters and disorder in the lattice. Since the former would also be expected to be present after irradiation in disordered alloys, the general conclusion was that the disordering effect of the irradiation was the primary reason for the changes observed. The present note reports some unusual effects of neutron irradiation on cold worked Nb(Cb)-Ti-V alloys. Specimens with nominal composition as given in Table I were prepared by repeated arc melting, taking the usual precautions, and homogenized at 1000°C for 24 hr prior to a reduction of 15:l by cold rolling. Samples of size 1 by 1 by 25 mm were sliced from the rolled strips. The sarmples chosen for irradiating were sealed in a cadmium tube (wall thickness 1 mm) prior to exposure in a fast neutron flux of 3.7 x 1019 neutrons per sq cm in the swimming pool reactor of Union Carbide Corporation&apos;s Sterling Forest Reactor. During irradiation the temperature of the samples was maintained below 70°C. The irradiated samples were stored at room temperature for approximately 1 year allowing the radiation levels to reach a safe handling value. Tc measurements were made employing a self-inductance coil technique. Jc data at 4.2"Kwere obtained using a four lead resistivity method. The latter measurements were made at the Francis Bitter National Magnet Laboratory by setting the perpendicularly applied field and increasing the direct current through the sample until a potential of about 10-6 v was detected. Tc values for the nonirradiated samples are listed in Table I together with the change in Tc after irradiation. A decrease in Tc of 0.25 or 0.3"K was measured for three of the four samples, the exception being the Nb-30 at. pct V-20 at. pct Ti alloy which showed a decrease of 0.5K. Jc vs applied field data are plotted in Fig. 1 for the Nb-20 at. pct V-30 at. pct Ti alloy. These curves are typical of the four alloys and surprisingly indicate that a degradation occurred in all the samples after irradiation. Hc2 values, obtained by extrapolation to zero current, are 8 to 15 pct lower than those measured for the nonirradiated specimens. The reasons for the degradation of Tc and jc ob-served in the alloys on exposure to neutron irradiation are not immediately apparent. The temperature at which the samples were maintained during bombardment (<70°C) and stored (20°C) make it highly unlikely that any significant change in the dislocation density of the materials through thermally activated migration and annihilation would occur. The explanation must lie in the type of damage produced by the bombardment. McEvoy et a1. reported no change in superconduct-

Jan 1, 1970

By P. O. Metz, E. L. Wiehe, J. S. Foster

BASED on their analyses of growth kinetics, ibbs&apos; and Harris et al.&apos; have recently proposed that the rate of subcritical graphitization of steel is controlled by the rate of diffusion of carbon in ferrite. illert has concluded that this is the rate controlling growth mechanism when the growing graphite particles are relatively large. Harris et al.&apos; also propose that the matrix surrounding a growing graphite particle is plastically deformed to make room for the growing particle. If this plastic deformation does occur, then preferred recrystallization may be observed around the larger particles. If diffusion of carbon in ferrite is rate controlling, then a graphite particle growing in a matrix of, say, tempered martensite or pearlite, should be surrounded by a carbide depleted zone. Thus, the simple metallographic experiments reported here are sufficient to determine whether or not diffusion of carbon in ferrite is the dominant mechanism for the growth of graphite in steel. The present investigation considers graphitization in five alloys. The compositions are given in Table I. Graphitization was carried out at three temperatures: 900°K (627"C), 950°K (677"C), and 990°K (&apos;717°C). Samples were reacted in both the martensitic and the pearlitic initial conditions. In the hypoeutectoid steels, Sets 7 and 8, .and in the hypereutectoid Fe-C binary, Set 6, graphitization was very limited after 1000 hr at temperature. In these samples graphite particles observed were less than 10 p in diam. In the hypereutectoid commercial steels, Sets 1 and 3, large graphite particles (100 p) did not appear until after 250 hr at temperature. The limited number of graphite particles observed and the range of sizes obtained prohibited measurement of growth rates in this investigation. However, in all samples showing graphitization, neithe a carbide depleted zone nor recrystallization was observed. illert&apos;s model predicts that diffusion of carbon in ferrite should become dominant for graphite particles larger than about 300 p in diam. The largest graphite particle observed in this investigation was approximately 450 p in diam. No evidence of carbide depletion was noted. Fig. 1 shows a typical graphite particle. This nodule is from an initially pearlitic sample from Set 3 after

Jan 1, 1970

By G. L. F. Powell

g samples of Cu-20 wt pct Ag alloy have been mdercooled to a maximum of 197°C by melting under a slag of commercial soda-lime glass in a vitreous silica crucible. No grain refinement of the primary copper was observed in samples undercooled to the maximum of 197°C. When the samples contained a small amount of oxygen, the copper dendrites were partially recrystallized at undercoolings greater than 97°C. In previous papers&apos;-3 reporting the grain structure of undercooled silver and copper, it was observed that grain refinement was dependent on both undercooling and oxygen content. Grain refinement occurred in undercooled silver when the degree of undercooling exceeded the range 153" to 175"C, while in Ag-0 alloys (0.12 wt pct) fine equiaxed grains were exhibited when undercooling was greater than 50°C. Similarly, copper samples undercooled as much as 208°C displayed fan-shaped growth from a single nucleation site, while the grain structure of Cu-O alloys (0.08 wt pct) was fine and equiaxed at undercoolings larger than 150°C. Thus the presence of oxygen greatly reduced the undercooling at which grain refinement occurred. It was also observed that the change in grain size resulted from recrystallization and was not due to an enhanced nucleation rate in the liquid-solid transformation. It is possible that the influence of oxygen on recrystallization is due primarily to its presence as a solute element. walker4,&apos; reported that, although a grain size change did not occur in pure nickel until the undercooling exceeded 150°C, small grains were observed in samples of Ni-Cu alloy solidified at small and large degrees of undercooling. Jackson et al.6 suggested that the fine grained structure of the Ni-Cu alloy resulted from the melting off of dendrite arms during recalescence. This remelting process may occur in alloys as a result of segregation during freezing which causes a variation in liquidus temperature from point to point within a dendrite. It was therefore decided to undercool copper with a metallic alloying element to ascertain whether the presence of a metallic solute would have a similar effect to oxygen in inducing grain refinement. A Cu-Ag alloy was chosen, since both metals had been shown to behave similarly on undercooling. The alloy Cu-20 wt pct Ag was selected since the eutectic constituent outlines the initial growth form of the primary copper, so that the as-frozen grain structure is not obscured if subsequent recrystallization occurs. This paper describes the results of undercooling experiments carried out with Cu-20 pct Ag samples undercooled to a maximum of 197°C and the effect of oxygen content on the grain structure of the undercooled samples. EXPERIMENTAL Melting was carried out in a small cylindrical resistance furnace using "fine" silver granulate and oxygen-free high conductivity copper. The procedure adopted was to melt the required quantity of silver in air in a clean vitreous silica crucible for approximately 15 min, freeze, and add granulated commercial soda-lime glass to form a complete surface slag cover, after which the sample was melted and frozen several times to reduce the oxygen content. The glass slag cover was approximately 3 in. thick. Pieces of copper (=50 g) were added to the crucible until the required quantity to make 350-g samples of alloy had been charged. Each piece was added quickly to the crucible which was held at a temperature slightly above the melting point of silver. The piece was quickly pushed beneath the glass to minimize oxidation and any oxide coating usually decomposed before the piece had settled down into the silver. After the full quantity of copper had been added, the melt was stirred with a silica rod to hasten homogenization and a Pt/Pt 13 pct Rh thermocouple enclosed in a vitreous silica sheath inserted for temperature measurement. Heating and cooling curves were recorded on a potentiometric chart recorder fitted with a zero suppression unit. The milli-voltage range of the recorder was adjusted so that temperatures could be read to 1°C. Heating and cooling curves were taken every hour until three consecutive readings gave the same solidus-liquidus range, consistent with the solidus-liquidus range for this alloy composition by reference to Hansen and Anderko.7 Metallographic examination of samples frozen at this stage, failed to show any variation in composition from bottom to top of the ingot. Consequently, it was considered that the melt was homogeneous at this stage and undercooling experiments were then car-

Jan 1, 1970

By T. Reed Scollon

With some uncertainty as to the type, rate and scope of adoption of air pollution control regulations, and with limited information on the availability of low-sulfur fuels, it is possible at best to make only a rough assessment of the impact of such regulations on coal. It is apparent from the following discussion that regulations of the type already effective or being considered, if adopted and enforced on a widespread basis over a short-term period and if substitute fuels were available, both quantitatively and economically, could result in a loss of more than one-half of coal&apos;s markets. Depending on the areal sequence of adoption of such regulations, some areas would be able to obtain limited amounts of low-sulfur coals. Others would not. Because of the usually greater distances from the large industrial coal consuming areas to new sources of coal supply, and because of the higher mine prices for low-sulfur coals, the delivered prices of such coals as might be available would be much higher than those for high-sulfur coals. A large and sudden demand for low-sulfur coals would exert further upward pressure on prices.

Jan 1, 1970

By Mark J. Klein

The creep of 043 was studied over the temperature range 1650" to 3200°F and over the stress range 3000 to 44,000 psi. The steady-state creep rate over this range of stress and temperature can be expressed by the equation where A is a constant, is the stress, and is -0.8 x 103 psi-&apos;. Over a narrow range of stress variations c0 a and for this proportionality n varies from 3 to 30 in accordance with the relation n = aB. Above about 2400° F, H, the apparent activation energy for creep, is 110,000 cal per mole, a value about equal to that estimated for self-diffusion in this alloy. Below 2400°F, H increases with decreasing temperature reaching a value of -125,000 cal per mole at 1700° F. In this temperature region, H appears to be a function of the interstitial concentration of the alloy. MOST of the detailed creep studies of dispersion strengthened metals have been concerned with metals having fcc structures. However, there are a number of important refractory alloys with bcc structures that derive part of their high temperature strength from an interstitial phase and whose creep behavior has not been well defined. This paper describes the creep behavior of the bcc alloy, D43, over the temperature range 1650" to 3200°F (0.4 to 0.7 Thm) and over the stress range 3000 to 44,000 psi. In addition to colum-bium, this alloy contains 10 pct W. 1 pct Zr, and sufficient carbon (-0.1 pct) to form a carbide dispersion throughout the matrix of the alloy. The effects of variations in temperature and stress on the steady-state creep rate of this alloy are presented in this paper. EXPERIMENTAL PROCEDURES Creep tests were made in a vacuum of 106 torr under constant tensile stress conditions using a Full-man-type lever arm.&apos; Creep specimens were machined from 0.020-in. D43 sheet (grain size -5 x l0-4 in.) processed in a duplex condition (solution annealed -2900°F, 40 pct reduction in area, aged 2600°F). The specimens were tested in this condition without further heat treatment. Specimen extensions over 1-in. gage lengths were continuously recorded using a high temperature strain gage extensometer. Differential temperature and stress measurements were used to determine temperature and stress dependencies of the creep rate. Activation energies were calculated from the changes in strain rate induced by abrupt shifts in the temperature during constant stress creep tests. The 100°F temperature shifts used in most of the activation energy determinations required 15 to 90 sec depending upon the temperature at which the shift was made. The dependence of strain rate on stress was determined by measuring the change in strain rate for incremental stress reductions during constant temperature tests. It has been shown that columbium-base alloys such as D43 are susceptible to contamination by gaseous interstitial elements during vacuum heat treatments.&apos; In this regard, it is unlikely that these alloys can be heat treated without some loss or gain of interstitial elements despite the precautions taken to control the heat treating environment. However, several factors suggest that changes in interstitial concentrations of the specimens during testing did not affect the results presented in this paper. First, the dependence of the creep rate on the stress or temperature determined during the course of a single creep test showed no variations with the duration of the test. A variation would be expected if a loss or gain in interstitial concentration during the course of the test affected results. In addition, precautions taken during this investigation to minimize interstitial contamination by wrapping the gage lengths of the specimens with various foils2 (Mo, Ta, W) did not produce a detectable change in the stress and temperature dependencies relative to the unwrapped specimens. The averages of duplicate analyses for carbon and oxygen in several specimens determined before and after creep testing are listed in Table I. The combined nitrogen and hydrogen concentrations which were ordinarily less than 50 ppm did not change in a detectable way with creep testing. The analyses show that only minor changes in carbon concentration occurred during creep testing except for specimen 4. This specimen which was tested at 3100°F lost a significant amount of its carbon concentration to the vacuum environment. Specimen 1 gained 100 ppm of O, while specimens 2, 3, and 4, which were tested at progressively higher temperatures, lost increasing portions of their initial oxygen concentrations during testing. RESULTS AND DISCUSSION The Temperature Dependence of the Creep Rate. The apparent activation energy for creep, H, was de-rived from creep curves similar to that shown in Fig. 1. Steady-state creep was rapidly attained at the beginning of the test and with each change in temperature. This behavior suggests that the alloy rapidly attains a stable structure with each shift in temperature or that the structure is constant throughout the test. Since the dispersion will tend to stabilize the structure, the latter is probably the case. The activation energy was found to be independent of the direction of the temperature shift and the magnitude of the shift (50" or 100°F). Although H was approximately independent of the strain, there was a tendency for it

Jan 1, 1970

Jan 1, 1970