Search Documents

By T. B. Massalski, C. S. Barett

THAT metals and alloys of the body-centered-cubic structure tend to become unstable at low temperatures is so nearly universal that any exceptions are worthy of special attention. Studies of the exceptions may disclose unrecognized examples of metastability or may shed light on the nature of some of the factors that stabilize the body-centered-cubic phase. Among the structurally analogous body-centered-cubic 3/2 electron compounds, the ß-phase in the Cu-Zn system deserves further attention for, although this phase is known to undergo partial mar-tensitic transformation on cooling below room temperature to the temperature of liquid nitrogen when the composition is less than about 42 wt pct Zn, it does not transform on cooling .when there is more than this amount of zinc present.1-4 The present authors recently found that an alloy containing 48.35 wt pct Zn could be transformed, however, if deformed severely at room temperature or lower temperatures.V his discovery led to the suggestion that the jerky flow of ß-brass single crystals6 and some of the unusual temperature variations of the yield strength and hardness are due to strain-induced transformation during the tensile tests. The present investigation is an extension of that work to cover the entire available range of compositions of Cu-Zn ß-phase alloys and certain interesting compositions of the ß-phase in the Cu-Zn-Ga system. Particular attention is given to the structure of the transformation product resulting from cold work, which is found to vary with alloy composition, and to the temperature range in which transformation can be induced in the various alloys. The composition limits of the body-centered-cubic phase in the Cu-Zn system, based on a recent review by Raynor7 and the precision determinations of Beck and Smith: are shown in Fig. la; in the ß region the alloys are disordered and in the ß region they are ordered. The rate of ordering is so rapid that a high degree of order is found in any of these alloys immediately after quenching. Experimental The nine binary and four ternary alloys of this investigation were prepared from cathode copper (99.99 pct purity), electrolytic zinc (99.99+ pct), and gallium (99.98+ pct). All alloys were cast in transparent Vycor tubing, homogenized at temperatures between 800° and 890°C, and quenched in water. Ternary alloys were cast in small amounts (-5 g) to predetermined compositions to fall along a line of constant electron concentration in the ternary system. Changes of weight during the casting operation of the ternary alloys were not greater than 1 part in 2500 from the intended nominal weights and therefore their compositions were accepted without chemical analyses. All binary alloys were cast in larger amounts (100 to 150 g), homogenized, and hot rolled after preheating to 800°C. The central portion of the rolled strip of each alloy was annealed for eight days and quenched in brine. All binary alloys were analyzed chemically. The compositions, heat treatments, and metallographic appearances of all alloys are tabulated in Table I. Various samples were subjected to heat treatment or cold work at several temperatures and were then examined metallographically or by means of X-rays. Specimens about 6x6x3 mm in size were examined metallographic ally to confirm the presence of the ß-phase after quenching. Cold work was introduced by fastening each sample between two pieces of steel, cooling the assembly to the desired temperature, and transfering it rapidly to a vise where moderate or severe compression was immediately applied. Cold working at liquid helium temperature was done by a chisel-shaped tool mounted on an X-ray spectrometer, described elsewhere.% A!! metal-

Jan 1, 1958

By D. L. McElroy, J. L. Schroder

Although the earliest use of rubber-tired haulage was in Illinois in 1936, the first unit of this type of equipment used in West Virginia was shipped into the state in 1938. All units placed in West Virginia mines have been of the shuttle-car type. The Katherine Coal Mining Co., at Lumberport, W. Va., was the second company in the United States to install shuttle cars and the Gay Coal and Coke Co., at Mt. Gay, W. Va., was the first company in the country to install the thin-bed type of shuttle car. On Oct. 8, 1940, nine companies in the state were using 41 shuttle cars, of which 10 cars at three mines were of the thin-bed type. All shuttle cars of the low type were in the southern fields and all of the high type were in northern West Virginia. At that time several units were in the process of delivery. The data in this paper were collected in the field by Mr. Schroder in connection with graduate work on mechanical loading of coal in West Virginia by mobile loaders. In the course of his field work he visited five mines in the state that were using shuttle cars. All of these mines were in high-volatile coal beds of strong structure. The more pertinent data obtained at the mines visited are given in Table I. The maximum distance of shuttle-car haul is the distance between the transfer station and the farthermost place. Tonnage per loading unit per foot of bed thickness is to some degree a measure of relative efficiency, especially if the different operations are working the same bed or have somewhat similar conditions. The thin-bed operation is well up in the higher bracket of the afore-mentioned column, which is as

Jan 1, 1942

By John L. Schroder, D. L. McElroy

ALTHOUGH the earliest use of rubber-tired haulage was in Illinois in 1936, the first unit of this type of equipment used in West Virginia was shipped into the state in 1938. All units placed in West Virginia mines have been of the shuttle-car type. The Katherine Coal Mining Co., at Lumberport, W. Va., was the second company in the United States to install shuttle cars and the Gay Coal and' Coke Co., at [ ] Mt. Gay, W. Va., was the first company in the country to install the thin-bed type of shuttle car. On Oct. 8, 1940, nine companies in the state were using 41 shuttle cars, of which 10 cars at three mines were of the thin-bed type. All shuttle cars of the low type were in the southern fields and all of the high type were in northern West Virginia. At that time several units were in the process of delivery. The data in this paper were collected in the field by Mr. Schroder in connection with graduate work on mechanical loading of coal in West Virginia by mobile loaders. In the course of his field work he visited five mines in the state that were using shuttle cars. All of these mines were in high-volatile coal beds of strong structure. [ ] The more pertinent data obtained at the mines visited are given in Table I. The maximum distance of shuttle-car haul is the distance between the transfer station and the farthermost place. Tonnage per loading unit per foot of bed thickness is to some degree a measure of relative efficiency, especially if the different operations are working the same bed or have somewhat similar conditions. The thin-bed operation is well up in the higher bracket of the afore-mentioned column, which is as

Jan 1, 1941

By J. L. Schroder, D. L. McElroy

Although the earliest use of rubber-tired haulage was in Illinois in 1936, the first unit of this type of equipment used in West Virginia was shipped into the state in 1938. All units placed in West Virginia mines have been of the shuttle-car type. The Katherine Coal Mining Co., at Lumberport, W. Va., was the second company in the United States to install shuttle cars and the Gay Coal and Coke Co., at Mt. Gay, W. Va., was the first company in the country to install the thin-bed type of shuttle car. On Oct. 8, 1940, nine companies in the state were using 41 shuttle cars, of which 10 cars at three mines were of the thin-bed type. All shuttle cars of the low type were in the southern fields and all of the high type were in northern West Virginia. At that time several units were in the process of delivery. The data in this paper were collected in the field by Mr. Schroder in connection with graduate work on mechanical loading of coal in West Virginia by mobile loaders. In the course of his field work he visited five mines in the state that were using shuttle cars. All of these mines were in high-volatile coal beds of strong structure. The more pertinent data obtained at the mines visited are given in Table I. The maximum distance of shuttle-car haul is the distance between the transfer station and the farthermost place. Tonnage per loading unit per foot of bed thickness is to some degree a measure of relative efficiency, especially if the different operations are working the same bed or have somewhat similar conditions. The thin-bed operation is well up in the higher bracket of the afore-mentioned column, which is as

