Search Documents

By H. F. Ramstad, F. D. Richardson

The equilibria (a) SiOz +Hz =SiO +H20 and (b) Si + SiO, = 2Si0 have beet1 studied at temperatures of 1425"to 1600°C ad 1310°to 1485°C respectively. The stattdard free energy changes for the tzrro reactions are given by the equatiotts Combination of the results for both equilibria leads to tiotz removes certain anomalies in existing high-terlzperature data for equilibria involving silica and silicon in iron. In many metallurgical processes and in many laboratory investigations silicon monoxide undoubtedly plays an important role. It is unfortunate therefore that wide differences exist between the results obtained by different investigators1-7 in their studies of such equilibria as In an attempt to put our knowledge of SiO on a surer basis, an exhaustive study has been made of equilibria [I] and [2] at temperatures ranging from 1300" to 1600°C. Reaction [I] was studied by measuring the amounts of silica which could be condensed from streams of Hz or Hz + HzO which had previously been brought into equilibrium with silica at temperatures ranging from 1425" to 1600°C. Reaction [2] was studied by measuring the material that could be condensed from streams of Hz or argon which had been brought into equilibrium with mixtures of silicon and silica at temperatures ranging from 1310" to 1485°C. EXPERIMENTAL Materials. The silicon was "superpure" grade and contained less than 0.1 pct impurities. The silica was prepared from pure mineral quartz; this was crushed and treated with concentrated hydrochloric acid to remove particles of iron, washed with water, and finally dried at 120°C. For the hydrogen + silica reaction, the silica was sized to —20+100 mesh. For the silicon + silica reaction, the two materials were ground to a fine powder in an agate mortar. The hydrogen and argon were commercial oxygen-free gases. The gas streams were controlled with capillary flow meters and the volumes were measured by wet gas meters. After passing through the meters, the gases were partially dried by silica gel. The hydrogen for the HZ + SiOz reaction was then passed through palladised asbestos at 300°C and dried with magnesium perchlorate. The efficiency of oxygen removal was checked throughout the experiments by passing the gas over an electrically heated strip of nichrome, used as an indicator as described by Rathman and de itt.' When mixtures of HZ + Hz0 were required, the partial pressures of water vapor (1.8 to 22 mm) were obtained by passing the hydrogen through oxalic acid dihydrateg' lo held at various controlled temperatures, O.l°C, by means of a water bath. The hydrogen for the Si + Si02 reaction was purified by passing it over a mixture of 3 parts of magnesium to 5 of lime heated to 600°."l1 u The argon for this reaction was passed through titanium powder (-3/16 in. + 100 mesh) heated to 900°C. The nitrogen used to prevent the reaction products escaping from the condenser (see later), was deoxidized by copper or iron at 600°C. All these gases were finally dried with magnesium perchlorate. Furnace, Temperature Contr01, and Measurement A molybdenum resistance furnace was used for both sets of experiments. The reactions were conducted inside a high-grade alumina tube, 36 in. long and 1 in. in diam as indicated in Fig. 1. With this arrangement an even temperature zone (2"C) 4 cm long was satisfactorily obtained. The temperatures were kept constant by means of a proportional controller actuated by a Pt-Pt 13 pct Rh thermocouple. This was placed between the two alumina tubes, so that the temperature at the junction was 1400" to 1450°C. Up to 1485"C, the temperatures were measured with Pt-Pt 13 pct Rh thermocouples. For higher temperatures an optical pyrometer was used, this being sighted (through the glass window 1 in Fig. 1) on the end of the alumina tube, that held the SiOz or Si +SiOz mixture, 10 in Fig. 1. The optical pyrometer was recalibrated whenever a change was made in any part of the apparatus situated in the hot zone. Successive readings with the optical pyrometer were reproducible to within 1"C. Equilibrium Apparatus and Procedure. Hydvogen and Silica Reaction. The apparatus is shown in Fig.

Jan 1, 1962

In the past few years the hydraulic aspects of rotary drilling have received considerable attention. It is generally recognized that accurate prediction of circulating pressures is desirable, particularly in those areas where a delicate pressure balance is necessary to prevent both blowouts and lost circulation. For calculation purposes the standard hydraulic equations for the pipeline flow of Newtonian fluids have been altered to facilitate similar calculations for non-Newtonian or plastic fluids.'.' This latter classification includes most colloidal types of drilling muds. Two properties: plastic viscosity, un, and Bingharn yield value, Y,, are commonly used to define the flow characteristics of a plastic fluid. These are commonly obtained with multispeed viscosimeters such as the Fann V-G meter.3 Normally these measurements are obtained at some surface temperature, and are not corrected to circulating temperatures for calculation purposes. This may result in considerable error, particularly in the laminar flow calculations which comnlonly apply to the drill pipe-borehole annulus. The purpose of this study was to investigate the effect of temperature on the flow properties of some water-base muds. A laboratory model Fann V-G meter was used for measurement of viscosities, yield values and gel strengths. An aluminum water jacket with O-ring seals was machined to fit around the mud cup of the Fann meter. All mud samples were prepared by thorough mixing in a Hobart mixer. Viscosity and gel strength measurements were made with the Fann meter using standard procedures. Flow properties of muds were studied at 80, 120, 160 and 180°F. The compositions of the test samples are given in Table 1. The typical effect of temperature on the flow properties of these drilling muds is illustrated by Fig. 1. It was found that plastic viscosity and apparent viscosity decreased with an increase of temperature. However, the curves were not linear, and did not appear to follow any definite trends or patterns. Yield point data showed much more scattering.' Since a usable method for predicting drilling mud flow behavior with temperature was desired, an attempt was made to find a simple relationship between these variables. The basic approach used by Havenaars for predicting the same effect on turbulent viscosity resulted in Fig. 2, which shows that

By C. D. Stahl, S. M. Farouq Ali, W. E. Culham

One- and two-dimensional mathematical models have been developed that simulate transient, two-phase flow of hydrocarbon mixtures in porous media in a manner that accounts for interphase mass transfer. Numerical simulations of one-dimensional depletion-drive experiments using a two-component hydrucarbon fluid were used to establish the validity of the mathematical models. In addition, the experimental and numerical data were used to demonstrate that production rate had a relatively insignificant effect on the recovery of individual hydrocarbon components from the experimental system, and that attainment of equilibrium between phases is possible for a wide range of liquid and vapor velocities in reservoirs containing light hydrocarbon fluids. Results of some two-dimensional numerical simulations are also presented. INTRODUCTION This study was undertaken to develop a mathematical model that would simulate transient, two-phase flow of hydrocarbon mixtures in porous media under conditions that result in interphase mass transfer and to test the validity of the assumptions used to set up the model. In addition, the study was designed to determine if production rate influences the recovery of individual hydrocarbon components from reservoirs producing by depletion drive. Two-phase flow in porous media, with interchange of components between the two phases, is important in many petroleum recovery processes. Studiesl-5 conducted within the last 3 years have outlined methods of solving multiphase flow problems incorporating mass transfer. Some of these studies have also indicated the importance of accounting for mass transfer under various producing conditions. An earlier work6 first demonstrated the importance of combining relative permeability data with equilibrium ratios in compositional balance methods. The mathematical model presented in this paper is formulated so that a phase behavior package, as described in previous papers,5,17 is not required as an integral part of the routine employed to solve for the primary dependent variables. The finite difference formulation is designed so that all the primary variables can be solved for simultaneously. This is accomplished by utilizing one basic set of equations. These innovations, which are in contrast to other models2,3,5 but are similar in some respects to the approach used by Taylor,l render the total problem computationally simpler than any of the previously referenced formulations. The mathematical model was developed by combining Darcy's law with a continuity equation for each hydrocarbon component. The principal assumptions invoked in the formulation were (1) that capillary forces and diffusional effects are negligible, and (2) that thermodynamic equilibrium exists in the reservoir at all times. No assumption as to the type of vaporization process was made in formulating this model. Experimental data were required to complete this study. These were generated by conducting several depletion drive experiments. The experimental apparatus consisted of a sandstone core enclosed in a pressurized casing. The apparatus was designed in such a manner that the core could be charged with a liquid hydrocarbon mixture and depleted at different production rates. The experimental tests were designed to determine the effect of production rate on component recovery. In addition, direct comparison of experimental and mathematically predicted

