Search Documents

INTRODUCTION The meaning of cash flow has already been defined as after-tax earnings available for use by the corporation or business. It is therefore apparent that a significant part of mineral economic evaluation is concerned with correct and accurate computation of taxable income under existing tax laws, and the proper treatment of tax allowances and deductions. The various methods of depreciation, depletion, and amortization are all means whereby the investor can earn back or recover, before taxes are imposed, the investments made for purposes of doing business or for the production of income. Each of these three types of tax deduction are applied to different investment items such as property, equipment, buildings, research, etc. The cost of buildings, equipment, machinery, and transportation equipment are representative of capitalized business costs that must be recovered by depreciation over the life of the asset. The property investment costs of acquiring natural resources such as minerals, oil and gas, and standing timber are examples of investment costs that may be recovered by depletion. Several other capitalized business costs, such as research expenses and trademark expenses, may be recovered by amortization. The term "capitalized" means that the in- vestment cost is not deducted as an operating expense in the year that it is incurred, but is spread out over a period of years through depreciation, depletion, or amortization deductions according to lawful schedules. Of course, any capitalized cost can only be deducted once. The other costs of doing business which are not capitalized, such as labor and operating costs, are said to be "expensed," i.e., they are written off for tax purposes in the year they occur. DEPRECIATION Depreciation, according to US tax law, is a deduction from taxable income permitted as a "reasonable allowance for the exhaustion, wear and tear of property out of its use or employment in a business." The original cost of the property is the usual basis for determining depreciation. Since depreciation write-off is applied over the useful life of the asset, the first step is to determine the estimated useful life or its allowable depreciable life. These two lives may or may not be the same. No average life is applicable to all types of depreciable assets under different conditions or in different types of business. It is the taxpayer's responsibility to determine the useful life of the asset according to his own operating conditions and experience, within the guidelines set by the US Internal Revenue Service (IRS). The Class Life Depreciation Range (ADR) System permits the taxpayer to choose from a range of depreciable lives at 20% below or above the present guideline life. The second step required in calculating depreciation is to determine an estimated salvage value for the asset. The salvage value is an estimate, made at the time of acquisition, of the amount of money which can be realized by the sale of the asset at the end of its useful life. The four depreciation methods most widely used in the US are: (1) the straight-line method, (2) the declining balance method (double declining balance), (3) sum of the years - digits method, and (4) units of production method. Straight-Line Depreciation Method The straight-line depreciation method is the simplest for calculating depreciation. The depreciation for each year is determined by subtracting the estimated salvage value from the first cost of the asset and dividing the result by the years of life. The method gives equal annual deduction over the depreciable life. Salvage value equal to or less than 10% of the first cost may be ignored. [first cost - salvage Straight line = years of life ] Declining Balance Depreciation Method Declining balance depreciation, or double declining balance (DDB) as most frequently applied, is used for accelerating depreciation in the early years of the asset and is the fastest depreciation method allowed by the US Internal Revenue Service. The declining balance

Jan 1, 1980

By G. K. Manning, W. E. Few

T has been known for some time that both'inter-granular carbide and intergranular oxide phases cause brittleness in molybdenum. Parke and Ham' indicated that 0.0025 pct 0 present in molybdenum after solidification from the melt was sufficient to induce intergranular brittleness during forging. However, material containing 0.06 pct C could still be forged. However, adding this amount of carbon to insure deoxidation causes the formation of carbides on freezing. These carbides are precipitated at grain boundaries, and Rengstorff and Fischer' have shown they have an unfavorable effect on the room-temperature ductility of molybdenum. Take? in 1928 proposed that the solubility of carbon was a constant value of 0.30 pct from room temperature up to 1800°C. In 1935, Sykes, Van Horn, and Tucker' obtained X-ray data suggesting that the solid solubility of carbon in molybdenum was between 0.07 and 0.09 pct at temperatures from 1500" to 2100°C. The a-solubility line for carbon was included in the Mo-C diagram given in the Metals Handbook as a tentative line based on the data of Sykes, Van Horn, and Tucker. In order to determine the partial phase diagrams for these systems, it was necessary to construct a high-temperature furnace suitable for heat treating molybdenum at temperatures up to 4000°F. A molybdenum resistance furnace was built for this purpose in which samples could be heat treated in a purified atmosphere and quenched from any desired temperature.' Experimental Work Solubility of Carbon: Arc-cast molybdenum supplied by the Climax Molybdenum Co. and 18-gage molybdenum wire samples produced by powder metallurgy at the Fansteel Metallurgical Co. were used in determining the a-solubility line for the MO-C system. The cast molybdenum contained 0.032 pct C initially. Two lots of 18-gage wire were used, one containing 0.011 pct C and the other 0.005 pct C initially. All material contained a carbide phase as received. Fig. 1 illustrates the carbides found in the untreated material. Heat treatment in a purified hydrogen atmosphere was found to reduce the carbon content in the as-cast molybdenum. Purified argonT also decarburized molybdenum samples at temperatures above 3500°F; however, the carbon content remained constant during heat treatment in purified argon at lower tem- peratures. The argon was purified by passing it over titanium at 750 °C and then through a magnesium perchlorate drying tower. The hydrogen was purified by passing it through a commercial hydrogen deoxidizer, followed by a drying tower. Specimens having different initial carbon contents were heat treated in argon at 3500°F for periods ranging from 5 min to 5 hr and were analyzed for carbon. Samples of the same final carbon content had the same microstructures regardless of time at temperature. Furthermore, a gradient between center and surface in the amount of undissolved carbide particles present was never observed. It was concluded that at these low carbon levels, and at the high temperatures involved, diffusion proceeds so rapidly that the specimen is always substantially homogeneous in carbon content, in spite of the gradual loss in carbon. Based on the above findings, the following procedure was used to determine the carbon-solubility limit at 3000°, 3500°, and 4000°F. At each of these temperatures, a number of samples were heat treated, as shown in Table I, and rapidly cooled to room temperature. Since the samples lost carbon in proportion to the length of time heated, a series of samples were obtained of varying carbon contents. By use of a rapid cooling rate, the phases present at the high temperature were retained at room temperature. All samples were examined microscopically. Those found to contain either no excess carbide phase or only a very small amount were then analyzed for carbon. The samples actually analyzed for carbon

Jan 1, 1953

By P. E. Arnold, G. S. Ansell

