Search Documents

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 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 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 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 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

INTRODUCTION Trade financing is a particularly important component of short-term financing for a minerals company since errors, extra costs, or payment delays/defaults can easily wipe out the profit margin on a particular shipment or can cause problems when arranging purchases of e.g. reagents. It is required to make any minerals operation function and does not have to be cross-border to require proper management. Sources of working capital are as follows: 1. Extended credit terms from suppliers 2. Inventory loans from banks, finance companies 3. Accounts receivable financing or factoring 4. Credit lines or overdraft facilities from banks. TRADE CREDITS These are provided under the familiar terms "net 10 days" or payment due by the "10th day of the following month." Thirty days for payment is usual but longer periods can be negotiated with regular suppliers. Shipment may be made on open account with specified periods for settlement. In the mining industry, a custom smelter gains trade credit by paying for concentrates three months after the month of shipment. INVENTORY LOANS Various forms of inventory finance are provided via "floor-plan" or warehouse financing. A set percentage, usually 60%-80%, of a good-quality, readily saleable inventory item (e.g. bullion, concentrates, ingots) will be provided by a bank or finance company which will usually have security over these products and has a right to inspect the goods from time to time. Legal security will he achieved by warehouse receipts, floating liens, chattel mortgages, trust receipts, or through collateral management service companies. The latter certify control over stock or distribution outlets and often carry insurance against fraud. Costs range around 1.5% - 2.0% p.a. or more plus some fees to cover out-of-pockets or other costs. ACCOUNTS RECEIVABLES Accounts receivable financing is based on the company's generation of satisfactory receivables which are continuously assigned to a bank or finance company ("factoring"). Because of the risks in being able to collect receivables, such as payments to a smelter for its blister copper shipments to a refinery, this form of financing costs as much as 3% - 5% of the debt amount purchased plus an additional charge for processing paperwork. Receivables can be purchased with or without recourse and on a notification or non-notification basis to the client's customer to, respectively, pay the invoice directly to the bank or not. LETTERS OF CREDIT Letters of Credit ("L/C's") are a common feature of international trade finance whereby a bank, at the request of the buyer, issues a written undertaking to an exporter to pay a specified sum of money within a set time on presentation of previously specified documents. This undertaking means that the exporter is no longer relying on the credit standing and integrity of the foreign importer. A dispute about shipment does not otherwise affect the obligation to pay upon presentation of a "clean" set of documents. These documents include the hill of lading and other documentation, often on standard "form" bases. There are a number of conventions and techniques in L/C's such as irrevocable/revocable L/C's, time L/C's, confirmation/advising L/C's and transferable L/C's. These are tools to adjust for the varying risks that the importer will not pay. Other routes where one is comfortable with an importer include documentary sight hills (D/P hill), documentary term bills (D/A bill), and, of course, open account.

Jan 1, 1985

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-

By B. C. H. Steele, C. B. Alcock

