Search Documents

By N. J. Olson, G. W. Toop

A series of equations previously derived for calculating ternary thermodynamic properties using binary data has been applied to the problem of predicting ternary phase diagrams and quaternary excess free energy. The methods are considered to be rigorous for regular ternary and quaternary systerns and empirical for nonregular systems. The equations have been used to predict ternary phase boundaries in the Pb-Sn-Zn system at 926°K and the Ag-Pd-Cu system at 1000ºK. Calculated quaternary excess free-energy values are presented for the Pb-Sn-Cd-Bi system at 773°K. A method for predicting the location of ternary phase boundaries would be a useful supplement to experimental measurements in ternary systems. This has been recognized with the considerable work that has been done to find models to predict or extend thermodynamic properties and phase diagrams in binary and ternary systems1-18 for which direct experimental measurements are limited. With the access to highspeed digital computers and mechanical plotting devices, it is currently rather easy to compare mathematical models with experimental data. The regular-solution model is consistent with systems which exhibit negative heats of mixing, positive heats of mixing, and miscibility gaps, and therefore it is applicable to simple phase diagrams. The purpose of this paper is to illustrate the use of regular-solution equations to predict, empirically, phase equilibria in some types of nonregular ternary systems. Corresponding equations for regular quaternary systems are given and used to calculate empirical quaternary excess free-energy data. METHOD FOR PREDICTING THE LOCATION OF TERNARY PHASE BOUNDARIES USING BINARY DATA Meijerin1,6 has used the regular-solution model to calculate common tangent points to ternary free energy of mixing surfaces and hence to determine phase boundaries in ternary systems involving miscibility gaps. He used the following equation to calculate ternary excess free energy of mixing values: stants characteristic of the binary solutions, and Ni is the mole fraction of component i. An alternate expression which gives for regular solutions as a function of binary values of along composition paths with constant N1/N2, N2lN3, and N1/N3 may also be derived:15 ternary r xs 1 ?c-*n.Ti*.U*. This expression for is more useful for the empirical calculation of ternary excess free-energy values for nonregular systems because actual binary AFXS data may be used in the expression rather than attempting to find suitable constants for Eq. [I]. The results of this feature of Eq. [2] are illustrated in Table I where calculated excess free-energy values for the Ni-Mn-Fe system at 1232°K are compared with experimental data of Smith, Paxton, and McCabe.19 Although regular-solution equations have been shown to give calculated thermodynamic quantities which agree quite well with experiment for single-phase nonregular ternary systems,14,15 care should be exercised in the use of the equations to predict thermo-dynamic properties of multiphase ternary systems in which strong compound formation is suspected. This precaution is consistent with the simple regular-solution model which for negative values of ai_j will indicate a tendency toward compound formation but even for very large negative values of ai-jwill not give multiphase binary or ternary systems involving a distinct stable compound. Hence, calculated ternary free-energy data using Eq. [2] might be expected to vary between being rigorous and poor, in the following order, for ternary systems which are: a) regular solutions, b) nonregular single-phase liquids in which random mixing is nearly realized, c) nonregular single-phase solids, d) nonregular multiphase systems exhibiting miscibility gaps, e) nonregular multiphase systems with binary compounds but no ternary compounds, f) nonregular multiphase systems with highly stable binary and ternary compounds. The calculated data will be expected to be least accurate for the last two cases. The general method adopted in this paper involves two-dimensional plots of ternary activity curves. The principle used is that tie lines indicating two-phase equilibria join conjugate phases a and B for example, for which a1(a) = a1(B), a2(a) = a2(B), and a3(a) = a3(B). Tie lines may be determined by plotting the ternary activities of two components along an isoactivity line for the third component and the unique points where the above equalities hold may be found graphically.

Jan 1, 1970

By AIME AIME

THE portrait of Herbert Hoover, President of the United States, recently painted by Philip A. de Laszlo at the request of the four Founder engineering societies, was unveiled at a special meeting in the Engineering Societies Building on Feb. 15. The portrait was delivered by Gano Dunn, chairman of the Hoover Medal Board of Award, who said, in part:,

Jan 1, 1932

By R. L. Perrine

A stability theory is developed for miscible liquid-liquid displacement within a porous medium. In the usual case considered, a high-density high-viscosity "oil" is displaced downdip by a low-density low-viscosity "solvent'! Perturbation methods are used to find the conditions under which the spreading mechanism changes from the stable dispersion process to unstable viscous fingering. We find that instability is conditional, that there is a dependence on the shape of a disturbance leading to a"diameter" effect and that very difficult experimental scaling problems may result. A useful consequence is the definition of a minimum "slug size" for stable miscible displacement. This should make possible optimum use of the solvent process for oil recovery. The results may apply to many situations in which one fluid displaces another of somewhat different fluid properties within a porous medium. INTRODUCTION The process of miscible liquid-liquid displacement looked very favorable when it first received serious consideration as a means to increase petroleum recovery. Some early laboratory results were interpreted to mean that only a small "slug" of solvent was needed, perhaps 2 to 3 per cent of a hydrocarbon pore volume. This "slug" could recover the oil in an entire reservoir provided the fluids remained miscible.1,2 Only a slow growth of the "mixing zone", or gradual transition from oil to solvent, should occur. This would follow naturally from mixing by a dispersion mechanism such as that described by Scheidegger.3 And, indeed, laboratory cores have shown this kind of behavior provided solvent viscosity was at least as great as that of the oi1.2,4-6 The same kind of behavior has also been observed with solvent viscosity lower than that of the oil after the distance traversed became large. Not all laboratory results have been this favorable, however.7-9 Mixing zone growth may become proportional to the distance traveled rather than the square root of this quantity. Such behavior accompanies high rates in short systems of large diameter, provided the viscosity ratio is adverse. If this more rapid spreading were to persist in a reservoir, the solvent requirement for miscible "slug" displacement could exceed 30 per cent of a pore volume. There is considerable economic importance in the difference between 3 and 30 per cent solvent. A 1ogical conclusion is that two different kinds of flow behavior are possible, with each kind leading to a different result. One kind gives efficient displacement, the other does not. In the efficient case, the solvent bank spreads out only by the mechanism we have termed dispersion. Under these conditions, displacement is as near piston-like as possible. In the second kind of flow, found only with adverse viscosity ratios, viscous fingering occurs. That is, permeability variations cause a small finger of low-viscosity solvent-rich fluid to move ahead of its average position within the mixing zone. This creates a path of low resistance to flow, and an even greater amount of solvent-rich fluid follows. Thus, the process is autocatalytic. Once started, the fingering mechanism rapidly becomes dominant. The question to be answered is this. Under what conditions will each of these two different kinds of flow occur? In particular, what conditions denote the transition from dispersion to viscous fingering as the solvent spreading mechanism? The answer may tell us whether or not the desirable, near piston-like miscible displacement process is practical. An answer to this problem can be obtained from theory by the use of perturbation methods. The procedure is as follows. We first formulate a mathematical representation of the system. Then a small disturbance (or perturbation) in solvent concentration profile is introduced and observed to see what happens. Unstable flow is indicated when a small disturbance will grow larger. This will lead to eventual viscous fingering. If on the other hand the disturbance dies out, the displacement is stable. In this latter case, with dispersion as the spreading mechanism, near piston-like displacement is possible. Thus, this paper presents the development of a stability theory for miscible liquid-liquid displace-

By J. Hudy