RELAXATION in metals has been studied in detail by many workers in recent years.1-5 These studies have shown that there is an energy-loss peak observed in a metal placed in mechanical resonance at low frequencies which may be correlated with various metallurgical processes. These experiments have been largely used to study single-phase materials. In order to study the effect of a finely dispersed second phase in a metal matrix upon this phenomenon, the relaxation behavior of an aluminum-aluminum oxide SAP-Type alloy in both the fine-grained as-extruded and coarse-grained recrystallized conditions was investigated by the use of internal-friction measurements. The alloys studied, designated MD 2100, consist of a matrix of commercial-purity aluminum containing a finely dispersed second phase of aluminum oxide. The aluminum oxide is present in this alloy in the form of fine platelets approximately 130A units thick and several microns on edge. The mean free path betwen dispersed platelets is approximately 0.5 p. The grain size of this alloy is extremely stable at temperatures as high as the melting temperature of the aluminum matrix. In the as-extruded condition the grain diameter is several microns, while in the recrystallized condition it is several millimeters in diameter. The damping characteristics of the alloy samples were determined using a torsional-pendulum device. The internal-friction constant was determined for this alloy in both the fine-grained as-extruded and coarse-grained recrystallized conditions over the temperature range from 343" to 923°K at frequencies of 0.9, 1.3, and 1.8 cycles per sec. An internal-friction peak was observed in testing the fine-grained alloys which was not observed for the coarse-grained alloys. This peak, which might be associated with a type of grain-boundary relaxation, was observed at each of the three different test frequencies for the fine-grained sample. The temperature at which the relaxation peak occurred for each of the three different frequencies is shown in Table I. From these data an activation energy of relaxation for each pair of frequencies was determined by comparing the temperature shift of the relaxation peak with frequency between each of the test frequencies, assuming that an Arrhenius-type rate equation is applicable. The activation energy of relaxation determined from these data is also shown in Table I. The mean value of the activation energy was about 6 kcal per mole. The activation energy for the relaxation process determined in this investigation indicates that the relaxation associated with the grain boundaries in this aluminum-aluminum oxide SAP-Type alloy is quite a different process than grain boundary relaxation in single-phase aluminum alloys. In these latter alloys, the activation energy for relaxation has been found to be the activation energy for self-diffusion, approximately 36 kcal per mole1,2 Where the activation energy is the same as the activation energy for self diffusion, it has been proposed that the relaxation phenomena is dependent upon vacancy formation and motion. It is apparent therefore, that, in this aluminum-aluminum oxide SAP-Type alloy, the relaxation must be due to another type of atom transport. It has been shown previously in studies of the high-temperature, low-stress, steady-state creep behavior of as-extruded aluminum SAP-Type alloys, that the activation energy for steady-state creep is one that may be characterized by the nucleation of dislocations from the grain boundaries. It is proposed that the relaxation behavior observed here is due to a similar effect. In this case, the internal-friction peak is a resonance behavior due to the time-dependent nucleation of dislocations from grain boundaries in the sample. A process of this type would have an activation energy which is undoubtedly stress-dependent and of the form Q = Q, + (dQ/du)u

Jan 1, 1962

By F. R. Scauzillo

Equilibrium data which should be useful in the design and/or evaluation of glycol dehydration units were prepared from an analysis of various published data and the correlation of these data by the use of the thermo-dynamic equilibrium ratio. The equilibrium ratios of water are used to solve the glycol absorber problem; such solutions are necessary to define the number of trays and the glycol circulation rate needed to meet drying requirements. Activity coefficients were obtained which relate directly to the equilibrium ratios of water in the water-TEG-natural gas system. These activity coefficients have been used to calculate the equilibrium dew points for aqueous TEG concentrations of 60 to 99.9 weight per cent for the temperature range of 40° to 120° F. They also provide a means of calculating equilibrium ratios for water in the water-TEG-natural gas system; this applies to any desired TEG concentration and to the temperature range from 40° to 120°F. INTRODUCTION During the last several years, the drying of natural gas with aqueous triethylene-glycol (TEG) solutions has become very prominent. Most users of TEG as a drying agent have been satisfied with the performance of TEG solutions at the conditions used; however, there always has been some discussion of the drying ability of TEG solutions at conditions not commonly encountered such as temperatures below 50° or 60°F, but more particularly temperatures above 100°F, and pressures above 1,000 psia. Today, when higher wellhead temperatures such as 120°F are more commonly encountered many are skeptical of TEG's ability to dry sufficiently well to provide dew points around 32°F, which is generally the maximum tolerable when attempting to dry a gas to contract specifications of 7 Ib/MMscf. These views probably have evolved to some extent from the days when suppliers would not guarantee dew-point depressions in excess of 65° to 75 °OF. Also, the feelings about TEG drying may have arisen from the lack of information about the equilibrium relation of water in the drying operations. Porter and Reid' have reported equilibrium data for a 95 per cent TEG solution, while Townsend2 reported equilibrium data for 95, 98 and 100 per cent TEG solutions. Townsend also presented a calculational method which obtains activity coefficients in the liquid phase and, subsequently, the equilibrium water content of a saturated gas over glycol from published3 atmospheric dew points. Wise, Puck and Failey' presented activity data for the water-TEG system at atmospheric pressure. In the past, the equilibrium ratios published for water in a 95 per cent by weight aqueous TEG solution have been used indiscriminately by many for all concentrations of TEG encountered in gas drying operations. Since essentially all such gas operations require the use of more concentrated TEG solutions, a study was undertaken to correlate the existing equilibrium data for aqueous TEG and gas systems and to provide some means of calculating the equilibrium ratio of water in the natural gas-water-TEG system where TEG concentrations were other than 95 per cent. The data presented herein consider only the equilibrium drying ability of the TEG solutions and do not consider the effect of temperature and pressure on the tray eficiency of contactors. In other words, this paper is concerned primarily with the development of the equilibrium relationship between water in the dried natural gas and the water in the lean TEG entering the top tray of the absorber. BASIC EQUILIBRIUM RELATIONS A relationship which relates the K value of water in the natural gas-water-TEG system to the vapor pressure of water at the system temperature, the total pressure of the system and the activities of the liquid and gas phases has been evaluated. It is well known that, for any component of a mixture at equilibrium between a gas phase and a liquid phase,

By W. Mckewan, W. M. Fassell