Jan 1, 1942

By R. L. Manegold

IN RECENT years there has been a decided trend toward bringing ores and coal out of open-pit and underground mines by long, sloping single-stage belt conveyors. Because the high investment cost of haulage conveyors dictates maximum speed and full loading for economy, conveyor belts are running faster than ever before, with a greater burden on the belts. In addition, larger blocks of ore are being mined and conveyed. Rougher use of large mining equipment generates more tramp iron; mechanical loading rejects nothing. Therefore the extraction of tramp iron and other deleterious materials for the protection of haulage conveyors at transfer points, chutes, and idlers is becoming an increasingly difficult job. In our experience, the best method of protecting long haulage conveyor belts is to extract the tramp iron, at least the large tramp iron, ahead of the main haulage conveyor or series of conveyors. It is most practical to collect the tramp iron at the earliest source -possible, usually directly after the primary car unloading dump or after the primary crushers, which are ordinarily not affected by tramp iron, see Fig. 1. The coal or ore from its primary source is fed by chutes or feeders onto a buffer or shock belt conveyor, a short, ruggedly built conveyor system that can be considered somewhat expendable. Its primary function is to bring the ore up to speed for transmission to the haulage conveyor. There is no assured method of protecting this conveyor from tramp iron damage at its feed end. Likewise, there is a possibility, although remote if the magnet applied to this belt is installed correctly, that the tramp iron extraction process will damage the shock belt, which is, however, short and comparatively inexpensive. Magnetic head pulleys with their maximum effective range of about 8 or 10. in. cannot be considered for large installations, even though they are desirable because of the continuous and, automatic tramp, iron removal feature. Further, the most dangerous form of tramp iron with respect to potential damage to a haulage conveyor is the long rod, roof bolt, or rail segment, which is not extracted satisfactorily because of the limited tangential contact of such an iron shape with the magnet. Magnetic head pulleys are also undesirable because the dribble or fines from the head pulley discharge at the same point as the tramp iron. Screens must be used to separate the fines from the magnetics, resulting in additional expense and loss of head room. Magnetic detectors have been used with success, but they only indicate the, presence of iron, and the conveyor belt must be stopped for its manual removal. The, stoppage of any conveyor means that all preceding equipment must be stopped or a cascade of material on the idle conveyor will result. Occasionally combinations of detectors and magnetic, head pulleys are used. The detector precedes the magnetic head pulley and detects only large pieces of tramp iron, while the magnetic pulley subsequently takes out the smaller ferrous metals. Suspended magnets overcome the disadvantages of magnetic head pulleys, detectors or a combination of both. There are several magnetic extraction problems involved in applying a suspended magnet at this

Jan 1, 1952

By L. C. Brown, C. St. John, C. C. Sanderson

The diffusion rate in Cu-Te couples is very sensitive to compressive stress, with a load of 20 psi making a significant difference to the width of the diffusion zone. At zero stress, two phases appear in the diffusion zone (Cu4Te3 and CuTe). Under compressive loading the third stable phase (Cuz Te) also appears, and its thickness increases progressively with increasing stress. The results are explained on the basis of an incipient Kirkendall porosity which restricts the transfer of atoms from the copper into the diffusion zone. DURING a study of the Kirkendall effect in Cu-Te couples prepared by clamping together the two components, it was found that the diffusion-zone width and shape in the plane of contact were not reproducible. Although the stresses involved in clamping are not normally sufficiently high to affect diffusion rates, preliminary tests established that the Cu-Te system is particularly stress-sensitive. The phase diagram for the system Cu-Te given in Hanssen1 shows that there is practically no solid solubility at either end of the phase diagram. Many areas of the diagram are not fully substantiated, but there appear to be three intermediate phases: Cu,Te—hexagonal in structure, having a grey luster; Cu4Te3—a tetragonal defect structure, having a red-purple luster; CuTe—orthorhombic in structure and having a golden-green luster. The existence of a fourth phase, the X phase at 37 at. pct Te, is considered doubtful. The composition ranges of the three stable phases are small, and are not accurately known. The phase diagram changes little with temperature up to 305°C, at which temperature a polymorphic transformation takes place in Cu2Te. The nature of the Cu-Te phase diagram indicates that the diffusion zone in a Cu-Te couple would consist of a series of layers of intermediate phases. The relative thickness of any one phase will depend on its diffusion coefficient and composition range.' In this type of diffusion couple it is often found experimentally that some phases are not visible at all in the diffusion zone due either to a small diffusion coefficient or to a restricted composition range.3 Since the composition ranges of the phases in Cu-Te are not known, it is not possible to determine diffusion coefficients in this system from a knowledge of the phase thicknesses. Several investigations have been carried out to determine the effect of compressive stress on diffusion rates in multiphase systems. Diffusion couples of Ni-A1 have been investigated by Storchheim et al.4 and by Castleman and Seigle.5 Two phases (ß and ?) appear in the diffusion zone under zero stress and the thickness of both phases is progressively reduced with increasing stress. According to Storchheim et al.4 a stress of 25,000 psi reduces the thickness of the diffusion zone by 50 pct. In a-brass—?-brass couples the thickness of the 0 phase formed in the diffusion zone was reduced by 20 pct at a stress of 20,000 psi.6 In other investigations the compressive load has been observed to increase the width of the diffusion zone. In A1-U, several investigators3,8 have found the width of the whase UA13 to increase with stress. According to casileman,8 the rate of formation of UA13 at 520°C is 75 pct faster at a stress of 20,000 psi as compared with a stress of 2500 psi. In Cu-Sb the effect of stress is greater than in the other systems described. According to Heumann9,10 only one phase (y) appears in the diffusion zone at a stress of 500 psi, but at a stress of 850 psi two phases (y and k) are present. If a diffusion couple containing both y and k phases is annealed at a low stress level, the y phase grows at the expense of the k phase. EXPERIMENTAL The diffusion couples were prepared from electrolytic copper bar stock with a nominal purity of 99.92 pct and from tellurium of 99.7 pct purity. The tellurium proved difficult to machine because of its brittleness and a technique was developed for casting the tellurium into a graphite slab mold and spark-machining specimens from this slab. Both the copper and tellurium were produced in the form of discs 2 in. diam by approximately 1/4 in. thick with surfaces ground flat to 3/0 emery paper. The diffusion apparatus is shown in Fig. 1. Auni-axial compressive stress was applied to the system through a simple lever system. A stainless-steel rod actuated by the lever arm lay inside a stainless-steel tube. The diffusion couple lay on top of the steel rod, and pressure was applied to the couple between the rod and a plug welded into the center of the tube. To ensure a uniform stress across the couple, a hemispherical boss and cup were used to transmit the load to the diffusion couple. A 400-w tube furnace with a uniform hot zone 3 in. long slid around the stainless-steel tube and maintained the assembly at temperature. A thermocouple situated 3 in. from the specimen operated a proportional temperature controller which maintained the specimen temperature constant to ±2°C. Most diffusion runs were carried out at 250C although a few tests were made at other temperatures in the range 235° to 300°C. The specimens were inserted and removed with the furnace at operating temperature, and took only 2 min to reach diffusion temperature—a time small compared with the total diffusion time. All the diffusion experiments were carried out in a hydrogen atmosphere, since consistent results were obtained in hydrogen and nitrogen atmospheres and in