An investigation has been made to determine the applicability of the Cendrex X-ray instrument for the measurement of the ash contents of washed bituminous coal products. The results obtained for selected products from five preparation plants in the northern Appalachian region are described. Emphasis of the study was placed on low-ash bituminous coals suit-able for utility or metallurgical markets. There was good agreement between the results obtained by conventional ash determination methods and by the Cendrex for all the products tested. A principal disadvantage of coal in competition with other fuels has been its inconsistent quality, especially with respect to ash content. Mine rnechanization and production from coalbeds of inferior quality have accentuated this failing. Present quality control methods in coal preparation plants are inadequate because (1) control samples taken over a production shift only serve to give the average quality of the product, (2) conventional sampling and analysis data are usually not available at the time a shipment is scheduled, and (3) the plant operator does not have the use of such data when the coal is being processed so that operating conditions may be adjusted for correcting variations in quality. Thus, the development and use of continuous-control monitoring instruments is of prime importance. The Federal Bureau of Mines arranged to investigate an instrument using x-rays for rapid continuous ash determinations. Developed by the Dutch State Mines and first described in 1958,' this instrument, which has found widespread acceptance in Europe, is designed to provide ash analysis of any desired coal feed or product stream with only a five-minute time lag between sample collection and analytical results. With limit controls, this unit could be used to feed back signals to preset the separating gravity of coal washers, control the ratio of various coals needed to produce an acceptable blend, and permit the sorting of a washed product into two qualities, one which is below and the other which is above a prescribed ash content. This unit can also be used discontinuously for rapid determination of coal samples taken in the usual manner. THEORY Di jkstra and Sieswerda of the Dutch State Mines studied the effect of chemical composition and ash content on the absorption coefficient of coal with the transmitted radiation method.' They found that with an X-ray wavelength of 1.79 A and values calculated for the absorption coefficients of ash and combustibles, the absorption coefficient of coal gave a linear relationship when plotted against the ash content. Thus, a method of measuring the absorption coefficient of coal could also be used to measure the ash content. However, practical attempts to measure the ash content of coal by comparing the intensities of the transmitted X-rays failed. The absorption coefficient of coal can also be obtained by irradiating a layer of coal with x-rays and measuring the intensity of the diffused reflection (backscatter method). The relationship between the diffused reflection and the ash content is actually hyperbolic but approximates a straight line over a narrow range. A measuring circuit that balanced the difference in intensity of the reflected radiation from a standard surface against that from a coal sample was used. By plotting this difference of reflected radiation against the ash content, a curve which approximated a straight line over a range of 8 percentage points of ash was obtained. DESCRIPTION OF EQUIPMENT The apparatus performs two distinct functions: (1) the preparation of the sample to a suitable size consist and moisture content for delivery to the x-ray element, and (2) the measurement of the ash content by the x-ray element. In a continuous operation the coal sample is delivered to a conditioner where it is dried to a maximum of 1% surface moisture and pulverized in a high speed mill to -48-mesh. This material is fed, at a rate of 250-300 gpm onto a rotating disk, shown in Fig. 1, where it is mixed with material already present. The mixture is rotated on the disk and

Jan 1, 1969

By J. T. Norton, H. Blumenthal, S. J. Sindeband

The interstitial phases formed by the transition elements with carbon, nitrogen and boron constitute a unique class of substances which are of considerable technical interest because of their well developed metallic properties and their high hardness and melting points. They are also of interest from the point of view of structure. It was pointed out by Hagg1 some years ago that the principal factor determining these structures was the relative sizes of metal and metalloid atoms. He predicted that if the ratio of the radii of metalloid to metal atoms was less than 0.59, the metal atoms would be arranged in a close-packed fashion of relatively simple type. This prediction has been amply confirmed. It appears, however, that this specification applies most generally to those structures in which the metalloid atoms occupy isolated positions in the lattice and there is no tendency for strong binding between the metalloid atoms themselves. The series of borides of the type MeB2 formed from the transition elements of the fourth (Ti, Zr) and fifth (V, Cb, Ta) groups of the periodic table provide an opportunity to investigate this situation further. If one takes as a basis the radius of 0.87 A for the boron radius and the Goldschmidt values for the metal atoms in 12 coordination, then the radius ratio varies from 0.54 for zirconium to 0.64 for vanadium, which latter is considerably greater than Hagg's limit of 0.59. It was of interest, therefore, to see if these five elements would form borides of the MeB2 type and have simple isomorphous structures. When this investigation was initiated, relatively little information had been published on the structure of this group of borides. McKenna2 reported the preparation of ZrB2 by a carbon reduction process. Its structure was hexagonal and the approximate lattice constants were given as a = 3.15 A, c = 3.53 A. Ehrlich3 reported the structure of TiB2 as being of the C 32 type with a = 3.02, c = 3.21 and c/a = 1.06 A. Recently, Kiessling4 published the results on ZrB2, showing it to be one in which the metal atoms are arranged in a simple hexagonal lattice with the boron atoms in flat sheets, midway between the layers of metal atoms. The structure is C 32 type with a = 3.169, c = 3.530 and an axial ratio =1.11 A. He also states that CbB2 and TaB2 are isomorphous with ZrB2. Preparation of Borides The borides of titanium, zirconium, columbium and tantalum were made by the electrolysis of fused salt baths after the method described by Andrieux.5 The bath, totalling about 6500 g in each case had the following general composition: Metal Oxide + B2O3 + CaO + CaF2. The actual proportions of the ingredients were determined for each individual case. Although Andrieux reports the preparation of vanadium boride by the same electrolytic method, several attempts proved unsuccessful. However, it was possible to prepare VB2 by the carbon reduction method used by McKenna to produce ZrB2. The mixture consisted of V2O5 + 6B2O3 + 11C and the reduction was carried out in an induction heated graphite crucible. This method proved satisfactory for ZrB2 and TiB2 and undoubtedly the others could be made in the same way. The borides, which were all in the form of fine gray metallic crystals from the electrolytic bath or gray powders from the reduction process, were carefully separated from other constituents and analyzed chemically for metal and boron content. The densities of the boride powders were measured using a standard 10 cc pyenometer at 18°C and ethyl benzene proved to be the most satisfactory liquid for this purpose. X Ray Technique The boride crystals, as prepared, were too small for single crystal examination so that powder techniques were employed. Patterns of the finely ground borides were prepared, using the Norelco recording X ray spectrometer. The lines on these patterns were indexed with the aid of Hull-Davy charts and the approximate dimensions

Jan 1, 1950

By D. Muraoka

Introduction Soda ash, anhydrous sodium carbonate, is produced from underground trona deposits occurring in the Green River Basin of southwestern Wyoming. Stauffer Chemical Co. of Wyoming, a jointly owned subsidiary of Stauffer Chemical Co. and Rocky Mountain Energy Co., mines the trona and processes it into high-purity soda ash. Stauffer Wyoming has been in operation since 1962 with a current production capacity of about 1.8 Mt (2 million st) per year of soda ash. The refinery has five parallel operating trains using the monohydrate process which is a series of chemical processing unit operations described below. Screening and Crushing The mined trona is brought to the surface in two counterbalancing ore skips. The ore is initially screened with refinery feed (less than peagravel size material) entering the process. Oversize is conveyed to an outdoor stockpile. Ore is recycled from under the stockpile via cone feeders as the refinery needs require. The ore is conveyed to hammer type crushers and screened. The refinery feed enters the process and the oversize is recycled to the crushers. Calcining The first step in refining trona to soda ash is calcining the trona ore. This is accomplished in concurrent flow, rotary kilns that are direct fired with natural gas. [2 (Na CO NaHCO 2H O) Heat 3Na CO + 5H20 + CO 2 3~3~ 2 ~~2 32] The trona ore is fed to the kiln at the burner end and is intimately mixed with the hot exhaust gases from the burner. Direct contact with the flame is avoided. The ore moves down the kiln by gravity and rotation aided by the kiln's internal lifters. The exhaust gases are cleaned to 99 + % particulate removal and discharged to the atmosphere. Dissolving the Calcined Trona The calcined trona is conveyed from the kiln discharge to a rotary drum dissolver where it is mixed with water and weak solution (unsaturated solution of calcined trona in water). The calcined trona is soluble whereas the impurities, primarily shale, are insoluble. The discharge from the dissolver is a slurry containing saturated solution (about 30% dissolved calcined trona in water), some undissolved calcined trona, insoluble impurities, and organic impurities. The dissolver discharge enters a rake classifier that makes the first solid/liquid separation in the process by removing the heavy undissolved material from the saturated solution. This undissolved material is conveyed to a second rotary drum dissolver and mixed with water to dissolve any remaining calcined trona. The secondary dissolver discharges to a second rake classifier that separates the solids from the liquid. The solids are now primarily insoluble impurities with very little unrecovered calcined trona. These wastes are conveyed to evaporation ponds for disposal. The liquid from the secondary rake is recycled to the primary dissolvers' weak solution system. Returning to the primary rake classifier, the saturated solution which contains suspended insoluble muds is pumped to a thickener. Flocculant is added to aid settling of the suspended muds. The flocculated muds are collected at the bottom of the thickener and pumped to vacuum cloth filters. The adhering saturated solution is separated from the muds that remain on the filter cloth. The recovered solution is returned to the dissolvers' weak solution system. The muds are reslurried and pumped to the evaporation ponds for disposal. Saturated solution, containing some suspended insoluble muds along with the organic impurities, overflows the thickener. Filtering the Saturated Solution The saturated solution overflowing the thickener enters a holding tank where it is pumped to a bank of pressure leaf filters. The suspended muds are captured on the filter leaves while the cleaned saturated solution is pumped to polishing filters. The polishing filters are also leaf pressure type, intended to remove any suspended solids that may pass the primary filters. Organic impurities are removed with activated carbon. The purified solution is now free of insoluble muds and has most of the organic impurities removed. Crystallizing and Drying the Monohydrate The purified saturated solution is pumped to evaporative crystallizers. The solution is circulated from the crystallizer vessel through a heat exchanger. Indirect contact with steam in the heat exchanger causes the saturated solution to boil, thus crystallizing the soda ash in its monohydrate form. Slurry crystal density is maintained in the crystallizer by simulta¬neous feed of fresh saturated solution and drawoff of the slurry.