The oxidation rates of copper have been determined at temperatures from 600" to 900°C in oxygen from 14.7 to 400 psi total oxygen pressure. The oxidation rate of copper is unchanged by oxygen pressures within the range studied. The observed data extend Feit-knecht's conclusion concerning the pressure independence of copper from atmospheric pressure to 400 psi. All samples studied in the above-mentioned temperature and pressure range have both Cu2O and CuO present in the oxide film. NUMEROUS papers have been published on the oxidation of copper, some of which have noted the effect of oxygen pressure on the oxidation rate of copper. It must be noted, however, that no data on the oxidation rates of copper at pressure in excess of 1 atm have been reported in the literature. (Le Chatelier¹ reported that a black oxide of silver was formed at 300°C and 15 atm pressure.) This is likewise true for the other metals. This investigation was initiated to obtain experimental oxidation rate constants for pure metals at elevated temperatures and high oxygen pressure. Copper was selected as the first metal since probably more reliable data are available on its oxidation rates and related physical and chemical properties than for any other metal or alloy. As a general summary, the following conclusions appear well established for the oxidation of copper: 1—The oxide coating consists of Cu2O with a superficial coating of CuO providing the ambient oxygen pressure exceeds the equilibrium pressure for the coexistence of Cu2O and CuO.² The following sequence of phases is then Cu/Cu2O/CuO/O2 (gas). The Cu2O consists of large crystals that have no relation to the orientation of the metal crystals, except for a very thin layer adjacent to the metal." The CuO consists of very fine crystals randomly oriented. If the oxygen pressure is below the equilibrium pressure for the coexistence of Cu2O and CuO, the coating consists of Cu,O only. 2—The rate-determining factor in the parabolic oxidation of copper is the diffusion of Cu+ ions through the Cu,O layer.4,6,7 The Cu+ ions are accompanied by electrons to maintain electrical neutrality. The mechanism of the reaction at the Cu2O/CuO interface and the method of growth of the CuO is in some doubt. 3—On copper, there are three reported pressure effects: (a) At pressures below 0.3 mm Hg oxygen pressure, the oxidation rate increases directly with pressure.', ' .Vxygen starvation at the surface may account for this effect. (b) At oxygen pressures from 0.3 to 63 mm, the oxidation rate increases as the 7th root of the oxygen pressure providing the temperature is such that only Cu2O is present in the film. (c) At pressures from above the equilibrium pressure for the coexistence of Cu2O and CuO (varies with temperature) to atmospheric pressure, the oxidation rate is independent of the oxygen pressure.' Equipment In order to determine the oxidation rates of metals in an oxidizing atmosphere at pressures up to 400 psi, the equipment is somewhat different than for low pressure work as is shown in Figs. 1 and 2. A Leeds and Northrup Micromax recording controller (A) is used to control the temperature of the Ni-chrome furnace winding in conjunction with a Micromax electric drive unit. The drive unit adjusts a 20 amp Variac (variable transformer) to control the power input into the Nichrome furnace winding. Due to the "self heating" of the sample by oxidation, frequently encountered in oxidation rate studies, it is essential that the sample temperature be measured independently. This is done by means of a series of six chromel-alumel calibrated thermocouples located spirally around the oxidizing metal sample. Each thermocouple is partially shielded from the direct radiation of the furnace wall by

Jan 1, 1954

By R. W. Heckel, R. A. Tanzilli

Treatment of the finite, diffusion-controlled, two-phase, moving-interface homogenization problem has been carried out by numerical methods for a wide range of input parameters. Solutions have been obtained for planar, cylindrical, and spherical diffusion geometries. The effects of parameters such as the interdiffusion coefficients (concentration-independent) of- both phases, the concentrations (time-independent) of each phase at the interlace, mean composition, initial size of the dissolving- phase, and initial composition of- the two phases are considered. The results are displayed to show both the position of the interjace and the degree of homogenization (expressed in terms of concentration at the matrix phase symmetry boundary) as a function of diffusion time. Typical exaw~ples to which the results may be applied are discussed. A limited number of exact, closed-form solutions to the multiphase, diffusion-controlled, moving-interface problem exists in the literature, e.g., Refs. 1 to 7. The exact, closed-form solutions generally consider the outermost phases to be infinite in extent and also assume an interfacial velocity function. It is through use of such boundary conditions that closed-form solutions are readily obtained. An extensive literature also exists on numerical solutions to the mathematically analogous solidification problem, e.g., Refs. 8 to 12. In addition, a numerical solution involving mass transfer with spherical interfaces and infinite boundaries has recently been carried out.13 A practically important class of diffusion problems exists, however, in which the boundaries of each phase away from the interface are finite. For example the finite boundary condition arises when one sets up a model for analyzing the homogenization behavior of two-component powder compacts in two-phase systems or for studying the solution behavior of precipitates. Although exact, closed-form solutions have as yet not been found to this class of problems, solutions are possible using numerical methods. The purpose of this paper is twofold: 1) to describe a numerical method for treating the two-phase, diffusion-controlled, moving-interface problem involving finite boundaries with planar, cylindrical, and spherical interfaces; 2) to present the results of a parametric study which identifies the role played by each variable appearing in the analytical model. R. A. TANZlLLI,Student MemberAIME,onleaveofabsence from the General Electric Co., Re-entry Systems Department and A. W. HECKEL, Member AIME, are Graduate Student in Materials Engineering, and Professor, Department of Metallurgical Engineering, respectively, Drexel Institute of Technology, Philadelphia, Pa. Manuscript submitted May 6. 1968. IMD MATHEMATICAL MODEL The three geometries analyzed are shown in Fig. 1. For each configuration a symmetry element is shown, across whose boundaries the mass transfer is zero. Mean composition, C, and the characteristic solute-rich phase half-thickness (or radius), 1/2, are the independent variables while the characteristic symmetry element half-thickness, L/2 (measured from solute-rich phase center to the matrix phase plane of zero mass transfer), is a dependent variable. The symmetry element for a planar interface accurately represents a real diffusion situation involving the stacking of symmetry elements of constant size. Because cylindrical and spherical symmetry elements of constant size cannot be stacked to fill space, the meaning of the dimension, L/2, for these situations may be considered in terms of the following approximations: and where NA is the number of parallel cylindrical elements per unit area on a surface normal to the cylinder axes and Nv is the number of spherical elements per unit volume. The model consists of a solute-rich phase matrix phase composite of limited mutual solubility. The concentrations at the interface are assumed constant with time; no assumption of local equilibrium at the interface is necessary except for the case where the

Jan 1, 1969

By E. C. W. Perryman

IN spite of the large amount of work which has been carried out on the recrystallization of aluminum and its alloys, there has been no complete investigation of the recrystallization characteristics of A1-Mg alloys. The most complete investigation is that reported by Bungardt and Osswald' who, working with commercial-purity materials, found that the temperature both of beginning and of complete recrystallization increased with magnesium content up to about 2 pct Mg, then decreased with a further increase in magnesium up to 5.0 pct, and thereafter increased again. Michaud and Segolz working with 99.93 pct A1 also found the temperature of beginning of recrystallization to decrease as the magnesium content increased from 2.13 pct to about 4.5 pct. Chossat3 using aluminum of 99.99 pct purity found that 0.12 pct Mg had no effect but that 0.4 pct Mg increased the recrystallization temperature by about 30°C. That no effect was found for 0.12 pct Mg does not agree with the later results reported by Demer and Beck4 who found that 0.025 pct Mg in 99.99 pct A1 increased the time for complete recrystallization and furthermore increased the activation energy for the recrystallization process. Kratz" using 99.997 pct A1 found that 0.5 pct Mg increased both the temperature of beginning and of complete recrystallization by 100°C. Chadwick and Hooper" working with commercial-purity material found that 1 pct Mg increased the temperature of complete recrystallization by about 40 °C for cold reductions in the range 17 to 72 pct. With smaller amounts of cold work produced by stretching (10 pct strain), they found no significant effect of magnesium. The latter result agrees with that obtained by Williams and Eborall' who found no effect of 2 pct Mg on the recrystallization temperature of commercial-purity aluminum strained 10 pct or less by stretching. Experimental Procedure Material Used: The chemical analyses of the materials used throughout this investigation are given in Table I. Fabrication: The alloys were cast in iron molds which had been washed thoroughly with alumina. After scalping, the ingots were annealed for 16 hr at 450°C, pressed, annealed at 450°C for 16 hr and hot rolled to 0.25 in. After hot rolling, the slabs were alternately annealed and cold rolled to 0.036 in. Cold-rolling reductions between intermediate anneals were kept at about 30 pct and a high intermediate annealing temperature of 450°C was used so as to obtain, as nearly as possible, similar starting grain sizes for the different materials. All intermediate anneals were done in a salt bath and were followed by cold-water quenching so as to retain all the magnesium in solid solution. The final cold reductions investigated were 20, 30, 40, 50, 60, and 80 pct. The initial grain sizes, that is, the grain size before the final cold reduction, were all very similar. Annealing Procedure: Specimens, 1 in. sq, were cut from the 0.036 in. sheet and annealed in a salt bath, the temperature of which was controlled to k2"C. Separate specimens were annealed for different times at 350°, 375", 400°, and in one or two cases at 415°C and after annealing were quenched in cold water. The time taken for a specimen to reach temperature was approximately 5 sec. Metal-lographic examination after annealing showed that no precipitation of Mg+413 occurred during cooling. Procedures Used for Determining Recrystallization Characteristics: All specimens were electrolyti-cally polished and etched by the method described previously1° and polishing was continued until about 0.002 in. had been removed from each side of the specimen. Using this technique, the smallest recrys-tallized grain which could be observed with cer-

