Search Documents

By J. L. Huitt, S. R. Darin

The use of a partial monolayer of propping agent to obtain a high flow capacity for a hydraulically induced fracture is discussed. From the results of laboratory work it was shown that a modified form of the Kozeny-Carman relation could be used to describe the flow in the partial monolayer propped fracture. With equations presented in the paper, the density pattern of the propping agent (number of particles per unit of fracture surface) that results in the maximum flow capacity for the fracture can be determined. The maxitnum flow capacity obtained with a partial monolayer is often an order of magnitude greater than the flow capacity obtained in greater width fractures containing multilayers of the propping agent. INTRODUCTION One of the predominant factors controlling the success of a hydraulic fracturing operation is the propping of the fracture. The trend in hydraulic fracturing recently has been an increase in the ratio of propping agent-to-fluid. The primary purpose of this increased ratio is to sustain a propped fracture of greater width. It was shown in a recent study' that for some formations low concentrations of propping agent would result in "closing" or "healing" of the fracture. In some cases, in an effort to insure sufficient propping agent concentration, the fracturing operation is designed to obtain a pack of the propping agent in the fracture. The placing of a pack of propping agent in the fracture provides a fracture of maximum width; however, the width alone does not control the flow capacity of the fracture. The flow capacity is dependent on the permeability of the fracture, as well as the fracture width. Thus, increasing the permeability to obtain larger flow capacities is as important as obtaining a greater width fracture. By increasing the concentration and particle diameter of the propping agent, fractures of sufficient flow capacities are obtainable for effective well stimulation in a number of formations. However, in some of the recovery techniques, the fracture capacity obtained by commonly used practices in fracturing is inadequate. One such recovery technique is gravity drainage via a horizontal fracture placed near the base of the productive zone. For this recovery technique, the fracture capacity desired is usually one to two orders of magnitude greater than the fracture capacity normally obtained in fractures propped with the size of sand commonly used (20-40 mesh). By the use of particles of larger diameter, the permeability of a pack2 is increased; however, even a multilayer pack results in a fracture flow capacity lower than that desired for some formations. For example, a fracture of 1/4-in. width, packed with 8-12 mesh sand of about 2,000 darcies permeability, would result in a fracture flow capacity (permeability X fracture width) of only 40,000 md-ft. If the fracture flow capacity needed is an order of magnitude greater than that cited in the example, an approach different from the formation of a pack in the fracture is needed. One approach is a partial monolayer of large-diameter propping agent. In working with this type of fracture system, it was found that the flow problems were unique to the system; therefore, the study discussed in this paper was made. THEORY AND RELATED STUDIES Previous studies, although related to that described in this paper, have not pertained directly to the flow in a fracture propped with a partial monolayer of propping agent. However, studies such as flow between parallel plates,' flow in packed columns,' fracture flow capacity' and flow in simulated fractures reduce the scope of this flow problem considerably. A fracture propped with a partial monolayer of dropping agent may be considered to be between two extremes — between an open fracture and a bed of particles packed In a conduit of fracture geometry. Thus, it seems appropriate to review the extreme cases.

By Robert E. Green, T. Hinton

Expressions are derived for tke orientation factors appearing in the Grana to -Lucke theory of dislocation damping, for fcc crystals. The factors are given as functions of the elastic constants, the direction cosines of the specimen axis, atld the direction cosines of the particle displacemeizts, for darnping due to dislocation motion on each of the twelve possible slip systems as ixell as lor simultaneous dislocation rnotion on all twelve systems. RECENTLY Green' has outlined the general procedure for calculation of the modes of wave propagation in a single crystal possessing a given crystallo-graphic orientation. In addition he has suggested how the results of such calculations along with calculations of the resolved shear-stress factors for slip enable the orientation-dependent factors (orientation factors) appearing in the Granato-Lucke theory of dislocation damping to be calculated. However, reen' assumed that the maximum stress on dislocations was in the direction of the particle displacement of the elastic wave giving rise to the dislocation damping. The present cal'culations show that this is an oversimplified approach appropriate only to special cases of high crystallographic symmetry. The exact calculation of stress on the dislocations can only be obtained by resolving the nine components of stress associated with the elastic wave into the appropriate slip direction and then summing. In this paper the general expressions for orientation factors are derived for fcc crystals in terms of the elastic constants, the direction cosines of the particle displacements, and the direction cosines of the specimen axis, bearing in mind the point mentioned above. Equations are presented which give the orientation factors for damping due to dislocation motion on each of the twelve possible slip systems as well as for simultaneous dislocation motion on all twelve systems arising from elastic-wave propagation. These orientation factors are appropriate to attenuation measurements made by the pulse-echo technique, where the elastic wave is a traveling wave. They will be compared with those for the case of a standing wave as derived by Granato and Lucke which are appropriate to measurements of internal friction by the resonance technique. The expressions for the orientation factors are compared with the five values given by insruch,' based on unpublished work of Granato and Liicke, for positions of highest symmetry in the stereographic triangle. The various equations derived for modes of wave propagation and orientation factors have been evaluated for the case of aluminum crystals at 1-deg intervals throughout the stereographic triangle. These results will be presented in the form of contour diagrams in a subsequent paper3 and their relevance to ultrasonic measurements discussed. ORIENTATION FACTORS The coordinate system used for the calculations is shown in Fig. 1 and all specimen axes are taken to lie in the shaded standard triangle. The twelve slip systems are indexed according to Table I. Following the suggestion of lers, we may define the orientation factor Ok in the Granato-Liicke theory5 by the equation

Jan 1, 1967

By James S. Browning