Jan 1, 1986

By Bob J. Rawicki

TYPES OF WET DUST COLLECTORS Basically, there are two types of wettype dust collectors. One is mechanical, incorporating pumps, motors, fans, sprays, filters, or flooded beds. These come in many forms, but their operating principle is basically the same. The screen or flooded bed is wetted by a series of sprays. Polluted air is drawn in through a fan; here the dust particles impinge on the screen and are flushed off into a settlement tank. The air then passes through spray eliminators to the atmosphere. A great deal of time and money has been spent developing this type of collector. In laboratory tests where clean water at constant pressure, air at constant volume, and dust at low but constant mass to volume were applied, these collectors gave high percentage dust collection efficiencies. Theoretically and technically, they are very good. In practice, however, under rugged mining conditions where there are no constants, and dust load, particle size, water pressure, and volume of air vary with time, their efficiency varies dramatically. Some of the most common problems are: 1. clogging of screens 2. blocking of nozzles 3. replacement. of damaged pump stators Any of the above result in loss of air flow and dust collection efficiency. Generally, there are so many mechanical moving parts, something is always breaking down. This type of collector is very expensive in parts and labor to maintain. The second type of dust collector, which I personally developed, with the help of the British National Coal Board of Great Britain, has none of the above mentioned problems and it works on a completely different principle. My collector, the Mark III Precipitaire Wet-Type Dust Collector, (an improved version of previously successful models) has no internal moving parts, no flooded beds, screens, spray nozzles, etc. It is highly efficient - 99.8% efficiency on total collected dust and 97% on respirable dust. The pressure drop through the collector is constant, hence, air flow remains constant. There is no maintenance expenditure, be it labor or parts. The only maintenance that there is, is desludging. This, of course, varies with the dust load. With a low dust load, desludqinq may only be required every 4 weeks, but this is a very simple operation that can be handled by unskilled men and it takes only a very short time. My collector is powered by a ventilation fan that is placed on the clean air side for a longer fan and motor life. Dust is extracted solely by water action. The collector's operating principle is: Dust laden air is drawn in along ducting into a tapered scrubber section, shaped to create a self-induced curtain of water. This action washes the dust particles from the air and the collected dust settles into the bottom of the tank. Cleaned air is then exhausted via a series of spray eliminators where it is completely dried and on into the atmosphere. These collectors have been sold in Great Britain, to the National Coal Board, and throughout Europe. They have been working very well for many years now. GENERAL APPLICATIONS Dust is very dangerous; silicosis and black lung kill and disable many workers each year. Many of the mine explosions can be attributed to coal dust. By the use and application of dust collectors, these deaths and disasters can be avoided. There really is no excuse. However good and efficient a dust collector may be, it's positioning and installation arrangement is critical. The collector may remove 99% of dust supplied to it, but it takes the skill and experience of a dust control engineer to design and fit ducting, hoods, etc. to remove polluted air at source and deliver it to the collector. In some cases this is very simple; in the case of longwall machines, it is very difficult. In general, dust collectors should be used to collect dust at every point of dust generation. I have noted a great reliance on dilution as a means of dust control. I personally think this is a dangerous and unreliable method. As mines become deeper and the use of larger and more productive equipment becomes more common, reliance on dilution will be impossible and dust collectors will be essential.

Jan 1, 1982

By John Tabor, R. B. Spivey

INTRODUCTION Shafts have been sunk in a number of ways. By and large, however, most have been sunk by drill, blast, muck, and slip form methods. Most methods have provided satisfactory means of completion. The major drawbacks to the conventional methods have been: (1) un- reliable cost estimates for completion; (2) unreliable rates of progress; and (3) unreliable integrity of the shaft lining. Recently , engineering design studies have been done that strongly indicate, that, by the use of existing proven technology, shafts can be completed faster, more cheaply, safer, and with a better lining that has been possible using conventional methods. THE METHOD The alternate method utilizes the techniques of big hole drilling, shield excavator tunnelling, and pre-cast concrete lining. Big Hole Drilling - Drilled shafts have been completed for several years and I don't intend to address the technology, only to point out the use of drilled shafts is an alternate method of shaft sinking. Ideally, in a multi-shaft mine, the first shaft would be blind bored and cased to a size, whereby it can handle the muck from the excavation of subsequent shafts. It should be an 8' to 10' completed diameter. The cost of such a shaft in sedimentary formations compares favorably with a conventional shaft. However, it is much faster. For additional shafts, the role of the drill is to complete small holes (3' to 4') to the desired horizon. This is merely a hole to allow slashing to a larger size with easier muck removal from an underground drift connected to the previously drilled shaft. The smaller holes, just as the larger one, should be cased and grouted into place so as to assure a smooth open hole. SHIELD EXCAVATOR Theory and History - Tunnelling Shields have been in common use for several decades, and steam powered shields were used in England over 100 years ago. The purpose of the shield is to provide temporary ground support, protect personnel during operation, and to house the excavating equipment. The excavating equipment approach has varied from attempts at a continuous or semi-continuous boring machine, the tunnel mole, to a mechanical or hydraulic digging arm. The shield may or may not be an active digging element. If not active, the boring tool or cutter is advanced ahead of the shield to full diameter and then the shield is moved up to the cutter, or the shield advances along with the cutter. If active, the cutter opens a conical free face, and the shield is either jacked forward so that the leading edge spalls rock to the free face or poling plates are thrust forward to spall to the free face. In either approach, various methods of muck removal are incorporated into the shield. These may be gathering arms, conveyors, or with moles, hydraulic transport, where the cuttings are sufficiently fine. The efficiency of the muck removal system is critical as excavation is rapid and the system could easily become muck-bound. Application to Shaft Sinking- Technically, the use of a shield to sink vertical shafts is no different than a tunnelling operation. It is, in fact, less complex in most respects. Muck and water removal are difficult in a blind-driven shaft. However, where an underground opening is available for handling muck and water, rapid shield sinking is possible, using a pilot hole for muck and water to fall through and be removed from below. The Tabor Mining Shield is basically a caisson approach, whereby, the ground is only seen at the face or bottom. It does have hydraulic jacks to force it downward and has segmented cutting edges which allows forepoling ahead of the shield. The hydraulic excavator allows mining in the center of the face around the previously drilled hole with the forepoling segments slabbing to this free face. First, the drilled hole is cased and grouted. The machine is a tube with an excavator. The work cycle is as follows: (1) The casing and grout are cut at predetermined intervals by shaped charges. (2) The casing is crushed by the excavator, removed from the hole and hoisted to the surface. (3) The excavator then rips out a face around the hole by digging the formation in tension and breaking the pieces

Jan 1, 1982

By R. E. Ogilvie, E. T. Peters