Jan 1, 1956

By Thomas E. Garnar

Zirconium and hafnium are curious elements because they are almost always found together in nature. Zirconium was discovered by Klaproth in 1789 and isolated 35 years later by Berzelius. Because hafnium's chemical properties are so close to those of zirconium, it was not discovered until 1923 when Coster and Von-Heresey detected it using x-ray spectrographic analysis. These elements occur most commonly in nature as the mineral zircon (ZrSiO4) and less commonly as the oxide baddeleyite (ZrO2). They also are found as a variety of other silicates. Zircon was identified as a component in alluvial and beach sands in 1895, but it was not produced in any quantity until 20 years later. Zircon is always a coproduct from TiO2 mining and processing. During World War I it was produced as a coproduct of beach sand mining for titanium minerals just south of Jacksonville Beach, FL, and was patented as a refractory. It was not until the 1930s when Zircon Rutile, Ltd. began mining at Byron Bay on the east coast of Australia that zircon was first used as a foundry sand. Later in the 1940s, NL Industries, Humphreys Gold Corp., and E. I. du Pont de Nemours & Co., Inc. began production of zircon sands from fossil beaches in northeast and northcentral Florida. Baddeleyite first became available as a commercial product in 1916, but never in the quantity of zircon. Zircon holds a unique position as an industrial mineral because it is used for both its physical and chemical properties as well as an ore of zirconium and hafnium metals. Australia is the major zircon producer followed by South Africa and the US. Zircon is also produced in India, Sierra Leone, Sri Lanka, Malaysia, China, Thailand, and Brazil. Undeveloped heavy mineral beach sand reserves containing zircon are known to occur in Egypt, Malawi, Senegal, and Tanzania. Locations are shown in Fig. 1. Mineralogy Zirconium and hafnium are always present together in naturally occurring compounds. Most commonly they form as the silicate sometimes containing iron, calcium, sodium, manganese, and other elements. Less commonly they are found as oxides in combination with titanium, thorium, calcium, and iron. Zircon is the most common form. It is extremely resistant to weathering, therefore, often it is present in ancient beach sands and placer deposits. Although resistant to alteration from external sources, it is vulnerable to internal alteration as the result of thorium and uranium substituting for zirconium either in the zircon lattice or in solid solution. Alteration to the metamict state takes place as radioactive emanations from these elements disorder the crystal lattice, accompanied by hydration, reduction in specific gravity, and changes in color. The mineralogy of zircon and baddeleyite and variation of physical and chemical properties are shown in Table 1. The effects of zircon's internal alteration and darkened color can be corrected by heating to 1000°C for 30 min or so. A certain time-temperature relationship makes it possible for decolorization and reordering of the crystal lattice to take place when heating at higher temperatures for shorter periods or at lower temperatures for longer periods. After heating, x-ray diffraction patterns are sharper, the grains may be harder, and except for particles of malacon, the grains will be white or colorless in the absence of grain surface coatings. The color change is only temporary and the particles slowly darken with time. Zircon fluoresces light yellow in shortwave ultraviolet light. If exposed to ultraviolet light for extended periods of time, white (calcined) zircon turns purple. With time the induced color begins to fade and eventually disappears (e.g., calcined Florida zircon exposed to ultraviolet light for 24 hr still retains faint, but discernible, purple color after six months).

Jan 1, 1983

By D. E. Baldwin