In choosing solid oxide electrolytes for use in the measurement of thermodynamic quantities at high temperatures, the two most important criteria are the values of the partial ionic and electronic conduc tizities. Measurements of these conducticities have been made for some Group IVA oxides and solid solutions of these oxides with CaO, Y2O3, and La2O3. The total conductivities as a function of temperature. oxygen partial pressure, and composition were determined, as well as the electromotive forces of cells in which the oxygen partial Pressures of the electrodes were between 10-7 and 10-29 atm at 1000°C. Additional factors that influence the experimental arrangements, such as the operating temperature, the appropviate cell atmosphere, and electrode kinetics, are also discussed. IT has long been recognized that galvanic cells incorporating solid electrolytes can possess many advantages for the measurement of high-temperature thermodynamic data. Unfortunately, the electrical-transport properties have only been established for a few solid electrolytes, although attempts have been made, and are still being made, to use such materials as porcelain' and magnesia.' Kuikkola and Wagner have discussed3 and demonstrated4 various applications of solid electrolytes, and, in particular, their work stimulated the electrochemical measurement of oxygen activities using cells incorporating the solid oxide electrolyte, Zr0.85Ca0.15O1.85 The excellent characteristics of this oxide electrolyte have been confirmed by the investigations of Kingery et al.5 (oxygen ion diffusion), Carter and Rhodes6 (zirconium and calcium ion diffusion), and Weissbart and Ruka7 (electronic transference number). However the observations of Peters and Mobius,8 later confirmed by Schmalzried,9 indicated that the electronic component of zirconia-based electrolytes could become significant at low oxygen chemical potentials. Under these conditions, thoria-based electrolytes have obvious advantages arising from their superior thermodynamic stability. It was known that certain thoria solid solutions possessed defect fluorite structures containing anionic vacancies (cf. Zro.85Cao.l50O1. 85), although data on their electrical properties were conflicting. Kiukkola and wagner4 concluded that thoria-based electrolytes exhibited appreciable electronic conductivity, whereas Peters and Mobius8 mentioned that the electronic conductivity only became important at high oxygen partial pressures, but gave no further details. In addition, both Kiukkola and wagner4 and Peters and Mann'" had determined the oxygen chemical potentials associated with the Fe-FeO, Co-COO, Ni-NiO, and Cu-Cu2O equilibria, and it was not known whether discrepancies in their results could be attributed to the different cell design, or the thoria-based electrolyte used by Peters and Mann. These uncertainties have been resolved during the present investigation into the measurement of low oxygen chemical potentials, and transport properties of the relevant oxide electrolytes are presented and discussed. Other experimental difficulties associated with low oxygen chemical potential measurements are also described, with particular reference to the Nb-O system. 1. THEORETICAL CONSIDERATIONS Oxide systems of potential use as solid electrolytes possess a large energy gap between the valence and conduction bands. Little is known about the defects present in these wide-gap metal oxides as a function of impurity content and deviations from stoichiom-etry, or the mechanisms by which electronic charges move in these solids. For these materials it is convenient to represent the measured conductivity as the the sum of three terms: where a, is the conductivity due to positive hole conduction (p type), and oe is the electronic conductivity due to electron conduction (n type). U there is a large concentration of ionic vacancies fixed by composition then the ionic contribution will not be a function of oxygen pressure. Although the possibility of a pressure-independent electronic component also cannot be ignored, it is likely that such a term would only become significant at very high temperatures, and in general the positive hole conduction will increase with increasing oxygen pressure, and the electron conduction will increase with decreasing oxygen pressure. It follows that the total conductivity is related to oxygen partial pressure as follows:

Jan 1, 1965

By C. J. Ramsburg

THE construction and operation of by-product coke-oven plants in America are essential to strong national defense and of the greatest importance to many widely diversified undertakings as well as to safety and health. CARBONIZATION PROCESS Coal used under a boiler or in a locomotive is consumed by combustion. Coal used in a by-product oven is distilled or carbonized, meaning that heat is applied to the coal, out of contact with air, causing a distillation that removes the volatile matter in the coal in the form of gas and vapor, leaving coke as a final product. From the gas and vapor distilled from the coal, many valuable substances can be secured, known as by-products. Hence the name "by-product oven." The simplest example of this principle can be shown by means of an ordinary clay pipe. If powdered coal is placed in the bowl of the pipe, the opening of the bowl is plugged with clay, and the bowl placed in a flame, the bowl will become hot, transferring the heat through the wall of the bowl to the coal within. When this coal attains a temperature of about 950°F., gas will begin to pour out of the pipe stem, and this gas may be lighted. A black liquid also will exude from the stem of the pipe. This is tar mixed with certain other products. What remains in the bowl after distillation is coke. A coke oven corresponds to the bowl of the pipe. A modern oven of the latest type is a refractory brick chamber approximately 42 ft. long and 14 ft. high, averaging 18 in. in width. The oven is filled with coal through closable charging holes in the oven top. Each end of the oven consists entirely of a self-sealing door, the removal of which permits the finished coke to be discharged by means of a power-operated pusher machine. The ovens are built in batteries, and the heating is effected by gas burned in vertical flues in the oven walls, each set of flues heating the walls of two adjacent ovens. Such ovens carbonize a charge of approximately 18 tons of coal in 18 hr., or have a daily throughput of 24 tons per oven, producing about I 7 tons of coke. In times of need, the output may be increased by raising the temperature of flues, thus shortening the coking time. There is one great essential difference in principle between distillation in a clay pipe and in a coke oven. In the 18-in. oven, the distillation proceeds from two parallel walls held at a high temperature (over 2000°F.) and instead of all the coal being heated at once, the distillation is a progressive one creeping inward. The average travel of this heating is 1/2 in. per hour from each vertical wall, and as the volatile products are formed progressively they must push away from the interior toward the wall through the previously formed hot coke, then up the walls and over the top of the coal charge, finally escaping through the outlet pipe. This results in a cracking of the components into many different and valuable substances not formed when the carbonizing temperature is low.

