Search Documents

By J. Phillip Neal, Kefton H. Teague, Carroll P. Rogers

Feldspars, the most abundant minerals of the igneous rocks, occur in numerous forms and mixtures. The feldspars of commercial significance are found in widely distributed pegmatites as large crystals usually free of iron-bearing impurities and thus suitable for hand-cobbing, and also in larger bodies where the ore contains various types of feldspar intermingled with quartz and relatively free of iron-bearing impurities, or at least readily unlocked from impurities present at relatively coarse mesh. The latter type of ore includes some pegmatites, and the "alaskite" bodies so common in the Spruce Pine district of North Carolina. On the west coast of the US and in Oklahoma there is-or has been-commercial recovery of natural feldspathic sand as a source of alumina in glass and ceramics. "Aplite" in commercial terminology is a light-colored igneous rock with a granitic composition and a fine sugary texture, often banded with readily removable iron impurities. Only one area in Virginia produces aplite for the glass industry; it is located near Montpelier in Hanover County. According to Deer et al. (1966) the word feldspar derives from the Swedish feldt or jolt (field) and "spath:" chunks of rock appearing in tilled ground overlying granite. Castle and Gillson (1960) indicate a Germanic origin, citing "spat," which is said to refer to any transparent or translucent material which is readily cleavable. The term "spar" has in the past been applied to minerals other than feldspar, such as barite, calcite, and fluorite. It is correct to call the latter "fluorspar." There is a true barium feldspar, but it is rare and of no present economic importance. What is recognized as feldspar or "spar" by present producers and consumers consists of three silicate minerals which, if pure, would have the formulas KA1Si3O8 (microcline or orthoclase), NaAlSi3O8 (albite), and CaAl2Si2O8 (anorthite). These are almost never found in pure form in nature, but occur together in great abundance in a three-component system. A considerable number of different combined ratios of these three exists, with the exception that isomorphism between potassium spar and calcium spar is very limited. User specifications for feldspar products are largely built around the prevalent naturally occurring ratios. Passages following will refer to "K-spar," "Naspar," or "Ca-spar." Nomenclature and Descriptive Terms Relating to both the mineralogy and the economics of feldspar are certain terms whose coverage it is useful to specify. Commercially, aplite is a feldspathic rock mined and beneficiated in one location in Virginia, in which both titanium and feldspar minerals are present-the latter now being the major, if not the only, economic product. Technically, the term aplite when applied to this rock is of questionable accuracy. But, in light of accepted commercial usage, the term will be here used to designate it. Regarding the Virginia rock, Ross (1941) is cited. Alaskite is another term requiring understanding if not precise use. A distinctive rock type customarily called by this name is found near Spruce Pine, Mitchell County, North Carolina. An important major feldspar source, it contains somewhat higher levels of plagioclase feldspar and quartz than alaskite as defined in the Glossary of Geology and Related Sciences, which states that `orthoclase, microcline, and subordinate quartz are the principal

Jan 1, 1983

By C. Land, R. Hultgren

Heat contents of four Cr-Ni alloys were determined in a diphenyl ethev calorimeter. Aside from ferromagnetic effects, Kopp's law of additivity of heat capacities is approximately fol-lowed, in disagreement with earlier results of A. W. Foster. In 1934 Foster1 reported heat-capacity measurements of dilute solutions of chromium in nickel, primarily as a study of the effect of alloying on the ferromagnetic Curie temperature. He found that, apart from the Curie temperature effect, the general level of heat capacities was reduced 10 pct by the addition of 1 pct Cr; 20 pct by the addition of 2 pct Cr. His measurements of pure nickel essentially agreed with accepted values. These results indicate that small amounts of dissolved chromium cause an increase in the Debye temperature of nickel, presupposing a change of frequency of interatomic bonding vibrations which is much larger than should be expected. If Foster's results were verified, therefore, they would have important implications for bonding theory. The only other heat-capacity measurements in this system seem to be those of Taylor and Hinton,' who studied anomalies which they attributed to super-lattice formation at 25 at. pct Cr. At temperatures above the anomaly, Taylor and Hinton found heat capacities at this composition were quite normal, approximately obeying Kopp's law of additivity. EXPERIMENTAL WORK In order to verify or disprove the remarkable results of Foster, the heat contents of four dilute solid solutions of chromiunn in nickel were measured in a diphenyl ether calorimeter. This instrument, described elsewhere,, operates on the principle of the Bunsen ice calorimeter, but uses diphenyl ether as the working substance. Diphenyl ether melts at the convenient temperature of 26.S°C and supplies a more sensitive measure of heat than ice. The hot sample, in an atmosphere of argon, is dropped from a furnace at a measured temperature into a chamber which is surrounded by a mixture of solid and liquid diphenyl ether. The heat from the hot specimen melts some of the solid without chang- ing the temperature of the calorimeter. The amount of solid melted, and therefore the number of calories released by the specimen, is measured from the increase of volume of the mixture, which pushes mercury down a graduated capillary tube. The apparatus has been standardized against the known heat content of pure platinum. Heat loss during the dropping of the specimen was estimated by comparing the heat content of an empty platinum capsule with a solid platinum ball of the same radius and degree of polish of the surface. Repeated tests demonstrate the precision of measurement to be closer than *0.3 pct. Smoothed experimental results are listed in Table I. The tabulated values for pure nickel are taken from kelley, several measurements were made which agreed closely with this tabulation. For chromium, however, no similar agreement was found. The values tabulated were taken from a thorough study of chromium with the diphenyl ether calorimeter which will be published elsewhere.= Suitable plots of the data showed the expected ferromagnetic anomalies in the low-chromium alloys, while the 11 pct alloy shows no ferromagnetic anomaly, but a different, small anomaly below 900°K. This is perhaps the same anomaly found in CrNi, which was once attributed to superlattice formation,' but is probably due to some other cause, such as antifer-

Jan 1, 1960

By J. S. Lasater

Resu1ts of experinmental measurernents of heat capacities and thermal conductivities of some typical porous rocks are presented. Measured heat capacities agree closely with va1ues calculated front known chenlical . compositions of the rocks. On the basis of this agreement, heat capacitir.7 of fluid-satttrated rocks were calc~ililted. Therrnal conductivities of the rock samples were measured under varioits conditions of fluid saturation. From these data therrnal dif-fvsivities were calculated. A significant variation of thermal diflusivity with trnlperature is indicated. Knowledge of the thermal characteristics of fluid-bearing porous rocks has become increasingly important with the development of thermal oil-recovery processes. Reliable thermal data for petroleum reservoir rocks are not available and approximate values are generally assumed for reservoir calculations. 'The effects of temperature, pressure and fluid saturation on thermal properties of porous rocks have not been carefully studied. The purpose of the present work was to measure and report the thermal characteristics of some typical porous rocks. A group of eight sedimentary rock samples, including limestone, shale, sandstone and silt-stone, was selected for the study. Properties reported include heat capacities, thermal conductivities and thermal diffusivities. The effect of fluid saturation on these properties was also considered. Determinations of thermal properties require precise measurements which are time consuming and difficult to reduce to a routine basis. Methods have not been developed for measurement of thermal proper-ties under conditions of pressure, temperature, and fluid saturation, which may be encountered in thermal recovery processes. As one approach to these problems, some preliminary work on methods of estimating thermal properties from other more readily measurable characteristics of the rock-fluid system is reported. The samples selected for the studies were oilfield cores from within, or closely associated with oil-productive zones. Descriptions of individual samples are given in Table 1. Clay minerals present in the samples were determined by differential thermal analysis (DTA). Quartz content was estimated from DTA cooling curves. For all but the thermal conductivity tests, a representative 200-gm sample of each core was disaggregated, thoroughly extracted in diox-ane and oven dried at 100°C. The cleaned and dried sample was care- fully split into appropriate portions for the several tests. Chemical analyses of the samples were obtained to permit calculations of heat capacities by methods which are outlined later. Resuits of chenlical analyses are given in Table 2. Heat contents of the test samples were measured by a Runsen-type calorimeter. The heat content of the unknown sample is measured relative to the known heat content of platinum. The apparatus is referred to as a relative error type because of the comparative measurements, and has the advantage of canceling systematic errors. Experimental results are reproducible within 0.5 per cent. A temperature base of 298.16 K is taken for heat content measurements. Measurements were made at five temperatures (260, 440, 620, 800 and 980° F). It was necessary to preheal the sample above the maximum experimental temperature to eliminate volatile constituents. Assuming the weight loss on pre-heat, ing was mostly water, a correction to the observed heat content was made to account for the water originally present in the sample. This latter procedure was followed in the determination of calculated heat contents as will be described later. The apparatus and the method of measure-

By J. M. Smith, D. Kunii

Effective thermal conductivities of sandstones filled with stagnant fluids were measured using a steady-state technique. Data were obtained for seven sandstone samples, taken from four different locations and ranging in permeability from 18 to 590 md. The measurements with gases (helium, nitrogen, air and carbon dioxide) covered a pressure range from 0.039 psia to 400 psig. Data were taken for four liquids — n-heptane, methyl alcohol, 79.8 weight per cent glycerol-water solution and pure water at atmospheric pressure. The experimental results were used to evaluate the theoretical equations for predicting stagnant conductivities developed earlier 4 The low-pressure measurements permitted evaluation of the consolidation parameter bpDp/ks (necessary to utilize the theory) for the various types of sandstones. Using these characteristic values, the tbeoretical equations correlated well with the experimental conductivity data for the several fluids and rock samples. INTRODUCTION An aspect of heat transfer in solid-fluid systems of considerable current interest is the effective thermal conductivity of porous media. The stimulus for study of the subject arises from the need for sound procedures for designing thermal methods of petroleum production. The general system occurs when there exists a flow of fluid through the pores of the solid material. However, a logical starting point in developing a theory for predicting the effective thermal conductivity in the general system is to attack the special case when the porous solid is filled with stagnant fluid. Since the flow rates anticipated in thermal production processes are very low, such stagnant conductivities k,O are also of practical significance. Kunii and Smith4 recently reviewed the existing data for stagnant conductivities and proposed a theory for heat transfer in porous rocks. This leads to equations for predicting the stagnant conductivity as a function of the properties of the fluid and solid phases. The existing experimental information1-3,5,7,8 covered a range of porosities and agreed well with the theory. However, the available data were insufficient to examine critically the effects of two important properties on koe, — the thermal conductivity of the fluid and the pressure. The objectives of the present study are twofold: (1) to measure the effect on koe of (a) the pressure for gases, and (L) the thermal conductivity of the fluid kg for bloth gases and liquids; and (2) to determine for rocks of different porosities and permeabilities the apparent solid conductivity ks and the consolidation parameter bpDp/ks required in applying the theoletical equations. Measurements extrapolated to zero pressure are helpful in evaluating ks and the consolidation parameter. However, the extrapolations must be based upon data obtained at low pressures because in small pores the molecular conductivity of the fluid may be unusually depressed. This occurs when the mean free path of the molecules is of the same order of magnitude as the pore diameter, that is, in small-diameter pores at low pressures. The phenomenon has been studied recently by Schotte6in beds of unconsolidated fine particles. With particles in the size range of 45 to 200 microns, a significant reduction of kg with pressure was found even at absolute pressures greater than 1 atm. In the sandstones employed in the present investigation, pore diameters lower than those experienced by Schotte are likely, so that the effect would presumably be greater. SCOPE OF MEASUREMENTS Stagnant conductivities were measured for four types of sandstones using a total of seven samples. The characteristics of the samples, including porosity and permeability values, are given in Table

By J. S. Crump, C. R. Hocott, A. E. Hoffman

Planning of the efficient operation of a gas-condensate reservoir requires a knowledge not only of the gross phase behavior of the system but also of the equilibrium distribution of the various components between the gas and condensate phases. This equilibrium distribution can be calculated with appropriate equilibrium constants. In this paper are presented equilibrium constants determined experimentally for the oil and gas phases initially present in the same reservoir and for the gas and condensate phases of the gas cap material at a series of pressures below the original reservoir pressure. Also presented is a method for the correlation of the experimentally determined equilibrium constants. The utility of the correlation is demonstrated further by an example of its use in adjusting the equilibrium data to permit their application to another gas-condensate system of similar composition. INTRODUCTION Planning of the efficient operation of a gas-condensate reservoir requires a thorouugh knowledge not only of the gross phase behavior of the particular hydrocarbon system but also of the equilibrium distribution of the various components between the gas and condensate phases. At the initial conditions of reservoir temperature and pressure, the original hydrocarbon materials in a gas-condensate reservoir or in the gas-cap of an associated reservoir exist in a single, homogeneous vapor phase. However, some condensation of hydrocarbons to a liquid phase usually occurs in the reservoir as pressure declines incident to production. Because of this condensation. the produced gas changes composition continuously. The composition of the produced gas, as well as that of the condensed liquid, can be calculated from the composition of the original reservoir material through the application of appropriate equilibrium constants or "K"-values provided these values are known. Equilibrium constants have been used for many years in problems of surface separation of gas and oil and natural gasoline recovery. However, no satisfactory equilibrium data at reservoir conditions have been available for the higher boiling hydrocarbons which acquire abnormal volatility at high pressure and are therefore present in the gas phase in condensate reservoirs. Early work in this field indicated that the behavior of the higher boiling hydrocarbons largely determines the behavior of gas-condensate systems at high pressure.' Consequently, a portable test unit was designed to permit determination of the distribution of the heavier hydrocarbons between the gas and liquid phases of a gas-condensate system. This paper describes the test equipment and presents the results of an investigation of the reservoir fluids from an associated oil and gas reservoir in the Frio formation in Southwest Texas. Also presented is a method for correlation of the experimentally determined equilibrium constants, together with an example of the use of the method in adjusting these data for application to another gas-condensate system of similar composition. FIELD TESTING AND SAMPLING In order to determine the compositions of the reservoir fluids, field samples were secured from two different wells, one completed in the oil zone, and one completed in the gas cap. The oil well was completed through perforations in 51/2-in. casing at a depth of 6,822-6,828 ft subsea. Pressure and temperature traverses made in the well indicated the pressure and temperature to be 3,831 psig and 201 °F, respectively, at a depth of 6.828 ft subsea. Subsurface samples were taken from this well after a four-hour shut-in period.

Jan 1, 1953

By S. Kado, G. Derge, L. Yang

Self-diffusion coefficients of silver in molten silver have been measured by means of the capillary-reservoir method in the temperature range 1002" to 1105°C. The results can be .represented by the equation D = (7.10 i 0.06) X 10-4 exp [(—8150 ± 1130)/RT]. The heat of activation agrees within the limit of experimental uncertainty with that for the viscous flow. The validity of the Eyring and the Stokes-Einstein equations are examined on the basis of the self-diffusion coefficients obtained and the known viscosity data of molten silver. MEASUREMENT of the self-diffusion coefficient of silver in molten silver was made in the temperature range 1275° to 1378°K, using the capillary-reservoir method and &"' as the radioactive tracer. The experimental arrangement is shown in Fig. 1. Silver of 99.99 pct purity (fine silver), was melted and degassed in vacuum in a graphite crucible situated inside a McDanel vacuum jacket. The crucible (4 in. tall and 1 in. ID) was protected both at the top and at the bottom with graphite shields and contained about 220 g of silver. The sample holder was a piece of graphite rod (1 in. long and 1/2 in. diam) with 4 capillaries drilled longitudinally into it. The holder was screwed onto a long graphite rod (12 in. long and 74 in. diam) which was connected to a steel disk hanging on a nylon string through a glass pulley. The steel disk kept the sample holder at the center of the McDanel tube and was heavy enough to overcome the buoyancy when the holder was lowered into the molten silver. By turning a handle connected to the pulley through a tapered glass joint, it was possible to raise or lower the sample holder without breaking the vacuum. Two such arrangements were constructed, sitting side by side in 2 Kanthal furnaces. One with non-radioactive silver in its crucible was used for the diffusion run and the other with radioactive silver in its crucible was used for filling the capillaries. To fill the capillaries with radioactive silver, the sample holder was degassed first and then lowered into the radioactive molten silver in vacuum. Purified argon gas was then admitted into the system to push the liquid silver into the capillaries. After cooling in argon, the filled capillaries were polished on emery paper to expose the top of the rod and used immediately as samples for the diffusion run which was carried out at the same temperature as that used for filling the capillaries. During the diffusion run, the samples were degassed in vacuum and lowered slowly into the nonradioactive molten silver before purified argon was admitted into the system. After a predetermined period of time, the samples were raised out of the melt and cooled in argon. The temperature of the melt was measured before and after the run with a chromel-alumel thermocouple inside a McDanel protection tube. Another chromel-alumel thermocouple attached to the outside wall of the McDanel vacuum jacket was used to indicate any fluctuation in the temperature of the furnace during the run, which was found to be less than ± 2°C for the longest run made (93 min). The temperature gradient inside the molten silver was found to be negligible. No significant difference was observed in the results obtained whether the melt was stagnant during the run or was stirred by raising and lowering the sample holder slowly a few times a minute for a distance of about 1/4 in. inside the melt.

Jan 1, 1959

By Julian Szekely

Through the simultaneous statement of Maxwell's equations. the turbulent Navier-Stokes equations and the differential thermal energy balance equation, a mathematical representation has been developed for heat transfer and fluid flow in electric arc furnaces. In broad outline the theoretical predictions appear to be consistent with experimental measurements. A discussion is presented of the general arc furnace control and optimization problem, together with the potential for the incorporation of the mathematical model given in the paper, into future control and optimization schemes. ELECTRIC ARC FURNACE STEELMAKING accounts for about 24% of the total steel produced in the U.S., consuming about 550kw-hr/ton of steel which amounts to about 1.4x1010 kw-hr annually. The percentage of steel produced by the electric arc furnace route has doubled since 1960 and is expected to increase further in the foreseeable future, particularly with the more extensive use of direct reduction technology. Table I shows, in terms of percentages, the principal items that contribute to the cost of running an electric arc furnace. It is seen that in addition to the fixed costs (labor, depreciation, taxes etc.) electrode costs and power consumption are the dominant inputs, and it follows that process changes that would ultimately reduce any of these cost components would have a significant impact on the efficiency of the operation. There are three main areas where progress in arc furnace design and operation must be directed in the next decade (Committee on Technology, 1981). 1. Increased furnace productivity 2. Reduction in processing cost 3. Improvement in the thermal efficiency of the furnace An estimate of the current operating efficiency of the electric arc furnaces may be obtained by comparing the theoretical minimum energy consumption per ton of steel with those achieved in operational practice. Theoretically 330-360 kw-hrs of electricity are required to melt one ton of cold scrap. For superheating, 8-12 kWh per 100°F are required. The actual energy requirements are given in Table II. It is seen that furnaces currently operate at 70-80% efficiency, with the major heat loss from the furnace resulting from the waste gases (Committee on Technology, 1981). While it is apparent that operational experience has produced a technology that is reasonably energy efficient this has been achieved in the almost complete absence of any detailed knowledge concerning the fundamental heat transfer mechanisms occuring in the furnace. The purpose of the present paper is to present a case for the usefulness of mathematical models in increasing the thermal efficiency of furnaces by investigating these heat transfer mechanisms and then combining this knowledge with existing electrode regulator mechanisms to produce truly efficient power cycles for melting and refining. In particular, existing mathematical models describing convection, radiation, and electro-magnetism in the furnace should be employed, and the development of optimization models which balance the cost of electrodes and energy against the thermal efficience of the furnace should be undertaken. Recent work in this laboratory (Szekely, McKelliget, 1981, 1982, 1983; Szekely, McKelliget, Choudhary, 1983) has made a useful start towards solving the problem of heat transfer within the electric arc furnace but the work still needs to be directly related to existing furnace power cycles and tested against plant scale measurements.

Jan 1, 1984

By I. Fatt

The chromalographic effect refers to the separation of constituents in a moving fluid phase which occurs when the phase is passed over a stationary phase, either solid or liquid, or large areal extent. In performing chemical analyses, the large surface contact areal is provided in the case of gas-solid chromatography by-passing the gas over a porous solid and for gas-liquid chronla tog rap hy by-passing the gas over a porous solid which has been coated with a nonvolatile liquid. The porous medium is usually placed in a cylindrical tube, although it may be a relatively flat sheet as in paper chroma-tography. According to one theory,' the separation results from the difference in time spent by various components in the stationary phahe. A gaseous component which is conlpletely insoluble in the stationary phase, or is not adsorbed on the solid, proceeds unimpeded through the column. The versatility and power of the various types of chromatogruphy for the analysis of complex mixtures are evident from their widespread use in performing chemical analyses and from the large number of symposia, books, and papers devoted to the subject in recent years. Associated developments are the application of chromatography to the determination of heats of solution, heats of adsorption, solubilities orK-values, etc. Gas-liquicl chromatography analyses may be characterized by the manner in which the samples are introduced and displaced from the column'. The three techniques used are known as (1) frontal analysis, (3) displacement analysis, and (3) elution analysis. The frontal analysis is performed by displacing an inert gas from the column by a continuous stream of the sample. As indicated in Fig. l(A), all but the first fraction leaving the column represent mixtures rather than pure components. The displacement analysis is conducted by filling the column with a sample, then displacing it with a more strongly adsorbed or absorbed vapor. The resulting chro-matogram is illustrated by Fig. 1(B). Finally, in the elution analysis a small sample of the mixture to be analyzed is introduced as a "plug" into a Rowing stream of carrier gas to give an elution-type chromato-gram. Fig. l(C). The three techniques may be used to provide models to which Row patterns in other physical situations may be compared. PURPOSE OF THE INVESTIGATION The purpose of this investigation was to determine whether natural reservoir rocks, when coated with crude oil, provided sufficient surface areas to give rise to a chromato-graphic separation of light hydrocarbon constituents. The tests were conducted with the light hydrocar-bons in the gaseous state. SANDS TESTED Two sands of widely different permeabilities have been studied: (1) Ventura core (16 md), and (2) Torpedo sandstone core (578 md).Table 1 gives a description of the cores that were tested. With the exception of one frontal-type flow arrangement which was performed on Core 3, the tests were of theelution-type in which a methane-propane mixture was eluted through the core by means of helium. The Ventura core, when impregnated with a weathered crude from the same field, became so impermeable to gas flow that it was necessary to conduct the tests on the granulated cores (Cores 2 and 3).

By D. J. DeLazaro, W. Rostoker, R. E. Riley, M. Hansen

KNOWLEDGE of the isothermal transformation behavior and the TTT chart method of graphically summarizing such information has been of invaluable aid to the ferrous metallurgist in understanding and developing heat treatments intended to provide specific mechanical properties. It is to be expected that similar lines of study directed to titanium-base alloys will be as profitable. This paper summarizes, on conventional TTT charts, the isothermal transformation products and pertinent reaction rate data for binary alloys of titanium containing 1, 3, 5, 7, 9, and 11 pct Mo, respectively. The results so presented were derived from metal-lographic examination. Some X-ray diffraction work was done to clarify certain points of question. Ti-Mo System To the authors' k.nowledge, TTT charts have only been used to describe the reaction rates of eutectoidal decompositions. It would seem that the kinetics of decomposition of any unstable solid solution might be conveniently described in this manner. The equilibrium diagram of the Ti-Mo system' is shown in Fig. 1. The similarity with the Fe-Ni system will be recognized. Additions of molybdenum serve to stabilize the /3 phase to increasingly lower temperatures. The solubility of molybdenum in the a phase is restricted to less than 1 pct at 600°C. Under normal conditions. a homogeneous ,8 structure quenched to a, temperature in the a f & field will reject a of invariant composition. The rate at which this occurs depends on the degree of undercooling and the degree of supersaturation of /3. The structural character of the rejected phase may also be expected to be related to the undercooling. Materials and Experimental Methods A sponge titanium was used which has a nominal impurity content as reported in Table I. The nominal composition of the molybdenum used is also included. The alloys were prepared in 150 g melts by arc melting in an evacuated water-cooled copper crucible of the type reported in an earlier paper.' The pancake-type ingots were forged to 1/4x1/4 in. bars and homogenized for 24 hr at 1000°C in evacuated Vycor bulbs. Specimens Vi in. thick were cut from these bars for heat treatment. The heat treatment cycles followed an invariant procedure. The specimens were soaked at 1000°C for 20 min under a helium protective atmosphere prior to quenching into a lead bath held at a temperature corresponding to a chosen degree of undercooling. After a prescribed period of time at this temperature, the specimen was removed and water quenched. A series of specimens so treated for progressively longer periods served to establish the times for initiation and completion of the decomposition of the /3 phase at a particular temperature. Although heat treatment of unprotected specimens in a lead bath does permit surface contamination by lead diffusion, this method is necessary to insure sufficiently rapid quenching rates. As long as the bath treatments were of less than 1 hr duration. the contaminated surface was not unduly deep and could be removed without excessive grinding. Metal-lographically, the contaminated zone could be readily distinguished from the pure binary alloy. Transformation Characteristics of Ti-Mo Alloys Metallographic examination showed that a supersaturated p phase may decompose by one of two processes. Alloys containing 1 to 9 pct Mo, when quenched in water from the homogeneous /3 field,

Jan 1, 1953

By D. Harker, E. Epremian

The constitution of the nickel-tungsten system has been studied by a number of investigators, the most recent of which are Ellinger and Sykes.1 On the basis of metallography, electrical resistivity and hardness measurements and some X ray diffraction work, they constructed a consitution diagram (Fig 1). Ellinger and Sykes report a0 for the alpha face-centered cubic phase (nickel saturated with tungsten) as 3.66?.‡ The gamma phase (tungsten saturated with nickel) is given as body-centered cubic, but no value of a. is recorded. At 43 pct tungsten by weight (approximately 20 at. pct) a beta phase is reported, but structure and lattice parameter are not known. Primarily on the basis of metallographic evidence, they conclude that beta is a new phase and rule out the possibility of formation of a superlattice since they- observe a new phase in the microstructure. Grube and Schlecht2 report an ordering reaction in a system similar to the Ni-W (Ni-Mo) and Harker3 has shown that NirMo forms a superlattice upon the pre-existing face-centered cubic lattice of the alpha phase. Elliniger arid Sykes observed a change in the diffraction pattern of a 4.3 pct W alloy upon aging and attributed it to the formation of a new phase. Thus, there appears to be some question as to the character of the beta phase reported by Ellinger and Sykes. The purpose of the present work is to determine precisely the structure and nature of the beta phase in the Ni-w system of X ray diffraction methods. Experimental Procedure One inch diameter ingots of 43.66 pct, 40.15 pct and 34.06 pct .tungsten by weight (19.8, 17.6 and 14.1 at. pct respectively) were melted and cast under vacuum from electrolytic nickel and pure wire filament grade tungsten. The ingots were soaked for 16 hr at 1300°C in a hydrogen atmosphere furnace, then swaged at this temperature and finally drawn to 0.020 in. diam wire. Some of the bar stock was retained in 1/2 in. diam rod for hardness and metallographic studies. Wire samples for diffraction work were sealed in evacuated quartz tubes and subjected to the various solution and aging heat treatments in a hydrogen furnace. X ray photograms were taken in a Debeye-Slierrer camera using copper K radiation filtered by nickel foil. Results X RAY DIFFRACTION The lattice parameters of the alpha and gamma phases at saturation were determined by an X ray photogram of the 19.8 at. pct W alloy after solution heat treatment at a temperature well within the alpha-gamma range above the peritectoid temperature (17 hr at 1150°C—oil quenched). The alpha phase was found to be face-centeted cubic with a. of 3.594 ± 0.001?, while the gamma phase was determined as body-centered cubic with a. of 3.158 ± 0.001?. Exactly the same phases and parameters were ohtained with the 17.6 at. pct W alloy for the same solution heat treatment. Thorough analysis was made of an X ray photogram of a 19.8 at. pct alloy sample in the quenched and aged condition (solution treated 1T hr at

Jan 1, 1950

By E. F. Sumner, C. K. Eilerts

A partial differential equation was developed and programmed to compute the cumulation and flow of the liquid phase of a gas-condensate fluid. It was found that such cumulation near the wellbore may attain saturations permitting liquid flow within periods of the order of 1 day. The radius at which this saturation obtains advances toward the external reservoir boundary with continued recovery, although it may move only 2 to 5 percent of that distance during the productive life of the well. This relatively high saturation can cause a serious dimunition in delivery capacity. The pressure dilferential required for flow of gas through a formation with such saturation may be as much as three times greater than the corresponding differential in a formation where saturation is less than 2 percent. Because of the radial geometry, total reservoir pressure differential required for flow with a zone of relatively high liquid saturation about the well was indicated to be two to three times greater than without the saturation. INTRODUCTION A computing program developed by the authors 1 to integrate the partial differential equation for transient radial flow of gas-condensate fluids required a prevalued function, W(Po) = kD(p0)/M(p0)z0(Po).......(1) related to cumulation of condensate in pores of the reservoir as pressure declines. Methods were available for relating mobility to saturation, saturation to liquid-vapor volume ratio, liquid-vapor volume ratio to pressure and, hence, mobility to pressure to obtain W(pD). However, certain assumptions had to be made in advance of program execution concerning movement of the liquid phase and change in radial distribution of liquid saturation as recovery progressed. Because reservoir fluid mobility and related pressure gradients are markedly dependent on liquid saturation, it was desirable to eliminate the necessity for these assumptions and in this respect improve the reliability of a result obtained by integration. Experimental research was begun to determine the relationship of variables that control the reservoir flow of low liquid-vapor ratio fluids. With the porous medium partly saturated, permeability to a two-phase flowing fluid could not be measured because the viscosity was unknown. Viscosity could not be measured without knowing permeability of the liquid-containing medium; but mobility, the ratio of permeability to viscosity, was measurable. This property, and liquid-vapor volume ratio of the flowing phases, have values determined by nature of the fluid and porous medium and, in addition, were found to be pressure, velocity and saturation dependent. A program developed by Kniazeff and Naville 2 for computing multicom-ponent, two-phase flow used polynomials to represent saturation dependence of permeabilities for gas and for liquid phases, and an effluent viscosity constant under reservoir conditions. The authors developed a program based on a partial differential equation of similar mass balance origin but designed to represent distribution of an unlimited number of components between gas and liquid phases utilizing the experimental information available as to variable dependence. Attention was given to liquid-vapor volume ratio of the flowing phases, a property which is essential as a liquid saturation boundary condition. A program for computing implicitly saturation increase at uniform intervals along the drainage radius was combined with an existing program for calculating pressures, and the two were used alternately in an iterative solution for both saturations and pressures. The function W(pD) needed for calculating pressures was replaced by the ratio of mobility to compressibility factor. BASIC EQUATIONS The partial differential equation for calculating

By W. D. Brown

Newell1 has shown that hydrogen is removed from steel in a vacuum at a temperature of 500° to 900° C. within 136 hr. Holm and Thompson2 also state that, especially when the hydrogen is high, the results obtained by extraction at 800°C. agree with those obtained by vacuum fusion. In this laboratory, an apparatus was designed (Oct. 1940) and used wherein the piece of steel was heated at 100°F. in a Torricellian vacuum (Fig. I). This vacuum was produced by raising and lowering a "leveling bottle or bulb" containing mercury. Extensive tests were made with this apparatus. The weakness of the procedure was the possibility of leakage of air into the tube around the rubber stopper at the bottom. The present plan is to raise and lower the mercury by applying vacuum from a mechanical pump, to the top of the tube and then to a bottle of mercury in which the bottom of the tube is immersed. The Apparatus The extraction is made in a Pyrex tube A, I in. i.d. It is drawn down at the top at N and an 8-mm. i.d. tube is welded on it. To this a Pyrex stopcock E, 4-mm. bore, is welded. The tube between N and E is preferably 8 in. long and is graduated in milliliters. The total length of A to the stopcock is 28 inches. Heat is supplied by a split-tube furnace. 0 in. long by 6 in. in diameter, hinged on the side. A short pyrometer couple P is inserted through the walls of the furnace. In use this is standardized by comparison with a standard couple in a I-in. tube, just opposite this short couple. In the I-in. tube A a glass rod R, flattened at one end, supports the sample S and the length of R is so adjusted that the sample is in proper place. (It may float in the mercury.) The tube A is inserted in a large rubber stopper to fit the bottle B and extends nearly to the bottom. It is so placed that the handle of the stopcock faces the operator. In this stopper also is a short tube D bent at right angles for the application of vacuum. Also in this stopper is an 8-mm. tube C reaching nearly to the bottom; this is bent so as to miss the heater and again bent to bring it in line with and as near as possible to the upper and reduced part of tube A. On this tube C is welded a 4-mm.-bore stopcock H. In order that the two tubes may be as near as possible, the handle of stopcock H faces the back, or is the reverse of stopcock E. The total height of II above the bottom is 28 in. As explained later, this is the pressure-measuring tube C. The bottle B is about half full of mercury. The mercury must cover the ends of A and C by about I in. when they are evacuated. By tees T and rubber tubing RT the tubes A, C and D, the latter through the three-way stopcock M, are connected to the vacuum pump. A sheet of asbestos at K prevents the stopcocks E and H from becoming so hot that they leak. A good grade of vacuum lubricant prevents leakage of E and H. In order to measure the pressure of the

Jan 1, 1945

By Walter J. Karplus

By means of a finite difference expansion a fluid flow field in cylindrical coordinates with axial symmetry, is simulated by a network of electrical resistors. A series of DC analog computing units, connected to the network bounrlary corresponding to the free surface effectively adjust the network parameterr to satisfy simultaneously Laplace's equation and the borrndary conditions. The system automatica1ly "re1axes" to the correct boundary shape in a fraction of a second. A detailed application of the method to the water coning problem in oil wells is presented. INTRODUCTION The water coning problem in petroleum production is a prime example of a class of problems known as free surface problems. A free surface arises in fluid flow regions when the shape of the boundary of a single phase flow region is affected in some manner by the pressure or potential gradients within the region. When there exists an interface of two fluids of different densities, the shape of this interface will be determined by the relative velocities of the two fluids and will therefore comprise a free surface. Other free surface problems of interest to petroleum engineers include the gas coning, water injection, and gravity drainage probIcms. In each case it is desired to predict the shape of the boundary of the oil zone as a function of the flow rate, the well penetration and the physical parameters and dimensions of the system. Such problems do not lend themselves to the classical analytical treatments of boundary value problems because in this case the boundary of the potential field and the potential distribution within the field are interdependent and must be calculated simultaneously. A further difficulty arises from geometric considerations. Since the fluid flow is of the radial type in most cases, such effective methods as conformal transformation and the hodograph are not applicable. It is the purpose of this paper to present a method whereby free surface problems of this type may be solved rapidly to the desired order of accuracy by means of a combination of a network of electrical resistors and conventional analog computer units. The water coning problem is treated in some detail. Then the procedure for extending the technique to other free surface problems is indicated. FORMULATION Most commercial oil reservoirs consist of an oil-saturated sand, bounded above by an impermeable strata and below by a water-saturated sand. Gravimetric separation results in a horizontal oil-water interface. When an oil well is drilled so that it partially penetrates the oil zone, and production is initiated at a constant rate, the oil-water interface will assume a new shape as a result of the pressure gradients in the oil zone. It should be emphasized that this is a steady phenomenon. After an initial transient period the water phase will be completely static and only the oil will be in motion. The new steady-state shape of the interface, as shown in Fig. 1, depends upon the depth of penetration, d, the thickness of the oil zone, I, the production rate, q, the hydraulic permeability zone of the sand, k, the viscosity of the oil, µ, and the densities of the oil and water, p.,. and pw.. The problem is to determine this relationship. Darcy's law + = - (P f g p z).....(1) it.

Jan 1, 1957

By Stanley L. Lopata, Eric B. Kula

The problem of preferred orientation has been considered in quantitative phase analysis by X-ray diffraction techniques. The average intensity of a diffraction peak can be obtained by integration over the whole pole figure, or by random rotation of the specimen in the X-ray beam, or by a combination of both. For rod or wire, where the pole figure shows rotational symmetry, the intensity is measured along a radius of the pole figure. Where no symmetry is present, rotational symmetry is imparted by rapid rotation of the sheet around its surface normal during measurement. These methods are applied to the determination of aus-tenite in an Fe-31 pct Ni rod and Fe-28 pct Ni sheet. X-RAY diffraction techniques are frequently used for quantitatively determining the relative amounts of two or more phases present in a structure.1,2 Briefly, the method requires a measurement of the intensity of a Bragg reflection for each phase. The method is based on the principle that the intensity of a Bragg reflection, IA, is directly proportional to the volume fraction of the given phase present, CA: T KRACA/2µ [1] where K is a constant independent of the diffracting phase, RA can be calculated and depends on the Bragg angle ? and the particular reflection measured, and µ is the linear absorption coefficient of the mixture.' For a mixture of two phases A and B, IA = RAcA [2] IA/IB=RAcA/Rbcb and CA+ CB = 1 [3] This allows the volume fraction of a phase to be calculated directly from the measured integrated intensities of the two phases present. This method is frequently used for quantitatively determining retained austenite in steels, but can be used for the quantitative determination of in titanium alloys, as well as phases in other alloy systems. The intensities are generally measured by diffracting only from planes whose normals are perpendicular to the specimen surface, with the implicit assumption being made that this intensity is typical of reflections from planes whose normals lie at other angles to the surface. If the specimen contains a texture or preferred orientation of the crys- tallographic planes by virtue of previous mechanical working, however, this measured intensity is correct only if the intensity level fortuitously is one times random. Examination of pole figures for actual rolled sheet shows that the intensity level at the center of the figure is seldom at the random level, and can very easily differ by a factor of five or more, leading to a commensurate error in the calculated amount of the phases. Textures have been commonly found in materials which are cold-rolled or recrystallized after cold rolling. This problem is becoming of greater importance since structural materials in use today, such as the high-strength stainless steels, frequently are further strengthened by cold rolling, and hence contain textures. Although high-strength carbon-containing alloy steels generally do not contain textures, some of the newer thermomechani-cal treatments will impart a texture even in these cases. In the case of hexagonal metals, such as titanium, probably every wrought sample will show some anisotropy. Since the standard X-ray techniques are inapplicable in all these cases, there is a clear need for some method of taking this anisotropy into account. The specific problem is to obtain the integrated intensity of a Bragg reflection for each phase, averaged over all orientations of the sample with respect to the X-ray beam. This can be done either by purely mechanical means or by a mathematical averaging. Mechanically, a small spherical speci-

Jan 1, 1965

By J. S. Aronofsky, J. P. Heller

This paper presents a mathematical analysis of the fluid mixing which occurs during flow through porous media. The analysis is based on the well-known diffusion equation with mass transfer term. It is pointed out that the use of this equation is justified by a single general assumption which does not specify any particular mechanism for the mixing process. Formulas are given for tracer fluid concentrations for two different boundary and initial conditions. Calculated numerical values compare closely with some published results of experiments in which no viscosity or density difference existed between displaced and displacing fluids. INTRODUCTION A phenomenological theory of the mixing and diffusion process is presented for the flow in a porous medium of two miscible liquids of equal viscosity and density. This is based on the classical diffusion equation which has been used to explain such processes as Brownian motion and heat conduction. More recently, this same differential equation has been derived (in a manner analogous to the work in this paper) for the problem of heat transfer during fluid flow in porous media.' The purpose of this paper is to show that the available experimental evidence justifies the use of this simple diffusion equation. The recent data of Koch and Slobod' and of von Rosenberg' are used to test this diffusion model by comparing concentration curve shapes and calculating "effective diffusion coefficients". The same differential equation for miscible displacement, as well as some experimental confirmation, has been published in Japan by K. Yuhara,4 and recently in this country by Day.5 Yuhara arrives at his con-clusions through analogy with turbulent diffusion, de- scribing the microscopic velocities in terms of "eddy motion" — although unlike the ordinary eddies of hydraulics. Day proceeds directly from Scheidegger's6 concept of the dispersion of velocities. After a very clear description of Scheidegger's theory, he applies it to calculate the motion first of a drop and then of a thin layer of brine into a bed of sand. For the latter case, Day derives and solves the same differential equation as Eq. 1. In our presentation no detailed assumptions are made as to the mechanism of displacement. In fact, several microscopic processes may well be involved—the only general assumption made is that their cumulative microscopic effect may be described by the addition of a second-order term to the differentirl equation of continuity. DERIVATION OF EQUATIONS Let us consider the displacement in one spatial dimension of Liquid A from a sample of porous media by Liquid B. A and B are assumed completely miscible and of the same viscosity and density. Both fluids are also assumed to be incompressible and it is specified that the total flow is the same for any cross section. It is evident that the rate at which Liquid B would be carried across such a section would be given by the product of the fluid velocity, u, and the concentration, C, of B in the fluid mixture assuming the constancy of u and C over the cross section. Since the volume of Fluid B is conserved, this leads to the one-way wave equation - u ?C/?x = ? ?C/?l where x increases in the direction of u—the "Darcy velocity", Q/A (or velocity flux per unit area), and 9 is the fractional porosity. Consider the case where fluid B is displacing Fluid A and the concentration distribution is represented by a step funcLon (with a vertical discontinuity). The solution of Eq. 1 predicts no mixing

Jan 1, 1958

By J. F. Evers, S. Sadiq, F. W. Preston, G. W. Swift

Proper scaling of unsteady-state flow in laboratory gas reservoir models can be achieved if the permeability/fractional porosity ratio is less than 0.5 md. Porous media meeting this specification can be obtained from mixtures of sand and Portland cement provided (1) the water/cement ratio is in the range of 0.35 to 0.40, (2) the sand/cement ratio is in the range of 1/1 to 2/1, (3) the sand grain size is in the range of 48-100 to 18-28 mesh, (4) the mortar curing time is 7 days or more and (5) the cured mortar is dried under vacuum to reduce the temperature level employed during the drying process. The areal homogeneity of the mortars tested is reasonably good, although significant trends in vertical permeability were noted. The Klinkenberg and non-Darcy coefficients of mortar are such that proper scaling of these effects in the laboratory model could be expected. INTRODUCTION A long term deliverability prediction for a gas well is generally based on data determined from flow tests that are of short duration relative to the producing life of the well. These data, together with an estimate of gas reserves in place, form the basis for establishing the parameters of the mathematical model used to predict performance at future times. Results obtained from a mathematical model represent only an interpretation of the physical process being modeled. Experimental verification of the production process from discovery to abandonment is needed where production occurs under carefully controlled conditions. Ideally, one might consider such an experiment to be conducted in the field. Even if the cost of such an experiment could be justified, the time required would make it impractical to consider. However, if the time scaling problem could be handled properly, an experimental model of the gas reservoir production process could provide the necessary data to study the mechanisms of the production process in a reasonable time period. An experimental model has been designed and constructed using the data of this study as a basis. The model reservoir system consists of one producing well located at the center of a right-circular cylinder of porous material. Ten to 15 observation wells are located at strategic offset positions from the producing well. These wells, and the producing well, are fitted with pressure transducers so that the pressure-time history can be recorded at various radial and azimuthal positions in the drainage area during the entire production life of the system. Time is scaled in the model so that the time required to produce gas from discovery pressure to abandonment pressure is about 2 hours. This compares with a 20-year producing life in the model's field counterpart and gives a time scale factor of approximately 1:100,000. DESIGN CRITERIA To establish design criteria for the experimental model, the reservoir is considered to be a finite right-circular cylinder of radius rr, and height h with a producing well of radius rw located at the center and completed through the producing formation. This reservoir consists of a porous medium which has a uniform permeability k and gas porosity and the gas contained in the reservoir has density p, viscosity µ and molecular weight M. Initially, the pressure at all points within the reservoir is uniform at po. and the reservoir temperature is T. Since the reservoir production process has been demonstrated to occur in two distinct transient modes for either Darcy or non-Darcy flow,l it is

By M. E. Wadsworth, J. R. Lewis, J. M. Quets

Samples of snythetic magnetite were reduced in hydrogen at various partial pressures and temperatures. The reaction mas found to be surface controlled and directly proportional to hydrogen partial pressure over the temperature range studied (400° to 1000°C). A transition in the mechanism occurred at approxirnately 590°C. This has been attributed to a change in the slow step of the reduction reaction resulting from the presence of stable wüstite at higher temperatures. Below 590°C the heat of activation for the reduction was found to be 14,700 cal per mole and it was 3,200 cal per mole above 590°C. THE reduction of the oxides of iron has been a subject of study and interest for many years. Edstrom1 has reviewed the work up to 1953 and his findings provide much basic information concerning details of the reduction of hematite (Fe2O3 ) and magnetite (Fe3O4) with both H2 and CO. He concluded that below 600°C the reduction of both hematite and magnetite was essentially identical, since the rate was controlled by a surface reaction at a magnetite gas interface and the iron formed was very porous. The decrease in the reduction rate for both hematite and magnetite immediately above 600°C was attributed to the formation of a stable wüstite phase across which diffusion was very slow. The wide difference in reduction kinetics for hematite and magnetite above 600°C was explained in terms of the character of the wüstite formed. The rate of magnetite reduction increased slowly with temperature beyond the minimum immediately above 600°C and rates equivalent to those observed at 550°C were not obtained until the temperature was increased beyond 850oC. Edstrom and Bitsianes2 studied the solid-state reduction of magnetite, the reaction of which can be written (4x-3)Fe + Fe3O4—4 FexO [1] where x is between 0.83 and 0.95. The rate was found to be controlled by the kinetics of transport of iron ions through the wüstite phase. Bogdandy and Riecke3 studied the kinetics of reduction of small grains of iron oxide in hydrogen and qualitatively attributed the slow step in the reduction process to the diffusion of hydrogen and water vapor through pores. McKewan4 has analysed the data of Joseph,5 Stalhane and Malmburg,6 and Lewis7 and has demonstrated that the kinetics of reduction of both hematite and magnetite are controlled by surface reactions and also that the variation in surface area during reduction must be accounted for. McKewan was able to correlate the reduction data by the equation rodo[1-(1-R)1/3] = Kt [2] where ro is the initial particle radius, do is the density of the iron oxide, R is the fraction reduced at time t, and K is a rate constant. The applicability of Eq. [2] for temperatures above 600°C where a stable wüstite phase is present, indicates the slow step in the over-all kinetics does not result from diffusion unless diffusion occurs laterally along the oxide metal interface or through a wüstite layer of constant thickness. Diffusional control through a layer of wüstite is unlikely, however, in view of the pressure effects demonstrated by Tenenbaum and Joseph,' Diepschlag,9 McKewan,10 and the results of the present study. McKewan has demonstrated a linear rate dependence on the hydrogen partial pressure over the temperature range 400" to 1050°C for the reduction of hematite.

Jan 1, 1961

By M. D. Hassialis

[Surface Areas of Concentrates and Collector Coatings (T.P. 2002, by A. M. GAUDIN and G. S. PRELLER, Min. Tech., May. Discussion by M. D. HASSIALIS and the authors) . I Activation of Minerals and Adsorption of Collectors (T.P. 2082 by J. ROGERS and K. 1'. SUTHERLAND, Min. Tech., Nov. Discussion by A. F. TAGGART and the authors). . . . 3 Correlation between Mineral Behavior in Cataphoresis and in Flotation (T.P. 2005 by A. If. GAUDIN and S. C. SUN, Min. Tech., May. Discussion by M. D. HASSIALIS). . . . . . 8] Surface Areas of Concentrates and Collector Coatings M. D. HASSIALIS.*-In the current flotation theories there exists a common gap spanned by a common assumption pertaining to the arrangement and orientation of the collector ions at the solid surface. Also common is the belief that contact angle and flotation performance are dependent upon the arrangement and orientation of collector ions. This stimulating and informative paper sheds some light upon this general problem. One of the indices of arrangement is the area available per collector ion; in this connection, Gaudin and Preller have applied the B.E.T. method for surface determinations. In a surprising number of cases, their results show a greater specific surface for tailing than for concentrate. This cannot be accounted for by the differences in specific gravity and is attributed to the presence of clay minerals of "primary" sluices in the tailing. This is not a satisfactory explanation, for experience shows that when free slimed clay minerals are present in a flotation cell they will be carried mechanically into the concentrate. A part of the answer to the decreased specific surface of the concentrate would appear to be its usual flocculated condition, yet this should not entirely account for the large spreads reported. The possibility occurs that the nitrogen adsorption isotherms for clean mineral and collector-coated mineral are different. If that were true, and if the coated mineral surface exhibited decreased adsorptive powers, the discrepancy might vanish. In discussing the disposition of collector ions at mineral surfaces, the authors adopt the more consistent viewpoint that there exists some agreement between the ions of the substrate and those of the superstrate. In the case of xanthate ions on galena, a possible arrangement is shown in Fig. 2. This arrangement shows every xanthate ion being shared by two lead ions, whereas the chemical composition calls for two xanthate ions for each lead ion (p. g). This whole discussion is misleading. In the first place, it is known that xanthate ions are not abstracted by unoxidized galena, hence Fig. 2 should show as the substrate not Pbs but PbS„ Om. It would be equally erroneous to use the lattice arrangement of PGSO4, since the surface properties of xanthate-conditioned anglesite differ from those of xanthated galena. It is highly probable that the state of the galena surface exploited in its flotation is ephemeral, since weathered galena is only poorly, if at all, collected by xanthate. Just how a substrate of oxidized galena is to be pictured depends upon one's conception of the transformation from galena to anglesite. This writer is inclined to the following sequence of events. Surface sulphide ions undergo progressive oxidation. As more oxygen is taken up by the sulphide ions they become displaced from their equilibrium positions in the galena lattice, thus accommodating for the larger volume of the S„ Om ions. Concomittant with this change the forces acting on the lead ions

Jan 1, 1946

By H. R. Thresh, A. F. Crawley

THIS paper reports the results of viscosity measurements on cadmium and indium undertaken as part of a continuing program of liquid metals research at the Mines Branch, Department of Energy, Mines and Resources, Ottawa, Canada. For several years this work has continued with support from the Canadian Zinc and Lead Research Committee and the International Lead and Zinc Research Organization. In the literature, existing viscosity data for cadmium are conflicting and for indium are scarce. To achieve reliable viscosity values, further work was deemed desirable. Viscosities were measured using an absolute technique applied to an oscillating closed right cylinder of liquid metal. The apparatus, constructed in this laboratory, and the techniques involved in its use have been described in an earlier paper.' In calculating viscosity coefficients from the experimental data. Roscoe's2 mathematical analysis was used. Because accurate viscosity determinations are dependent on accurate density data, and since the available data in the literature are conflicting, the densities of liquid cadmium and indium were also determined in this laboratory by a pycnometric technique. This increased the control exercised over all the variables contributing to the accuracy of viscosity measurements. For both density and viscosity determinations, cadmium and indium of 99.999 pct purity were used. These materials were obtained from Cominco Ltd. In order to suppress excessive evaporation losses and help prevent oxidation, a technique developed in measuring the viscosity of zinc' was used for cadmium. This involved maintaining an atmosphere of purified hydrogen in the viscosimeter during the rela- tively long periods while the experimental temperature was changed and stabilized between individual determinations of the logarithmic decrement. Then, immediately prior to a decrement measurement. the viscosimeter was rapidly evacuated to a pressure of about 2 X 10-6 torr. On the other hand. with indium. no difficulties arise from vaporization or oxidation. In this case. therefore. the apparatus was held under a constant vacuum of about 2 x 10-6 torr. When all sources of random and systematic error have been evaluated, the technique is considered to yield data accurate within 10.5 pct. The viscosity of cadmium was measured to a maximum temperature of 451.8°C and that of indium to 346.5°C. These temperatures are. respectively. 130° and 190°C above the melting points. Two separate runs were made on indium. The results are given in Table I and, when plotted. yield the smooth curves evident in Figs. 1 and 2. Values of A and C to fit the well-known Andrade1 equation were determined by a regression analysis. The equation is: where n = coefficient of viscosity. i = specific volume. T = temperature in °K Despite its empirical nature. Andrade's equation has been shown by many workers to describe accurately the temperature dependence of viscosity. During the current research program, this equation has been found to fit the data precisely for all the metals and alloys that have been studied.5,6 The A and C' values for cadmium and indium are given in Table 11. Comparing the magnitude of the present data with those of previous workers,7-10 as has been done in Fig. 1, it is evident that the results of Fisher and Phillips10 deviate markedly from the rest. With the exception of one point from Gering and Sauerwald,7 all other data agree within ±2.5 pct. However, the agreement between the present results and the work of Menz and sauerwald9 is especially good over the entire experimental temperature range. Such agreement may be expected since Menz and Sauerwald's

Jan 1, 1970

By W. R. Foster, J. M. McMillen, A. S. Odeh

The complete equations of average linear momentum balance for a single-phase fluid in an incompreisible, homogeneous, porous medium are derived. The derivation begins with Euler's equation of motion for a continuum and uses an integral transform recently developed by Slattery. 7 For steady flow of a compressible, Newtonian fluid, the usual equations of motion result. For transient flow, the space-time description of the pressure is determined in the lowest approximation by the telegrapher's equation. From the analysis a new phenomenological coefficient results which connects the viscous traction to the deriv ative of the linear momentum density. The magnitude of this coefficient determines the velocity of sound through the pore structure in this approximation to the pressure field. INTRODUCTION The modification of Darcy's law of momentum balance for steady, single-phase flow through porous media has been discussed for many years. The first modification was suggested by Forcheimerl who added terms of higher order in the velocity. These can be expected to appear because the underlying microscopic equations of momentum balance are themselves nonlinear in the point velocity field. The Reynolds tensor pvv, which represents the convective flux of momentum density, appears in the momentum balance equation. Only in rectilinear flow (parallel stream lines) does the divergence of this tensor vanish. Since the steady flow stream lines in most porous media are not parallel, nonlinear dependence of the pressure gradient on the velocity should naturally appear. This nonlinearity has nothing to do with turbulence in the ordinary sense of random fluctuations in the pressure and velocity fields. It arises simply because the stream lines converge and diverge, even for steady flows. Klinkenberg2 demonstrated that the permeability coefficient in Darcy's law depends on the absolute pressure or, alternatively, on the density field. However, because he neglected inertial terms of the Forcheimer type, his correction coefficient could not be represented by a constant but tended toward a constant as the velocity decreased. Forcheimer's and Klinkenberg's modifications can be combined in a rigorous way to account for both inertia and slip during steady flow. This will be shown in a future paper. The transient change of pressure in porous media has been described by the diffusion equation.3 This form results from eliminating the velocity and density fields from a combination of the equations of motion in the form of Darcy's law, the continuity equation and an equation of state. Fatt4 suggested that the cause of deviations from the prediction of the diffusion equation for pressure transients lies not in the choice of Darcy's law as the equation of motion but on the existence of dead-end pores which might invalidate the averaged equation of continuity. On the other hand, Oroveanu and Pascal5 noted chat the time derivative of the momentum density must be included in the equations of motion since this measures the local rate of change of momentum density. Their differential equation for pressure is the telegrapher's equation (neglecting gravity). However, their form of this equation predicts that the speed of pressure propagation through the pore structure is the same as that through the bulk fluid. M. K. Hubbert6 attempted a derivation of Darcy's law by volume averaging the Navier-Stokes equations. Since these equations represent momentum balance at a point within an open set of points containing the fluid itself, Hubbert's volume averaging cannot lead to tenns involving transfer of momentum between the fluid and the walls of the pores. Once these viscous tractions are lost by choosing a control volume containing only fluid, they cannot