A method has been devetoped for predicting the produced concentration profile for a miscible slug in a five-spot pattern. The technique consists of dividing the five-spot into radial etements and applying an approximate radial solution of the dispersion equation to each element. To test the method, laboratory experiments were conducted in five-spot models. The models were 2-ft X 2-ft X 2-in. Berea sandstone slabs simulating one quarter of a five-spot. An aqueous system was used with tritiated water as the tracer fluid. Effiuent concentration profiles were obtained for slug sizes ranging from 0.7 to 10.2 per cent of the pore volume. Agreement between predicted and experimental profiles was excellent. A prediction was also made for a field tracer test. For this case the predictive technique was modified to account for stratification by use of a layer cake model. INTRODUCTION Dispersion in porous media is of growing interest to the petroleum industry because of the increasing importance of secondary recovery operations. One of the controlling factors in the recovery of oil by miscible displacement or by the use of waterflood additives is the degree of mixing between the fluids of interest. Consequently, one of the prerequisites of any recovery prediction is an adequate method of predicting dispersion. Also, any quantitative interpretation of waterflood tracer tests is dependent upon the ability to determine the dispersion of the tracer as it flows through the reservoir. During the past few years there have been reported in the literature the results of a large number of investigations, both theoretical and experimental, concerning dispersion in porous media. (An excellent review of these results has been presented by Perkins and Johnston." Only a few of these investigations, however, deal with non-lin-ear flow systems. Raimondi et al.9 have presented an approximate solution of the dispersion equation for radial, diverging flow. Lau et al. 7 and Bentsen 1 have presented laboratory data for radial systems and have shown that the dispersion can be adequately predicted by Raimondi's solution. Bentsen has also presented a small amount of data for a radial, converging system. Brigham and Smith' have applied Raimondi's solution, with approximating assumptions, to a five-spot pattern to predict the behavior of a field tracer test. The afore-mentioned investigations have been concerned almost entirely with radial diverging flow and have neglected the effect of converging flow. For this reason the present investigation was undertaken to develop a method of predicting dispersion in a five-spot pattern including the effects of both diverging and converging flow. The method which was developed consisted of dividing the five-spot into diverging-converging radial elements and applying an approximate solution of the disperson equation to each element. DEVELOPMENT OF PREDICTION METHOD The development of the prediction method consisted of two phases: construction of a radial model of a five-spot pattern and derivation of equations to describe the dispersion in the radial model. FIVE-SPOT MODEL A radial model which simulated the flow characteristics of a five-spot pattern was constructed by dividing the five-spot into diverging-converging radial elements (Fig. 1). Because of symmetry the model represented only an octant of a complete five-spot. Each element had an included angle of l° (making a total of 45 elements) and a radius r, such that the sum of the individual element areas was equal to the area of the five-spot octant. The normal five-spot fractional flow behavior was determined by averaging the experimental five-spot data of Caudle and Witte,' Dyes, Caudle and Erickson 5 and Fay and Prats. 4 The resulting curve is shown as a solid line in Fig. 2. To match this curve with the radial model it was

Jan 1, 1967

By J. S. Levine, M. Prats

Several methods are available for calculating the performance of solution-gas-drive reservoirs from the PVT properties of the oil and from the relative permeability and other properties of the formation. These methods require a number of simplifying assumptions. The present method of computation has made use of a high-speed computer to solve simultaneously the nonlinear partial differential equations that describe two-phase flow by solution-gas drive in order to calculate the performance of a reservoir. Some of the results obtained by the nonlinear partial differential equation solution are compared with those obtained with an approximate method, which has been called the semisteady-state solution. The pressure and saturation profiles from the wellbore to outer boundary calculated by the two methods are compared for one constant-terminal-rate case and two constant-teminal-pressure cases. The agreement in these profiles, as well as in the values of average reservoir pressure and cumulative recovery, leads to the conclusion that, for most engineering calculations, the semisteady-state method will give a reasonable approximation to the numerical solution of the differential equations describing solution-gas drive. An unfavorable (as regards ultimate oil production) set of relative permeability curves was used in the calculations in the belief that the effect of the parameters which were studied would be emphasized to a greater degree. Furthermore, the reservoir was assumed to be completely homogeneous, and these results should not be considered applicable to any other type of reservoir. Gravity effects are not considered. The absolute permeability was varied from 25 to 0.5 md. At an economic limit of 2 B/D, the recovery for a 25-md reservoir is about 1.8 times as great as that for a 0.5-md reservoir. The effect of permeability on the producing gas-oil ratio is minor. Once PVT properties of the oil and the relative penneability properties of the reservoir are fixed, the pro- during gas-oil ratio is found to be a function of the fraction of oil recovered. Well spacings of 10, 40 and 80 acres were considered. For the assumed homogeneous-reservoir properties, the effect of spacing on recovery at an economic limit of 2 B/D was very slight. Certain dimensionless groups can be used to extend the results to other fields having different penneabilities, spacings, reservoir thicknesses, well radii and porosities, so long as the PVT and relative permeability properties are similar to those used in this paper. INTRODUCTION An important aspect of reservoir engineering is the prediction of the performance of the reservoir based on the limited information normally available early in the life of a field. Usually, soon after a field has been discovered, it is necessary to decide upon the spacing to be employed and the production program to be used for most efficient utilization of reservoir energy. Where these decisions have not been dictated by legal or political considerations, they have frequently been based on experience factors which indicate that two fields with similar characteristics will probably have similar performances. To the extent that the industry has been successful in economically exploiting oil reservoirs, this rule-of-thumb method has been found to have merit. A large amount of theory has been developed with which it is possible to predict the performance of a reservoir by using certain known properties of the oil and the formation. Because of the difficulties of the mathematics involved for a drainage area having square or rectangular boundaries, the problem usually has been simplified and idealized by assuming each drainage area to have radial symmetry and the field to be represented by a number of these drainage areas. Muskatl was one of the first to formulate the theory for two-phase flow and to solve the equations for a few cases by numerical integration, using the relatively slow means of computing available at that time. His equations describing

By E. H. Lee, F. J. Fayers

Problems in reservoir analysis can usually he cupressed in terms of a system of nonlinear partial difjerential equations. A rnethod for .setting physically reasonable boundary conditions for these systerrls by ayplication of the theory of characteristics is discussed. A short description i.c. given of the technique for finding chnmcteristics of general first order and second orrler equations. Specific examples are quoted jrorn eqrrations occurring in theoretical reservoir studies. The princip1e. ossocicred wit11 the use of chnracteristic to set physically realistic boundary conditions me developed, and thr method is discussed in terms of boundary condition for the noncupillary, compressible, two-phase flow problen: jor the cupillary, ir~cowlpressible, two-phase flow problems for the tlrtergent flooding prohlenl; as for the compressible. noncapillary, two-phase flow pro11letn. INTRODUCTION In most problcms in reservoir analysis the physical situation is described by a system of nonlinear partial differential equations. The general solution to the set of equations usually comprises a wide range of functions ancl a particular solution is normally selectetl by applying suitable boundary conditions to the problem. In some cases it is clear from the physical situation what form the boundary conditions should take, but it is usually desirable to have an independent mathematical check of the suitability of the chosen conditions. From the mathematical standpoint, a system of partial differential equations with boundary conditions is said to be properly set if the combination determines a unique solution which is continuously dependent on the equation and the boundary conditions. Such a solution is sometimes termed "reasonable". In investigating the properties of the mathematical representation, it is useful to determine the structure of the system by the application of the theory of characteristics. This theory determines in the space of the independent variables, lines or surfaces which play a special role in thc system's behavior. The characteristics determine the possibility and nature of the propagation of discontinuities in the derivatives of the dependent variables. Along the characteristics certain differential relations hold which for some problems allow a so- lution to be obtained analytically. In the case of hyperbolic equations the characteristics can bc used as 3 basis for a numerical solution. If a complicated nonlinear system is to be solved, a study of behavior of the characteristics will indicate what boundary conditions should be chosen. By using the method of characteristics, several advances have rccently been made in solving reservoir engineering problems. Sheldon, Zondek and Cardwell' have discussed how the the method may be used grap-phically to solve the incompressible two-phase displaccmerit problem when the gravity term is included in the equations. Fayers and Perrine' have applied the method to study the equations representing detergent flooding. When compressibility and solution-gas effects are introduced into the two-phase flow problem, the equations and their characteristics become more complicated. It will be seen that in this problem the characteristics are useful for indicating a reasonable set of boundary conditions but cannot readily be used to compute a numerical solution. This problem has been solved on the IBM 701 using normal finite difference methods hy West, Garvin and Sheldon." The material to be prcscntcd in this report is divided into three sections. The first discusses the determination of the characteristics of first and second order partial differential equations; the second describes the selection of suitable boundary conditions by the application of the principle that the future cannot affect the past; the third discusses the use of characteristics for setting the boundary conditions of the two-phase solution-gas flow problem. A useful introduction to the method of characteristics has been given by Hildebrand.' A summary of the method will be given in this section of the paper. The First Order Equation The quasi-linear equation of first order is It is linear in the partial derivations of z. The properly linear equation has the form The characteristics of are curves along which the

By L. A. Weitzenkamp, N. M. Bashara

A study of the electrical conductivity of gold films less than 200 in thickness indicates a negative temperature coefficient of resistance and a thermal actiuatlon energy of less than 0.25 ev. The films consist of discrete metal islands of linear dimensions from a few angstroms to 1000A. Films were able to withstand applied fields up to 60,000 v per cm before breakdown occurred. A model for the conduction mechanism is proposed which involtes a Potential barrier due to work junction, irriage forces , and an electrostatic potential arising from separating an electron from a neutral island leaving behind a positively charged metal island. Activation energies of less than 0.25 ev ave consistent with this model, which includes a tunneling andor Thermal mechanism for conduction between the separate metal islands. An apparent satuvation effect in conductivity at high fields is observed experimentally. The saturation effect may be explained by the influence of the electrostatic contribution to the potential barriers for the transition through several islands adjacent to the original island which was activated to give up an electron. The proposed model gives Poor quantitative agyeement for the intermal electrical field between islands. ThIS paper will relate the structure of discontinuous metal films to the electrical properties and will emphasize the following topics: 1) the relation of the conduction mechanism to particle size and a comparison of particle size from electrical measurements with electron-microscope studies; 2) the saturation of the conductivity at high electric fields; and finally, 3) a comparison between the external field applied and the apparent field between particles. che- latter comparison points to a need for a modification in the model for the conduction mechanism first proposed by orter' and expanded by Neugebauer and webb2 which will be dealt with in the discussion. The paper is based on results on gold films on glass microscope slides. The clms were vacuum-evaporated and were about lOOA thick, thickness being measured by interferometry. It has been observed that metal films of this thickness do not exhibit bulk metallic conduction. Films thinner than some critical thiclmess, depending upon substrate temperature and other factors, have conductivities which increase with increasing temperature with thermal activation energies of less than 0.25 ev. Electron-microscope studies3 have shown that these films consist of discrete metal islands having linear dimnsions from a few angstroms to approximately 1000A. When the metal is being vacuum-deposited on the substrate, islands on the order of a few angstroms appear first. As the deposition continues the islands may tend to agglomerate into larger particles and may eventually grow together into a continuous film,3''' In our tests the conductivity of the gold films was found to be about ten orders of magnitude less than the bulk conductivity of gold. When oxygen was introduced into the test chamber the conductivity decreased by a factor of about one-third. This can be related to an increase in work function for gold caused by adsorption of oxygen, as reported by Oullet and Rideal.s Exposure to an atmosphere of water vapor caused an increase in conductivity by a factor of 2 or 3 when a low field was applied. When the films were exposed to water vapor under a high field the conductivity increased by as much as two orders of magnitude. The activation energies measured were much less than the work function for gold, 0.25 ev or less, with the smaller activation energies observed for thicker films. It will be shown that the activation energy is a function of particle size and spacing. In general the model for the potential between metal islands would include the effective work function, the image force, the applied field, and an electrostatic term. The electrostatic contribution arises from the potential required to separate an electron from a conducting sphere. According to Neugebauer and webb2 the maximum potential due to the field and the electrostatic term only is

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'chenko' 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' 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'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.' The presence of silicon in iron is said to reduce the activation energy for diffusion,' 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' 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 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'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?' 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.' 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,