For the past several years, USBM has conducted laboratory and continuous process development work on the weathered mica pegmatites ores of Alabama, Georgia, and North Carolina to determine the feasibility of recovering commercial-grade mica concentrates by flotation. The research has resulted in the development of two processes for flotation of mica: 1) an acid-cationic method for recovery of coarse mica and 2) an alkaline-anionic-cationic method for recovery of fine mica. The research indicated that the two processes may be effectively utilized separately or in combination for flotation of mica from weathered pegmatite ores. This report summarizes the process development work, and demonstrates the feasibility of producing commercial-grade concentrates with a high recovery of the contained mica. Mica mining in the United States is confined largely to pegmatites and schists in a few well-defined areas. The largest region extends from central Virginia southward through western North and South Carolina, northern Georgia, and east-central Alabama.'" Whereas in western North Carolina, these pegmatite dikes are weathered to only a shallow depth, they are deeply weathered in the Carolina Piedmont and in Georgia and Alabama. Particularly in Georgia and Alabama, mica is found in the decomposed slates and schists adjacent to the pegmatites, and in deeply weathered granites. Most flake mica is obtained as the major marketable product by crushing and milling pegmatites4 To a lesser extent, mica is produced as a byproduct of feldspar" and spodumene. Past research by U.S. Bureau of Mines (USBM) has developed two methods for mica flotation. One method utilizes acid-cationic flotation for recovery of mica.' This method requires thorough desliming of the ore with consequent fine mica losses. The other utilizes an alkaline-anionic-cationic process for recovering fine-size mica from pegmatite ores after desliming sufficiently to remove clay materials but not so drastically as to remove the fine mica.8-18 The latter process has also been successfully applied to a micaceous schist ore from California.' This report summarizes the results of these research studies. Both of the flotation processes developed by USBM have been adopted by the mica industry and are being commercially utilized separately and in combination for flotation of mica from weathered pegmatite ores. Description of Ores The weathered mica pegmatite ores used in the investigations were obtained from commercial mica operations in Alabama, Georgia, and North Carolina. The

Jan 1, 1971

By P. H. Wendland

Calculations are presented for the design of a silicon photodiode in which the incident light beam makes multiple passes between the detector surfaces. Total internal reflection is used for this "light-trapping" effect. By this means, the optical path length can be extended to several millimeters, while the electrode separation remains less than 102 cm, as required for nanosecond response time. Data are presented for a Schottky barrier photodiode constructed on a multiply reflecting silicon base wafer. It is shown that the long-wavelength response is considerably extended in such structures without a corresponding sacrifice in high-speed response. The development of efficient and powerful lasers at 1.06 p has stimulated interest in detectors which operate at this wavelength. In typical silicon photodiodes, for detecting 1.06 p radiation, the requirements for high speed and high sensitivity are mutually exclusive. Since the absorption coefficient is only 25 cm-', a lo-'-cm path length is required to absorb 92 pct of the incident 1.06 p radiation. If the electrode separation is greater than 10 cm, however, the carrier transit time will be greater than 1 nsec. This problem can be solved by allowing the incident light beam to make multiple passes between the electrodes. The optical path length can then be extended to several millimeters, as required for complete absorption, while the electrode separation remains less than 10' cm, as required for nanosecond response time. In a typical photodiode geometry, one ohmic contact and one rectifying contact are formed on the two opposite surfaces of a base wafer, and the wafer thickness determines the electrode separation. The objective of the multiple reflection design is to allow all 1.06 p radiation to enter the detector front surface and to form the back detector surface so that no 1.06 radiation can exit. Total internal reflection at the back detector surface is well-suited for light trapping of 1.06 p radiation because the relatively large dielectric constant of silicon leads to a critical angle of 16.5 deg for total internal reflection. LIGHT TRAPPING It is well-known that, as light passes from one medium such as air into another medium such as glass or silicon, the angle of refraction is always less than the angle of incidence. In the limiting case, where the incident rays approach an angle of 90 deg with the normal, the refracted rays approach a fixed angle +, beyond which no refraction is possible: this is called the critical angle. It follows from Snell's law that where = critical angle, n - index of refraction of air, n' - index of refraction of the medium. Applying the principle of reversibility of light rays, all internal angles of incidence greater than +, will produce total internal reflection and "light trapping". The index of refraction of silicon at 1.06 p is 3.5,' and the critical angle is thus 16.5 deg. Fig. 1 shows these relationships for silicon. This very small critical angle in silicon is significant because all incident angles between 16.5 and 90 deg will produce total internal reflection and "light trapping". This effect can be implemented with a "prismlike" geometry, so that incident light can be introduced into the sample without loss and "trapped". PHOTOSIGNALS A precise knowledge of the absorption coefficient at 1.06 in silicon is of critical importance to the design of fast and efficient silicon photodiodes for 1.06 radiation. Dash and newman2 show a value of 25 cm-l, and our measurements have corroborated this value. Assuming that the collection of photoinduced minority carriers is perfect, the quantum efficiency of a photodiode is dependent only on the absorption coefficient. It then follows from Lambert's law that where QE is the quantum efficiency in pct, a is the absorption coefficient, d is the optical path length, and the reflectivity at the surface is assumed to be completely suppressed by an optical interference layer. Fig. 2 gives the maximum quantum efficiency for 1.06p radiation of a silicon photodiode with optical path length d, using Eq. [2]. The ultimate response time of a fully depleted photodiode to an incident light pulse can be considered to be the arrival times of all photoinduced carriers at the contacts, i.e., the minority carriers at the junction interface and the majority carriers at the oppo-

Jan 1, 1968

By Harrison Everett Ashley

(Pittsburg Meeting, March, 1910.) Slimes are usually defined as all material passing a certain-sized sieve, which is invariably the finest sieve employed by each metallurgist in his tests; 100-mesh and 200-mesh have been taken as the limits by different writers in recent volumes of the Transactions ; 150-mesh is common practice. The reader of the first two volumes of Richards 1 is impressed by the inflexibility with which all grinding-arrangements are praised or condemned according as they make little or much slime. Until recently, one of the principal aims of the mill-man was to keep down the slimes. He was obliged to have separate apparatus for sand- and slime-treatment. The latter was sometimes puzzling and erratic. Often it was found most .economical to omit treatment of the slimes entirely, or to store them in immense basins or heaps, waiting for possible future improvements in methods of treatment. Recently, improved processes of agitation and filtration have made possible the successful extraction of values from many slimes. The tendency, when the metal is exceedingly finely disseminated, is to cut out all sand-treatment and to slime everything. These processes have been almost entirely mechanical in their conception, and it is the object of the present paper to show that certain chemical considerations may greatly hinder or assist their successful operation. It has long been a source of dissatisfaction to chemists that their analyses fail to be of much help to the clay-worker; that it is a matter of common occurrence for clays of closely identical ultimate chemical constitution to show the widest possible divergence in physical properties. * Published by permission of the Director of the U. S. Geological Survey. Ore Dressing (1903).

Aug 1, 1910

By B. J. Reid, J. N. Greenwood