The X-ray Kossel-line method has been used preaioz~sly for measuring lattice parameters to accuracies of 1 part in 100,000.5 A second application of this method is described for determining the crystallographic orientation of a randomly positioned single-crysta1 spherical volume that can he as small as 50 µ in diameter, within accuracy limits of ±1/2 deg. The theory, experimental procedure, and interpretation of Kossel-line patterns and an experimenta1 verification of the predicted orientation relationship of the cph k and fee a Cu-Si phases are presented. This orientation relationship can be described as (111)a, (00.1)k; [110]a [11.0]k. In addition, the lattice parameter of the a phase was found to he a. = 3.62154 ± 0.00014Å. The Kossel-line method when used in conjunction with electron microanalysis is shown to he capable of providing a complete chemical and structural analysis of a given crystal. THE most easily accomplished methods for determining the orientation of a single crystal are variations of the Laue X-ray diffraction method. Although certain materials can be oriented to within ± 1 deg accuracy by observation of exterior macroscopic features, such as etch pits, cleavage faces, or growth features, the Laue method is generally preferred for routine laboratory application. Both back-reflection and transmission patterns are coordinated by appropriate reference charts and are plotted in terms of reflection plane poles (normals) on a stereographic projection. Orientation is deduced by relating the pole distribution (which is fixed by the crystallographic symmetry of the specimen) to two specified external reference directions. The Laue method has several limitations. 1) Orientation can rarely be determined to better than ±l deg accuracy. Principal errors involve inaccuracy in specimen-to-film distance, measurement confidence of individual Laue spots, and stereographic plotting. 2) Because of the indeterminacy of the X-ray wave length diffracted to a given spot, it is not possible to determine supplementary crystallographic data, such as interplanar spacings and lattice parameters. 3) The method is generally limited to specimens of cubic symmetry or to specimens of high symmetry and known structure. 4) The relatively large cross-sectional area of the incident X-ray beam generally precludes the measurement of relative grain orientation in a polycrystalline material.* _____ Several of these limitations can be overcome by application of the Kossel-line method, which has been previously employed for precision lattice-parameter determinations.'-= For this method, a point source of divergent monochromatic X-rays is generated within the crystal by means of an incident electron beam. The divergent X-rays which fulfill the Bragg law are diffracted by the specimen and are recorded on a film placed either in transmission or back reflection. Analysis of the film yields a direct measurement of orientation and lattice-spacing values to an accuracy of ±1/2 pct. As analyses can be obtained from spherical specimen volumes as small as 50 µ in diameter, the method provides a means for structural analysis of second-phase or impurity precipitates within a given matrix. The primary limitation of the Kossel-line method is the requirement for an electron microanalyzer or similar apparatus capable of producing a finely focused electron beam. This paper is designed to present the theory, experimental procedure, and geometrical interpretation of Kossel patterns. The experimentally determined orientation relationship between the k and a phases occurring in the Cu-Si system and a precision measurement of the a lattice parameter are presented as a practical application of the method. THEORY OF KOSSEL LINES Divergent X-ray beam photography utilizes an effective point source of characteristic X-rays which, when diffracted from a single crystal, form numerous diffraction and absorption cones that are recorded on film.7 The cones generated from a source lying within the crystal are called Kossel lines.' Although the X-ray scattering from a divergent point source contained within a crystal is described in terms of Laue dynamical theory,9 the directions of the diffracted spectra can be ade-

Jan 1, 1965

By W. D. Ludemann, J. E. Dor, L. A. Shepard

Recovery of the creep resistance of 99.99 pct pure Al was studied at temperatures 540°, 573°, 600°, and 611°K. Poly-crystalline specimens crept under a stress of 950 psi to a strain of 5.5 pct were allowed to recover for periods of from 1 min to 16 days under a residual stress of 4.4 psi. Increased creep rates upon reapplication of the 950 psi stress evidenced softening of the material. The activation energy for the recovery process was found to be 64,000 cal-per mol. THREE major observations have clearly revealed that the creep of polycrystalline aluminum above about 510°K is controlled by the rate of climb of jogged edge components of dislocations past the barriers impeding their motion: 1) Extensive polygonization, characteristic of climb, takes place during creep.'9' 2) The activation energy for creep equals the estimated activation for self-diffusion. 3) The stress law for creep coincides with theoretical predictions based on the climb mechanism.4,5 The decreasing rate during the primary stage of creep must be ascribed to the introduction of additional barriers to the glide of dislocations. During secondary creep the density of barriers must remain constant, the rate at which new barriers are formed being equal to the rate at which they recover. For this reason the nature of the barriers, and their rates of formation and recovery, are significant to a complete understanding of creep. It is proposed here to study the kinetics of the recovery of barriers to creep of high-purity polycrystalline aluminum in the climb range. The technique will consist of creeping a series of tensile specimens under a prescribed stress to given state following which the specimens will be recovered for various times and temperatures under approximately zero stress. The amount of recovery will be determined by comparing creep curves following the recovery. EXPERIMENTAL PROCEDURE The material used for this investigation was a 99.99 pct pure Al supplied by the Aluminum Co. of America in sheets of 0.001-in. thickness having an H-18 temper. The spectroscopic analysis of impurities gave Cu-0.004 wt pct, Fe-0.002 pct, Si-0.001 pct, others-0.000 pct. Creep specimens were machined with their tensile axes in the rolling direction of the sheets. Annealing in a molten potassium nitrate-potassium nitrite bath for 1 hr at 686K, followed by air-cooling, resulted in a grain diameter of from 0.22 to 0.27 mm. Creep machines were equipped with constant stress lever arms of the type suggested by Fullman, Carreker, and Fishher,' which maintained the applied stress to within 0.04 pct of the reported value of 950 psi. Creep strains were calculated from extensions over a 6-in. gage section to the nearest 3 X 10"5. The temperatures of testing, 540°, 573", 600°, and 611°K, were obtained by immersing the specimen and extensometer assembly in a molten KNO2-KNO3 bath. Temperatures during the periods of creeping could be maintained constant to within better than ± 1°K. Correlation of creep data at different temperatures of testing demanded the use of an appropriate temperature-compensated time 8. Previous investigations have shown that the creep of pure metals at temperatures above about 0.55 of their melting temperatures can be correlated by the functional relationship3" e=f(0), s= const. [1] where e = total plastic strain during creep f = a function that depends on stress a = stress 6 = te-Q/RT = a temperature-compensated time t = time under test Q = 35,500 cal per mol = activation energy for creep R = gas constant, and T = absolute temperature In the present series of tests all of the creep was conducted under the same stress of 950 psi. To provide identical initial conditions for the recovery all specimens were first crept to a temperature-compensated time of = 47 X 10-14 min which resulted in a creep strain of 0.055 + 0.005 regardless of the creep temperature. During recovery, a small stress (4.4 psi) was permitted to remain upon the specimens to maintain tautness in the pulling assembly and to prevent mechanical damage to the soft specimen. After each recovery period, the full stress was reapplied and

Jan 1, 1961

By Carl Altstetter, Roy Adams