By L. F. Mondolfo, B. E. Sundquist

The crystallographic orientation relationships resulting when lead is nucleated from the liquid by Ni, Cu, Ag, and Ge were determined. For each nucleating agent several definite orientatioz relationships were found. These relationships seemed to be controlled by good symmetry relations and low crystallographic disregistry between mating planes. For any given nucleating agent the under colling for nucleation was found fairly constant and independent of the orientation relationship and consequent disregistry. It was also found that, upon re melting and refreezing the Pb, the orientation relationship was changed. These findings prove that crystallographic disregistry is not the controlling factor in heterogeneous nucleation from the liquid. The results of this investigation tend to confirm the theory presented in a preceding paper that heterogeneous nucleation starts with the formation of an adsorbed layer of nucleated metal on the nucleat-ing impurity. Evidence is given that cavities in the nucleating agent act as centers of nucleation. IT has long been known' that solid extraneous particles are active in catalyzing phase transformations that occur in a system, particularly condensation and crystallization. It is well established that these heterogeneities act as catalysts by providing surfaces upon which nuclei of the precipitating phase can form with activation energies smaller than those required for homogeneous nucleation. Numerous investigations have shown that in this process of heterogeneous nucleation: a) the nucleus forms with one, or several, definite crystallographic orientation relationships with the nucleating phase2-4 and b) that there is a small range of undercoolings or super saturations characteristic of the nucleation of a given solid on a given Substrate.5-10 Turnbull and vonnegut11 have developed a theory based on theories developed by Volmer12 and Turn-bull and Fisher1= for heterogeneous nucleation from gases and liquids, that relates the super saturation or undercooling required for nucleation to the dis-registry between the lattices of the nucleus and the nucleating agent. This theory predicts that nucleation should occur with the orientation relationship between the nucleus and nucleating agent that minimizes the disregistry. Further, it predicts that the undercooling or super saturation necessary for nucleation should be a function of the disregistry. Numerous investigations have dealt with the orientation relationships resulting from the condensation of vapors onto crystalline solid substrates2,3 and a few with the nucleation of one phase by a second phase in solidification4,14. Others have dealt with the supersaturation8-10 and undercooling5-7 associated with nucleation in condensation and solidification respectively. However, there is virtually no report that gives both of these factors for the same system. In this investigation a study was made of the undercoolings and orientation relationships resulting when Pb is nucleated from the liquid by Ni, Cu, Ag, and Ge. It was the purpose of this investigation to check the Turnbull-Vonnegut theory, i.e., the importance of crystallographic disregistry between nucleating catalyst and nucleated metal. The results indicate that disregistry is not an important factor in nucleation and that the nucleation process is probably somewhat more complex than current theories suggest. EXPERMENTAL PROCEDURE Small single crystals of nickel, copper, silver, and germanium were prepared from materials of four to five nines purity, and the Pb used was also 99.999+ pet pure. Cu and Ag single crystals were prepared by sealing small chips of Cu or Ag in an evacuated quartz capsule and heating the capsule at 2000°F for 1 hr before cooling. Nickel crystals of 200 diam were also prepared in evacuated quartz capsules, but melting was done by heating the capsules in an oxy-acetylene flame for a few minutes. These spheres were invariably polycrystalline so

Jan 1, 1962

By F. P. Bullen