It has been shown1 that cavities are formed in the grain-boundaries of copper and 70:30 brass (as well as in magnesium) by the application of tensile stresses at elevated temperatures. For a given rate of strain, cavities are formed in brass at a lower temperature than in copper, and they are larger; it was thought that this earlier onset of cavitation in the alloy might be caused by the distortion produced in the metal lattice by the alloying element. If the cavitation is caused by the vacancy-migration mechanism, ' the distortion of the lattice would reduce the energy required for vacancy formation, and thus provide vacancies more readily to assist in the cavitation process. It was decided to study a range of alloys to find out whether increasing amounts of an alloying element would cause any progressive increase in, or modification of, the cavitation in copper. Nickel was chosen for this purpose, since it is soluble in copper throughout the whole range of alloys. A further advantage was its low vapor pressure, which enabled the preparation of the alloys in vacuum and the consequent reduction of their hydrogen content. It has been suggested that cavities may be nucleated by the migration of hydrogen atoms to the less regular atomic spacing at the boundary, and their combination to form molecules. The latter being too large to diffuse back into the lattice, form stable cavities into which vacancies could condense. PREPARATION OF ALLOYS All alloys, ranging from pure OFHC copper to 43.7 at. pct of nickel (carbonyl), were prepared in a vacuum furnace employing high-frequency induction heating. When a charge was melted, vigorous gas evolution occurred initially, after which the pressure was reduced to 10-4 mm Hg, using an oil-diffusion pump. After remaining molten for a further 1/2 to 1 hr to aid degassing, the melt was allowed to solidify in the crucible. In this way contamination by gases evolved from the surface of a separate mold was avoided. The list of alloys prepared is shown in Table I. The percentage of hydrogen present in the OFHC copper was reduced from 0.0054 at. pct to 0.0012 at. pct, as measured by a vacuum-fusion method. The alloys were cold-rolled to 0.03-in.-thick strip, and were homogenized in vacuum for 1 hr at 900°C after every 30 pct reduction. At a thickness of 0.05 in., the strip was held at 450°C for 24 hr under high vacuum (10-5 mm Hg) in an attempt to reduce the hydrogen content still further. Tensile specimens were prepared with a gage length of 1 in., and were finally annealed for 2 hr at 800°C to give an average grain diameter of 0.1 mm. TENSILE DEFORMATION The specimens were deformed in a tensile straining machine operated by a synchronous motor at a constant rate, equivalent to 1.7 pct per hr on a 1-in. gage length. The specimen and grips were contained in a vacuum jacket. A pressure of 10-2 mm Hg was maintained, which was sufficient to prevent surface tarnishing of a specimen held at 650°C for 24 hr. Tests were made between 20. and 650°C and the temperature was controlled to within ±3° over the gage length. The load was measured by the deflection of a cantilever spring. a) Mechanical Properties-The ductility curves shown in Fig. 1 fall into two groups. Figs. 1(a) and (b) for copper and the 5 pct nickel alloy, respectively, show a continuous fall in ductility with a steep fall between 200° and 400°C. Figs. l(c) to (f) show the same sudden loss of ductility between 200° and 400°C. With higher temperatures, however, there is a recovery in ductility values to a maxi-

Jan 1, 1959

By L. Kaufman