Jan 1, 1967

Theoretically mobile primary crushers at the mine working face can be considered the ideal comminution system. A crusher which follows the ore can reduce materials handling costs because it eliminates about one ton of truck deadweight per ton of ore hauled to the crusher. It also eliminates the backhaul of this same deadweight ton. Development of mobile primary crushers has centered in Germany where many soft iron ore and limestone deposits are flat-lying and thus readily adaptable to the use of shiftable conveyors which easily followed the working face. These ores were easily crushed in impact- type mills which had extremely large receiving openings and therefore large reduction ratios. The feeding of such a unit was done with a front-end loader which dumped its load to an inclined apron feeder, which in turn fed the impactor. Under the crusher was an inclined conveyor which transferred the crushed ore to a belt wagon for delivery to the shiftable overland conveyor. Low headroom overhead eccentric single-toggle jaw crushers and short-shaft gyratories were also used on these units, which generally travelled on tracks. While these early mobile primary crushers were versatile and served their purpose, they generally suffered from low capacities. Manufacturers have quickly responded with larger crusher, feeders, and moveable supports, however; and today mobile crusher s have capacities of up to 1500 metric tons per hour, a figure comparable to the capacity of some 54 in. by 74 in. gyratories. An example of a quarry portable crushing plant is at the Bath cement plant of Canada Cement Lafarge, Limited. Two 12 cu yd front-end loaders are digging the stone and hauling it to the portable crushing plant which feeds a field conveyor system. The crushing plant is moved and the field conveyors extended whenever the length of the haul exceeds 300 ft. At the crushing plant, the loaders discharge the stone onto a pan feeder which feeds into the hopper of a Weserhutte-Hazemag compound crusher which combines primary and secondary rotors in the same housing. These are driven at 1200 rpm by 800 hp and 1,000 hp motors (600 and 750 kw). The maximum 3 ft feed (914 mm) is reduced to 3/4 in. (19 mm) at 800 tph. The Weserhutte steered walking mechanism, operated by hydraulic cylinders, takes five -foot strides. Another example is at Alcoa of Australia's operation at Huntley in Western Australia. At the mine face, front-end loaders load 50 mt rear-dump trucks that transport the bauxite ore to the mobile crusher which follows the mine operation. The crusher can turn 360° within its own length, walk at a speed of 34. 5 ft per minute (10.5 m per

Jan 1, 1978

By Wei-Kao Lu

A general rate equation is derived for gas -solid reactions in metallurgical processes by considering the contributions of chemical reaction at inter-phase boundaries and diffusion through the solid product layer simultaneously. This equation takes care of the whole range continuously from purely diffusion control to purely chemical reaction control. The parabolic law,' McKewan's equationa derived with the assumption of interface area control and the equation derived by Kawasaki et al. with the assumption of purely diffusion control are found as the limiting cases of the corresponding forms of the general rate equation. Many metallurgical processes involving gas-solid reactions have been extensively studied. The common examples of these processes include the reduction of oxides to metallic products, the roasting of sulfides to oxides or chlorides and the oxidation of metals. Gas-solid reactions consist of at least two kinds of steps; one involves the transport of reactant and product to and from the reaction zone and the other is the chemical reaction taking place at the interphase boundaries. In the limiting cases one or the other of the above mentioned steps will be much slower than the other and is then rate controlling. For example, the oxidation of metals is diffusion controlled and obeys the parabolic law.' On the other hand, if the diffusion rate is much larger than the chemical reaction rate then it will be controlled by the chemical reaction rate. This latter case occurs for the reduction of iron oxides studied by McKewan. If the rate of diffusion and of the chemical reaction have the same order of magnitude, a more detailed analysis is necessary. Thus at the beginning of the reaction when the solid and gaseous reactants are in direct contact, the reaction should be chemical reaction rate controlled as no diffusion process is involved. At the other extreme as the reaction zone moves further into the solid, the diffusion distance becomes greater and the diffusion rate continuously decreases. Eventually, therefore, for a sufficiently large sample the rate must become diffusion controlled. As the transition from one controlling step to the other must be continuous, we cannot describe this transition period by either of the two limiting cases. In this paper the author attempts to derive a general rate equation for gas-solid reaction systems. This is effected by considering the rates of diffusion and chemical reaction simultaneously from the state of pure diffusion control to pure chemical reaction control. I) DERIVATION The system under consideration is made of at least three phases, namely the gaseous one, the solid reactant, and the solid product as shown in Fig. 1. The notations are listed at the end of the paper. The mass transfer process involves the diffusion of reactants and products across the solid product layer. It is reasonable to consider that no chemical reaction will take place in the solid product layer with respect to both gaseous reactant and product. For the characteristics of local equilibrium of the system, the diffusion process is really in a quasi-steady-state. It will be proved in the following paragraphs. Then we can invoke the law of conservation of mass valid for all chemical species passing through the layer. The solid reactant sphere has an initial radius, ro, after a fraction of completeness of the reaction, the unreacted core has a radius ri, where r, is in the solid product region. Let Cg be the activity of gaseous reactant in the

Jan 1, 1963

By Seymour G. Epstein

Using a modified capillary-reservoir technique, electromigration of dilute solutions of silver in liquid bismuth has been measured with relatively good precision. The effects of experiment duration, temperature. current density, and silver concentration have been investigated. Chemical diffusion coefficients of silver in liquid bismuth were also measured at temperatures ranging from 300°to 600°C. The effective valence of silver in bismuth was calculated as a function of temperature from a comparison of electric mobilities and diffusion coefficients. In an alloy containing 1 wt pct Ag, the electric mobility of silver increased from 1.65 x X sq cm per v-sec at 300°C to 2.38 x 10-4at 600°C; the diffusion coefficient in-creased from 3.0 to 17.5 x 10-5 sq cm per sec over the same temperature range. Calculated values of effective valence decreased from 0.28 at 300°C to 0.10 at 600°c. ALTHOUGH there has been considerable progress in the theoretical treatment of electron-transport properties in liquid metals, the mechanisms of mass-transport phenomena in liquid metals are still not clearly understood. This is especially true for the case of electro migration, alternatively termed elec-trotransport or electrodiffusion, which is mass transport induced by an applied electric field. A thorough knowledge of the mechanisms governing electromigration will considerably aid in understanding the electronic state of an impurity atom in a liquid metal. Electromigration has been long observed and studied in liquid metals and alloys. However, no satisfactory theory to explain the observed effects has been developed. The disparity in the present theoretical treatments, summarized by Verhoeven,1 has been attributed to the difficulty in describing the momentum transfer between electrons and ions and to the lack of precise measurements in liquid metals, especially at elevated temperatures. Prior to the last decade most of the electromigration studies were qualitative, in that only the direction of migration of alloy constituents was determined with any degree of reliability. These results have been summarized by schwarz,2 who was able to measure the electric mobility of several elements in mercury. More recently electromigration in liquid metals has been measured by a variety of techniques, with a varying degree of success. Mangelsdorf3,4 made precise measurements on dilute amalgams at room temperature, using a technique relating changes in composition to changes in alloy resistivity, but was unable to make measurements at elevated temperatures. Using a modification of this technique, Verhoeven and Hucke5 measured electromigration and resistivities in liquid Bi-Sn alloys. However, the resistivity technique is not sensitive to small composition changes and is impractical for very dilute alloys. Verhoeven and Hucke6 also measured electromigration of several solutes in liquid bismuth by adapting the capillary-reservoir method for determining diffusion coefficients in liquids. The reported precision of these measurements was poor, however, primarily due to the large uncertainties in the chemical analyses to determine concentration changes. Belashchenko7 has summarized the recent Russian investigations of electromigration, many of which, however, are of questionable accuracy. The present study was initiated to develop an experimental technique which would be both versatile and sufficiently sensitive to yield precise measurements of electromigration of solutes in dilute metal alloys at elevated temperatures. Using this technique, a comprehensive study of silver in liquid bismuth was made to determine the effects of experimental variables on the apparent electric mobilities. Chemical diffusion coefficients of silver in liquid bismuth were also measured over the same temperature range for comparison with the electromigration data. EXPERIMENTAL TECHNIQUE Electromigration Measurements. Electromigration of silver in liquid bismuth was measured by the capillary-reservoir technique, similar to that used by Verhoeven and Hucke,6 but modified to yield more