Measurements of enthalpy changes and transformation temperatures me reported for the fcc -hcp marten-sitic transformation in pure cobalt single and multi-variant crystals andpolycrystalline specimens. From these measurements the free-energy difference between the two phases was calculated for both heating and cooling and for repeated cycles through the transformation. For heating of single crystals the average value of enthalpy change was 113 cal per mole but on cooling it was 84 cal per mole. This difference is interpreted in terms of lattice imperfections which are introduced during transformation. A calorimeter suitable for measurement of enthalpy changes of athertnal transformations is described. COBALT exhibits an allotropic phase transformation at around 417°C' where an hcp phase (a = 2.507, c = 4.0686, c/a = 1.6228)~ stable at low temperatures transforms to an fcc phase stable up to 1495°C. The transformation has a hysteresis of about 30°C— the exact magnitude of the hysteresis depending on the number of transformation cycles, grain size, impurity content, and prior treatment.'-6 The athermal character of the transformation and surface relief due to transformation shear indicate that the transformation is martensitic. Pole mechanisms for the transformation have been discussed by Sebilleau and Bibrin, Bilby, Seeger,' and Basinski and Christian. These authors have proposed that the cooling transformation is effected by rotation of a Shockley partial dislocation around a pole dislocation with a screw component of 2a/3(111). Such a rotation would change the stacking of close-packed planes from the abcabc fcc stacking to the ababab hcp stacking sequence. This mechanism would predict a lattice correspondence of {111}fcc~~(0002~cp, have been observed. One would expect the heating transformation to occur similarly, involving dislocation movements on the (0002)hcp. However, a high density of pole dislocations and high dislocation ve- locity are required to account for microscopic observations. Furthermore, the high back stress developed on the Shockley partial dislocation as it rotates toward the sessile dislocation associated with it makes it doubtful that this is in fact the transformation mechanism. A mechanism for spontaneous nucleation of partial dislocations1*l2 seems more reasonable. If a single crystal of cobalt is cycled through the transformation but never heated above 600°C, it will remain single."13 However, if the specimen is heated through the transformation and held for a short time at or above 1000°C before cooling, the cooling transformation usually involves all four {lll}fcc habits, resulting in neighboring orientations of hcp differing by 70.5 deg. This condition is termed multivariance and has been observed by Bibring and sebilleau5 and Nelson and ~ltstetter.' The latter authors used single crystals of cobalt grown in an electron beam zone refiner. They proposed that the strongly directional cooling caused the initial cooling transformation to occur on the {lll)fcc most nearly perpendicular to the temperature gradient. Transformation shear on one set of {lll}fcc planes immobilized dislocations on the other {lll}fcc planes involved in the initial transformation. Since there is only one (0002)hcp plane orientation, the product of subsequent heating transformations could have no more than two variants—the original fcc orientation or its twin. If, however, the specimen was annealed at a high temperature, such as 1000°C, there would be relaxation of dislocations into low-energy configurations, and in the absence of a sharp temperature gradient transformation on all four (111 Ifcc planes could occur. If a dislocation mechanism is responsible for the transformation, the perfection of the lattice should affect the reversibility of the transformation. Mutual constraints of neighboring grains in polycrystalline material would inhibit dislocation movement resulting in increasingly sluggish transformation. Indeed, Nelson and Altstetter found the M, to decrease and the A, to increase in going from single crystals to multi-variant specimens to polycrystals. The effect is enhanced by decreasing the grain size of a solid specimen or a powder.6'14 Data on the free-energy difference necessary to initiate the transformation gives information about the driving force of the transformation. This is an important consideration in the understanding of the transformation mechanism. The free-energy change associated with the transformation of the hcp phase to the fcc phase at some temperature T different from the equi-

Jan 1, 1969

By John C. Calhoun, R. C. Newman

A report is made of experimental work performed on the capillary retention of water within porous solid systems, the displacement being accomplished with air and various organic liquids. A portion of the experiments were designed to measure the lowering of vapor pressure of water within a porous solid with subsequent conversion of such vapor pressure data to capillary pressure values. The form of the capillary pressure curve at high capillary pressures has been elucidated from this data. Certain theoretical approaches are presented to indicate the correlation between work done by previous investigators. The present work is correlated also with theory. Surface area values are calculated for the various core systems studied. INTRODUCTION Leverettl in 1941 presented a paper which gave the essential concepts of capillary behavior in porous solids. His presentation included both theoretical and experimental aspects of the problem. He defined the term "capillary pressure" and applied it to an ideal porous system. His work included the definition of a dimensionless quantity which was a function of fluid saturation and which could be correlated with porosity and permeability for clean, unconsolidated sands. His experimental work was performed by the drainage of water from packed columns of unconsolidated sands. Subsequently, other authors have presented experimental techniques which permit the determination of capillary pressure data on small core samples. Data have been presented on porous systems other than unconsolidated sand and with the use of liquids other than water. Notable among these are the works of Hassler, Brunner and Deahl,2 Bruce and Welge,³ Amyx and Yuster, and Purcell. Others have indicated the applicability of such data to field problems without regard to the theoretical aspects of the capillary pressure saturation function. Within the recent published literature there have also been reported studies of the correlations between capillary pressure data and other fundamental properties of porous solids. It was the intent of the present reported work to investigate the behavior noted by water and various other liquids within porous systems in an attempt to amplify existing correlations between capillary properties, surface properties, and other fundamental characteristics of porous systems. The experimental work described here was performed as two different investigations but due to its nature it is reported as separate phases of the same general topic. CAPILLARY PRESSURE AND VAPOR PRESSURE LOWERING It has long been recognized that the vapor pressure above the curved surface of a liquid is a function of the curvature of the liquid surface." The capillary pressure is also a function of the curvature of the liquid surface. Both capillary pressure values and vapor pressure lowering are, therefore, functions of the liquid saturation of a porous solid. A quantitative relationship between the vapor pressure lowering and the capillary pressure can be simply developed by considering a porous solid system containing water in equilibrium with its vapor. The capillary pressure at any height within the system is defined as Pc= Pv -P1.......(1) where PI is the pressure in the vapor phase and PI is the pressure in the liquid phase. Consider that both the water and its vapor are continuous phases within the system and that at any point the two phases are in equilibrium. The gradients in the liquid and vapor columns within the porous system are: PI = plgdh dPv = pvgdh By integrating these gradients from the point where h = 0 to the height h, values of PI and PI at the point h are found. Let the point where h = 0 be that where Pc = 0. At this point also Pv is equal to Pve, the equilibrium vapor pressure of the liquid above a flat surface. From the first of these integrations since p is independent of height: P = p1gh + P, From the second, by substituting pv = PvM/RT, where M is the molecular RT weight: Pv M Pvo- RT gh By combining the two last equations to eliminate h, one obtains: P1 RT Pv Substituting, values for Pv and Pi in equation (1) gives PyRT RT Pv which is approximately Pc=RT/Mp1 In Pvo/Pv .......(2) This equation relates the vapor pressure above a curved surface with the capillary pressure across the curved surface in terms of measurable quantities. The equivalent of this equation is given by Freundlich."

Jan 1, 1949

By N. De Nevers

A calculation method has been developed for carbonated water flooding. This method takes into account the effects of oil viscosity reduction and oil swelling, due to carbon dioxide transferred to the reservoir oil from the carbonated flood water. It shows the effect of changes in carbonation pressure and carbonated water slug size. The method is based on a Buckley-Leverett-type linear flow model. The mathematical approach is similar to that developed by Welge, et al, for enriched gas drives. However, it contains a very important advance over the previous papers in its solution to the problem of injection of a "slug" of carbonated water rather than continuous injection of fluid of constant composition. The method also shows the chromatographic nature of the transport of Co 2 in a carbonated water flood. Sample calculations indicate that viscosity reduction is the most important source of incremental oil recovery. Swelling contributes to a lesser extent. Raising the carbonation pressure increases recovery, but this effect diminishes as pressure is increased. increasing the slug size also increases recovery, but again this effect diminishes as slug size is increased. INTRODUCTION In recent years there has been considerable interest in carbonated water flooding as a secondary recovery method. Several carbonated waterflood field tests are now in progress.1'3 Yet there is no published method for predicting the results of a carbonated water flood. This paper presents a method for predicting the oil recovery for ; carbonated water flood. In a carbonated water flood. a carbonated water slug is injected into the reservoir. Then this slug is driven forward by injecting plain water behind it. The field operator must choose the carbonated water slug size and the CO2:water ratio in the slug. The calculation method presented here shows the effects of both slug size and CO2:water ratio on oil recovery. The calculation assumes that the addition of carbon dioxide to flood water increases the oil recovery by two mechanisms: 1. Carbon dioxide transfers from the flood water to the oil, thereby reducing the viscosity of the oil. At the lower viscosity, the oil flows more readily and is more easily swept out by flood water. 2. The carbon dioxide which transfers to the oil phase causes the oil phase to swell. Assume that the residual oil saturation (volume per cent of pore space) is the same whether the oil is carbonated or not. Then the same volume will contain less stock tank oil if it is carbonated, due to oil swelling. The plain water driving the carbonated water slug forward will strip the CO2 out of this residual oil., causing this oil to shrink. This swelling-followed-by-shrinking leaves behind a lower residual oil saturation than would have been left behind by a plain water flood. Since some oils swell to almost twice their original volume at reasonable CO2 pressures, this effect may be large. Two additional benefits sometimes claimed for carbonated water flooding are: 1. Adding CO2 to flood water greatly increases the water injectivity in some cases.1 (Injectivity = BPD of water injected/ft of sand x differential bottom-hole pressure, psi.) CO2 forms a weak acid in water; injecting carbonated water is, therefore, somewhat like acid treating. If a formation responds to acid treating, it will probably have a higher injectivity for carbonated water than for plain water. Although the increased injectivity increases the profitability of a carbonated water flood over a plain water flood, it does not change the production curve based on total fluid injected. Therefore, the method presented here does not show any effect of increased injectivity. 2. It is claimed in one patent4 that carbon dioxide will react with some constituents of the oil to form detergents. These detergents are expected to increase the oil recovery. There is little data available to support this idea. The calculation method presented here takes no credit for extra oil recovery due to the formation of any such detergents.

Jan 1, 1965

By U. N. Bhrany, J. H. Brown

The specifications of particle size and the size analyses of fine particulate materials are commonly presented without reference to the method of analysis. A review of the various sizing methods showed that not only numerous sources of error but also different definitions of size are inherent in the different methods of analysis. To demonstrate the magnitude of the differences, the size distribution of finely crushed quartz was determined by several methods and a critical comparison of the data was made. The experimental results show that appreciably different sizes may be indicated for one material depending on the method of analysis used. However, because the shape of the size-distribution curves determined by almost all methods was constant, a simple method of correlating the data can be adopted. Conversion factors, by which size measurements obtained by the various methods can be expressed in terms of screen sizes, are presented as one method of correlating the results. The size analysis and size specification for lump rock and particulate materials are commonly determined with sieves. With this technique, a readily reproducible measure of size can be obtained merely by duplicating the sieve and the method of sizing. For coarse materials, such as lumps of rock, this technique is simple and effective; and the effects of properties of the material such as shape can be readily understood. When very small particles are measured, however, problems arise. Small sieves are difficult to construct and maintain, and size separations become sensitive to the analytical procedure. For these difficult sieve sizes and for still smaller particles, other sizing techniques employing, for example, sedimentation, light extinction, and microscopy are available. Because these various techniques employ different size-dependent responses to establish size, it is appropriate to examine them and the results they yield to determine their relationships and their limitations. The 'size' of a lump of rock is not an absolute quantity, but rather a defined one. With sieves, for example, the size is defined by the screen aperture and no reference is made to the actual shape of the particle, yet only in the case of a regular body, such as a sphere, can the sieve size be related to the shape. With sieves only two dimensions of a particle can be suggested, for the sieve makes no allowance for the length of the particle passing through an aperture. Similarly, the size of a particle as measured by a microscopic technique need bear no relationship to the size as measured by sieves, although the two measurements may be related statistically. In fact, all size measurements are not only defined quantities, but the measurement actually used is an average that is controlled by the technique employed and, hopefully, reflects the property in question of the material of interest. Many papers have been presented that illustrate the accuracy of the various sizing methods, present the advantages of certain methods for certain applications, and describe the principles and technical backgrounds of the various methods. No attempt will be made in the present paper to review this work in detail. The purpose of the present paper is to examine the limitations of the various sizing methods, and to demonstrate experimentally the relationships between the sizes as determined by the various methods of analysis. It is hoped that by recognition of these limitations and relationships, errors in size specification and size interpretation will be avoided. REVIEW OF SIZING METHODS In specifying the size of a material, not only the data but also the method of measurement must be stated. For example, the size of the lump might be specified as its total volume, as its weight, as its maximum dimension, or as its minimum cross section; any one definition might be completely satisfactory for one problem, but not for another. In practice, size analyses follow four basic approaches and fortunately, the conditions that restrict each approach are known. These major methods include: 1) screen analysis, 2) direct measurement of particle dimensions, 3) determination of equivalent spherical size in response to fluid flow, and 4) specific surface determination. Each of these techniques has several variations and some new techniques are also available, but the limitations of the basic techniques are sufficient to permit at least a preliminary