An equation is derived relating AF a", the difference in free energy between austenite and martensite, to temperature and composition in the iron-chrmnium and iron-chromium -nickel systems. This equation is used to calculate at M,for low-carbon stainless steels. A method for utilizing the energetic equations to calculate M, in stainless steels is suggested. In addition, an explanation is given for the anomalous effect of chromium in forming a y loop concurrent with the lowering of M, in iron-base alloys. It has been demonstrated that the effect of alloying elements on the Ms temperature is strongly dependent upon the influence of the alloying element on the thermodynamic properties of the y and a phases. Notwithstanding the fact that alloying elements affect the Ms temperature by changing the lattice parameters, elastic constants, and inter-facial energies of the a and y phases, their primary role is a thermodynamic one. Thermodynamic analyses have been presented for the iron-carbon'-4 and iron-nickele systems in which equations are derived for Fa and FY, the free energies of the b.c.c. (or b.c.t.) and f.c.c. phases. These equations, which can be obtained with the aid of the equilibrium diagrams, can be used to help predict the Ms vs composition curves. This paper deals with a calculation of the thermodynamic properties necessary to compute AFffLY for iron-chromium alloys. In view of the lack of direct information concerning the activity of chromium in the a and y phases of iron, the regular solution approximationa,g was utilized to bridge the gap. After the necessary equations were derived for the iron-chromium system they were used to calculate AFff— Y for stainless steels by combining them with the previous results obtained for iron-nickel alloys. The results of these calculations offer a means of rationalizing the anomalous behavior of chromium in iron-base alloys. The latter depends on the regular solution approximation in contrast to the explanation offered by zener , which is based upon the "magnetic specific heat effect." THERMODYNAMICS OF THE IRON-CHROMIUM SYSTEM If y represents the atomic fraction of chromium in an iron-chromium alloy, then the free energy of the b.c.c. phase Fa, the f.c.c. phase FY, and the difference in free energy between the f.c.c. and b.c.c. phases of the same composition are given through [3] and Fmy represent the free energies of mixing of the a and y phases; while @era -' and uFea 4Y represent the difference in free energy between f.c.c. and b.c.c. chromium and iron respectively. The values of AEFe0ff4Y in the temperature range between 0" and 900K ( .

Jan 1, 1960

By C. E. Brukl, E. Rudy, St. Windisch

The ternary-alloy system Nb-Mo-C was investigated by means of X-ray, melting point, DTA, and metallo-graphic techniques; a complete phase diagram for temperatures above 1500°C was established. Above 1960°C, niobium monocarbide and the cubic (Bl) high-tenzperature phase in the Mo-C system form an uninterrupted series of solid solutions. The ternary range of the pseudocubic q MoCl-, is very restricted. Dimolyb-denum carbide dissolves up to 44 mol pct Nb2C (2240°C), whereas the maximum solid solubility of MO2C in Nb2C does not exceed 5 mol pct. The order-disorder transformation temperatures in Mo2C and Nb2C are lowered by the mutual metal exchanges. Six invariant (p = const) reactzons occur in the ternary system; three correspond to class 11-type four-phase reactions involving a liquid phase, one to a class I (eutectoid)-type, and two further isotherms are associated with limiting tie lines. The results of the Phase diagram investigation are discussed, and the thermodynamic interpretation identifies the low relative stability of the binary sub-carbides in conjunction with the large stability diflerences between niobium and molybdenum carbides as the cause for the formation of a stable equilibrium between the monocarbide and the metal phase in the ternary reson. Due to their refractoriness, the carbides of the high-melting transition metals have received increased interest in recent years as base materials in composite structures for aerospace applications at high temperatures and for the development of self-bonded cutting tool materials; other novel fields of application include power reactors, where operation at high temperatures becomes essential for attaining high power efficiencies. In these applications, the increased reaction rates at high temperatures require a close consideration of the chemical interactions between the alloy constituents. As a consequence, a detailed knowledge of the phase relationships in the alloy systems is required in order to provide a sound basis for developmental -type work. Partly as a result of the considerable experimental difficulties associated with the investigation of this high-melting alloy class, no complete studies of ternary metal-carbon systems have been performed until recently. Even the high-temperature phase relationships in the binary transition metal-carbon systems have been delineated only during the past few years to a degree of accuracy required for a more detailed study of ternary or higher-order alloys. In recent investigations of binary and ternary systems of refractory transition metals with carbon, boron, and silicon,&apos; alloys from the ternary systems Nb-Mo-C became of interest because of the demonstrated possibility2,3 of obtaining compatible composites based on metal + monocarbide combinations. In the meantime, however, studies in other, but related, ternary metal-carbon systems, such as Ta-W-C, have shown that the solid-state equilibria may change significantly toward higher temperatures (>2000°C), and that extrapolations based on low-temperature equilibrium data are, in general, not very reliable. Although the lower-temperature (<2000°C) phase relationships in the Nb-Mo-C system are similar to those found in Ta-W-C, a cursory thermodynamic analysis of the equilibria indicated4 that complete solid-solution formation between Mo2C and Nb2C should not occur at higher temperatures. The present work was conducted in order to experimentally verify these expectations and, in addition, to provide phase equilibrium data in the melting range of the alloys. In the boundary systems, niobium and molybdenum are known to form a continuous series of solid solutions.576 The continuous solubility was also confirmed by Kornilov and Polyakova,7,8 who also observed a minimum melting point at 22 at. pct Mo and 2345°C. The phase diagram investigations of the Nb-C system by Storms and Krikorian9 and Kimura and Sasaki10 were recently supplemented by Rudy et al.11,12 The system contains a high-melting monocarbide with the B1 structure, Table I, and a subcarbide, Nb2C, which exists in at least two different states of sublattice order at low temperatures47"-&apos;3 and a disordered state above approximately 2500°C.11,12 The melting-point measurements by Rudy et al .11,14 are in close confirmation of the data by Kimura and Sasaki.10 The rather complex phase relationships in the Mo-C system were only recently Clarified.15,18 The system is characterized by three congruently melting, intermediate phases, MozC, ? MoC1-x and a Mol-,, Table I, of which only Mo2C is stable at temperatures below 1650°C. Substoichiometric MozC exists in several states of sublattice order which interconvert in homogeneous phase transitions. Hyperstoichiometric compositions cannot exist in the ordered state. Upon cooling through a critical temperature range, the

Jan 1, 1968

By P. Datta, L. L. Handy

For a wetting phase displacing a nonwetting phase from a porous medium the distribution of the residual fluid may depend on displacement conditions. Although this subject has been debated in the literature, only a few experiments have been cited to support the various conclusions. Experimental results presented in this paper show that fluid distributions are dependent on imbibition pmcedures. Results agree qualitatively with predictions from the pore doublet model. If the rate of water imbibition is restricted, the nonwetting phase is trapped preferentially in .the larger pores as expected. But if the rate of water imbibition is unrestricted, trapping occurs somewhat more in the smaller pores. These conclusions were deduced primarily from relative permeability measurements. INTRODUCTION Relative permeabilities are known to depend on the saturation history of the porous medium. For either continuously increasing or continuously decreasing wetting phase saturations, however, they have normally been assumed to be single-valued functions of saturation. Severat investigators have compared relative permeabilities measured by different methods and have reported acceptable agreement.1-5 Studies of simple models of the displacement process suggest that the distribution of residual fluids can be influenced by the displacement method. If this is so, relative permeabilities also should depend on the method of displacement. These predictions have not been supported by experimental data. The object of this paper is to investigate experimentally some of the predictions of these model studies. The particular question of importance is whether steady-state and unsteady-state methods should be expected to give the same values of relative permeability. Much of the reported data showing agreement between steady-state and unsteady-state methods have been obtained on un consolidated sand. Johnson, Bossler and Naumann, however, compared steady-state and unsteady-state results for Weiler sandstone and found no significant differences.5 Levine made a detailed study of pressure and saturation distributions for a laboratory waterflood in a consolidated system but made no attempt to calculate performance from steady-state relative permeability data Bail and Marsden in a somewhat similar set of experiments did attempt a comparison of observations with predictions but the results were inconclusive.7 Excluding gravity, two forces affect the distribution of wetting and nonwetting phases during an unsteady-state, immiscible displacement: viscous and capillary forces. If relative permeabilities are indeed the same when measured by steady- or unsteady-state methods, the fluid distributions at the same fluid saturations must be similar for both methods. Furthermore, the implication is that the relation between capillary and viscous forces is the same for both methods insofar as the effect on microscopic fluid distributions is concerned. Unsteady-state displacements are normally run at high rates to maximize the ratio of the viscous to capillary forces. The objective is to reduce the effect of capillary forces on macroscopic fluid distributions. For floods in which the phase which wets the porous medium displaces the nonwetting phase, capillary forces compete with viscous forces in determining fluid distributions.8 The residual nonwetting phase is discontinuous. Competitive aspects of viscous and capillary forces and the resulting effect on water-oil distribution in porous media have been illustrated in an elementary way using a pore doublet model. This model and predictions from it have been discussed at considerable length by Rose and Witherspoon, Rose and Cleary and by Moore and Slobod.9-13 Rose and Witherspoon show that so long as water invades the model at a rate equal to or greater than the free imbibition rate, the water will move through the larger capillary of a doublet first and trap oil in the smaller capillary. On the other hand, if a negative pressure is applied to the water

Jan 1, 1967

By H. W. Mead, C. E. Birchenall

The rates of growth of graphite nodules in cast irons are calculated for a model of a growing graphite sphere surrounded by a shell of austenite through which carbon and iron are diffusing. The carbon is supplied by a cementite-austenite mixture at the outer edge of the shell. The calculated rates from carbon diffusion agree fairly well with the measured values, but the rates based on iron self-diffusion are much too low. Because the iron must be removed for the process to continue, o plastic deformation mechanism is suggested. THERE has been an unfortunate tendency in the literature of growth and dissolution reactions to identify the slope of a plot of the logarithm of the growth or dissolution rate vs reciprocal of temperature with the activation energies determined by diffusion measurements in chemically similar systems. If the solubility limits at the growing or dissolving phase interfaces change appreciably with temperature, such procedures may lead to serious errors in specifying the mechanisms controlling the rates of these processes. A reaction involving the simultaneous growth of graphite from and dissolution of cementite in austenite will be discussed here to show that substantial oversimplifications have been made in earlier studies of this process. From a dilatometric study of the overall y-range graphitization process, Bunin and Danil&apos;chenko&apos; found an apparent activation energy of 90 kcal per g-atom. They concluded that, since this value was much nearer the approximately 68 kcal per g-atom observed for self-diffusion in pure iron than it was to 32 kcal per g-atom for diffusion of carbon in iron, the rate-controlling process must be the diffusion of iron away from the growth interface to provide space for graphite to grow. Bunin and Salli&apos; bolstered this agreement with the observation that graphite forms preferentially at surfaces, cracks, and pores where free volume is present initially. Ziegler, Meinhart, and DeaconQ ad found an apparent activation energy for graphitization of 98 kcal per g-atom. This value, like that of Bunin and Danil&apos;chenko, applies to nucleation and growth occurring simultaneously. It is not generally assumed that nucleation is entirely diffusion-controlled. Treatments such as prior cold-working and prequenching, which strongly affect the rate of graphitization, would be expected to have little effect on diffusion processes occurring at the graphitization temperature. In any case 90 kcal per g-atom is not particularly close to the 54 kcal or lower to be expected for diffusion of iron which is carbon-saturated.&apos; The presence of silicon in iron is said to reduce the activation energy for diffusion,&apos; so its presence in cast irons would be expected to have a similar effect. The growth rates of nearly spherical graphite nodules have been studied by Brown and Hawkes." The initial rate of growth before impingement of adjacent carbon-depleted areas was found to be parabolic, and rate constants were measured at several temperatures. Their apparent activation energy for the growth process only was roughly 45 kcal per mol. Burke and Owen&apos; made growth-rate measurements on several alloys, recording appreciable induction periods, initial nucleation rates, and average radial growth rates based on an assumed linear rate of radial growth. They found that increased silicon content increases the ovel-all reaction rate and the initial rate of nucleation, and decreases the induction period. Assuming 1) spherical nodules. 2) constant volume rate of nucleation, 3) austenite in equilibrium with cementite at positions where the cementite is not depleted, and 4) a constant diffusion distance during the growth of a nodule, they calculated an apparent overall activation energy of 68 kcal per mol, which was related to an activation energy for nucleation of 79 kcal per mol and activation energy

Jan 1, 1957

By D. C. Reynolds

Recent experiments with II-VI compounds have shown that they hazle considerable potential for laser applications over a broad region of the optical spectrum. It may be possible to cover the spectrum continuously from 3200A (ZnS) to the far infrared (CdHg:Te) since HgTe is a semimetal. At this writing laser action has been observed in ZnS, Zn0, CdS, CdSe, CdS:Se, CdTe, and some of the CdHg:Te alloys. Of particular interest are those lasers operating in the zlisible and near ultraciolet regzons of the spectrum where detectors of high sensitivity are available. The lasing transitions in II-VI compounds are bound exciton transitions some of which have been identified in auxiliary experiments. High efficiencies and low thresholds for lasing hare been achieved almost exc1usively in plutelet-type crystals. The greater crystalline quality exhibited by the phtelet-type material is shown to result from the crystal growth habit. Phonon scattering- of conduction electrons to the ground-state exciton is discussed ill relution to Lou thresholds and high efficiencies for lasing- observed in the CdS:Se solid solutions. The first successful semiconductor laser operation was achieved in the III-V compounds. It is possible to choose a material in this group that will operate between approximately 0.65 and 8.5 . There are at least two reasons why one would like to have a laser operating at shorter wavelengths. First, it would be easier to experiment with a laser operating in the visible region of the spectrum, and also more desirable to have high-in tensity visible light sources. Second, the most sensitive photomul-tiplier detectors are available in the visible and near ultraviolet regions of the spectrum. It is known that II-VI compounds are direct-band-gap semiconductors and as such offer the potential of operating at any specified wavelength between 3200 (ZnS) and 7772A (CdTe). Light emission from II-VI compounds has been the subject of numerous investigations for many years. These investigations were all primarily concerned with noncoherent emission. It has been only recently that coherent emission from these compounds has been observed. To date, laser operation has been demonstrated in CdS, CdSe, and the solid solutions of CdS:Se, ZnS, ZnO, and CdTe. These compounds cover an appreciable portion of the optical spectrum from the ultraviolet to the near infrared. In considering laser applications, the use of lasers in communication&apos;s systems offers many desirable features. In any operation of this type one must consider the losses in transmitting the radiation from the source to the detector. Atmospheric absorption in the visible and near ultraviolet is variable and greater than in certain regions of the infrared. It might be concluded that for long-range communication systems an infrared laser operating in a spectral region that is coincident with a transmission window in the atmosphere would be preferable. However, one cannot overlook the possibility of operating a system in the sensitive region of a highly sensitive photomultiplier detector or other light-amplifying system. LASER CONSIDERATIONS To produce a source of coherent radiation it is necessary to achieve a population inversion. In the case of semiconducting materials it is necessary to raise the electrons from one energy state to a higher-energy state relative to it. In semiconductors, this population inversion can be achieved by three different techniques: 1) Current Injection. This technique uses a p-n junction biased in the forward direction. Large numbers of electrons are injected from the n region into the p region, and recombination occurs close to the junction. An inverted population is obtained in this region and the recombination radiation propagates parallel to the junction. This type of pumping has been used in the GaAs junction-type lasers but has not been successfully employed in the II-VI compounds. 2) Optical Pumping. In this case, one uses photons to obtain a population inversion by exciting electrons to higher-energy states. The pump sources are flash lamps or arc lamps and, occasionally, other laser sources when such sources have the appropriate energy for exciting the electrons. The disadvantage of this type of pumping is that flash lamps put out a rather broad spectrum of radiation, whereas the laser material has a rather narrow region of absorption. This results in an inefficient process. Laser sources provide efficient pump sources but the number of usable wavelengths is limited. 3) Electron Beam Pumping. In this technique, the laser sample cavity is bombarded with electrons having energies in the range from approximately 10 to a few hundred kv. The bombarding radiation excites electrons from valence to conduction band states in the semiconductor, giving rise to an inverted population. This type of pumping has been used successfully in several 11-VI compounds.