Jan 1, 1967

By James H. Hance

Discussion of the paper of JAMES H. HANCE, presented at the New York meeting, .February, 1916, and printed in Bulletin No. 110, February, 1916, pp. 299 to 326. FREDERIC P. DEWEY, Washington, D. C. (communication to the Secretary*).-I wish to congratulate Professor Hance and to thank him for his extremely valuable addition to the fund of data upon this subject. It presents clearly conclusive and incontrovertible proof of the unreliability of drill samples of this class of metal for determining the value of such bullion. I am afraid, however, that his hope that some one in the mint service might continue such an investigation upon a more elaborate and exhaustive scale is, for the present at least, doomed to disappointment. As a practical matter, the mint service is not now particularly interested in segregation in crude bullion. So much trouble and friction arose, -both in the original purchase and in the transfer from one institution to another, in handling segregated bullion, that it became necessary to adopt drastic measures and now it is seldom attempted to determine the value of such bullion in the service. When bullion known or even suspected of being segregated is presented for purchase, it is at once strongly refined in the pot and the disturbing elements so far removed that concordant assays may certainly be obtained from the dip samples and a reasonably fair agreement obtained from drill samples. The mint service prefers that the owner should refine his bullion himself, but if he does not we will do it at his expense and return the slag to him when desired. Much work has been done in the mint service trying to arrive at the value of deposits of segregated bullion without strong refining. It was no unusual thing for such deposits to be melted three times at the purchasing office and to have 30 to 40 assays on the metal made before buying it. Yet when shipped to a mint it was found necessary to melt two or three times more and to make 15 to 25 assays before the mint would accept the bar, and there were considerable differences between the valuations of the two offices. In what was undoubtedly the worst case ever investigated, the purchasing office melted the deposit, 953 oz., three times, with a loss of 105.5 oz., and made 46 assays; the mint melted three times, with a loss of 6.6 oz., and made 28 assays. The mint allowed the purchasing office 1.178 oz. of fine gold more than the purchasing office claimed, while an extensive investigation of the mint samples by the Bureau, in which over 100 independent assays were made in various service laboratories, indicated that the metal carried 0.7 to 0.9 oz. more of fine gold than the mint allowed. Preliminary assays on the original metal

Jan 5, 1916

By R. C. McFarlane, T. D. Mueller, J. E. Walstrom

In evaluating uncertainty, experiments are usually performed repeatedly and then conclusions are drawn from the distribution of results. With the advent of high-speed electronic computers, it is possible to perform experiments using mathematical models constructed to simulate complex experiments or operations. Statistical methods are then applied to the results of the simulated experiments. This procedure forms the busis of this paper. Demonstrated is the need for properly accounting for uncertainty in petroleum engineering problems. How uncertainty affects solutions is evaluated in three example illustrations. The method used to evaluate uncertainty in petroleum engineering studies is the Monte Carlo simulation procedure.'-" INTRODUCTION The solution to most technical problems may be derived from interrelationships among several quantities called variables or parameters. There may be only a few variables or several hundred. Interrelationships among parameters may be explicit or implicit, well established or only approximate. Some variables that fully or partially depend on the magnitude of others are called dependent variables. Input variables for most practical problems are not precisely known; there is usually an uncertainty in their value. The degree of uncertainty may vary from one variable to another. Variables that are known accurately are called determinates.' For instance, the gravity of crude obtained from a particular pool may be known precisely, and therefore is a determinate. The degree of precision with which a quantity can be determined increases as data describing the pool are accumulated during the development of the field and the producing life of the pool. The uncertainty of a parameter may result from difficulty in directly and accurately measuring the quantity. This is particularly true of the physical reservoir parameters which, at best, can only be sampled at various points, and which are subject to errors caused by presence of the borehole and borehole fluid or by changes that occur during the transfer of rock and its fluids to laboratory temperature and pressure conditions. Uncertainty may also result in attempting to predict future parameter values. This type of uncertainty is particularly evident in investment analyses involving future costs, prices, sales volumes and product demand. Uncertainty in the solution to investment problems is often called risk, and its study is called risk analysis.' Uncertainty also enters into biological and sociological analyses in which indeterminate factors are often important due to limited control of the experimental material. It is customary, in evaluating uncertainty, to perform repeated experiments and to draw conclusions from the distribution of the results of these experiments. With the advent of the high-speed electronic computer, it is possible to construct mathematical models which simulate complex experiments or operations and to perform the experiments repeatedly, utilizing the models. Statistical methods are then applied to the results of the simulated experiments This method forms the basis of the investigation reported here. PROBABILITY DISTRIBUTIONS FOR VARIABLES The uncertainty in the value of a variable may be indicated by a probabilistic description accomplished by expressing the quantity by a probability distribution. Many recognized probability distributions can be used to describe physical quantities. Recent studies used various types of distributions to describe core analysis data.',' However, for the examples in this paper, the uniform and triangular distributions are believed to reasonably approximate the data used (Fig. 1). The uniform distribution confines the variable between an upper and a lower limit. The variable may lie anywhere between the two limits. This distribution is used when no one range of values for a variable is more probable than any other, but information or intuitive reasoning indicates the variable will lie somewhere between the chosen limits. The triangular distribution is used for a variable when more data are available to indicate a central tendency of distribution. This allows postulating a "most likely" value to the distribution and upper and lower limits. In this case, as for the uniform distribution, the variable is not expected to assume a value less than the lower limit or greater than the upper limit. However, with improved quality of data it can be postulated that the variable will tend to assume a value close to the most likely value, and that there will be a decreasing probability for values away from the most likely value. The area under either of these probability distributions is equal to unity since it is assumed that there is a 100 percent probability that the variable will lie somewhere under the curve. An ordinate erected at any particular value of the variable divides the area under the curve into two parts: the area to the left of the ordinate represents the probability that the value of the variable will be equal to or less than the value of the variable at the position of the ordinate, and vice versa. The probability is zero that the variable will have any specific deterministic value. If two ordinates are drawn for any two values of the variable, the probability that the variables will have a value lying between these ordinates is equal to the area under the curve lying between the ordinates.