Jan 1, 1962

By John C. Calhoun, R. C. Newman

A report is made of experimental work performed on the capillary retention of water within porous solid systems, the displacement being accomplished with air and various organic liquids. A portion of the experiments were designed to measure the lowering of vapor pressure of water within a porous solid with subsequent conversion of such vapor pressure data to capillary pressure values. The form of the capillary pressure curve at high capillary pressures has been elucidated from this data. Certain theoretical approaches are presented to indicate the correlation between work done by previous investigators. The present work is correlated also with theory. Surface area values are calculated for the various core systems studied. INTRODUCTION Leverettl in 1941 presented a paper which gave the essential concepts of capillary behavior in porous solids. His presentation included both theoretical and experimental aspects of the problem. He defined the term "capillary pressure" and applied it to an ideal porous system. His work included the definition of a dimensionless quantity which was a function of fluid saturation and which could be correlated with porosity and permeability for clean, unconsolidated sands. His experimental work was performed by the drainage of water from packed columns of unconsolidated sands. Subsequently, other authors have presented experimental techniques which permit the determination of capillary pressure data on small core samples. Data have been presented on porous systems other than unconsolidated sand and with the use of liquids other than water. Notable among these are the works of Hassler, Brunner and Deahl,2 Bruce and Welge,³ Amyx and Yuster, and Purcell. Others have indicated the applicability of such data to field problems without regard to the theoretical aspects of the capillary pressure saturation function. Within the recent published literature there have also been reported studies of the correlations between capillary pressure data and other fundamental properties of porous solids. It was the intent of the present reported work to investigate the behavior noted by water and various other liquids within porous systems in an attempt to amplify existing correlations between capillary properties, surface properties, and other fundamental characteristics of porous systems. The experimental work described here was performed as two different investigations but due to its nature it is reported as separate phases of the same general topic. CAPILLARY PRESSURE AND VAPOR PRESSURE LOWERING It has long been recognized that the vapor pressure above the curved surface of a liquid is a function of the curvature of the liquid surface." The capillary pressure is also a function of the curvature of the liquid surface. Both capillary pressure values and vapor pressure lowering are, therefore, functions of the liquid saturation of a porous solid. A quantitative relationship between the vapor pressure lowering and the capillary pressure can be simply developed by considering a porous solid system containing water in equilibrium with its vapor. The capillary pressure at any height within the system is defined as Pc= Pv -P1.......(1) where PI is the pressure in the vapor phase and PI is the pressure in the liquid phase. Consider that both the water and its vapor are continuous phases within the system and that at any point the two phases are in equilibrium. The gradients in the liquid and vapor columns within the porous system are: PI = plgdh dPv = pvgdh By integrating these gradients from the point where h = 0 to the height h, values of PI and PI at the point h are found. Let the point where h = 0 be that where Pc = 0. At this point also Pv is equal to Pve, the equilibrium vapor pressure of the liquid above a flat surface. From the first of these integrations since p is independent of height: P = p1gh + P, From the second, by substituting pv = PvM/RT, where M is the molecular RT weight: Pv M Pvo- RT gh By combining the two last equations to eliminate h, one obtains: P1 RT Pv Substituting, values for Pv and Pi in equation (1) gives PyRT RT Pv which is approximately Pc=RT/Mp1 In Pvo/Pv .......(2) This equation relates the vapor pressure above a curved surface with the capillary pressure across the curved surface in terms of measurable quantities. The equivalent of this equation is given by Freundlich."

Jan 1, 1949

By D. M. James, J. C. Martin