Jan 1, 1942

By J. W. Dougherty

THE Hollinger gold mine is in the Porcupine mining division, District of Cochrane, in the northern part of the Province of Ontario. The adjacent town of Timmins, with a population of 25,119, lies 490 miles by railway north of the city of Toronto. GEOLOGY The property is situated on the northwest limb of a major syncline, which pitches northeast at approximately 50°. The principal rocks of the Hollinger area are of pre-Cambrian age and consist of altered, schistose, Keewatin lava flows intruded by stocks of schistose, quartz-feldspar porphyry generally assigned to the Algoman period. The structure of the Keewatin series is extremely complex as a result of close folding and some faulting. The schistosity, trending north 70° east and dipping steeply to the southeast, is most highly developed along a subsidiary anticline in which the major porphyry bodies occur. The ore bodies, which exhibit three principal structures, consist of quartz veins and associated pyritized wall rock. Some lodes are com-posed of en echelon veins, others exhibit a main quartz body with branch-ing stringers and a third consists of a single definite lead. Gold, in a coarse visible form, may be confined to quartz in some ore bodies but in others it is finely divided and is found not only in quartz but also asso-ciated with pyrite in the adjacent wall rock. In many veins the dis-tribution of gold is such that there are several ore sections, separated by areas of low-grade material, in the same vein although its physical appearance remains unchanged.

Jan 1, 1940

LIBERIAN ore has been called the richest iron ore mined in the world, and thus far, the output from the Liberian Mining Co.'s Bomi Hills mine has lived up to its reputation. Iron content has been variously given as 68.47 pct to 68.825 pct. A recent United Nations survey gave it first place, surplanting the previously top-rated Swedish product. The first shipment of Liberian ore to the United States reached Baltimore in June, 1952 and went into blast furnaces at Republic Steel's Cleveland plant, where excellent results were reported. The physical make-up of the ore was hard, dense lumps, free of fines and with ideal size distribution. The Bomi hills mines were scheduled to be in full operation by the end of 1952. At the close of 1951, 200,000 tons of ore had been delivered to Baltimore. Liberian Mining Co. states that production is about 90,000 tons of ore per month at present. Between 1200 and 1500 men are employed at the mine site, with about 25 pct of them Americans. The Bomi Hill project was financed with the aid of the Export Import bank, which loaned Liberia Mining part of the sum needed to initiate and put the development into production. No definite expansion plans are in the making but extensive diamond drilling is being pursued. Original estimates of annual production for Bomi Hills ran in the neighborhood of 1.5 million tons, but the figure is no longer applicable. Liberia's potential as an iron producing country has been known for a long time, but the physical obstacles in the way of large scale production were such as to prohibit even contemplation of the project. Until recent years, the only route from the mine site to Monrovia was a foot trail. Once to the coast, over swamps and unbridged rivers, the ore was still useless because of the lack of loading facilities. Transportation Solved One problem was solved with the completion of the artificial harbor at Monrovia. Harbor construction was well under way when the mine project began in 1948. The island location of Monrovia also made the building of a bridge across the St. Paul River necessary. The 980 ft reinforced concrete structure, accommodating railroad and vehicular traffic, was opened in 1949. Meanwhile, preliminary work on railroad construction started. Grading was completed quickly. In April, 1951, the first train made the trip from Monrovia to the Bomi Hills. About half of the cost of the entire project went into railroad construction. More than 2200 men, often forced to resort to primitive methods, were engaged in clearing a path through the heavy forests between the mine and port. Swamps were drained and filled, and bridges, including the William V. S. Tubman, named for the Liberian President, were constructed. At the port, a Robbins Conveyor Belt Co., orehandling plant was installed. Ships can be loaded from the stockpile at the rate of 3000 tons per hr and vessels need spend only one day in port. Liberian Mining is operating under an 80 year lease, giving the company exclusive mining rights to all minerals, except gold, platinum, and diamonds, within a 40 mile radius of Bomi Hills. The contract provides a fixed royalty to the Liberian Government of five cents per ton of ore extracted, and a premium -based on the New York market price for Bessemer-grade pig iron. The entire project started because an American,