By F. F. Poland

The purpose of this paper is to describe the recent applications and improvements made in the process and equipment for the recovery of metallic zinc from secondary metals by means of high temperature electric resistor furnace distillation. The process and high temperature furnace, as developed by Revere Copper and Brass Incorporated, Research and Development Department, was described in a previous paper. † Since the publication of that paper the process and equipment have been licensed to others. Consequently its field of application has been expanded, particularly for the treatment of scrap brass, galvanizers dross and scrap zinc. Other applications are being considered, particularly the melting of metallic titanium and special alloys requiring a controlled atmosphere. Consideration is also being given to modification in design that would allow the melting of copper cathodes for the production of oxygen-free and electrolytic tough pitch copper shapes without the necessity of blowing and poling. In the design of a furnace for melting copper cathodes it is intended that the fundamental principle of charging the solid material into a large pool of molten metal through a double door charging vestibule as shown in Fig 3 of the previous paper, † would be adhered to and, in addition, either dehydrated and desulphurized producer gas or commercial nitrogen, that is, approximately 96 pct nitrogen, 2 pct hydrogen 2 pct CO, also dehydrated, would be used to provide a protective atmosphere for the resistors and to flush out more soluble gases from the molten copper. The operation of a small unit producing approximately 1000 lb of copper per hr has shown that the commercial nitrogen will allow the production of tough pitch electrolytic copper when the metal is transferred from the meltdown unit contemplated to a special type induction unit where the oxygen absorption and temperature can be controlled. Standard practice in treating low grade scrap brass has been to oxidize the zinc by blowing air through the molten metal. This operation results in the recovery of approximately two-thirds or less of the zinc content of the brass in the form of an impure zinc oxide having only a fraction of the value of the zinc in metallic form. The balance of the zinc is lost to slag, fume and drosses. Fig 1 shows the revised flow sheet of the process for the treatment of scrap brass. For this application and as shown on the flow sheet a low frequency induction furnace has been substituted for the Wilkins-Poland type as a melting unit. This type of furnace has been found the most suitable for the purpose as will be shown by the following discussion. Serap Brass The melting of refinery scrap brass preparatory to distillation presents a melting problem quite different from the usual melting of brass mill scrap for recasting because refinery scrap is a heterogeneous mixture of alloys each having different melting and boiling points, and coated with a variety of corrosion products and foreign materials. In addition, brass and copper plated pieces of sheet iron, stainless steel, ceramic materials, and others are present that do not melt except at temperatures above the boiling point of the liquid bath. In the previous paper it is noted that difficulty was encountered in the melting of refinery grade scrap brass without excessive loss of zinc by volatilization and that 5 to 10 pct razorite flux was used to alleviate this condition. While this flux was successful in preventing the excessive loss of zinc by volatilization it interposed a low heat conductivity layer between the source of heat and liquid charge and resulted in lowering the melting rate of the furnace. This same problem presents itself in all types of furnaces in which the metal is heated by radiation and conduction from above the charge. It has been found in melting this type material that the melting rate per hour or per square foot of hearth area

Jan 1, 1950

By A. S. Odeh

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

By Andrew Cosgarea, William R. Upthegrove, Morteza Mirshamsi

The capillary-reservoir technique was used with lead-210 and cadmium-115m to determine the self-diffiLsion coefficients of both cadmium and lead in the liquid binary Cd-Pb system. The self-diffusion coefficients of pure cadmium and pure lead were obtained and were compared with the theoretical predictions. Good to excellent agrement between the experimental and predicted values was obtained. The self-diffusion coefficients of cadmium were tneasuved in alloys containing 2.50, 9.13, 17.40, 31.00, 45.00, 69.00, and 97.00 lot pct Cd by determining- the amount of cadniiutn-115m which diffused out of a small-bore capillavy into an infinite reservoir during- a given time peviod. Sinzila7-measurements were made with lead-210 to determine the self-diffusion coefficients of lead in these identical alloys. Diffusivities were determined from measurenzents performed in the temperature interval of 290" to 480°C. The results were correlated with the Ar-vhenius equation, and the maximum variation of the equation parameters (Q and Do) was also inrestigated . THE theory of diffusion in liquids, particularly liquid metals, is relatively undeveloped in contrast to that for the gaseous and solid states. Although the practical application of liquid metals as heat-transfer media has become increasingly important, few liquid-metals systems have been investigated. Experimental data of fundamental significance in this field are not readily obtained, which may explain but not justify the present lack of knowledge. What work has been completed is primarily restricted to liquid diffusion of pure metals; little work has been done in liquid-metal diffusion of binary mixtures. A review of liquid-metal diffusion theory and research is available elsewhere.1-4 In an effort to add to the knowledge of liquid-metal systems and to increase the basic understanding of the diffusion process in liquids, a study of diffusion in the binary-liquid system, Cd-Pb, was undertaken. The capillary-reservoir technique5 was employed to measure the self-diffusion coefficients of cadmium and lead in molten binary alloys. Measurements were made with seven selected compositions and over a temperature range from 290° to 480°C. The experimental apparatus consisted essentially of the following items: constant-temperature bath, diffusion cells, capillaries, capillary-filling device, and a radioactive tracer counting system. EXPERIMENTAL APPARATUS Constant-Temperature Bath. A cylindrical steel vessel, 8 in. in diam and 15 in. deep, surrounded by an insulated heating coil was used with a sodium-potassium nitrate salt mixture heating medium. The bath was maintained slightly below the desired control temperature by the furnace-heating element; and a 250-w heater, actuated by a Bayley proportional temperature controller, was utilized for the final control of the temperature. A constant-speed mixer stirred the salt to insure a uniform temperature within the bath. Four calibrated Chromel-Alumel thermocouples were placed at various positions in the salt bath to verify the absence of temperature gradients. The observed temperature variation during any diffusion run was less than 0.l°C. The entire furnace assembly was mounted on four shock absorbers to exclude building vibrations and the stirrer propeller blades were adjusted so not to induce vibrations within the reservoir. A schematic diagram of the furnace and the constant-temperature bath is shown in Fig. 1. Diffusion Cell. The diffusion cells and associated parts were the same, except for slight modification, as the one used by walls1 in this laboratory, and are shown in detail elsewhere.' A graphite crucible, 4 in. long and 40 mm (1-1/2 in.) ID, enclosed in a 60-mm (2-1/4 in.) Pyrex tube cell about 18 in. long, was used as a container for the melt. The reservoir (molten alloy in the graphite crucible) was usually about 2 to 2-1/2 in. deep. Graphite was used because of its satisfactory nature as a refractory material and the low solubility of carbon in molten Cd-Pb alloy.677 The Pyrex cell was closed at the bottom and fitted at the top (open end) with a 2-in. Dresser coupling. A brass flange was welded to the top of the coupling. The upper part of the diffusion assembly was bolted to this flange with an O-ring seal. The lower part of the diffusion cell was supported in a 3-in. brass cylinder which was open to allow for circulation of salt around the cell. The top assembly consisted of two synchronous motors, a drive shaft, a thermocouple well, and controlled-atmosphere inlets and outlets. One motor was used for rotation of the capillaries at a rate of 1/2 rpm in the reservoir during the diffusion run. The other motor was used for the vertical positioning of the capillaries and the capillary holder by means of a simple screw drive. The capillary holder and drive assembly were lowered into the reservoir for the run and raised after the desired diffusion time at a rate of approximately 0.4 in. per min. Capillary holders were made of graphite. These