Jan 1, 1970

By C. C. Harris

The characteristics of some common size distribution equations are critically discussed. A generalized form of several well-known size distribution equations is obtained from a differential equation describing statistical distributions. The equation contains three parameters and can describe the major features of size distributions in the fine and the coarse size regions. A graphical method for its implementation is provided. The application of this and other equations to sets of data are compared both for the quality of fit and from a comminution kinetics viewpoint. If a narrow size range of a brittle material is broken sufficiently to obliterate the feed sizes, but not so severely that excessive secondary breakage occurs, a plot of cumulative fraction by weight undersize (Y) vs. sieve size (X) on a log log grid (sometimes called the Gates-Gaudin-Schuhmann plot)1-3 gives a straight line of slope ~ 1 over much of its range. On closer inspection, deviations in the extreme coarse range covering perhaps 10 to 30% of the total sample may be apparent (frequently, a change in slope from ~ 1 to a value different from unity) together with over-all slight irregular departures from a smooth line. One model of breakage postulates that the same fracture pattern persists throughout all size ranges.4 The fracture pattern is characterized by the slope of the line on a log log grid. Accordingly, a distribution of sizes broken in the manner specified earlier is expected to produce a size distribution having the same slope as that of a broken narrow size range. Additionally, the slight irregularities mentioned above should even out in the summation process, giving a smooth straight line of slope ~1 on a log log plot. This idealized state of affairs does not, however, describe the distribution of a multi-event process such as that for a tumbling mill product. On the average, these curves (plotted again log log) tend to have a slope in the fine region different from and usually less than unity,5 while the coarse size region curves with increasing slope for cases of mild reduction, and with decreasing slope when moderate or severe size reduction has occurred. In addition, there are usually slight irregular deviations from a smooth curve. Most curves display two distinct regions — fine and coarse — and a few curves show one or more intermediate regions. Feed size studies6,7 show an effect which theory is required to explain. For the same material, mill loading and other operating conditions, and the same time of grinding, plotting the dimensionless ratios Y vs. X/(feed size) does not reduce the size distribution data to the same curve. Y vs. X does not correlate the coarse region, but it can provide a crude correlation of the fine region which improves as size diminishes and as grinding time increases. The size distribution in the coarse region is somewhat more dependent on the feed size than is the distribution in the fine region wherein the degree of dependence diminishes as time proceeds; the distribution of sizes in the fine region are determined largely by the nature of the material and the comminution conditions. A comprehensive model of comminution must recognize that several different patterns or modes of breakage can occur in a mill; 5,8,9 that there can be some selection in that some size ranges are broken more than others; and that, while some particles may be the product of a single fracture event, or may even remain unbroken, others result from multiple rebreak-age. Thus, breakage in a tumbling mill is more complex than the Schuhmann4 model admits in that several different types of comminution micro-events* occur rather than just one type, and these events, though different in size scale, are of the same overall pattern as that which is visualized by the model. The concept of an average comminution micro-event has, therefore, a mathematical rather than a physical connotation, at least for tumbling milling. Equations with which to describe particle size distributions have been sought for over half a century.10-13 No equation presently in general use was first derived from an analysis of the statistical mechanics of breakage; whatever theoretical basis the existing

Jan 1, 1969

By Carl F. Lutz, Robert F. Mehl

A layer of brass was formed on 0 brass using a vapor-solid reaction technique. The variation in composition with distance within the phase layer and across the a -ß interface was determined by an electron probe. The diffusivities were calculated as a function of concentration; the diffusion coefficient in brass was found to change by a factor of twenty-five over a compositional range of 8 at. pct. An explanation is proposed to account for the large change in the diffusivity, based on possible stresses in the diffusion zone; it is suggested that the concentration dependence of the thermodynamic factor might decrease the activation energy with increasing zinc content. ALTHOUGH the body of data on diffusion coefficients (D- values) in terminal solid solutions and in isomorphous systems is quite large, there is but little quantitative information for intermediate alloy phases. This paper recounts measurements of D-values for the ? phase in the Cu-Zn system. The experimental method employed consists in the forming of alloy layers and determining and analyzing the concentration-distance (c-x) curve. The experimental method employed could have been applied to several systems; the Cu-Zn system was chosen because the D-values in the a and ß phases are well known, and can thus be used for purposes of comparison. The rates of growth, of intermediate alloy layers are known for a number of systems. In nearly all cases the thickness varies parabolically with time, as to be expected if the rate of growth is controlled by the D-values. In a few cases, nonparabolic behavior has been noted.1 Nonparabolic growth may be the result of a) a variation of the composition at the phase interface with time, b) interface-controlled kinetics, c) short diffusion times coupled with long vacancy lifetimes (in vacancy diffusion),2 and d) crystal reorientation during diffusion in anisotropic systems.3 In the Al-Ni system4 and in the Al-u5 system parabolic growth kinetics are slowly approached after an initial transient period. In general those phases stable at a given temperature in a system A-B appear as separate phase layers when A and B are put in contact and given a diffusion he at-treatment. In most cases the compositions at the interface of adjacent phase layers are those read from the phase diagrams for the termini of the respective phases. This discontinuity in concentration is taken as independent of time, and the growth of one phase at the expense of another is assumed to be independent of interface kinetics, i.e., the rate of interface movement is controlled by volume diffusion in the phases concerned. wagner6 has given many solutions for multi-phase diffusion processes, assuming that the chemical diffusion coefficient, D, is independent of concentration, as have others. Jost7 has pointed out that the familiar Boltzmann-Matano solution is as valid for a (c-x) curve exhibiting concentration discontinuities at phase boundaries as for the usual single-phase case. Heumann and associates8,9 have applied this solution to the multi-phase case, but lacking concentration data within the several layers were forced to assume linear concentration gradients, thus calculating only average D-values. The purposes of the present study were: 1) to determine the dependence of D on concentration, 2) to calculate the intrinsic diffusion coefficients D?Cu and D?Zn, 3) to establish the time-law for the movement of the interface, 4) to determine the concentration limits at the ß-? interface, and, 5) to attempt to clarify the mechanism of diffusion. EXPERIMENTAL PROCEDURE The ? phase is exceedingly brittle; conventional solid couples and conventional sectioning methods were thus inapplicable. Vapor-solid couples were accordingly employed. Unfortunately the equilibrium vapor pressures of Zn for the ? phase are unknown; to escape this awkwardness, samples were enclosed in a capsule containing powder of an alloy of known composition, sufficiently large in quantity as to be effectively an infinite source, maintaining the concentration of Zn at the sample surface constant with time. The sources of Zn-vapor which can be employed in reaction with Cu, or a brass, or ß brass, to form a layer of ? brass satisfactorily wide in concentration range, are alloys of the phases ?, or ? + E or ? + d (depending on the temperature).