Empirical relationships between fracture stress, orientation angle, and diameter of crystal have been determined at 77°K. Orientation ranges of markedly different behavior were found—a law of constant normal stress&apos; of value a (diameter)-1/2 for the fracture of ductile crystals, and a condition of shear stress (or strain) for more brittle crystals. The observations are not consistent with current theories. An interpretation is advanced which is also applicable to observations on the effect of prestrain at room temperature on the subsequent fracture stress at 77°K and to the effect of cyclic stressing on the cleavage strength.&apos;&apos; The law of constunt normal stress&apos; and the brittle ductile transition are also explained. The interpretation is more consistent with the initiation of cracks by intersecting dislocations than with theories based on stress -concentration by dislocation arrays. ZINC single crystals are particularly suited to the study of cleavage because fracture occurs on the basal plane over a wide range of crystal orientations. Analysis of the conditions of stress and strain at fracture in crystals of different orientations should indicate which parameters control the cleavage process. Unfortunately, controversy has arisen over the correct empirical relationship between tensile fracture stress and orientation. Schmid&apos;s observations1,2 favored a &apos;law of constant normal stress&apos;, as observed in other materials.2 For zinc, however, the observed values are far below the theoretical strength and cannot represent the true limit of cohesion between neighboring atomic planes. Hence, the interpretation of such a &apos;law&apos; is not straightforward. Deruyttere and Greenough3,4 found a complex variation between tensile fracture stress and orientation; this variation did not agree with a &apos;law of constant normal stress&apos;. Two theories have been advanced to account for their observations: a) the propagation of cracks from low-angle boundaries,5 and b) the release of energy from piled-up dislocations during crack-propagation.4 The present work resolves the apparent discrepancy between the observations and shows that neither of the above theories are applicable to the tensile fracture of zinc single crystals. A phenomenological explanation, along the lines suggested by Gilman,&apos; is advanced and successfully applied to previously unexplained effects. EXPERIMENTAL DETAILS &apos;Crown Special&apos; redistilled zinc was used, except for one comparison series op tests using &apos;Tadanac&apos; electrolytic zinc. Crystals of 1 mm diam, subsequently called &apos;1 mm crystals&apos;, were grown from the melt in vacuo in precision-bore Pyrex tubes internally coated with graphite. Several specimens 1 in. long were cut from each crystal and chemically polished. Jigs were used to minimize handling strains, and crystals were mounted in the Polanyi machine the day prior to testing to allow recovery from any such strains. One-mm crystals were chemically polished for long periods to obtain 0.1 mm (approx) crystals. One-mm crystals were cemented into miniature gimbals by &apos;Araldite&apos; casting resin. The Appendix gives the reasons for using gimbals and the results obtained by other methods. More complete details of all techniques are given elsewhere.7 The symbols and terminology used are as follows: X = orientation angle (angle between tensile axis and line of greatest slope in the basal plane). T = tensile stress (on true cross-section) S,N = shear and normal stress (components of T with respect to the basal plane) ? = shear strain D = crystal diameter. The subscript &apos;f&apos; will be used to denote values at fracture. PART I-ANNEALED CRYSTALS EXPERIMENTAL OBSERVATIONS One-mm crystals were used to establish the variation of fracture stress at 77°K with orientation at fracture (Xf), Fig. 1. For 18 deg = Xf = 55 deg, a &apos;law of constant normal stress&apos; was observed. For Xf > 55 deg, the fracture condition approximated to a constant shear stress. At Xf< 18 deg, twinning occurred before fracture so that the results were not typical of homogeneous single crystals,4,8— such specimens will not be considered herein. The dependences of fracture stress upon Xf were of similar type for 6 mm,* 1 mm, and 0.1 mm crystals,

Jan 1, 1963

By O. D. Gonzalez, R. A. Oriani

A thermo-osmosis technique has been used to measure the heat of transport, Q* , of hydrogen and of deuterium dissolved in a iron and in nickel, and of hydrogen in Feo.6Nio.4 in the tempevature range 400° to 600°C. For all these systems, Q* is negative and has a large temperature coefficient; an isotope effect can be established for the solutes only in nickel. The magnitude of Q* is considerably larger than the activation energy, E , for migration in the case of these isotopes in a iron, so that all variants of the Wirtz model must be rejected. The phenomenological definition of Q* is developed to show that Q* is related to the mechanisms by which the activation energy is dissipated back into the lattice, and that Q*/E is a function of the ratios of the mean free paths of electrons and of phonons to the distance of jump of the diffusing atom. A correlation is shown to exist between the sign of Q* in thermal diffusion and that of the effective charge in electro migration, and may be understood as due to the large ratio of the mean free path of electrons to that of phonons. THE study of nonisothermal diffusion in the solid state is of interest because certain aspects of the mechanism of diffusion are made manifest which remain concealed though still present in the isothermal case. The thermodynamics of irreversible processes characterizes nonisothermal diffusion, known as thermal diffusion or as the Ludwig-Soret effect, by a quantity, Q*, called the heat of transport, which essentially measures the sign and the magnitude of the steady-state concentration gradient produced by the imposed temperature gradient. The problem then is to understand the quantity Q* at the atomistic level. There does not exist a theory for the atomistic interpretation of the heat of transport. There are, however, two classes of simple models which seem to afford some physical insight into Q*. Wirtz&apos;s model1, z and its subsequent generalizations partition the activation energy, E, for the migration of the atom spatially about the region of the migrating atom, so that Q* is the difference between the portion of the activation energy centered about the original site of the jumping atom and that portion about the arrival site. Clearly, the magnitude of Q* cannot on this model be larger than the activation energy for migration. The model of Oriani3 also considers the spatial distribution of E not only before but also after the atomic jump, and Q* is related to the variation of the spatial distribution of E that is produced by the jump. The simplest kind of system with which to compare the consequences of these ideas is the interstitial solid solution, since one may easily choose a temperature at which only the interstitial component moves, the solvent lattice serving as a completely satisfactory frame of reference. Thus, one avoids complications arising when both species in a binary system move. In addition, one avoids the necessity of knowing the energy for the formation of a vacancy, something which is needed for the analysis of thermal-diffusion data on metals which diffuse by the vacancy-exchange mechanism. However, the number of well-measured interstitial systems is extremely small, and, in particular, there are almost no data in the published literature for the temperature dependence of the heat of transport. If one works with solutes which are gases in the pure state at ordinary temperatures, one may use the technique of thermo-osmosis4 which permits thermal diffusion to be measured over a small temperature difference so that the temperature dependence of Q* can be more easily determined. The choice of dissolved hydrogen in iron and nickel was based on these considerations as well as on the desire to look for an isotope effect associated with localized

Jan 1, 1965

By J. A. Herbst, D. W. Fuerstenau

This paper examines the zero order production phenomenon in the context of the size discretized batch grinding model. A restrictive interrelationship between the selection and breakage parameters of the model, which is mathematically sufficient to ensure zero order behavior, is delineated. Based on this interrelationship, a scheme is developed for predicting the values of the selection and breakage parameters for all size fractions from a minimum of experimental data. To test this scheme, dolomite was ground in a laboratory batch ball mill. Successful simulation of the comminution behavior of this system was achieved using the batch grinding model with the parameter values obtained from the proposed scheme. Historically, a preponderance of comminution research has centered on attempts to relate breakage energy to the performance of comminution machines. Concomitantly, grinding data were interpreted almost exclusively in terms of inherently empirical energy-size reduction relationships1-3 or "laws of comminution"4- 6 which were based on highly oversimplified descriptions of the fracture process. In some instances, these relationships provide a crude basis for the correlation of experimental data, but, invariably, this approach is inadequate for meaningful process simulation. The control and optimum design of comminution circuits require a mathematical model capable of depicting the size reduction behavior of every size fraction for grinding conditions of technological importance. Energy-size relationships do not provide this detailed information. For example, an energy-size analysis for batch ball-milled dolomite has recently been completed in the authors&apos; laboratories.7 This analysis provided a satisfactory correlation between the hypothetical size modulus, which characterizes only the finer portion of the size distribution, and the input energy. However, this relation did not constitute an adequate description of the complete system. In general, one of two distinct viewpoints toward simulation has been adopted in recent comminution research. The first approach focuses on the fracture of single mineral specimens, with the essential aim of representing the over-all process in terms of the breakage characteristics of individual particles and the characteristics of the stress field which the particles experience within the particular size reduction device. As illustrated by Harris&apos;s review,&apos; the actual incorporation of single-specimen fracture information into a description of the behavior of a multiparticle comminution system has seldom been attempted. The only model to accomplish this incorporation was the one developed by Schönert.8 The Schonert model treats single-passage grinding machines in terms of a distribution of effective loads acting on a particle, a distribution of required particle breakage energies, and the breakage product distribution. However, the use of single-particle fracture behavior to derive models of the size reduction process is at present limited to machines in which the particle residence time is approximately equal to the time required to apply stress. As a result, Schönert&apos;s analysis cannot be used to advantage to simulate the complex environment which prevails within a tumbling mill. For multiple-passage systems a second approach to the description of comminution processes has been more fruitful. This approach entails the formulation of a mathematical model which is phenomenological in nature in that it lumps together the entire spectrum of stress-application events which prevail in a system under a given set of operating conditions. The appropriately defined average of these individual events is then considered to characterize the over-all breakage properties of the device. Thus, to analyze the performance of a tumbling mill, the manner in which the particles of a particular size (or size fraction) are stressed need not be distinguished. Instead, a single parameter is assumed to represent the resistance of that size to fracture, given the average grinding environment which exists in the mill. The isolation of such a parameter and a related set of quantities, which constitute the breakage product size distribution for the average event in this size fraction, allows the formulation of physically meaningful descriptive equations capable of yielding precise and detailed information for simulation.