Jan 1, 1967

By M. R. J. Wyllie

Eight laboratory-prepared aqueous base drilling muds representing common mud types, and 15 aqueous base drilling muds sampled in the field, have been used in an experimental investigation of the relationship between the streaming potential developed across a mud filter cake and the differential pressure causing filtrate flow. It has been established that the relationship is of the form, ES = kPy, where Es is the streaming potential, P the differential pressure and k and y are constants for any particular mud system at constant temperature. Beyond a certain minimum thickness this relationship is independent of filter cake thickness. Experiment has shown that both k and y are temperature dependent parameters and that in certain mud systems there is a critical value of P above which the simple relationship Es = kPY ceases to hold. The constant k tends to be a function of the mud resistivity at 77ºF, but the exponent y has not been found to be related to any commonly measured mud property. Certain field experiments are described which serve to demonstrate that the total SP (self potential) developed opposite a formation may be quantitatively accounted for as the algebraic sum of an electrochemical potential and an electrokinetic potential. A chart to permit the simple use of mean solution activities in the interpretation of electrochemical SP data is given and its use described. Various field methods for computing the electrokinetic SP component of the total SP are considered and rejected as being unduly complicated and inaccurate. It is shown that the employment of low resistivity drilling muds is essential if optimum accuracy is desired in the use of the S.P. curve for the purpose of computing connate water salinities by the electrochemical SP method. The use of such low resistivity drilling muds is accordingly recommended. INTRODUCTION In a previous paper' a quantitative analysis of the electrochemical compovent of the SP (self potential) curve was made and it was shown that the electrochemical* emf, E, could be expressed as E = ZtRT/F In a,/a,, where t is the transference number of the chloride ion, R is the gas constant, T the absolute temperature of the formation water and mud, a, the mean activity of the connate water (assumed to be sodium chloride) and a, the mean activity of the mud (as sodium chloride). The essential accuracy of this expression was confirmed by an investigation 2 of the SP kicks opposite sands of known water salinity. The logs examined had been routine runs in oil wells in a number of Mid-Continent states. Further confirmation of the original analysis of the electrochemical SP is implicit in the paper by Tixier3 wherein a relationship between the SP and the ratio of connate water and mud resistivities, based only on field data, is given. Allowing for the fact that mud and connate water resistivities were used instead of activities, and that electrokinetic effects were ignored, the constant reported by Tixier is essentially 2tRT/F. The error involved in the universal substitution of resistivities for activities will receive further reference below. Possible reasons for the fact that the electrochemical properties of shales do not appear to be significantly dependent upon either the nature or the quantity of cation exchange materials, such as clays, they contain has been considered at some length by Wyllie and Patnode.4 In the same paper, the probable effects of shale thickness and compaction on electrochemical efficiency, and hence the inherent difficulty of obtaining reliable laboratory-scale data, were also pointed out. However, it is

Jan 1, 1951

By M. R. J. Wyllie

Eight laboratory-prepared aqueous base drilling muds representing common mud types, and 15 aqueous base drilling muds sampled in the field, have been used in an experimental investigation of the relationship between the streaming potential developed across a mud filter cake and the differential pressure causing filtrate flow. It has been established that the relationship is of the form, ES = kPy, where Es is the streaming potential, P the differential pressure and k and y are constants for any particular mud system at constant temperature. Beyond a certain minimum thickness this relationship is independent of filter cake thickness. Experiment has shown that both k and y are temperature dependent parameters and that in certain mud systems there is a critical value of P above which the simple relationship Es = kPY ceases to hold. The constant k tends to be a function of the mud resistivity at 77ºF, but the exponent y has not been found to be related to any commonly measured mud property. Certain field experiments are described which serve to demonstrate that the total SP (self potential) developed opposite a formation may be quantitatively accounted for as the algebraic sum of an electrochemical potential and an electrokinetic potential. A chart to permit the simple use of mean solution activities in the interpretation of electrochemical SP data is given and its use described. Various field methods for computing the electrokinetic SP component of the total SP are considered and rejected as being unduly complicated and inaccurate. It is shown that the employment of low resistivity drilling muds is essential if optimum accuracy is desired in the use of the S.P. curve for the purpose of computing connate water salinities by the electrochemical SP method. The use of such low resistivity drilling muds is accordingly recommended. INTRODUCTION In a previous paper' a quantitative analysis of the electrochemical compovent of the SP (self potential) curve was made and it was shown that the electrochemical* emf, E, could be expressed as E = ZtRT/F In a,/a,, where t is the transference number of the chloride ion, R is the gas constant, T the absolute temperature of the formation water and mud, a, the mean activity of the connate water (assumed to be sodium chloride) and a, the mean activity of the mud (as sodium chloride). The essential accuracy of this expression was confirmed by an investigation 2 of the SP kicks opposite sands of known water salinity. The logs examined had been routine runs in oil wells in a number of Mid-Continent states. Further confirmation of the original analysis of the electrochemical SP is implicit in the paper by Tixier3 wherein a relationship between the SP and the ratio of connate water and mud resistivities, based only on field data, is given. Allowing for the fact that mud and connate water resistivities were used instead of activities, and that electrokinetic effects were ignored, the constant reported by Tixier is essentially 2tRT/F. The error involved in the universal substitution of resistivities for activities will receive further reference below. Possible reasons for the fact that the electrochemical properties of shales do not appear to be significantly dependent upon either the nature or the quantity of cation exchange materials, such as clays, they contain has been considered at some length by Wyllie and Patnode.4 In the same paper, the probable effects of shale thickness and compaction on electrochemical efficiency, and hence the inherent difficulty of obtaining reliable laboratory-scale data, were also pointed out. However, it is

Jan 1, 1951

By A. E. Jenkins, O. C. Roberts, D. G. C. Robertson