Jan 1, 1962

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 Craig R. Barreft, Oleg D. Sherby

The steady-state creep characteristics of pure polycrystalline copper were studied in the temperature range 400" to 950°C and in the stress range 400 to 7000 psi. Tests were conducted in dry deoxidized hydrogen and all specimens exhibited high ductility before fracture. Two activation energies for creep were established as a function of temperature. Above about T/T, = 0.65, the activation energy for creep was 48,000 cal per mole whereas between 0.5 and 0.65 T, it was 28,000 cal per mole. The observed temperature dependence of the activation energy for creep follows the same trend as exhibited by the self-diffusion coefficient of copper in copper single crystals and silver in silver single crystals. It is therefore believed that the rate-cmtrolling mechanism in creep of copper in the stress and temperature range studied is diffusion-dependent and is probably associated with a dislocation-climb process. The low activation-energy region is attributed to enhanced volume difFusion from the presence of dislocations which act as short-circuit paths for diffusion. THERE is extensive experimental evidence1,2 that the activation energy for creep, Q,, for pure poly-crystalline metals tested at temperatures above one half the absolute melting temperature (0.5 Tm) is a constant and equal to the activation energy for volume self-diffusion, Qsd This equality is generally interpreted as meaning that creep in this temperature range is controlled by a dislocation-climb process. The existence of high-temperature creep activation energies not equal to the respective self-diffusion values has been suggested both theoretically3"5 and experimentally.6-8 Two examples of apparent nonequality between Q, and Qsd have been reported for the cases of silver7 and copper,' where activation energies for creep much lower than those for self-diffusion persist to rela- tively high temperatures (up to 0.6 T,). These observations suggest that dislocation climb may not be the only important process in the high-temperature deformation of these metals. That dislocation climb does indeed occur relatively slowly in these two metals due to their low stacking-fault energy has been suggested by many investigators.9-10 Although the creep behavior of copper has been studied extensively at low temperatures, the high-temperature activation energy for creep is not well-defined. The data of Feltham and Meakin8 indicated that Q, is temperature-dependent above one half the melting temperature, equaling a value of -30,000 cal per mole at temperatures up to 550°C and increasing to a value of -49,000 cal per mole at 700°C. No data has been reported above 700°C. The value of 49,000 cal per mole is very near the value of Qsd for copper reported by Kuper et al." as 47,100 cal per mole. In view of the lack of high-temperature creep data and the apparent temperature dependence of Q, for copper the present investigation was undertaken to determine the value of Q, over the temperature range of 400" to 950°C. It was hoped that such a study would clearly document the temperature dependence of Q, and yield significant information as to the nature of the rate-controlling creep mechanisms. MATERIALS AND PROCEDURES The material used in this investigation was high-purity OFHC copper. After rolling and recrystal-lization treatments spectrographic and chemical analyses showed only 0.0008 pct Ag, 0.0002 pctCa, 0.0002 pct H2, and 0.002 pct O2 as trace impurities indicating a purity of greater then 99.995 pct. Specimens were machined from cold-rolled sheet and the gage sections polished with 2/0 emery paper prior to annealing treatments. Annealing was done in both vacuum and dry deoxidized hydrogen atmospheres. No variation in creep strength was found as a result of the different annealing atmospheres. The final annealing temperature was either 700°, 800°, or 1000°C depending on the temperature range to be studied. These treatments resulted in average grain diameters of 0.03, 0.4, and 1.0 mm, respectively. The precise mechanical-thermal history of these

Jan 1, 1964

By P. N. Nikiforuk, M. C. Zoerb

An analog computer model has been developed of a large, walking dragline. This model permits changes in the configuration of the dragline, or changes in its digging cycle, to be readily investigated on a computer. As a result, it is possible to determine ways of improving the dragline performance without having to carry out many expensive experiments on the actual machine. This model can also be used as a simulator for dragline operator training. INTRODUCTION The movement of large volumes of material is always a costly operation and the increasing trend towards the use of large capacity equipment does not necessarily produce the economy desired. This problem is of particular concern to the mining industry in its open-pit operations where large draglines are used to remove the overburden. Since this operation often forms a very large part of the total cost of open-pit mining it is particularly important that it be kept as economical as possible. This requires that the dragline be operated at its maximum possible productivity. In order to optimize the operation of a dragline it is desirable that a suitable model of the dragline be available for study. Once such a model is available the effects of changes in the configuration of the dragline, or changes in its digging cycle, could be readily investigated on a computer and methods of improving its performance could be determined without having to carry out expensive experiments on the actual machine. Such models would be of use not only to dragline users, but also to the manufacturers of draglines for design purposes. Presumably studies of this type could eventually lead to computer-con trolled stripping operations. As far as can be determined very little work1'2 has been done on the derivation of dragline models. Such models must be derived before work of the nature indicated in the previous paragraph can be carried out. As a small contribution to this field this paper describes work that has been carried out by the senior author during the past few years on the derivation of one such model. This model describes the drag, hoist and swing operations of a 35 cu yd, commercially available, walking dragline. This model was derived for use with an analog computer and can be used either for optimization studies, or for the training of dragline operators. DESCRIPTION OF DRAGLINE General: The particular dragline for which an analog computer model was developed is owned by the Great West Coal Company of Estevan, Saskatchewan. It is a walking dragline and its function is to remove 40

Jan 1, 1967

By E. T. Turkdogan, J. H. Swisher