Jan 1, 1968

By R. D. Burlingame, T. L. Joseph, Gust Bitsianes

DESPITE almost fifty years of commercial practice, the sintering of iron ore has received little fundamental study. Much of the theoretical work1-&apos;has dealt with the constitution of sinter produced under widely varying conditions. While these studies have broadened our knowledge of the changes that occur in the sintering zone and in the freshly formed sinter during the early stages of cooling, they provide little insight into the changes that precede the formation of sinter. These preliminary changes merit study as a part of the overall process. Hessle. working with beds of Swedish magnetite concentrates, was one of the first investigators to study the sintering process in its entirety. On the basis of temperatures observed at various levels of the bed during sintering, he postulated a number of distinct reaction zones to account for the chemical changes leading to the formation of sinter. A more direct method of attack is that of arresting the sintering zone after it has progressed part way through the bed. A study of a vertical cross section through such a quenched bed provides direct information on the changes taking place at various levels. This method was used by McBriar et al.&apos; to show that several well-defined zones of chemical change existed within beds that were typical of British sintering practice. The same general method of attack was developed independently in the present investigation to study partially sintered beds typical of American practice. Experimental Sintering Equipment The sintering operation was carried out on an experimental scale with the equipment shown in Fig. 1. The refractory-walled sintering chamber A was 11 in. deep and averaged 9 in. in diameter. Air was introduced through a tapered flow section B, which contained the orifice C for accurate metering of the incoming air. This section was located directly above the square ignition housing D, which in turn rested upon the sintering chamber A. The bed was ignited with burner E. The required suction for the operation was furnished by a fan F, which had an air capacity of 500 cfm (stp). Hot exhaust gases from the sintering chamber were cleaned in the dustcatcher G before entering the exhaust fan. In the study of partially sintered beds, it was essential to find some technique for removing the entire charge from the sintering pot without disarranging the unsintered bottom portion. This problem was finally solved by sintering the charge in a removable basket, which snugly fitted the sintering chamber. This basket was constructed of two thicknesses of window screen and was lined with a 3/16-in. layer of asbestos paper. The bottom of the basket consisted of two thicknesses of wire screen, which were fastened to the basket wall. For high fuel mixtures, additional insulation was provided by a somewhat thicker layer of asbestos cement. Preparation of Partially Sintered Mixtures The moist feed was carefully placed in the sintering basket, to prevent segregation of the particles, which varied widely in size and composition. A thermocouple was placed in the center of the basket with the hot junction halfway down, and the mixture was evenly distributed around it. During ignition and throughout the sintering of the upper half of the bed, the hot junction temperature increased very little. When the sintering zone reached the halfway point, as indicated by the sudden increase in the hot junction temperature, the charge was quenched. During quenching the suction was turned off and the orifice was tightly stoppered to prevent further influx of air. At the same time, nitrogen was admitted to the sintering chamber through the orifice tap. As soon as the nitrogen had displaced the air and products of combustion, the charge was removed from the sintering pot for immediate dissection. It is impossible to preserve the exact zone structure of the bed at the instant that combustion is arrested unless the downward transmission of heat is also immediately stopped. Fortunately, heat transfer is very slow in beds containing a stationary fluid, especially if the particle size is small. It follows that the minimum quantity of nitrogen should be used to displace the air and that static conditions be established as soon as possible. A very steep temperature gradient across the combustion zone for some time after the quench was evidence of in-