Jan 1, 1952

By F. E. Streigert, J. H. Mallery, A. H. Geisler

The existing data on precipitation in Cu-Be alloys have been reconciled with some new data to present a more complete and rational analysis of the atomic rearrangements involved in the process. Two different coherent precipitates were proposed. Each is responsible for a type of diffuse X-ray diffraction effect. THE nature of the atomic rearrangement involved in the precipitation of the compound y-Cu-Be from supersaturated solid solutions of the copper-base alloys has been the subject of two extensive investigations.'-" he precipitate in this alloy system has a rather simple structure; at equilibrium it has a body-centered cubic unit cell with copper atoms on the corner sites and beryllium atoms on the center site. While a simple mechanism would be expected to obtain for the precipitation process, those which have been proposed do not adequately explain all of the reported data. The purpose of the present report is to rationalize some of the previous observations and to contribute new data on the subject obtained by the usual techniques of structure analyses. The Cu-Be system of alloys represents one of the most extensively investigated precipitation systems, so a detailed review of the pertinent literature on structure is in order. Striations: Since one of the important features of the mechanism—the habit plane—is determined by microscopic examination, microstructural features that have been observed should be considered at this time in order to assess their role in the precipitation process in Cu-Be alloys. One of the most prominent features appears as regularly oriented markings or stripings (ripples or striations) across the crystals on a suitably polished and etched sample of the aged alloy. These have been observed and interpreted by many workers. Masing and Dah14 observed striations in both the a and /3 phases of a 2.5 pct Be alloy, and attributed them to the beginning of the precipitation of very fine y crystals along favored planes of the a matrix and to the local decomposition of the £ crystals into a plus y. These markings were apparent in both the a and /3 phases of the micrographs given by Smith5 for a 2.75 pct Be alloy, while Moore" detected them in the a matrix phase using polarized light. Desch7 attributed the markings to fine precipitate particles lined up on a family of matrix planes like Masing and Dahl had done, and he identified the planes as those of the (111) family, the slip planes of the matrix. The striations were presumably the structure which Northcott- eferred to as "martensitic." He found that this structure was present in samples that had been aged under widely different conditions with hardness values ranging from 160 to 400 VHN. Later Guinier and Jacquet' also found striations in both the a matrix and secondary /3 phase. Their analyses showed that the striations were along (110) planes of the matrix phase, apparently in contradiction of the results of Desch. They found that the striations were developed during aging on the surface of a specimen that had been polished prior to aging.' On this basis they concluded that the striations represented mechanical deformation by slip along (110) planes of the matrix resulting from precipitation stresses.'" Similar but fainter striations were observed in the same manner on an A1-Cu alloy sample, in which the change in volume was much less than for Cu-Be alloys and

Jan 1, 1953

By Thomas Fraser, David R. Mitchell

THE primary purpose of the loading plant is to transfer the finished product from the preparation machines to the railroad car, truck, or barge in which it is to go to market. Secondary purposes of the loading plant are: to prevent excessive degradation of sizes during loading, to reduce segregation as much as possible, and to prevent excessive quantities of very good or off-grade coal from entering any one car. Nonuniformity of consignments is one of the chief sources of complaints by consumers and coal merchants. A properly designed loading plant will do much to reduce complaints and insure uniformity, not only from shipment to shipment but also in each individual consignment along the length of the car and in cross section. DEGRADATION IN STORAGE, LOADING AND SHIPMENT The object in loading is not only to handle the coal without breakage and segregation, but also to produce a well-trimmed load that will ride to its destination without spillage o~ breakage en route. The small sizes present in a car-of prepared coal as received by the purchaser consist of: missed small sizes not screened out in the preparation plant, breakage during loading, and breakage caused by handling en route. For domestic sizes of anthracite shipped from the Pennsylvania anthracite field to Atlantic coast cities, the dealers' loss because of breakage in loading and shipment amounts to from 2 to 5 per cent and for anthracite shipped to Middle West cities, the loss is 8 to 9 per cent. This difference is in part because dealers in the Middle West rescreen their coal more thoroughly than is usual in the East. One dealer handling Illinois No. 6 coal shipped from southern Illinois mines to Chicago reports an average size loss of 10 per cent when lump coal is unloaded with a fork. This runs as low as 5 per cent for carefully loaded coal and as high as 15 per cent for carelessly loaded coal. An Ohio dealer estimates the degradation in loading and shipment of West Virginia coals consigned to Ohio cities as follows: Pocahontas lump