bi the fivst part of the paper results ave given on the rate of decarburization of Fe-C melts ln CO2-CO atmospheres at 1580°C. The rate -controlling step is believed to he that irvlloluing dissociation of curbotz dioxide on the suvfuce of the melt. 4 genevral reaction mechanistm is poslnlated jor gels-t11eta1 veactions oc-curit~g on the surface of iron coutcotamncited with chemi-sovbed osygesL. Oxygen the present work on decavbuvization of liquid iron and previous studies on the kinetics of nitrogen absorption and desorplion are discussed in terms of the postulated mechanism, ManY of the early studies of rate of decarburization of liquid steel were of an exploratory nature and laboratory exppriments carried out pertained to open-hearth or oxygen steelmaking processes. References to previous work on this subject may be found in a literature survey made by Ward. Using more sophisticated experimental techniques, several investigators have recently studied the kinetics of decarburization of molten Fe-C alloys in oxygen-bearing gases. For example, Baker et al2.' reported their findings on the rate of decarburization of liquid iron, levitated by an electromagnetic field, in carbon dioxide-carbon monoxide-helium atmospheres. In these levitation experiments the samples used were small in size, e.g., -0.6-cm-diam spheres weighing -0.7 g, and the rates were measured for decarburization from about 5 to 1 pct C at 1660°C. The rates obtained under their experimental conditions were considered to be controlled primarily by gaseous diffusion through the boundary layer at the surface of the levitated melt. Parlee and coworkers3 measured the rate of absorption of carbon monoxide in liquid iron. The rates were found to follow first-order reaction kinetics, yielding a reaction velocity or a mass transfer coefficient in the range 0.2 to 0.4 cm per min. The coefficient was found to decrease with increasing carbon content of the melt. These investigators attributed the observed rates to the transfer of carbon or oxygen through the diffusion boundary layer adjacent to the surface of the melt. In the work to be reported in this paper, an attempt has been made to study the kinetics of gas-metal surface reactions involved in the decarburization of liquid iron. EXPERIMENTAL The experiments consisted of melting 80-g samples from an Fe-1 pct C master alloy in an induction furnace and decarburizing in controlled CO2-CO mixtures at 1 atm pressure and 1580°C. The master alloy was prepared by adding graphite to electrolytic "Plastiron" melted in racuo. None of the impurities in the master alloy exceeded 0.005 pct. The reacting gases were dried by passage through columns of anhydrone; in addition, CO2 impurity in carbon monoxide was removed by passage through a column of ascarite. A schematic diagram of the apparatus is shown in Fig. 1. A 1.25-in.-diam recrys-tallized alumina crucible containing the sample was placed inside a 3-in.-diam quartz reaction tube, all of which was surrounded by an induction coil. A 450-kcps induction generator was used as the power source. Water-cooled brass flanges, which contained the gas inlet, gas exit, and sight port, were sealed to the top of the reaction tube with epoxy resin. The reacting gases were metered with capillary flowmeters and passed through a platinum wire-wound alumina preheating tube, 0.25 in. ID and 11 in. long. The gases were preheated to about 1300°C. A disappearing-filament optical pyrometer was used to measure the melt temperature. The pyrometer was initially calibrated against a Pt-6 pct Rh/Pt-30 pct Rh thermocouple. The temperature was controlled to within +10°C by manually adjusting the power input to the induction coil. In a typical experiment, an 80-g sample of the master alloy was melted in a CO2-CO atmosphere having pcO2/pco = 0.02 and flowing at 1 liter per min. A negligible amount of carbon was lost and no significant reduction of alumina from the crucible occurred during melting, e.g., 0.005 pct Al in the metal. After reaching the experimental temperature of 1580°C, the gas composition was changed to that desired for a particular series of decarburization experiments. The duration of the transient period for obtaining the desired gas composition at the surface of the melt was about 20 sec . The flow rate of the reacting gas was maintained at 1 liter per min. After a predetermined reaction time, the power to the furnace was turned off. During freezing, which took about 10 sec, the amount of gas evolution was not sufficient to result in a significant loss of carbon. The samples were analyzed for carbon by combustion and in a few cases they were analyzed for oxygen by the vacuum-fusion method. RESULTS A marked increase in the rate of decarburization of iron with increasing pcO2/pco ratio in the gas stream is demonstrated by the experimental results given in Figs. 2 and 3 for pco2/pco ratios from 0.033 to 4.0. In one series of experiments, denoted by filled triangles in Fig. 2, the reacting gas was diluted with argon (48 vol pct) resulting in a slower rate of decarburization. Samples from two series of experiments with pco2/pco = 0.033 and pco2/pco = 0.10 (with argon dilufion) were analyzed for oxygen. In these Samples the oxygen content increased with reaction time

Jan 1, 1968

By P. E. Rider, K. A. Gschneidner, O. D. McMasters

The solid solubilities for thirteen rare-earth metals in gold were determined by using the X-ray parametric method. Solubilities ranged from 0.1 at. pct for lanthanum in gold up to 8.8 at. pct for scandium in gold. The solubilities from lanthanum to gadolinium were very small and essentially constant, but a sharp increase occurred from gadolinium to scandium. The large solubilities for the heavy rare-earth metals were not expected because of the large size and electrochemical differences between rare-earth atoms and the gold atom. Contributions from first- and second-order elasticity theory plus an electronic contribution were found to reasonably account for a more favorable size factor. Electron transfer from the rare-earth metal to the gold Is thought to occur such that the resultant rare-earth and gold electronegativities are favorable for solid-solution formation. It was also found that this mutual adjustment of size and electronegutivity does not occur if the pure-metal size factors are greater than a critical value of 25 pet. The eutectic temperatures for ten systems were determined and these remained fairly constant at approximately 809 "C for the lighter lanthanide metal-gold systems until the Er-Au system was reached, at which Point the eutectic temperature successively increased reaching a maximum of 1040°C in the Sc-Au system. This rise was correlated to the size factor becoming more favorable for solid-solution formation at erbium. The valence state of ytterbium was found to change from two in the pure metal to three when ytterbium is dissolved in the gold matrix. RECENT results1 reported concerning the solubility of holmium in copper, silver, and gold, showed that the solubility of holmium in gold was quite large, 4.0 at. pct, compared with 1.6 in silver and 0.02 in copper. The small solubilities of holmium in silver and copper are quite reasonable in view of the large size difference (22.2 and 38.2 pct, respectively), large electronegativity difference (0.59 for both systems), and possible unfavorable valency factor (assuming one for silver and copper and three for holmium). The large solubility in gold, however, is unexpected because these same factors are also unfavorable for holmium and gold (22.5 pct size difference and 0.69 electronegativity difference), and because the light rare-earth metals, lanthanum, cerium, and praseodymium, have negligible solid solubilities in gold.2 In view of this unexpected behavior, it was felt that a study of the solid solubilities of most of the rare-earth metals in gold would be desirable to better understand the factors involved in the formation of solid solutions. Of the rare-earth metals added to gold in this study, only ytterbium is divalent in the pure metallic state (the other rare-earth metals are all trivalent) and many of its physical properties (such as the metallic radius, electronegativity, and so forth) are much different from those of the normal trivalent rare-earth metals.' The properties of ytterbium are such that one would expect solid-solution formation to be less favorable for ytterbium in gold than for any of the normal trivalent rare-earth metals. But chemically ytterbium is known to possess a stable trivalent state, and it is quite possible that ytterbium may alloy as a trivalent metal under certain conditions rather than as a divalent metal. Because of the dual valency nature and because so little is known about the alloying behavior of ytterbium, the gold-rich ytterbium-gold alloys are of special interest. EXPERIMENTAL PROCEDURE Materials. The gold used in this investigation had a purity of 99.99 pct with respect to nongaseous impurities. In general the rare-earth metals were prepared by reduction of the corresponding fluoride by calcium metal.3 The impurity contents of the metals used in this study are given in Table I. Preparation of Alloys. Two- or 3-g alloy samples were prepared by arc melting. The samples, with the exception of some of the Er-Au alloys, had weight losses of 0.5 pct or less. All alloy concentrations noted in this paper are nominal compositions. After arc melting, the alloys were wrapped in tantalum foil, sealed off in quartz tubing under a partial atmosphere of argon, homogenized for approximately 200 hr at 780°C, and then quenched in cold water. X-Ray Methods. The X-ray parametric method was used in determining the solubility of the rare-earth metals in gold. filings were sealed in small tantalum tubes by welding under a helium atmosphere. The tantalum tubes were then sealed in quartz tubing under a partial argon atmosphere, and annealed for times ranging from 1/2 to 3 hr (the length of time was inversely proportional to the an-