Jan 1, 1957

By T. W. Guy

IN approaching the problems of sampling coal, a brief statement of certain facts that are more or less taken for granted may be helpful: 1. The coal-mine operator needs reliable data as to the physical and chemical characteristics of his coal in order to plan and control intelligently its mining, preparation, and marketing. The consumer needs such data in order to select the equipment and the coal best suited to meet his requirements, to check the quality of fuel delivered, and the efficiency of its use. Various government control agencies require such data as a basis for classifications and establishment of limits, for use in carrying out the requirements of federal, state or municipal laws. 91. All such data must be obtained by analyses, measurements, or tests made on samples, which are taken to represent the coal tested. It is now becoming generally recognized that, in accordance with statistical laws, test results from samples are not expected to represent the exact qualities of the lot or lots from which they are taken, except within statistical, or "expected limits. Such limits are definitely specified by some present standards of the American Society for Testing Materials and the British Standards Institution. 3. Samples that give ash within satisfactory limits are usually assumed to be sufficiently accurate for other analyses, except under special condi¬tions such as meeting close specifications for coals showing wide variability in sulphur, fusibility of ash, or other quality. 4. In modern coal sampling the increment is recognized as the basis for taking samples. For a given accuracy the number and weight of increments specified are increased with: (a) the expected ash, and (b) the top size designation of the lot to be sampled. 5. It is extremely important that every sampler recognize that: the specified number and weight of increments are minimums, and that the precautions to be used in order to take them properly are also minimums; that in the reduction of the gross sample the sizes specified for the coal after each stage of the crushing are maximum; that the quantities specified to be taken after each crushing stage are minimum; that any short cuts or failures to meet the requirements specified will reduce the statistical accuracy of the results, although it is possible that one or two checks on a few results may seem to show otherwise.

Jan 1, 1943

By N. Albon, J. F. Miller, R. W. Coutant

The deposition of nickel from nickel carbonyl onto amorphous substrates has been studied, with attention being paid to the specific effects of the physical and chemical nature of the substrate and to the effect of oxygen on the crystal morphology. With substrates having microscratches or etch lines, nu-cleation was observed to occur preferentially at these surface features. Flame polishing of the substrate surfaces was found to remove any regularity in preferred nucleation sites. Nucleation was observed to occur less readily on substrates having relatively low surface energies, e.g., Teflon and carbon, than on materials having relatively high surface energies such as fused silica. At low oxygen concentrations [Ni(CO)4/O2 > 54 x 1031, apparently pure and well-formed nickel crystals were observed to nucleate and grow irrespective of the presence of oxygen. At higher concentrations of oxygen, a secondary growth process, dependent upon oxygen concentration and resulting in a gradual change in crystal morphology, was observed. It is concluded that the influence of oxygen on the nucleation and growth of nickel is determined by the relative kinetics of the various possible reactions, and that under the experimental conditions employed the 0-Ni interaction occurs at some nucleation or prenuclea-tion stage. THE purpose of the research on crystal nucleation in chemical vapor deposition as described in this paper has been twofold: 1) to examine some effects arising from the use of amorphous substrates, and 2) to study the specific effects of oxygen on the nucleation and growth of nickel in a nickel carbonyl-carbon monoxide system. The nickel-nickel carbonyl system was chosen for study for several reasons: 1) nickel films are of potential interest in electronic technology, and 2) the deposition can be carried out at moderate temperatures, thus making possible detailed study of the processes occurring and the use of a variety of substrate materials (including some which melt at relatively low temperatures). Fairly well-characterized chemical kinetics and a broad background of information are available on the py-rolysis reactions and on the mechanisms of Ni-O interactions. NUCLEATION AND GROWTH THEORY The theoretical background for the analysis of nucleation and growth by vapor deposition for a one-component system (physical vapor deposition) and the problems encountered therein have been reviewed by Hirth and pound.&apos; More recently Mande12 has indicated the difficulties of applying such theory to chemical systems. One of the major problems involved is the possible multiplicity of chemical reactions and deposition routes. For example, those reactions which are probably of significance in a discussion of the pyrolytic deposition of nickel from nickel carbonyl under the conditions of the present experimental work are illustrated in Fig. 1. As is true with most cases of chemical deposition, many of the quantitative data necessary for a complete analysis of the system are either not available or not directly measurable. However, with the present system, sufficient data are available for a discussion of nucleation and growth on a qualitative basis. Classically, the rate of nucleation of a monatomic species being deposited from a vapor phase is given by