Jan 1, 1969

By H. A. Shaw, H. G. G. Whitton

This paper describes the use of the electric-arc furnace for the production of tough-pitch, horizontal cast copper shapes and the production of copper anodes from tank house anode scrap. This installation is the only one of its kind in the United States. SINCE the initial use of electric-arc furnaces for refined copper melting on a commercial scale at Copper Cliff in 1936, this type of melting has been of extreme interest to everyone in any way connected with copper casting and fabricating. Early in 1948, Kennecott Copper Corp. decided to build an electrolytic copper refinery at Garfield, Utah. One of the major decisions to be made at that time was the type of furnace to be used for melting cathodes for casting into refined copper shapes and for melting anode scrap for the production of anodes. Based on the results obtained from arc furnaces by the Tnternational Nickel Co. of Canada, at Copper Cliff, and by the American Smelting and Refining Co., at Baltimore, coupled with certain important economic factors such as having an adequate amount of relatively cheap power available at the Magna Central Power Plant, high cost of poles and oil in the Salt Lake District, and the lack of natural gas at that time, the final decision was to install electric-arc furnaces rather than the standard type of rever-beratory furnaces. Experience has shown certain parts of the original installation to have been inadequate or in need of modification. Although the present status of arc-furnace melting at the Utah Refinery is far from perfect, nevertheless substantial progress has been made. At the Utah refinery, copper is melted in three 15-ft diam, size "M" Lectromelt furnaces, two of which are used for melting cathodes and one for remelting anode scrap. Each furnace is energized from a 13,800-225 v, 3-phase, 60-cycle, 6000-kva (at 40°C rise) oil-immersed water-cooled transformer. Auxiliary equipment for each refined copper furnace includes a charging machine, launders, ladle, 40 ft diam Walker-type casting wheel, bosh, and inspection conveyor. The anode furnace differs in that it is charged by hand from an elevated platform and the anodes are picked off the 24-ft diam casting wheel by a takeoff hoist and loaded directly into specially designed trailers, Figs. 1 and 2. Electrical Equipment Electric power to operate the arc furnaces is generated at the Kennecott-owned Magna Central Power Station two miles distant and transmitted to the Refinery outdoor substation at 44,000 v. At this point it is stepped down to 13,800 v and carried through concrete-lined underground tunnels to the switchgear and then to the primary side of the furnace transformers. The power is conducted from the Y-connected secondary bus bars of the transformer to the electrode arms by means of 22 flexible copper cable conductors per phase. Each cable has a cross-sectional area of 1,750,000 circular mils, giving a total area of 38,500,000 circular mils per phase. The power is carried from these cables to the three individual electrode holders by means of water-cooled copper bus tubing. The electrode clamps are incorporated in the electrode holders and are released by air pressure working against a spring in the clamp control cylinder. The movement of the electrode arms is either automatic or manual. The automatic control is used during melting and is regulated by three Westing-

Jan 1, 1954

By Charles L. Guettel

The temperature coefficients of electrical resistance of 31 alloys in the nickel-rich corner of the Ni-Cr-Fe system were determined. The results indicate that a range of binary Ni-Cr alloys has lower temperature coefficients than a specific region of Ni-Cr-Fe alloys which has been reported and quoted by others as providing the lowest coefficients in the ternary system. DESIRABLE properties of resistivity and temperature coefficients of resistance have been the basis for the successful commercial application of various Ni-Cr and Ni-Cr-Fe alloys. In particular, low or negligible temperature coefficients of resistance are necessary for wire used in precision wire wound resistors and potentiometers involved in closely controlled electrical and electronic circuitry. Alloys most successfully applied to this purpose basically have not been ternary Ni-Cr-Fe alloys but rather heat-treated Ni-Cr alloys in the region of approximately 73 to 79 pct Ni and 20 pct Cr with significant additions of other elements such as aluminum, silicon, iron, copper, manganese, and cobalt. While the addition elements play a vital role in the nature and control of the temperature-resistance curve of any one of these alloys, the inherent relatively low temperature coefficient of the binary Ni-Cr base composition itself in the approximate region of 80 pct Ni and 20 pct Cr contributes fundamentally to the eventual negligible temperature coefficient obtained with the applied alloys. W. A. Dean 1 has reported mean temperature coefficients (20&apos; to 100oC) of as-cast samples for the Ni-Cr-Fe system, the nickel corner of which is given in Fig. 1. Dean had to use as-cast samples because alloys from a considerable region of the system were unworkable. Dean refers to the work of Hunter and sebast2 who used annealed wire samples to show that, for chromium contents of 10 to 25 pct, the temperature coefficient increases as the iron content is increased at the expense of nickel. Although the data of Hunter and Sebast were limited in the nickel-rich areas to eight samples of 60 to 91 pct Ni with only four samples containing iron, the possibility that no Ni-Cr-Fe alloy has a lower temperature coefficient than certain binary Ni-Cr alloys should have been given serious consideration. However, Dean&apos;s work has been referred to by J. S. Marsh,3 who, quoting Dean, states that for the temperature coefficients "the smallest values-about 0.0001 per oC (100 ppm)—were found for alloys containing about 55 pct Ni and 20 pct Cr" (balance iron). A region such as that reported by Dean would be significant to those trying to develop alloys for the applications described. Therefore, a series of three-pound laboratory open-air induction melts of Ni-Cr and Ni-Cr-Fe alloys were melted, conventionally processed to 0.005 in. diam wire, and continuously strand annealed at 1900° F. The basic Ni-Cr compositions, with balances essentially iron other than usual minor impurities, and the corresponding temperature coefficients of resistance and resistivities are given in Table I. The resistivities are approximate and are not to be considered as a basis for relating them to their position of composition in the ternary system. The wire diameters were determined with a micrometer, and not by weight and density measurements which would be required for accurate resistivity values. However,

Jan 1, 1963