Jan 1, 1965

By H. K. Birnbaum

The effect of 1ow-temperature stabilization anneals on the structure of the 0 phase Au-Cd alloys and on the diffusionless transformations observed in these alloys was examined by X-yay diffraction techniques. A phase separation in the ß-phase region was proposed to account for the experimenta1 results. The effects of quenching from elevated temperatures on the transformation behavior of these alloys were shown to be consistent with the proposed mechanism. IT has been shown that the high-temperature ß phase (CsCl structure) of the Au-Cd alloy system transforms to a phase having an orthorhombic (D2) ß&apos; structurel1-3 for compositions near 47.5 at. pct Cd and a tetragonal (4/m, m, m) ß" structure* in the vicinity of 50.0 at. pct Cd. Both transformations are diffusionless, crystallographically reversible, and occur on cooling at about 60° and 30°C respectively. The temperature interval from the beginning to the end of the transformation is of the order of 5°C in each case. Although the transformations are normally athermal, some of them have been reported to occur isothermally.= wechsler6,7 has shown that the effects of quenching a 49.0 at. pct Cd alloy from elevated temperatures are consistent with the retention of a nonequilibrium number of lattice vacancies. Annealing of these quench effects results in a broadening of the X-ray reflections.8 After a suitable quench, the 47.5 at. pct Cd alloy transforms to a phase having not the p&apos; orthorhombic structure but another structure which has properties similar to that of the ß" tetragonal structure.5.9 This change in the type of transformation has also been obtained after long anneals in the ß-phase region at about 70oC10 The present investigation was primarily concerned with the structural changes accompanying the above transformation phenomena. The change in transformation product and accompanying physical changes during an anneal in the ß phase have been termed stabilization effects. Experimental Procedure —The results reported in this investigation were obtained with the use of a Norelco diffractometer fitted with a temperature-controlled cryostat. The specimen temperature was controlled to better than ± 0.l°C during the measurements. CrKa radiation monochromated electronically with the use of a scintillation counter and pulse height analyzer was utilized. Specimens containing 47.5 and 50.0 at. pct Cd were prepared by sintering filings obtained from homogenized ingots of the proper alloy composition. (Gold of 99.999 pct purity and cadmium of 99.98 pct purity were used). All heat treatments were carried out with the specimens capsulated in vacuum ( < 10 % mm Hg) or in a He-H gas mixture. The quenching technique used in these experiments was to drop the pyrex capsule which contained the specimen from the annealing furnace into water, the temperature of which was controlled. The pyrex capsule shattered on contacting the water resulting in a relatively rapid quench. After the heat treatment, the specimens were mounted in the diffractometer and were left undisturbed in the diffractometer specimen holder during each sequence of measurements. EXPERIMENTAL RESULTS A) Low-Temperature Annealing—The transformations which were considered "normal" for these alloys were those obtained athermally during furnace cooling at approximately 50°C per hr after an elevated temperature anneal. Under these experimental conditions, the specimens were observed to transform to phases having structures whose diffraction patterns could be indexed as the ß&apos; orthorhombic structure for the 47.5 at. pct Cd and as the 0" tetragonal structure for the 50.0 at. pct Cd alloys. The transformation temperatures on cooling were approximately 60" and 30°C, respectively. Under the "normal" conditions both transformations were observed to go to completion, i.e., the entire volume of the ß phase was transformed to the product phase. In some specimens an extremely weak ß 110 reflection was observed at 20°C indicating that a small amount of retained ß was present. The effect of low-temperature annealing on the nature of the diffusionless transformations was examined for the 47.5 and 50.0 at, pct Cd alloy. The specimens were annealed in evacuated capsules at temperatures in the vicinity of 600°C (as specified in Table I) for 24 hr and were then cooled to 100°C at a rate of 50°C per hr. The specimens were then removed from the capsules and mounted in the diffractometer without allowing the specimen temperature to drop below 80°C. Annealing at the low temperatures was accomplished in the diffractometer by means of the cryostat which was mounted around the specimen. During the low-temperature anneals the lattice parameter, integral breadth of the reflections, and ratios of the integrated intensities of the fundamental and super lattice reflections for the 0 cubic phase were periodically determined. After annealing for the required time, the specimens were slowly cooled in the diffractometer and the diffraction patterns were recorded as a function of temperature. The specimens were cooled until the phase transformations were completed, following which the specimens were heated and diffraction

Jan 1, 1960

By R. C. Hassinger, D. V. von Rosenberg

Transverse dispersion has received considerably less treatment in the literature than has longitudinal dispersion. Different methods for determining transverse dispersion coefficients have been used in different investigations, and the results obtained have not been consistent enough to permit accurate generalizations as to the effect of various physical parameters on the magnitude of these coefficients. A numerical solution to the differential equation describing transverse dispersion in the absence of longitudinal dispersion was obtained to enable one to calculate the dispersion coefficient from experimental results. The more general dispersion equation including longitudinal dispersion also was solved numerically to give quantitative limits of a dimensionless group within which the assumption of negligible longitudinal dispersion is justified. Possible experimental procedures were examined, and one utilizing a cylindrical packed column was chosen for the determination of transverse dispersion coefficients. Values of these coefficients were determined for a system of two miscible organic fluids of equal density and viscosity, for two sizes of packing material over a wide range of flow rates in the laminar regime. The dispersion coefficient was found to decrease, for a constant value of the product of packing size and interstitial velocity, as the size of the packing material particles increased. INTRODUCTION Longitudinal dispersion has received extensive treatment in the literature, and consequently is better understood than its orthogonal counterpart, transverse dispersion. Many mathematical models of dispersion processes assume that transverse dispersion is rapid enough to damp out any radial** concentration gradients and therefore may be neglected. Laboratory and production results, however, indicate that this is a poor assumption.l,2 Various experimental procedures for determining transverse dispersion coefficients have been used in previous investigations, but the results have generally been expressed by similar correlations. The transverse dispersion coefficients obtained, however, have often varied considerably for given values of the correlation parameters. We feel that further experimental determinations of transverse dispersion coefficients will help alleviate some of the inconsistencies in these empirical correlations. One assumption implicit in all previous investigations is that of negligible longitudinal dispersion in the experimental system. An attempt to justify this assumption often is made using intuitive reasoning, but it is apparent that this reasoning must break down as the condition of zero flow rate is approached. A mathematical examination of the equations describing the system yields physical limits outside of which the assumption of negligible longitudinal dispersion is invalid. BACKGROUND In a porous medium, the "effective molecular diffusivity" De is less than the molecular diffusivity D measured in the absence of a porous medium, due to the tortuous path which a diffusing molecule must travel. Various authors3-5 have reported values of the ratio De/D in the range of 0.6 to 0.7. When there is fluid flow within the porous medium, mass transfer occurs by convective dispersion as well as by molecular diffusion. These are separate phenomena and can be treated as such on a microscopic scale. However, the mathematical complexity is such that only extremely simple geometries could be considered, and the results hardly would be applicable to the complex geometries existent in actual porous media. A more