The results are presented of a study of the application of analytical methods to the solution of two-phase flow into single wells. Approximate analytical expressions for the pressure distribution in two-phase flow are found ;;; a number of conditions. The results obtained from the analytical solutions me found to be in good agreement with results obtained by finite difference techniques using a high speed digital computer. Mathematical solutions for four sets of boundary conditions are presented. All of these solutions are composed of a short-term transient plus a steady or quasi-steady state. The rates of decay of the short-lived transients are analyzed. It is found that the durations of the short-term transients may be characterized by a parameter defined as the time constant which can be determined from simple relations. It is shown also that if the outer radius is much greater than the radius of the well. the short term transients decay at rates which are proportional to the square of the exterior radius, and the rates of decay are only slightly dependent upon the radius ratio. Numerical solutions based on finite difference techniques are presented for a number of conditions. The numerical solutions are in good agreement with the predictions based on the theoretical analysis for small and moderate drawdowns. Examples involving large drawdowns indicate that the nonlinearities in the equations of flow do not appreciably alter the longevity of the short-term transients. In all cases the time required for the short-tern transients to disappear is predicted satisfactorily. INTRODUCTION The mechanism by which oil and gas flow into a single well is of vital interest to the petroleum industry. The fundamental equations of two-phase flow which describe this mechanism are nonlinear partial differential equations. Numerical solutions of these equations describing pressure transients have been obtained with the aid of electronic computers. Although solutions obtained in this manner take into account a large number of effects, the reduction of this information to useful generalities is difficult. One method of obtaining generalities is the use of linearized approximations of the nonlinear equations. Since it is possible to obtain explicit solutions of the linearized equation, general properties of the role of pressure in the flow mechanism may be ascertained. Results obtained from this approach are limited to some extent by the linearizing assumptions. The severity of these limitations may be evaluated by comparing solutions of the linear equation with numerical solutions of the more exact nonlinear equations of two-phase flow. In the past considerable amount of work has been devoted to studying pressure build-up using the single-phase flow theory. Unfortunate1y, most prep sure build-up tests involve multiphase flow. A small amount of work has been done studying pressure build-up where the flow is two-phase. The encouraging results of these studies suggest that useful results may be found from additional studies of not only pressure build-up, but also the rapid transients associated with placing a well on production. This paper presents the basic theory of the pressure transients associated with placing a well on production and with closing it in. The paper is concerned chiefly with two-phase compressible flow; however, the results also apply to single-phase flow. The results are based on analytical solutions of the flow equations, and they are verified by numerical solutions using finite difference techniques. Much of the previous work on compressible flow into wells has been confined to single-phase flow. Some work has been done on compressible two-phase steady state flow, and solutions of the equations of flow have been found by finite difference techniques using high-speed computers. Muskat1 presents some rather general solutions to the equations of single-phase compressible flow into wells. Much work has been done on pressure build-up in wells (see, for example, Refs. 2-6). Almost all of the work on pressure build-up concerns single-phase flow with the exception of Ref. 5 and part of Ref. 2. Some work has been done on pressure fall-off in injection wells.718 Muskat9 presents the solution of the equations for radial steady state two-phase compressible flow. Handyl0 extends Muskat's solution to obtain

By R. G. Wheeler, J. W. Goffard

The Larson-Miller parameter was used to correlate time, temperature, and indentation creep of magnesium, aluminum, and some of their alloys. In the temperature range 300" to 450°C, the short-time Meyer hardness of pure magnesium was less than that of the magnesium alloys tested, but for long times the pure magnesium has greater indentation creep resistance. Aluminum (1100 alloy) had 1.5 to 2.5 times more indentation creep resistance than magnesium at 300" and 450oC, respectively. Hardening of aluminum with a dispersion of Al2O3 was effective in the time and temperature ranges studied. New technologies have required the development of new materials and the utilization of the more familiar materials for new and unusual applications. The use of magnesium and aluminum and some of their alloys, because of their desirable nuclear characteristics, light weight, low cost, and ready availability, has been extended to the 300" to 450°C temperature range. In this temperature range the basic consideration of these materials must be their rate of plastic flow rather than offset yield strengths. The indentation testing reported here arose from a need for design data for the load-holding ability of supports made of these materials. Test Procedure—Hardness indents were made with a 0.275-in.-diam quartz indentor and a 10.65-lb load. The indentor was made by fire-polishing a spherical surface on the end of a fused quartz rod. The samples were held at temperature in a graphite crucible controlled to ±2°C. A thermocouple was attached to the sample and test temperatures were recorded. The diameter of the spherical indentation was measured at the end of a test period and the compression stress (Meyer Hardness) was determined by: H___________load__________ m = projected area of indent Samples were 1 in. in diam and at least 1/4 in. thick. It was observed that at the higher temperatures and longer times, the quartz indentor would stick to the magnesium sample. The quartz indentor was, therefore, frequently inspected and fire-polishing repeated when necessary. The area of sticking was always a small fraction of the area of indent and was therefore considered to have an insignificant effect on results. Correlation of Hot-Indentation Test Data with Time-Temperature Parameter—Sherby and Dorn' have correlated creep or tensile data of a' solid solutions of aluminum with a temperature and strain-rate parameter suggested by Zener and Holloman. underwood2 used this parameter to correlate creep properties of some steels with hot hardness, and upon the basis of this correlation a means of obtaining creep properties from short-time (and inexpensive) hot hardness tests has been demonstrated. Since the validity of the correlation of creep properties with a time-temperature parameter and the correlation of creep properties with hot hardness have been shown, it follows that hot hardness may correlate with the time-temperature parameter. The hot-indentation data obtained was expressed as Meyer hardness, and was shown to be time and temperature dependent. Correlation of Meyer hardness, time, and temperature with the parameter was made using the relationship: Hm = Meyer hardness t = time, hours T = absolute temperature, OK K = constant A value for the constant K was calculated by equating In l/t + K/T at different temperatures and times but at the same hardness. The correlation was tested by plotting Hm vs the parameter, In 1/t +K/T. Since materials are being sought which have high hardness at low indentation creep, i.e., a high Meyer hardness for long time at high temperatures, low values of the parameter are ofthe most interest. TEST RESULTS Magnesium—Pure magnesium (99.98 pct) cut from extruded rod was indentation tested perpendicular to the rod axis at temperatures of 300°, 350°, 400°, and 450°C for times ranging from 6 sec to 112 hr. Fig. 1 shows the time dependency of Meyer hardness at the four constant temperatures. Fig. 2 shows the correlation of the Meyer hardness of pure magnesium with the time-temperature parameter using a K of 22,720 in Eq. [I]. At the bottom of Fig. 2, the effect of doubling the time of indentation t2 = 2(t1), on the abscissa for any time is shown graphically. This effect is of constant magnitude. Also shown graphically are the magnitudes of the effects on the

Jan 1, 1960

By Linda Lee, R. E. Reed-Hill, C. R. Smeal

Four fcc metal surfaces, etched to make them responsive to polarized light, have been studied with an electron microscope. Jones'prediction that these surfaces are grooved has been verified. Optical-goniometer measurements made on commercially pure nickel indicate that the groove walls are poorly defined (100) facets. Surface mientations close to either 4001 or {ill).- show no extinctions on the polam'zid-light microscope. An explanation is offered for this orientation dependence. It has also been deduced that polarized-light extinctions on a grooved cubic-metal surface should not be used directly in crystallographic-orientation determinations. The nature of the etching solutions that produce these surfaces is considered. A cubic-metal surface may be made responsive to polarized light in two basic ways."' In one, an anodic film, believed anisotropic, is formed on the surface. A plane-polarized-light beam, at normal incidence, should be reflected from this type of film with an ellipticity that varies with grain orientation. Under crossed polarizers of a polarized-light microscope, each grain may be distinguished from its neighbors by a difference in reflected-light intensity. An alternate treatment,' more generally applicable to cubic metals and of principal concern here, involves etching of grain surfaces. The characteristic features of this surface were first deduced by Jones3 from light-microscope observations. She concluded that, in general, grain surfaces were furrowed so that light was reflected from two parallel sets of etched facets. A grooved surface produces elliptical polarization of a plane-polarized beam because the beam does not strike groove walls at normal incidence. Jones also observed that the furrows must have faces inclined to each other by approximately 90 deg, since the light returns along the incident path, and that a grooved surface should show four maxima and minima of reflected-light intensity during a 360-deg stage rotation of the polarizing microscope. Positions of maximum extinction were predicted to occur when groove axes were parallel to either polarizer or analyzer vibration directions. Because of the limited resolving power of Jones' optical microscope, her deductions were primarily indirect. Proof of the correctness of her conclu-sions, as demonstrated by the electron microscope, will be given as well as evidence concerning the crystallographic nature of the etch furrows and the types of etching solutions that produce them. EXPERIMENTAL PROCEDURE During a comprehensive study of hot plastic deformation in nickel and nickel alloys, an etch was evolved that produced a surface with an excellent polarized-light response on Nickel 200 ("A" Nickel). The etching solution and associated metal-lographic procedures are given in Table I. In evaluating this etch, a study was made of the topological features of the etched surface and their relation to the underlying crystalline structure. As part of this investigation, large crystals (2 mm average diameter) were grown in a 1-cm-sq Nickel 200 specimen. After the surface was etched, the crystallographic orientations of ten grains were determined by the standard back-reflection Laue technique of Gren-inger.4 Maximum extinction positions during a microscope stage rotation were also measured for the ten grains. Groove-wall positions on the surfaces of the ten grains were measured with a two-circle optical goniometer. The technique was essentially that of Barrett and Levenson.5 Several grain surfaces were photographed with a Philips Model 100A electron microscope. All specimens were replicated with collodion, or collodion backed with formvar and chromium shadowed at 18 to 20 deg. The basic material was Nickel 200. However, an electron-micrographic study was also made of surfaces developed by polarized-light etches on other fcc metals (90-10 a brass, Monel 400, and lead). All etching procedures are given in Tables I and 11. EXPERIMENTAL RESULTS Fig. 1 shows typical electron micrographs of three different fcc metal specimens and an optical micrograph of a fourth. All photographs show a grooved structure corresponding closely to Jones'3 predictions. Also, as noted by Jones, extinctions were always observed when groove axes were either parallel or perpendicular to the microscopes' vertical cross hair. The symmetry of the grooves, with respect to the twin boundaries in Fig. 1, implies that furrows have a crystallographic basis. The coarse-grained Nickel 200 specimen was used to study this basis. Facet Orientation. The poles of the ten Nickel 200