Using an inverted-cone coil at 450 kHz, it has been possible to levitate iron (FeS), cobalt (CoS), and nickel (NiS) sulfides. Important nontransition metal sulfides such as ZnS, PbS, and Cu2S have proven impossible to levitate although Cu-Fe-S ternary alloys containing 30 wt pct S and up to 10 wt pct Cu, and Cu-Co-S and Cu-Ni-S ternary alloys containing 30 wt pct Cu have been levitated. The levitation technique has been used in preliminary experiments on the vaporization from liquid sulfides and the reaction of liquid metal-sulfur alloys with oxidizing atmospheres. The course of the reactions with pure oxygen were followed using highspeed photography and two-color pyrometry. ELECTROMAGNETIC levitation is now established as a basic laboratory technique in high-temperature research but its application has been restricted mainly to metals and alloys. Applications have included alloy preparation,' metal purification,2'3 determination of liquid metal densities and emissivities,4,5 and studies of metal supercooling,4 alloy thermodynamics,6 and vaporization phenomena.7-9 The application of the technique to compounds has not been considered previously. The successful investigation of the reactions between dilute iron alloys and oxidizing atmospheres10'1 has prompted the current physico-chemical studies involving levitated metal sulfide drops and flowing inert or oxidizing atmospheres. This paper presents the results of such a study and provides a basis for future studies involving a wide range of other compounds of metallurgical interest. The successful levitation of many metal sulfides and mattes provides a method of studying the oxidation reactions fundamental to flash-smelting and similar pyrometallurgi-cal operations under closely controlled laboratory conditions. In addition the system allows the use of a controlled atmosphere (e.g., a gas stream of a certain H2/H2S ratio) with a particular chemical potential to study the relevant thermodynamic equilibria or the mass transfer processes between the atmosphere and the levitated drop under conditions where the hydrodynamics of the system can be closely defined. The optimum frequency for the levitation melting of metals in an inverted-cone coil type inductor is within the radio frequency range 400 to 500 kHz. At frequencies lower than 10 kHz the rate of heat generation is usually insufficient to melt the levitated charge' or where melting is achieved, "dripping" from the charge is encountered.'' At frequencies above 2 mHz the levitation force decreases. Metals, alloys and preheated elemental semiconductors such as germanium and silicon, have been levitated but the levitation of only a few metal compounds has been reported. Jostsons13 and the authors have levitated liquid titanium-oxygen alloys containing 50 at. pct 0 while clark14 has reported the levitation of mixtures of FeS and MnS for short periods. With a "cold crucible" inductor sterling15 has melted ferrites by preheating them by induction in a 4 mHz field and melting at a lower frequency. However this second type of inductor has been designed purely for the melting of materials without contamination; there is only a small gas film between the charge and the inductor and the electromagnetic levitation effect is of secondary importance. For this reason further discussion will be restricted to the use of the coil type inductor. The assessment of the suitability of a particular metal compound for levitation is based upon the following two criteria: i) thermal stability, and ii) physical "levitability". In this paper these two criteria will be considered separately. The thermal stability of a solid or liquid metal compound with respect to a gaseous environment depends upon its chemical reactivity with that environment or, in the case of an inert atmosphere considered here, its volatility. The physical criterion as to whether or not a particular compound can be levitated is based upon a comparison between those physical properties of the compound determining "levitability" which are defined by the fundamental equations of levitation theory as developed by Okress et a1.,16 and the properties of the metals. Since it is not practical to cover the vast field of metal compounds, further discussion will concentrate on the metal sulfides but the treatment would be applicable to any metal compound. THE THERMAL STABILITY OF METAL SULFIDES The temperatures usually encountered during levitation in inert atmospheres cover the range 1400" to 2000°C. The stabilities of the condensed states of the sulfides under these conditions are considered in relation to the periodic classification by reference to Table I. Two general classes of sulfides emerge. The solid sulfides of elements of group IIB and of groups further to the right are volatile while those sulfides of group IB and of groups further to the left are nonvolatile solids. The sulfides described as volatile may be dismissed as unsuitable for levitation. The stabilities of the more favorable nonvolatile sulfides under the anticipated conditions must be studied more closely From Table I it is seen that the alkali metal sulfides exist as liquids in the temperature range of in-

Jan 1, 1970

By H. A. Deans

A three-parameter mathematical model for one-dimensional flow in porous media is developed. The objective of the model is to predict accurately the longitudinal dispersion associated with the flow of either gases or liquids over a wide range of Reynolds number. A qualitative analysis of the model parameters is given. Published experimental pulse-response curves are compared with results predicted by the model. Several general types of problems are given for which the model can be used. INTRODUCTION The problem of multicomponent single-phase flow through porous media is encountered in the study of petroleum reservoirs, gas chromatographic and ion-exchange columns, industrial fixed-bed con-tacters and elsewhere. These particular examples span the wide range of flow conditions possible; Reynolds numbers of less than l0-4 are not unusual in oil-production problems, while values in excess of 103 are common in large fixed-bed operations. The study of flow-dependent transport phenomena is complicated both by the changes in the character of the flow over this range, and by the irregularity of the flow boundaries inherent to porous media. Dispersion is one of the important phenomena known to depend fundamentally on flow conditions as well as on fluid and medium properties. As used in this paper, the term "dispersion" refers to the observed mixing of fluid elements of different composition which occurs in flow systems. The actual mechanism may be one or more of a number listed below. Only dispersion in the direction of the mean flow (referred to as axial or longitudinal dispersion, or mixing) is considered here, although lateral dispersion arises as part of certain coupled mechanisms. MECHANISMS FOR LONGITUDINAL DISPERSION A number of distinct mechanisms are known3,5,7,8 to contribute to the phenomenon of longitudinal dispersion. The more important of these are as follows. 1. Molecular diffusion in the flow direction. 2. Turbulent (cell) mixing. 3. Lateral transport processes coupled with velocity and/or residence time distributions, including: (a) "Taylor" diffusion caused by the interaction of velocity profiles in individual voids with lateral molecular diffusion; (b) separation and remixing or interdiffusion of streams having different velocities around particles; and (c) the couplipg of gross velocity profiles, caused by viscous instability or inhomogeneous porosity, with lateral dispersion. 4. Finite mass-transfer rate between a porous matrix and the flowing phase, and finite diffusion rate inside elements of the porous matrix. It is a well documented fact that these mechanisms, whether acting individually or in combinations, produce essentially similar integral effects. The observed dispersion can be described approximately by solutions of the diffusion equation with a properly chosen dispersion coefficient. Necessary conditions for equivalence between dispersion and diffusion are discussed by Klinkenberg and Sjenitzer.7 Apparently these conditions are not satisfied exactly in many real systems. "Best fit" solutions of the diffusion equation often depart appreciably from measured breakthrough curves for the usual pulse- or step-forced experimental systems. In particular, pulse response curves usually show some degree of asymmetry and "tailing" which cannot be reproduced by the one-parameter diffusion model. The three-parameter model proposed herein attempts to produce these effects in a realistic manner. Asymmetry and tailing are predicted for certain sets of parameter values, while for other sets the new model agrees with the diffusion model. VELOCITY DEPENDENCE OF DISPERSION Molecular diffusion apparently controls at sufficiently low Reynolds number in both gas and liquid systems.*l3 At Re > 20, turbulent cell mixing is dominant for gases flowing at low pressure; the data of McHenry and Wilhelm8 show that the mixing cell is approximately 1 particle diameter long in random sphere packs. on the other hand, the data

By O. R. Morris, G. L. Hawkes