Jan 1, 1969

By S. Shapiro, G. Krauss

The strutural and cr~stallo~aphic features of the plates of austenite produced by the martensite to aus-tenite or reverse martensitic tramformation have been determined in an Fe-33 wt pct Ni alloy. Micro-focus X-ray techniques and single-surface trace analysis 072 bulk samples yielded two distinct habit planes, (0.174, 0.309, 0.~35)~ and (0.375, 0.545, 0.749jM. The former plane uas the one predorninantly observed and its existence was verified by transn~ission electron microscopy. The orientation relation-ship between the reversed austenite and the parent martensite was approximately the same as the Nishiyamna and other relationships reported for the direct tramformation. Replica and thin-joil obser-vations show that both high densities of tangled dis-locations and occasional twins constitute the fine structure of the reversed austenite. Application of the phenomenological theory to a variant of the predominant habit plane defines an irrational plane and direction for the second shear in accordance with the comnplexity of the fine structure. The shear accompanying the reverse martensitic transformation is at least 0.51, the maximnum value of- the tangent oJ the tilt angle measured on surface replicas. A mechanism relating the fine transformation twins in martensite to the nucleation of reversed austenite of the predominant habit is proposed. The reversal of Fe-Ni martensite takes place both at the edges of martensite plates and in a piecewise fashion within them.&apos; The shear-type nature of the reverse transformation is verified by the surface relief which accompanies both edge-type reversal2 and the fragmentation of plates of martensite&apos; as a result of rapid heating above the A, (austenite start temperature). The orientation relationship between the edge-type reversal product and the parent martensite, as determined by transmission electron microscopy, is reported to be within 4 deg of either the Kurdyumov-Sachs or Nishiyama relationships,&apos; but there is no work at present in the literature relating to the crystallography of the platelike reversed austenite. The fine structure of reversed austenite after heating to 50°C above the Af is reported4 to be composed primarily of tangled and jogged dislocations in concentrations up to 10" per sq cm. A replica investigation of partially reversed Fe-Ni martensites5 corroborates the increase in dislocation density following reversal and presents indications of other possible modes of fine structure. This paper reports on an investigation performed to examine in detail the morphology and crystallography of the plates of reversed austenite and the shear which accompanies their formation in a matrix of Fe-Ni martensite. EXPERIMENTAL PROCEDURE Discs of a high-purity Fe-32.9 wt pct Ni single crystal served as the starting material. The single crystal had been produced in the course of an earlier investigation6 and the carbon content was determined at the time to be 0.006 wt pct. The M,, and the A,, were respectively -120" and 300°C. Partial transformation to martensite was performed at -125°C and reversal of some of the martensite was accomplished by heating in the temperature range 340" to 355°C. Most samples were heated to the reversal temperature by immersion&apos; in a salt bath for 2 min. For surface-relief studies some polished and etched Samples were heated in a tube furnace for ten min in a hydrogen atmosphere maintained over the samples throughout the heating and cooling cycles. Samples were prepared for metallographic examination by electropolishing and etching with an HC1-HNOs-H20 et~hant.~ On one of the surface-relief samples two sets of fiducial scratches were placed on the etched sample by drawing it over a slurry of 0.06 p alumina on a "microclothJ&apos;. The orientations of individual plates of martensite were determined by X-ray analysis. A Rigaku-Denki microbeam X-ray generator in conjunction with a Micro-Laue camera with facility for precision location of the sample in front of the beam was employed for this purpose. The collimator size was 30 p and the specimen to film distance was 5 mm. The Laue photograms were enlarged to an equivalent 3-cm specimen to film distance for analysis. The orientation of each of the plates of martensite was compared to that of the parent austenite and the relationship be -tween the two phases was, in all cases, within a few degrees of those predicted for the direct transformation. The orientation relationship between one of the larger plates of reversed austenite and its parent martensite was determined in a similar manner. The habit plane of the islands of reversed austenite in the X-rayed plates of martensite was determined by a single-surface trace analysis. Each reversal island had with the parent phase one comparatively straight boundary which was presumed to represent the habit plane trace. The pole locus technique7 was applied to traces from six different plates of martensite to determine the indices of the habit plane. A 40-cm stereographic net was employed for this analysis The morphology and fine structure of the reversed austenite were studied by electron microscopy of pre-shadowed direct carbon replicas,5 and the macroscopic shear was evaluated by a two-stage replica technique similar to that employed in electron fractog-

Jan 1, 1968

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&apos; 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,&apos; flow in packed columns,&apos; fracture flow capacity&apos; 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 R. D. Taylor, Jim Douglas Jr., H. H. Rachford Jr., P. M. Dyke

A major obstacle to the use of wetting agents in .secondary recovery by water flooding is the adsorption of the agents on the sand. As a result of adsorption, the surfactant always lags behind the floodwater front. Consideration of the chromatographic theory of adsorption indicates that the detergents will not lag as much if used in very high concentrations. An investigation was made of the possibility of using high concentrations economically by flowing slugs of wetting agents followed by normal flood water. The experiments consisted of adsorption studies on Alundum powder and Berea sandstone. Flow rests on a 12-in. Alundum core and 22-in. Berea core were used to determine rate of detergent movement. The results of the flow experiments indicate that the relative rate of surfactant advance is, indeed, sensitive to the concentration of the agent. A 10 per cent slug moved with a rate that war 78 to 95 per cent as fast as the rate of advance of the flood water. By contrast, one with 25 ppm (the number of parts of commercial detergent in a million parts of water on a weight basis) concentration moved less than one-fourth as fast as the flood water, and calculations indicate that in very long porous systems the rate of movement of the lower concentrations will be a small fraction of the rate of advance of the flood front. The results. of the adsorption studies were utilized to calculate the rate of advance of the detergent when only the initial concentration was known. The calculated rates showed substantial agreement with the experimental flow tests in the high concentration ranges. The adsorption results were also used to estimate the cost of the materials for a slug-type surfactant flood in the field. In addition to the faster rates of movement, the concentrated detergent slugs removed much more oil than the dilute solutions. However, the effectiveness of the slug process depends on many variables. The quantity of oil removed is increased markedly by increasing the flooding rate. The efficiency is also influenced by the type of crude, type of reservoir rock and initial water saturation. Therefore, a careful analysis of each reservoir system is required before the economic value of the process can be determined. INTRODUCTION It is well known that the displacement of oil by invading water during water flooding is far from complete. It is generally agreed that the unrecovered oil is retained in the porous medium by the capillary forces which may be relatively large compared to the forces generated by the flowing water. Therefore, it was logical that some early workers should turn to surface-active materials to reduce the capillary forces to facilitate the release of oil. As early as 1927,&apos; a patent was granted for the use of surface-active materials in water flooding. In 1932, when soap solutions were passed through Bradford and Venango sands, it was reported that the results were inconclusive, erratic and that "further investigation is needed to determine exactly the function of the solution and to obtain a clearer insight into the phenomena involved."&apos; Some of the modern scientific reports conclude with a similar statement,&apos; showing that the lack of agreement on the mechanism of oil removal by wetting agents is still very widespread even though several comprehensive studies have been reported.&apos;." Although there is a lack of agreement as to the general effectiveness of the detergents for water flooding, most investigators do agree that all of the common detergents are strongly adsorbed onto the solid surfaces of the reservoir. In the early calculations it appeared that all additives would be lost before reaching much of the formation area which contained the additional oil to be removed. Experiments indicated that if the usual small waterflood concentrations of wetting agents were used, the rate of advance of detergent through the formation would be only a small fraction of the rate of advance of the flood front. Indeed, some investigators4 felt that the use of wetting agents would never be economically feasible because of their adsorption. For example, DunningG estimated that the wetting agent in concentrations of 25 ppm, would advance only 0.05 times as fast as the flood front. Ojeda, et al,&apos; found that a surfactant in a concentration of 10 ppm moved less than 0.01 times as fast as the flood front. It is significant, however, that both investigators found that increased concentrations of wetting agents moved faster, relative to the flood front, than solutions at the lower concentrations. Ojeda showed that an extrapolation of his data indicated a relative rate of 0.5 at 1 per cent concentration, while Dunning6 estimated a relative rate of 0.46 for a 1 per cent concentration. It was obvious that these concentrations could not be used for continuous injection because the cost of the injected detergent would far exceed the value of additional oil produced. Traditionally, detergents are used in very low concentrations for they show good

By Charles R. Wraith

The property of the Roan Antelope Copper Mines, Ltd., is situated at Luanshya, north central part of Northern Rhodesia, and is connected with the main line of the Rhodesia Railways, Ltd., by a 24-mile branch line from N&apos;dola. Most of the material and supplies that are not available in Africa are shipped into the country via Beira, on the East Coast of Africa, a distance of 1473 miles from Luanshya and the nearest ocean port. The output of copper is shipped by rail to Beira and thence by boat to the market. Passenger travel is mostly via Cape Town, a distance of 2162 miles from Luanshya. Excellent passenger trains operate on regular schedules between Cape Town and Northern Rhodesia. Luanshya is 13" south of the equator at an elevation of 4000 ft. above sea level, and has a mild healthful climate. During the rainy season, from October to March, the rainfall averages 53 in. The area is profusely covered with vegetation. An ample supply of water for industrial and domestic purposes is available throughout the year. The more important supplies for smelting purposcs, such as wood, coal, limestone, silica, quartzite, and clay for "dolly" mud, are conveniently available. Excellent grades of coal, coke and firebrick are obtained from the Wankie Colliery Co., Ltd., Southern Rhodesia. The Wankie colliery has just completed the erection of equipment for the manufacture of silica brick. An ample supply of quartzite of good grade is available for this purpose. Silica brick is also obtained from the Belgian Congo and from India. Limestone is quarried from a large deposit a few miles from Bwana Mkubwa, which is 30 miles by rail from the smelter. An abundant supply of clay for " dolly " mud and for lining launders is available within a short distance of the smelter. Metallurgy Copper is contained in the ore as chalcocite, chalcopyrite and bornite. Chalcocite predominates in the eastern end of the orebody almost to the exclusion of other copper minerals. The amount of chalcocite decreases with a corresponding increase of chalcopyrite and bornite toward the western extension, where the latter predominate.

Jan 1, 1934

By Charles Wraith

THE property of the Roan Antelope Copper Mines, Ltd., is situated at Luanshya, north-central part of Northern Rhodesia, and is connected with the main line of the Rhodesia Railways, Ltd., by a 24-mile branch line from N&apos;dola. Most of the material and supplies that are not available in Africa are shipped into the country via Beira, on the East Coast of Africa, a distance of 1473 miles from Luanshya and the nearest ocean port. The output of copper is shipped by rail to Beira and thence by boat to the market. Passenger travel is mostly via Cape Town, a distance of 2162 miles from Luanshya. Excellent passenger trains operate on regular schedules between Cape Town and Northern Rhodesia. Luanshya is 13° south of the equator at an elevation of 4000 ft. above sea level, and has a mild healthful climate. During the rainy season, from October to March, the rainfall averages 53 in. The area is pro-fusely covered with vegetation. An ample supply of water for industrial and domestic purposes is available throughout the year. The more important supplies for smelting purposes, such as wood, coal, limestone, silica, quartzite, and clay for "dolly" mud, are con-veniently available. Excellent grades of coal, coke and firebrick are obtained from the Wankie Colliery Co., Ltd., Southern Rhodesia. The Wankie colliery has just completed the erection of equipment for the manufacture of silica brick. An ample supply of quartzite of good grade is available for this purpose. Silica brick is also obtained from the Belgian Congo and from India. Limestone is quarried from a large deposit a few miles from Bwana Mkubwa, which is 30 miles by rail from the smelter. An abundant supply of clay for "dolly" mud and for lining launders is available within a short distance of the smelter.

Jan 1, 1933

By Thomas O. Ziebold, Robert E. Ogilvie

Experimental analysis of diffusion samples covering the entire Cu-Ag-Au system at 725°C has been carried out. Experimental coefficients are computed at twenty-eight points in the lermary system. It is found that the direct coefficienl for copper is nearly independent of the silver, content for 1ow-silver alloys and that the cross coefficient for silver correlates with the ther-nzodynatnic properlies. This lalter observalion leads to the conclusion that the mobililies of copper and gold are nearly equal. Thermodynarmic actiuilies for Cu-Ag-Au were computed from binary data and adjusted to be consistent with neasured tie lines across the two-phase region which extends into the ternary diagram at 725°C. Application of these therknodynarnic calculaliorls to the dijftcsion data terifies the Onsager reciprocal relations to the extent that experimental uncertainlies will allow. It is jound that diffusion neav the critical point of the 1200-phase field in the isothermal section at 725°C causes the composition gradients for copper and silver to approach infinity while the direct and cross coeficients become equal. From the standpoint of experimental Procedures, it has been shown that microanalysis of ternary samples may be carried out using a simple expression for the conversion of X-ray zntensily to composilion. The accuracy of&apos; this method has been demonstrated by the analysis of standard alloys and by the reproducibility of composition determinations using differed X-ray lines. DIFFUSION in a ternary alloy system is basically and significantly different from diffusion in binary alloys. When only two elements are present there is just one independent composition variable, and the diffusion profile of composition vs distance across a sample must change monotonically through all solid solutions bracketed by the terminal alloys. If two-phase regions are stable for compositions between the terminal values, these must appear as planar interfaces in the diffused sample, with the composition profile exhibiting sharp discontinuities. Adding a third element may alter this picture completely. Because of the additional degree of freedom in the composition, we may see composition profiles which do not change monotonically through solid-solution ranges, and stable, nonplanar phase interfaces. Even though two-phase regions may lie between the terminal compositions of a diffusion sample on the constitution diagram, they will not necessarily appear in the diffusion profiles. In binary diffusion only one coefficient is needed to describe the interdiffusion process. From Fick&apos;s law this coefficient relates the diffusive flux of one component to its own composition gradient. In ternary diffusion we must allow for cross interaction between the two independent species. By a linear extension of Fick&apos;s law we introduce "cross" or "off-diagonal" diffusion coefficients which relate the flux of one component to the composition gradient of the other independent constituent. Thus, we require four interdiffusion coefficients to describe the transport process in a ternary system. In its most general form,1,2 the rate equation for mass transport by isothermal diffusion in a system of s components states that the diffusive flux of component i, Ji, is a linear summation of all driving forces, Xk, each force being multiplied by an appropriate compliance coefficient:

Jan 1, 1968

By E. P. Klier, Jane Jellison

The course of the transformations on cooling in a series of Cu-A1 alloys has been followed by means of thermal analysis, X-ray diffraction, and optical metallogrAphy. Specimen size was found to have a pronounced effect on the time of completion of the pearlite transformation due to modification in the unmber of nuclei. The ordeving transformation zvas observed for all cooling rates for all alloys studied, and is considered insuppvessible. At intermediate cooling&apos; rates, a transitional struchire for 011 alloy compositions and based ox the bainite reaction was observed. The alloys of aluminum and copper possess a eu-tectoidal region at from about 9.4 to 16 pct Al. Above the critical temperature, the disordered bcc phase, 0, is stable. At subcritical temperatures, the 0 phase decomposes under equilibrium conditions into the stable phases, a and y,. At the eutectoid composition, these phases develop a lamellar pearlite structure, in appearance much like that found in steels. During rapid quenching from supercritical temperatures, the diffusion processes required for the formation of the equilibrium phases are suppressed, and the 0 transforms to an acicular martensitic phase, designated B or y&apos; depending on the aluminum content, or to a mixture of the two. Early investigators of this system concentrated on determining the equilibrium phase diagram, &apos;&apos;&apos; the structure of the equilibrium phases,3 and later the isothermal decomposition of the 0 phase at subcritical temperatures, or the decomposition of the martensite on tempering. Other work which treated the nonequilibrium transformations concentrated largely on the martensite transformations and the structure of the b and y&apos; phases.15"17 Since the pearlite and martensite transformations in this system resemble those in steels, the one being formed by nucleation and growth of adjacent lamellae of the equilibrium phases, and the other by a diffusionless, shear-type transformation, the present study has been undertaken to investigate the structures formed on continuous cooling, with the objective of constructing cooling transformation diagrams analogous to those derived for alloy steel systems. In particular it was desired to determine whether parallels between the A1-Cu and iron-iron carbide systems on continuous cooling extended to encompass the intermediate phase, namely the bainite structures, which form by a combined martensitic-diffusion mechanism,&apos;&apos; but which are supersaturated with respect to the equilibrium phase upon which their crystal structure is based. To this end the transformations in representative alloys of the Cu-A1 B eutectoid system have been followed, using the methods of thermal analysis, X-ray diffraction, and optical microscopy. EXPERIMENTAL PROCEDURE Materials. The materials tested were those used by Klier and Grmko 5 in their study of the isothermal decomposition of these alloys. The alloys were melted in graphite crucibles and chill-cast into 5 by 2 by 1/4 in. ingots. The ingots were then hot-rolled to sheets 1/8, 1/16, and 1/32 in. thick. The compositions are given in Table I. Specimen Preparation. Specimens for thermal analysis were in the form of small squares approximately 1/4 in. on a side. For the 11.9 and 13.5 Al-Cu alloys, the sheet thickness was 1/16 in. and for the 10.5 A1-Cu alloy the 1/32-in. sheet was used. One lead of a chromel-alumel thermocouple was spot-welded to either face of the square with the specimen itself thus forming the hot junction of the thermocouple. Specimens for X-ray analysis were prepared as filings from martensitic sheet stock and were sealed in glass ampoules under 1/2 atm of argon. Metallographic examination of the powder specimens was carried out on a portion of the heat-treated filings mounted in lucite. All specimens were examined metallographic ally after etching

Jan 1, 1965

By Kenneth E. Gray

The turbulent flow drag coefficients, or friction factors, have been experimentally determined for the cut-tings normally encountered in drilling operations. The gas law and average drag coefficients for characteristically flat particles (limestones and Shales) and for angular to sub-rounded particles (sandstones) were used to extend Newton&apos;s equation to a more useful form. The resulting equations are expressions for the slip velocity of a particle as a function of the particle size and shape, the bottom-hole injection Pressure, and the Bottom-hole temperature. The velocities necessary to lift actual rock cuttings as observed in the laboratory were compared with velocities predicted from the derived equations. Results indicate that there is no single correct circulating velocity and that the commonly used 3,000 ft/min linear velocity may he sufficient to lift only very mall particles. The advantage, in terms of horsepower requirements, of cycling high pressure air to obtain lower compression ratios was shown. INTRODUCTION A promising departure from standard rotary hydraulic drilling is air or gas drilling. In many instances the rate of penetration and bit life increases resulting from this method have been very substantial. Further the application of air drilling in areas of lost circulation or water sensitive pay zones has been highly successful. Purpose For the drilling contractor, the practical question remains, "How much air pressure and volume should I have?&apos; lt is necessary to have a reasonable knowledge of the air velocities and pressures required in air drilling operations for the proper design and most advantageous use of surface equipment and for adequate removal of cuttings from the borehole. Eased on experience, most operators agree that at least 3,000 ft/min linear velocity is necessary for satisfactory hole cleaning.&apos; It was the purpose of this investigation to evaluate the drag coefficients, or friction factors, for the cuttings normally encountered in drilling operations, i.e., sandstones, limestones and shales, and to utilize the values obtained to develop an expression for the minimum velocities and pressures (hence the volumes) necessary to lift these cuttings. THEORY Several investigators have done extensive work on the ability of drilling mud to lift bit cuttings. They conclude that the ability of a drilling mud to lift cuttings is dependent on the density difference between the rock being drilled and the drilling mud, the cutting size and shape, the mud flow constants, and the flow state (laminar, transitional or turbulent) of the mud. However, drilling muds are not true fluids but have a variation of viscosity with velocity in laminar flow and a constant apparent viscosity in turbulent flow. A survey of the literature revealed that no experimental work had been done to determine the velocities necessary to lift actual rock cuttings using a compressible fluid such as air or gas. In 1953, Nicolson6 stated that in air drilling fluid mechanics the particles could be assumed spherical in shape and that ; constant drag coefficient of 0.5 could be used due to the highly turbulent flow of the circulating air. By equating the turbulent resistance on the particle to the gravitational force on the particle, he obtained the expression, ^ = 2.67^y ,.......(1) where V, is slip velocity of spherical particle, ft/sec; D is particle diameter, in.; p, is particle density, 1b/ft8; and p, is fluid density, lb/ft3. However, as stated earlier, most operators seem to agree that adequate hole cleaning can be obtained by circulating a sufficient volume of air to give a linear velocity of 3,000 ft/rnin in the annular space. Basic Equation The velocity with which a solid particle freely falls through a fluid will increase until the accelerating force, gravity, is equal to the resisting forces, buoyancy due to the fluid displaced by the particle and friction due to the relative motion of the particle through the fluid. When the accelerating and resisting forces become equal,