Jan 1, 1965

By A. E. Jenkins, L. A. Baker, N. A. Warner

The electromagnetic levitation technique has been successfully applied to rate studies of the de-carburization of liquid Fe-C alloys from 5.5 to zero pct C at 1660°C using gas mixtures containing 1 to 100 pct CO2 with CO or helium as diluents. The rate was shown to be independent of carbon concen -tration in the melt and only slightly affected by total pressure. The results were shown to be consistent with control either wholly or predominantly by gaseous diffusion. As well as avoiding side reactions with the crucible, levitation allowed precise calculation of reaction rates since the geometry of the system was accurately known. It also enabled interpretation using the established mass-transfer correlations for a gas flowing past a single sphere, The existence of gaseous-diffusion control and the extremely high rates of reaction reported for gas flowing past small spheres are of great practical significance in relation to current investigations into the use of spray-refining techniques in continzcous steelmaking. PREVIOUS work on the kinetics of carbon removal from liquid iron has generally involved the use of oxygen as the gaseous oxidant. Investigations have been made on commercial steelmaking furnaces and in addition laboratory-scale experiments have been attempted with a slag phase present. The general lack of agreement and the obvious difficulty in interpreting the results is at least partly due to the complexity of the systems used. This work has been reviewed by carter1 and more recently by Bods-worth2 and ward.3 The gas-liquid metal reactions have also been studied with slag-free melts by the conventional technique of blowing oxygen onto the melt contained in a crucible. The reaction between the oxidant and carbon could be followed by sampling either the melt or the gas phase. Side reactions with the crucible material have precluded direct interpretation of the results. Fujii, 5-7 using O-A mixtures, attempted to correct for the crucible reaction and reported that the reaction between carban and oxygen occurred at the melt surface and was controlled by gaseous diffusion for carbon contents greater than a critical value and by carbon-diffusion control in the metal for carbon contents less than the critical. The critical value was given as approximately 0.15 pct C. However, whereas this was the case for oxygen contents in the gas stream up to 10 pct, between 10 and 20 pct the rate was found to be approximately independent of oxygen content, and Fujii attempted to explain this observation by postulating the presence of an oxide "film".7 Turkdogan et a1.&apos; have studied iron fuming during decarburizing and suggested gaseous-diffusion control of the reaction above a critical carbon concentration, given as approximately 2 pct C. Filippov9-12 reports in a similar manner to Fujii but makes no mention of the crucible reaction with the carbon in the melt, even though this would be expected with the temperature, crucible material (MgO), and carbon contents used. Both CO2 and O2 were used as oxidant gases. In the present work the decarburization of liquid Fe-C alloys was studied using an adaptation of the electromagnetic levitation technique, and thus the complexities introduced by crucible side reactions were eliminated. In this technique a single sphere of the Fe-C alloy is freely suspended and heated by electromagnetic induction in a flowing stream containing the oxidant gas. The full range of carbon levels up to saturation was studied, the temperature of the melts in this instance being held constant at 1660°C. The oxidant was CO2 which was used either pure or in association with a diluent gas, the total gas pressure being maintained either at 1 atm or at a reduced pressure. Thus the reaction under study was:

Jan 1, 1964

By John Drew Ridge

For the early conservationists before the first decade of the twentieth century, conservation meant largely the planned preservation of water, forest, soil, and wildlife resources-renewable natural resources-through efficient use. Water was primarily thought of as a resource that must be kept available for humans to enjoy in a recreational sense and as a habitat for wild animals, birds, and fish. Little thought was given to conserving water for actual human or industrial needs. If industry entered the picture it was as an actual or potential polluter of water, which was the basic necessity for drinking and bathing and for outdoor recreation. Even the forests were only held worthy of conservation largely because they were areas in which wildlife could, or should, abound and in which human beings might find relief from the cares and problems of everyday life. Finally, the soil was to be preserved in its natural fertility through proper fertilization, cultivation, and crop rotation; in this the conservation program did enter the realm of the practical. For through soil conservation, the world was to be able to raise increasingly large crops of constantly improved organic foods which the steadily growing population of the earth was to need to survive and improve. Relatively little time, however, was spent on considering what profit the farmers of the world were to make from the scientific farming on which they were to embark. In short, the appeal of the early conservationist program was much the same as that of the entire era of Theodore Roosevelt of which it was a part; it was to the heart rather than to the mind; it was to emotion rather than to reason. Because of the way in which the problem was approached, conservation took on the character of a crusade rather than that of a practical plan for the betterment of man through economically sound methods of improving his environment and the uses he makes of it. Conservation had, then, all the advantages of a stirring appeal to the emotions, but it did not have the lasting results that would have been achieved by a program with a sound economic basis.

Jan 1, 1959

By J. Szekely, V. Stanek

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

Jan 1, 1970

By Douglas W. Fuerstenau, Kalanadh V. S. Sastry