Jan 1, 1950

By Nicholas J. Grant, Ulf Kalling, John Chipman

THE operation of a blast furnace is dependent to an important extent upon the sulphur content of materials charged and the desired limit of sulphur in the product. It has long been known that the blast furnace is the most efficient tool for desulphurization in common use and that this efficiency is associated with the strongly reducing conditions of the hearth and is enhanced by increased basicity and fluidity of the slag. The chemical reactions of desulphurization may be studied from the viewpoint of the ratio of the process or of the final equilibrium conditions. Both kinds of studies contribute to an understanding of the process and both are included here. A simple measure of the desulphurization power of a slag is given by the ratio: Pct sulphur in slag (Pet S) Pct sulphur in metal [Pct S] This ratio was used by Holbrook and Joseph',' to measure relative desulphurizing powers under controlled laboratory conditions. It was also used by Hatch and Chipman3 as a measure of the equilibrium distribution. For the latter purpose it would be preferable to employ thermodynamic activities rather than percentages, but until very recently this has been impossible for lack of data. Now, thanks to the work of Morris and Williams and Morris and Buehl," the effects of carbon and silicon upon the activity of sulphur in the metal are known. The confirmation of this work and its extension to include the effects of other elements by Sherman and Chipman and by Rosenqvist and Cox' make it possible to calculate the activity of sulphur in pig iron of any composition. Hence it is now possible to use data on the equilibrium distribution of sulphur to find its activity in the liquid slag and to approach an ultimate solution of the thermodynamic aspects of the problem. The rate of transfer of sulphur from metal to slag is the problem of major industrial importance and indeed the principal need for equilibrium data has been as a necessary adjunct to the kinetic studies. The rate of approach to equilibrium under laboratory conditions seems slow compared to the requirements of industrial practice, and it is clear that further laboratory studies of rates are needed. In the research reported below, the items which were investigated were the following: I—The role of mechanical stirring on the approach to equilibrium. 2—The role of MgO in desulphurization as compared to CaO. 3—The role of MnO in desulphurization. 4— The limiting reactions which constitute the slow steps in desulphurization. Experimental Procedure The experimental set-up and procedure previously described by Hatch and Chipman" were essentially followed with several small modifications. The graphite crucible containing the slag and metal charge was altered to provide considerably more active stirring and mixing of the slag and metal in the carbon monoxide atmosphere. For this purpose the crucible was machined to provide two deep cylindrical wells which were interconnected at top and bottom as shown in Fig. 1. A graphite screw with a flat thread and of shallow pitch (4 threads per in.) spinning at 600 to 800 rpm was used to lift the slag and metal over the partition between the two wells and throw them over into the second well, where the metal settled through the slag into the reservoir at the bottom. It was possible to see actual particles of slag and metal being thrown over the partition. In this respect, the stirring was more vigorous than used in the work of Hatch and Chipman. A charge of 400 g of wash metal was first melted, and 20 g of FeS was then added to yield a bath containing 1.65 pct S. Immediately 400 g of slag (as pure mixed oxides) was added and fused. The slag was generally fused in 1 hr * 10 min. Within 30 to 45 min after melting, the temperature was adjusted to 1525"C, and the first slag and metal samples were taken. The slag was picked up on the end of a cold Armco iron rod, whereas the metal was sucked into a silica tube. The wash metal composition was (in percent): 4.29 C; 0.022 S; 0.021 P; 0.38 Si. The slags used were of four fixed starting compositions covering a wide range of acid-base ratios shown in Table I. Deliberate variations in MgO were made in these slags to check the role of MgO in blast-furnace desulphurization. Changes due to additions and reactions were followed by analysis of samples. Additions of Mn and MnO were made to most of the heats to note the role of Mn and MnO on desulphurization. Three heats (62 through 64) were made in an open pot induction crucible (graphite) using a