A galvanic cell, using as electrolyte a fused salt solution of calcium carbide and as electrodes carbon and a Fe-C alloy of known composition, has been set up to study the thermodynamics of Fe-C alloys in the temperature rmzge 800" to 1000°C. Time independence and reproducibility of the cell electromotive force were taken as evidence of the reversible behavior of the cell. Carbon was believed to be present in the electrolyte as the so-called acetylide ion, C;-. The plots of the cell electromotive force us temperature for a specific alloy composition were straight lines within the limits of experimental error. The average Partial molar enthalpy of carbon in iron relative to pure carbon was found to be +10,610 i 93 cal per g-atom C. Thermodynamic analysis of the data has led to the following equation for the carbon activity, ac, based upon pure carbon as the standard state: In ac = In Zc + 10,560/RT + (10.02 + 77O/T)ZC - 2.350 where ZC is the lattice ratio [nC/(nF, - nc )] and T is the absolute temperature. This equalion gives carbon activity values generally slightly lower than those from gas equilibration studies reported in the literature. METAL LOGRAPHIC examination of a polished cross section of the steel anode used in the electrolysis studies of fused salt solutions of calcium carbide by Morris and Harry revealed extensive carburization of the steel by the electrodeposited carbon. This carburization was reflected in the variability, with time, of the applied potential to the electrolysis cell, necessary to maintain a constant current density at the electrodes. This observation suggested the setting up of a galvanic cell of the "alloy concentration" type to study the thermodynamics of some metal-carbon alloys. Cells of this general type have been widely used for the study of alloy systems.2 In view of the availability of published data in respect of the austenite phase of the Fe-C system, it was decided to carry out measurements upon these alloys before proceeding to studies of less well documented systems. The galvanic cell may be written: where [C] is carbon dissolved in iron. The electrolyte was a fused salt solution of calcium carbide, containing some 5 to 10 mol pct of carbide. The cell reaction is believed to be: C(s)-[CI [I1 Carbon forms an interstitial solid solution in iron, with the atoms located in the octahedral interstices. In the fcc crystal structure of austenite there is one octahedral interstice per iron atom. Thus, the lattice ratio, ZC, shown by Gurney3 to be the fundamental concentration parameter in the context of interstitial solutions, is given by: where nc and nFe are the number of carbon and iron atoms, respectively. chipman4 has recently shown empirically the advantages of using this concentration parameter instead of the more usual atom ratio or atom fraction. The cell electromotive force, E, assuming reversible behavior, is related to the carbon potential or the partial molar free energy of carbon in the solid solution relative to pure carbon at the same temperature and pressure, GP at the composition ZC, by the equation: where z is the carbide ion valency and F is the Faraday constant. An activity of carbon, ac, in the solution relative to the value of unity assigned to pure carbon, and an activity coefficient, qC , are defined such that: where R is the gas constant and T the absolute temperature. GF is further related to the relative partial molar enthalpy Hm, and the temperature coefficient of the cell electromotive force, (aE/aT)Zc, by the equations: Measurement of the cell electromotive force thus enables calculation of the relative partial molar thermodynamic properties of carbon in iron, if z is known. At E = 0, the solid solution is in equilibrium with pure carbon. More convenient for many purposes is the standard state based upon the infinitely dilute solution, Henry's law. The relationship between the activity coefficient of carbon based upon this standard state, and that based upon the pure carbon standard state, qC , may be obtained by considering the free energy of transfer of carbon from the latter standard state to the former. The relationship is: where +:H is the activity coefficient of carbon in the hypothetical standard state based on a reference of

Jan 1, 1969

By D. B. Robinson, C. A. Macrygeorgos, G. W. Govier

Experimental data have been obtained on the volurrletric behavior of ternary mixtures of methane, hydrogen sulfide and carbon dioxide at temperalures of 40°, 100" and 160°F up to pressures of 3,000 psia. The results indicate that the compressibility factors for this system do not agree with compressibility factors for sweet natural gases at the same pseudo-reduced conditions. The deviation increases as the temperature and methane content decrease. Discrepancies of up to 35 per cent were observed. A careful analysis has been made of the existing pUrblished data on compressibility factors for binary systems containing light hydrocnrbons and hydrogen sulfide or carbon dioxide. It has been found that the deviation of actual from predicted compressibility factors for methane-acid gas mixtures is a function of the methane content and the pseudo-critical properties,.v of the mixture. The ratio between actual compressibility factors for methane-acid gas mixtures and compressibility factors for sweet natrlral gases at the same pseudo-reduced conditions has been currelated over a range of pP,, from 0 to at least 7 arid a range of pT, from about 1.15 to at 1east 2 0 with an error not exceeding 3 per cent and over most of the range within I per cent. The validity of the correlation for mixtures containing appreciable hearvier hydrocorbons has not been fully established, but it is shown to be preferable than the use of a corretation based only on hydrocarbons. INTRODUCTION Although a relatively accurate method for predicting compressibility factors of pure materials is provided by charts based on reduced properties and the assumption that the compressibility factor is a unique function of T P and z the determination of the correct values of compressibility factors for gas mixtures is somewhat difficult. Two general methods of dealing with gaseous mixtures have been proposed. The first assumes a direct or modified additivity of certain properties of the mixture in terms of the properties of the individual components. Examples of this method are based on the familiar laws of Dalton and Amagat. The second method averages the constants of an equation of state applicable to the pure components. Both of these methods are of limited value in engineering calculations because the first usually provides reliable answers only over narrow ranges of pressure and temperature and the second is cumbersome to handle. In petroleum engineering practice accurate estimations of the volumetric behavior of natural gases arc frequently required. To fulfill this need, several generalized compressibility charts have been developed.' ' Of these, the one prepared by Standing, el al is widely used at present. In the construction of charts of this type a third method for dealing with mixtures has been followed. It is based on correlation of pseudo-critical properties as outlined by Kay and calculated from the critical properties of the individual components in a mixture. Although these charts provide relatively accurate information on the compressibility of dry or wet sweet natural gases, they are less reliable when used for gases containing high concentrations of hydrogen sulfide or carbon dioxide or both. Thus, an experimental program, although time consuming, is the best means now available for the determination of the volumetric behavior of sour or acid gas mixtures. An increased interest in the behavior of these gas mixtures, particularly in connection with some of the fields in Western Canada where the acid gas concentration of the reservoirs may be as high as 55 per cent and where hydrogen sulfide alone may be as high as 36 per cent, provided the incentive for this study. It was the purpose of the investigation to determine the volumetric behavior of selected mixtures of methane, hydrogen sulfide and carbon dioxide over a range of temperature from 40" to 160°F and at pressures up to 3,000 psi. EXPERIMENTAL METHOD The apparatus used in this investigation was basically the same as that described by Lorenzo.'" The amount of each pure component used in preparing the gas mixtures was measured over mercury in a glass-windowed pressure vessel. The pure components were then transferred individually in the desired amounts to a second glass-windowed pressure vessel where the volumetric behavior of the mixture was determined. Volume was varied by mercury injection or withdrawal. The capacity of the cell was about 125 cc. Temperatures in the cells were measured with copper-constantan thermocouples and a Leeds Northrup semi-