A mathematical model is developed for flotation in a countercur-rent column where continuously generated air bubbles rise through a downward flowing pulp. The model is based on the assumption of axially dispressed plug flow of both the liquid and air-bubble phases. An appropriate flotation-rate expression, which takes into account the changing concentration of solid particles on the bubble surface, is formulated and applied to the column flotation process. Although a complete solution of the resulting equations is not possible for the general case, solutions for specific limiting cases have been obtained. From these, it is possible to predict theoretically the dependence of the column performance on certain variables. Specifically, concentration profiles as a function of height in the column have been calculated for various values of the characteristic dimensionless parameters. The effect of axial dispersion on the perforance of the equipment is pointed out. In recent years, there has been considerable progress in delineating the fundamental physicochemical principles underlying the unit operation of flotation and in advancing methods for kinetic analysis and process modeling of flotation.&apos; In general, nearly all such process research has been concerned with the analysis of conventional flotation equipment and circuits. From these investigations, the possibility of developing and designing flotation apparatus for more efficient separation appears limited with conventional equipment. How- ever, recently there has been an attempt to apply coun-tercurrent flow techniques to flotation in the context of a flotation column.- some of the results from pilot-plant studies2 indicate possibilities for success of this process at least in obtaining higher-grade concentrates though not higher recoveries. A flotation column operates on countercurrent principle in which bubbles continuously generated at the bottom of the column rise through the downward flowing slurry. A schematic diagram of the flotation column is shown in Fig. 1. The physicochemical principles of bubble-mineral attachment, using this flotation column, are essentially the same as those which prevail in conventional equipment, the only differences being in the flow patterns of the various phases. The salient operating feature of the column as compared to those of a conventional cell is that neither mechanical nor turbulent pneumatic agitation is employed. The column is

Jan 1, 1971

By George A. Colligan, John M. Lo

Electrical resistivity as a function of temperature was investigated in undercooled gallium and under-cooled Ga-In alloy of eutectic compositian to examine structural changes in undercooled melts. A continuous linear extension of the resistivity was determined for 99.999 pct pure gallium when undercooled below the stable Ga(l)-liquid, and the metastable Gafl)-liquid equilibria. In icndercooled Ga-In, a significant continuous negative departure from linear resistivity behavior was observed below the eutectic temperature. All resistivity measurements were made by the four-point probe technique. ThE temperature range for investigations into the physical properties of liquids before freezing has been greatly extended with the undercooling of massive melts of pure metals and alloys. Experimental determinations of thermal and transport properties and diffraction patterns in undercooled systems should throw further light on the structural behavior of liquids. Many measurements showing prefreezing anomalies in viscosity1&apos;5 have been regarded677 as evidence for structural groupings in the liquid state. However, the majority of viscosity results,~)~-~~ conforming to the law 17 =A exp(~/RT), do not support theories of atomic aggregation. For the system Cd-Sb, which exhibits a metastable solid phase, Fisher and Phillipsg reported that, in undercooled alloys, the activation energies for viscous flow below the stable liquidus are much higher than those above it. Similarly, Miller et a1." observed a doubling of the activation energies for electrical conductivity in undercooled Cd-ib alloys. These effects were suggested as indicative of atomic association and clustering in the liquids at lower temperatures. With the recent establishment of metastable phases in pure gallium22 and its alloys,23 attempts were made during this investigation to undercool bulk samples and to determine the temperature dependence of electrical resistivitv. Pure gallium, when undercooled below -17°C at atmospheric pressure, can crystallize to either a stable Ga(1) phase (mp, 29.78"C) or a metastable phase identi-

Jan 1, 1967

By L. L. Handy

The development of a theory for miscible liquid displacement requires evaluation of the variables which affect growth of the mixing zone between solvent and displaced oil. Factors which appear to be important are individual fluid viscosities, viscosity ratios, flood rate, fluid densities, flow characteristics of the porous medium and molecular diffusion coefficients of the fluid components. The primary purpose of this paper is to evaluate diffusion effects. Theoretical treatments to date have been limited to floods for which the viscosity ratio is one. Two principal theories have been proposed. Von Rosenberg adapted for porous media a theory derived for capillary tubes by G. Taylor.&apos; ,&apos; In this theory molecular diffusion perpendicular to the direction of flow is a primary factor in sharpening the flood front. Slow floods give sharper fronts for a given distance traveled than fast floods. An alternative theory considers miscible liquid displacement as a statistical problem.3,4,5 Diffusion is not an important factor in this theory, but it leads to the same general type of equation as von Rosenberg&apos;s. Both theories predict S-shaped concentration profiles with the same dependence on distance traveled. The statistical or "dispersion" theory predicts rate independence, however. To supplement rate studies a direct measurement of a diffusion effect would be helpful in evaluating which of the two proposed mechanisms best describes miscible liquid displacement for one-to-one viscosity ratio systems. No quantitative theory has been proposed for floods in which a low-viscosity fluid displaces a high-viscosity fluid. It might be anticipated, however, that the extensive fingering observed in floods with adverse viscosity ratios would increase opportunities for an exchange of components between displaced and displacing liquids by a diffusional process. Even if molecular diffusion were not an important mixing mechanism for one-to-one viscosity ratio systems, it could be significant in systems with adverse viscosity ratios. METHOD FOR EVALUATING DIFFUSIONAL MIXING In solvent flooding the displacing liquid can become mixed with the displaced liquid by a number of processes; in general, either mechanical or diffusional in character. To evaluate the contribution of diffusion, a method was sought that would distinguish diffusional mixing from mechanical mixing. Suppose two soluble tracers are added to the displacing liquid and that one tracer has a diffusion coefficient much greater than the other. Then, if diffusion is important during miscible liquid displacement and if diffusional transport is primarily transverse to the direction of flow, substances with high diffusion coefficients will have sharper concentration profiles than those with comparatively low molecular diffusion coefficients. In order to have material balance the profiles for the two tracers must intersect. For one-to-one viscosity ratios the intersection would occur at the 50 per cent concentration point. It is generally assumed that diffusion in the direction of flow is negligible. This assumption appears reasonable because flow velocities are ordinarily orders of magnitude greater than diffusional velocities. If this is not a valid assumption, however, the shapes of the fronts will be further degraded by longitudinal diffusion. The effects of rate and molecular diffusion coefficients are the opposite for longitudinal diffusion to those for lateral diffusion. Higher rates or lower diffusion coefficients would tend to give sharper fronts. The double tracer method offers several unique advantages in evaluating effect of molecular diffusion. First, the method gives a direct measure of the effect of diffusion on mixing zone size. Second, the effect of a difference in diffusion coefficients is determined from a single experiment. Comparison of several experiments at different rates is not required. Third, the method is applicable for any viscosity ratio. Fourth, the method is unaffected by density differences between the two fluids which might result in rate dependent gravity effects. EXPERIMENTAL PROCEDURE Suitable tracers are substances whose diffusion coefficients differ as widely as possible, but whose adsorptions on rock surfaces is negligible. The two tracers selected were methanol and sucrose in water solutions. The approximate diffusion coefficients for methanol and sucrose, as given in the International Critical Tables, are 1.3 X 10-5 and 0.3 X 10-5 m2/sec, respectively. Concentrations in the water solutions were 200 gm/liter of solution for both methanol and sucrose. These high