Jan 1, 1952

Report of the Secretary GENTLEMEN Herewith are transmitted reports from the Treasurer and of the principal standing committees of the Institute. To these special reports members are referred for details as to the year's work. Despite the unfavorable business conditions obtaining throughout the year, the Institute was able to carry on its various activities without curtailment and without drawing on the surplus accumulated through previous years. In fact, small additions were made to the Life Membership and various other funds. It is highly desirable that such additions should continue to be made, especially to the Life Membership fund, for reasons reported last year. At present, money must be used from current receipts in order to enable the Institute to fulfill its obligations to those who on joining the organization commuted their dues by a single payment. In addition, with the passage of years the number increases of those who have paid dues forty years or more and who, by reason of age or disability, are unable longer to practise their profession. These men have supported the Institute through their active professional life and it would be a pity to lose them as members through resignation or for failure to keep up their dues. In accordance with the By-laws the Board has authority to advance such members, on application, to the Honor Roll of Senior Members and to relieve them from further payment of dues. It has been the policy of the Board to do this liberally, but to keep up this, as well as to meet other obligations, it is necessary for money to be provided. It has been estimated by the Treasurer that nearly $50,000 additional should be placed in the Life Membership fund as soon as possible to meet these obligations without using current income and to that end it is necessary, not only to avoid running behind currently but to earn, if possible, some surplus each year. The small amount transferred to the Life Membership fund from surplus this year was, as a matter of fact, made possible by certain members paying personally bills that were properly chargeable to the Treasury. Such generosity is common among those active in Institute work and probably no year would have closed with a surplus except for it.

Jan 1, 1929

By John A. Putnam, Morrough P. O’Brien

DATA on the physical and thermodynamic properties of hydrocarbons have been made available in recent years but the formal method of applying these data to flow in porous media appears not to have been fully developed. By making usual thermodynamic assumptions and employing existing data on permeability to two-phase flow, Darcy's law may be applied to obtain the flow characteristics for representative combinations of flowing medium and geometry. The numerical examples presented herein are based on hydrocarbons for which complete physical and thermodynamic data are available under equilibrium conditions. INTRODUCTION Thermodynamic data are obtained in the laboratory under equilibrium conditions and their application requires the assumption that equilibrium also exists during flow. This assumption is at least usual in engineering analyses, though not easily justified in many instances, and will be followed here. In regard to hydrocarbons, laboratory experience in determining the pressure-volume-temperature relationships indicates that gas-liquid mixtures must often be stirred vigorously for periods of 15 to 20 min. before reaching equilibrium. On the other hand, boiling and condensation are known to occur on solid surfaces and the large surface exposed to the flowing substance by a porous medium may bring about rapid approach to equilibrium. For the purpose of providing at least a guide to experimentation and a standard for later comparison with experimental results, it will be assumed that the phases present in cross sections perpendicular to the direction of flow are in thermodynamic equilibrium. Another assumption usually made in the solution of problems of flow in porous media is that the fluid remains at a constant temperature; i.e., isothermal flow. This assumption is arrived at by considering the porous medium and surroundings to be an infinite heat source and the rate of heat flow to be infinite at a boiling front. Its validity depends in large measure upon the character of flow, the number of components, and the physical and thermodynamic properties of the fluid. However, the assumption of isothermal flow represents a limiting case. The other limit can be represented by the assumption that the fluid is thermally isolated; i.e., irreversible adiabatic. The thermodynamic process characterizing the actual flow of a given fluid system probably lies somewhere between these two limiting cases and might be thought of as a flow system in which the confining walls are maintained at a constant temperature. For isothermal flow the solid material of the porous media is considered as the infinite heat source capable, without loss in temperature, of giving up heat to the fluid as the latter tends to cool by expansion. For the thermally isolated system, the solid material is considered to give up no heat to the fluid. In the intermediate case the impervious walls confining the porous media are at constant temperature from which heat is transferred along various paths, solid and fluid, throughout the porous media. It would appear that an actual petroleum-

Jan 1, 1941