Jan 1, 1964

By Ralph Hultgren, Prodyot Roy

Manganese vapor pressures from 1250° to 1500°K were measured by conventional Knudsen and torsion-effusion methods in twelve Fe-Mn alloys with compositions from 9 to 80 at. pct Mn. The Knudsen re-sults agveed approximately with previous measurements found in the literature which indicated the solutions were nearly ideal. However, except for the higher manganese compositions, the initial torsion readings indicated much (up to 50 pct) higher vapor presszcres than the Knudsen method. These high readings decreased steadily with time. The results are interpreted as due to depletion of surface concentration of manganese during evaporation. Thus the initial torsion readings are the most nearly correct and Knudsen methods on alloys are to be regarded with suspicion unless it can be demonstrated that diffusion rates are rapid enough to replenish surface concentrations depleted by evaporation. For the lowest manganese contents and highest temperatures, depletion causes initial torsion readings to be too low. FEW thermodynamic data are available for alloys of the high-melting transition metals due primarily to experimental difficulties at the high temperatures involved. Measurement of vapor pressure is one of the most promising techniques to be applied to this problem. From the vapor pressure of a component of an alloy phase, its activity or its partial molar Gibbs energy can be directly calculated. From measurements over a sufficient range of temperatures and compositions the partial and integral Gibbs energies, enthalpies, and entropies may be determined through Gibbs-Duhem integration. For solid phases of variable composition, which commonly are found in alloys, measured vapor pressures are low because of depletion of the surface concentration of the more volatile component as it is selectively vaporized. Diffusion from the interior tends to restore the depleted concentration. The seriousness of the effect depends on the relative rates of diffusion and vaporization. Although depletion has been recognized as a factor,' direct measurements of its rate are lacking and quantitative estimates of its effect on measurements are uncertain. In the present work the vapor pressures above a series of Fe-Mn alloys have been measured by the torsion-effusion method. This method permits a measurement of the vapor pressure of manganese as soon as the sample comes to temperature and continuously throughout the run. The rate of decline of the apparent vapor pressure measures the rate of depletion and the torsion reading at zero time should be correct. At high manganese concentrations (70 and 80 at. pct Mn), the torsion reading remained nearly constant; surface depletion was negligible. There was a slow drop in measured pressure due to bulk loss of manganese from the sample. However, for lower manganese contents, especially at high temperatures, the depletion effect was considerable and the apparent vapor pressure decreased steadily with time. For these alloys the Knudsen method should give pressures which are too low. The measurement of pressure in the Knudsen method cannot be made until enough manganese has been vaporized to be weighable. The Knudsen result is therefore an integrated average between the initial and final vapor pressures. To verify this, conventional Knudsen measurements were also made with results in approximate agreement with previously published Knudsen measurements for these alloys. As expected, for the lower manganese contents the initial torsion readings of pressure were as much as 50 pct higher than the Knudsen. THE APPARATUS The apparatus, shown in Fig. 1, is capable of operating to temperatures up to 2000°K. A vacuum of 5 x 10"6 mm Hg can be maintained. Temperatures were measured to +3° by a W-Re thermocouple placed in the dummy cell, N. The torsion cell consisted of two alumina crucibles with holes in their covers. They were placed in the molybdenum holder shown suspended from the torsion wire. Rectangular torsion wires (1 by 4 mm) were found to have superior sensitivity and less residual distortion than circular (2-mm-diameter) wires. The precision of angular measurement was found to be *0.025 deg; the variation is no doubt due to temperature fluctuation and vibration. For the Knudsen experiments the alumina cell was set on the support, N. SAMPLE PREPARATION 800-g samples were prepared by melting together electrolytic iron and electrolytic manganese and pouring the melt into a water-cooled copper mold. Melting and pouring were done under a helium at-

Jan 1, 1965

By R. A. Wattenbarger, H. J. Ramey

A systematic study has been made of the application of the real gas pseudo-pressure m(p) to short-time gas well testing. The m(p) function can be used in real gas flow problems to account for the variation of viscosity and z-factor with pressure. A mathematical model was solved numerically to generate solutions of various real gas flow problems. The effects of turbulence, formation damage and wellbore storage were included in the model. The analysis of simulated well tests showed that the interpretation methods normally used for liquid flow are generally accurate when the m(p) is used. For practical rates without turbulence, the solutions are not rate-sensifive, and the flow capacity.. kh and skin effect can be determined accurately from either a buildup or a drawdown test. However, the kh calculated from a drawdown test can be significantly low when turbulence is present. A case was simulated in which this error wa.s 36 percent. Turbulence does not affect the determination of kh from buildup tests. The proper determination of kh from p and pa buildup and drawdown plots is developed by analyzing their relationship to the m(p) method. A simple equation is given that can be used for long-range gas well performance forecasting. This expression is compared with a method pre-sended by Russell et al." Introduction The formation flow capacity and wellbore damage condition can be determined for liquid-producing tests by means of buildup and drawdown tests. These tests make use of short-time pressure transient data rather than stabilized flow tests that are often used for gas well testing. Tracy' presented a method of gas well testing that was based on ideal gas equations and utilized p2 buildup plots. This method was shown to be good for low pressure wells. MatthewsZ suggested plotting p and using an average slope and average gas properties. This method was more successful on high pressure wells. Al-Hussainy, Ramey and Crawfordh howed that variation in gas properties could be simplified by using the real gas pseu do-pressure m(p). Some cases showed that the use of the m(p) function Original manuscript received in Society of Petroleum Engineers office Aug. 8. 1967. Revised manuscript received June 11, 1968. Paper (SPE 1835) was presented at SPE 42nd Annual Fall Meeting held in Houston Tex., Oct. 1-4. 1967, and at SPE 4th Annual Eastern Regional ~eetin$ held in Pittsburgh. Pa., Nov. 2-3. 1967. 0 Copyright 1968 American Institute of Mining, Metallurgical, and Petroleum Engineers, Inc. 'References given at end of paper. provided an accurate method of interpreting buildup and drawdown tests. Because the work of Al-Hussainy et al. included only a few actual cases, there was a need to explore the use of the m(p) method for a greater variety of flow conditions. The occurrence of turbulent flow around the wellbore often is an important factor in gas well testing. Swift and Kiel,' and Carter et a1.' treated turbulence in gas well testing, but did not consider the variation of viscosity and z-factor with pressure. This paper presents the results of an investigation of the application of the m(p) method to buildup and drawdown testing. Vhe investigation included the effects of turbulence, formation damage and wellbore storage. Gas Flow Equations To formulate the mathematical model, many of the assumptions usually used in well testing theory are applied. The system has radial geometry with a closed outer boundary and is composed of a horizontal porous formation that has uniform and isotropic rock properties and uniform thickness. Allowance is made, however, for a radial region of reduced permeability near the wellbore. This region represents formation damage. The geometry of the system is shown in Fig. 1. Darcy's law does not always apply to gas flow. A more general expression is needed for non-Darcy, or "turbulent" REGION OF DAMAGED PERMEABILITY SEALED BOUNDARIES ClRCULAR WELL BORE Fig. 1—Radial flow model. AUGUST, BT7

Jan 1, 1969