Search Documents

By P. M. Kelly, E. P. Butler

The stability of a 70/30 Mn-Cu alloy aged to peak damping has been investigated using electron microscopy, X-ray diffraction, and torsional pendulum measurements. Storage at room temperature or at 100" C causes a reduction in damping capacity. This is shown to be due to the formation of subdomains which "lock" the main domain structure and so hinder their movement. After storage the damping capacity can be restored by a retransformation treatment, which generates structures identical to those found in newly aged specimens. Plastic deformation by rolling also causes a decrease in damping capacity. In this case the damping capacity can only be restored by retransformation if the amount of deformation is less than about 5 pct. THE general metallography and structure of these high damping capacity alloys has been discussed in Part I of this work.1 Part II deals with some of the variables that affect the damping capacity. In particular, the stability of the damping capacity of a 70/30 alloy aged for 2 hr at 425° C and quenched has been investigated. The temperature sensitivity of the damping capacity of a 70/30 alloy aged for 2 hr at 450° C has been reported by Birchon.2 The damping capacity begins to decrease with increasing temperature, and is very low above the M, temperature of 125°C. The cubic — tetragonal transformation is completely reversible and there is no apparent thermal hysteresis associated with the change in damping capacity. Another important effect was first reported by Dean and his coworkers.3 They noticed a decrease in the damping capacity with time for specimens of a quenched tetragonal alloy. More recently workers at the Admiralty Materials Laboratory 4 have studied the stability of an aged 70/30 alloy. They found that the damping capacity decreased to half its original value after about 1000 hr at room temperature. EXPERIMENTAL The techniques used for electron microscopy and X-ray diffraction have been described in Part I. The damping capacity was measured by means of a torsional pendulum, using rectangular specimens, 3 by 0.3 mm with a 6 cm gage length. The specific damping capacity, S.D.C., was calculated as a percent amplitude decay in free vibrations over 10 cycles, starting at an amplitude which corresponded to a 700 psi surface shear stress. Because of the stress dependence of the damping,' the values obtained were about a third those reported for a 5000 psi surface shear stress.' RESULTS a) Effect of Aging Temperature. The effect of aging temperature on the damping capacity is shown in Fig. 1. The damping capacity reaches a peak of 12 pct after aging at 400° to 425°C) then begins to decrease as the volume fraction of a Mn increases. At higher aging temperatures the structure becomes cubic again with a low damping capacity. b) Stability of the Aged Structure. The damping capacity of the aged alloy decreases with time at room temperature. This process could be accelerated by storing at 100°C, when about 100 hr were required for the damping to fall to half its original value. The results of room temperature and 100°C storage, for times up to 1000 hr, are shown in Fig. 2. The damping capacity starts to decrease after about 24 hr at 18° C and 1 to 2 hr at 100°C. It was found, however, that the damping could be almost completely restored by retransformation. The specimens were held at a temperature above the M, for a short time, 250° C for 5 min, and then quenched. The values of damping capacity for specimens that were aged, stored, and then retransformed are shown by the full points in Fig. 2. X-ray measurements of the c/a ratio and hardness readings on the specimens did

Jan 1, 1969

By D. de Fontaine

A convolution procedure is described for obtaining an analytical expression for the composition profiles of solute concentration:; in binary solid solutions. The free energy of solutions containing such compositional inhomogeneities can then be evaluated and the stability of spherically symmetric clusters (G.P. zones, for short) us a function of their radii and central compositions can be studied accordingly. Two characteristic zone radii are defined: the minimum radius (RM) corresponding to the smallest zone for which growth is energetically favored over dissolution, and the critical radius RK) corresponding to zones associated with the minimum work of formation. The values of RM and RK coincide at the two classical limits of nucleation and growth and spinodal decomposition. As the spinodal is approached, however, RM and RK diverge increasingly and it is proposed that the initial] decomposition morphology is governed more closely by the minimum, rather than by the critical radius since the former possesses no singularity within the m eta stable phase boundary. Reversion is briefly discussed in the light of the foregoing and it is shown that a set of reversion temperatures located between the coherent miscibility gap and the coherent spinodal must correspond to a given distribution of low temperature zones. SINCE the early works of Guinier1 and preston,2 the kinetics and morphology of GP zones have been studied extensively, see Ref. 3, for example, for a recent review. Nevertheless, a completely satisfying theoretical treatment of the decomposition of meta-stable or unstable binary solid solutions is still lacking due to the almost unsurmountable mathematical difficulties that one encounters in attempting to solve the pertinent diffusion equation which, in principle, describes the whole process. This equation, first derived by Cahn for the initial stages of spinodal decomposition,4 is a nonlinear partial differential equation of the fourth order in anisotropic space, a general analytical solution of which is completely out of the question. From this equation, some prediction can be made, however, concerning the later stages of spinodal decomposition5 and a recent numerical solution6 of the equation in the case of one- and two-dimensional composition modulations of a solid solution inside the spinodal has led to the following considerations: the sinusoidal morphology7 is to be expected only for solutions of average composition c, located near the center of the coherent miscibility gap.8 For solutions with co closer to the spinodal decomposition, depleted shell zones9 should constitute the first recognizable decomposition product, later evolving to iso- lated, enriched preprecipitates in a uniformly depleted matrix as this must, in all cases, constitute the final metastable morphology, barring, of course, loss of coherency and change in crystal structure. In agreement with the qualitative description of Bonfiglioli and Guinier,l0 it thus appears that the various models discussed by previous authors are, in fact, idealized successive stages of the general decomposition process. Unfortunately, corresponding calculations for a solid solution outside the spinodal are not yet available because of two additional difficulties: it is essential that the diffusion equation be solved in three-dimensional space, and furthermore, one must include in the equation an appropriate "source term" describing the statistics of the creation of composition fluctuations required to trigger the nucleation events. At present, a general treatment of this problem does not appear to be tractable. Nevertheless, by using a modelistic approach, i.e., by selecting a priori certain composition variation profiles and by evaluating their associated free energies, one can draw certain conclusions relative to the stability of these isolated composition modulations or "zones", both outside and inside the spinodal. This is the object of the present investigation. THEORY We shall now describe a standard four-parameter spherically symmetric GP zone, evaluate its free energy according to the Cahn and Hilliard11 formula (Eq. [5], below), and plot the corresponding energy surfaces as a function of zone amplitude and radius as

Jan 1, 1970

By G. W. Healy, W. Craft, D. C. Hilty

By means of a theoretical extension of the Cr-C temperature relation in molten chromium steels to low chromium contents and by a correlation of the ratios of chromium to iron in the slag and metal, a method has been developed for estimating the amount of metallic oxidation during the oxidizing period of a chromium steel heat. Application of this method has indicated that due to temperature limitations metallic oxidation may be excessive for very low carbon steels charged with more than a small amount of chromium. IN the melting and refining of chromium steels the oxidation of carbon to a low level is accompanied by the oxidation of chromium and iron in considerable amounts. For the subsequent recovery of chromium and associated iron from the slag, the amounts of metal oxidized must be known for efficient reduction and control of composition. Further, in order to arrive at an accurate understanding of the process, so that chromium-bearing scrap can be used effectively, information is required on the relation between composition of the charge and weight of metal oxidized to reach the desired carbon content. Crafts and Rassbach1 recently developed an empirical relation between chromium, iron, and manganese oxidized per ton of steel charged, and final carbon content and temperature. This relation did not show any dependence on amount of chromium charged. However, it was felt that such a dependence might be present, so a further analysis of the data was undertaken. Cr-C Relation at Low Chromium Levels Crafts and Rassbach's results were based on temperatures estimated from the Cr-C relation previously established experimentally by Hilty,2 and on a single temperature observation (immersion thermocouple) on a heat made at low temperature from a charge of carbon steel scrap containing only a small amount of residual chromium. The temperatures estimated from the Cr-C relation were for heats charged with chromium-bearing scrap and which contained substantial amounts of chromium (2 to 10 pct) at the end of the oxidizing period. No data, other than that from the low temperature heat just mentioned, were available for the low and intermediate chromium ranges. In order to extend the range in which temperature may be estimated, the Cr-C relation requires further consideration. The Cr-C relation, expressed as suggests that carbon is zero at zero chromium in the bath. This is obviously absurd, because as chromium approaches zero the carbon content must approach the limit established by the Fe-C-O equilibrium. It is evident, therefore, that the Cr-C relation as defined by Eq. 1 is not valid below some minimum chromium content that may vary with temperature. On the other hand, it is probable that any such limiting chromium content is less than 4 pct at 3200°F, since the experimental data from which the Cr-C relation was derived included several observations at that level with no deviation. In any event, it is apparent that further evaluation of the results presented by Crafts and Rassbach is dependent upon a means for estimating the Cr-C temperature relation in the low chromium region. From the observations of Chen and Chipman9 and thermodynamic data given in Basic Open Hearth Steelmaking,4 the following relation can be calculated for melts containing small amounts of chromium and carbon assuming a carbon monoxide pressure of 1 atm: % Cr -21,250 K' = (%C)2; logK1 == —t— +13.88 [2] Moreover, by plotting Eq. 2 for any given temperature on cartesian coordinates, it can be extrapolated graphically to the limiting carbon content at zero chromium for the specified temperature without serious complications. By means of Eq. 2 and its extrapolation, the Cr-C temperature relations originally established experi-

Jan 1, 1954

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 Robert L. Davidoff

INTRODUCTION The impact of taxation differs substantially from state to state in a complex fashion that depends on both the physical characteristics of a mineral deposit and on the price of the mineral commodity relative to its costs of production. A comparative analysis of the tax structures of eight states (Arizona, Colorado, Idaho, Montana, Nevada, New Mexico, Utah and Wisconsin) showed wide variations in rates of return for any given type of mineral property, with no state having either the most burdensome or least burdensome tax structure under all circumstances. Three hypothetical mineral operations were analyzed: a producing open-pit copper property, a non-producing underground platinum operation, and a non-producing underground gold operation. None of them corresponds exactly to an existing property. Each of these property types is further defined to have ore feed grades at three different levels, leading to three potential profitability levels, referred to in this paper as "economic", "marginal" and "subeconomic". Thus the differential impact of state tax structures is measured for mineral properties of different types, as well as for different profitability levels. Montana is used as the base case. This study does not address the topic of corporate structures or tax consolidation. There are no costs assumed external to the operation of the deposit (e.g. exploration to locate additional deposits). Two types of analyses were performed for this study. The first calculates the total cumulative tax payments in each state using all applicable taxes (including U.S. federal income taxes and examines the resultant differences in discounted cash-flows rate of returns DCFROR's amoung the eight states. The second is a sensitivity analysis that examines each type of state tax separately (i.e. state income, property, and severance taxes) in order to isolate and examine the differences in tax payments that are due to the difference in rates and methods of computing specific types of taxes in each state. The results of this study serve to illustrate the complex nature of interactions among various taxes. For example, property tax is a deduction in calculating severance tax in some state. Changes in property-tax and severance-tax payments both enter into the calculation of state and federal income-tax liability. This, in turn, could change the timing and usage of tax-loss carry forward amounts and investment tax credits. Additionally, the physical characteristics of a property, the method and costs at which the mineral will be extracted, and the profitability level at which each property operates will affect tax treatment. Each mineral operation was evaluated at three different predetermined profitability levels, referred to using the terms economic, marginal and sub-economic. These three profitability levels are defined by the ore feed grades of the primary commodity that yield in constant terms, approximately an "economic" 18% DCFROR, "a marginal." 10% DCFROR, and a "sub-economic" 5% DCFROR, respectively, under the Montana tax structure. In the analyses, the mining and processing parameters for each of the three types of properties were assumed to be identical for each of the eight states. It was assumed that each operation produces at full operating capacity throughout its life. This implies the level of mineral commodity demand will support the inferred commodity price, such that each operation will be able to sell all its output. The analyses performed for this study were done using a comprehensive economic evaluation simulator that is used in the Bureau of Mines, Division of Minerals Availability MAS program to perform DCFROR analyses. All economic evaluations performed upon these data (see appendix 2) were in constant January 1981 dollars. By-product prices used for the analyses were: [Gold $425.00/troy ounce Palladium$200.00/troy ounce Silver $ 10.00/troy ounce] Typical cash-flow calculation layout is given in Appendix 3. The initial step for each evaluation is to define an "economic" profitability level scenario. An 18% DCFROR rate in Montana was considered the

Jan 1, 1985

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. C. Knopf, I. M. LeBaron

In the Carlsbad area potash is dry-mined and wet-concentrated. Wet concentration involves recircu-lation of saturated brines, with resultant difficulties of brine disposal and inherent losses in recovery through solution of the KCl. A dry beneficiation process for potash minerals would offer many advantages over present methods. International Minerals & Chemical Corp. has been working on development of a dry process for some time. This article is a progress report. Potash ores available in this country are largely mixtures of halite and sylvite contaminated with slimes and often with many sulfate and magnesium minerals. Electrostatic concentration of KCl when the contaminating minerals are in low percentage has been accomplished. Continuous production of better than 60 pct K2O concentrate from high potassium sulfate ores has not been demonstrated at the present time. Laboratory Procedure: Laboratory work on electrostatic beneficiation of potash minerals has been carried out using both pure crystals and mined ore. Electrostatic separation of the two main components of potash ore, sylvite (KCl) and halite (NaCl), is accomplished by first heat-treating the ore to allow for electric charging of the particles. Preferential exchange of' electrical charges between different mineral particles through surface contact is the principal mechanism involved. After proper conditioning of the particle surfaces, the oppositely charged minerals are dropped vertically into a horizontal electric field having a field gradient from about 5000 to 15,000 v per in. The particles are free-falling in the field. A suitable power supply is used to produce the proper total voltage for the electric field. Total required voltage is, of course, a function of electrode spacing. A typical electrode arrangement with power supply is shown in Fig. 1. In general, the total voltage on any system is 75,000 to 90,000 v. Power consumption is negligible, as leakage at insulators and corona discharge are the main losses in the system. Optimum preheat temperatures for pure material components generally are not the same as for mined minerals. Fig. 2 shows clearly that maximum recovery on a 50 pct pure KCl, 50 pct pure NaCl mixture occurs at a preheat temperature of about 600°F. In this test the preheat mixture was allowed to cool to about 225°F before dropping through the electric field. The optimum temperature for sep- aration is practically never the optimum temperature of preheat. Optimum separation for mined Carlsbad potash ores occurs when a preheat temperature closer to 900°F is used, as shown by Fig. 3, which plots percent sylvite (KCl) in the concentrate and tails as a function of preheat temperatures. Here again the ore was allowed to cool to about 225°F before separation. It should be noted that efficiency of separation increases continuously with increasing preheat temperatures until incipient fusion of the ore takes place at about 950°F. The increase in optimum preheat temperature for Carlsbad ores over the pure minerals can be traced to the effect of the clay slimes intimately associated with the mined potash. If the slimes associated with the ore are removed by washing with alcohol or other non-aqueous material, the necessity of preheating to the high temperature is removed. This temperature range has been proved critical within these limits during continuous pilot plant operation. Concentrate and tailings grade and optimum recovery also are functions of the temperature of the potash mineral particles during separation in the electric field. Fig. 4 presents data indicating that feed to the rougher pass should range between 150°F and 275°F for best results. It has been found by subsequent pilot plant work that these temperatures can be maintained in actual operation.

Jan 1, 1959

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 B. H. Bergstrom, J. J. Gilvarry

Previously the authors showed that the Gilvarry equation correctly describes the distribution of fragment size in single fracture, provided the exoclastic particles showing original surface of the specimen are removed from the distribution. In the present work this conclusion has been strengthened in two respects. First, it has been shown that the existence of two peaks in the differential distribution for the endoclastic particles, corresponding to the effect of surface and volume flaws in the interior of the specimen, can be demonstrated from the average of data for a large number of specimens. Second, it has been shown by means of a Coulter counter that data obtained from specimens fractured in gelatin display the behavior at fragment sizes approaching 1 p which is demanded by the limiting form of the Gilvarry equation in this case. In previous work, a rigorous derivation of the proper distribution function for fragment size in single fracture of a brittle solid has been given by Gilvarry,' based on a closely defined physical model and deduced strictly by the laws of probability. Experimental work has been reported by Gilvarry and Berg-str~m,~ confirming the theory in its broad outlines. Subsequently, Gilvarry and Bergstrom extended the theoretical discussion to include comminution and presented experimental data bearing on this point. They have described further experimental results for the case of single fracture recently. The experimental evidence regarding single fracture was rather convincing. However, it was based on only five specimens for which fragment sizes were determined by sieve analyses and three for which measurements in the finer sizes were made on a Coulter counter. The purpose of the present paper is to extend the experimental results as obtained by means of both sieve analyses and the Coulter technique. In the former case, sufficient specimens are included so that the average data should vield ac- curately the detailed behavior of the distribution function for the coarser fragment sizes. In the latter case, measurements for the finer sizes are reported which correspond to single fracture more closely than prior results and fill a lacuna left in the experimental verification. In addition, an extension of the theory given recently bGilvarry for the case of fracture of a two-dimensional solid will be outlined, and experimental results obtained by him will be discussed.6 THREE-DIMENSIONAL CASE Gilvarry' has shown that the Schuhmann size modulus is a measure, in a sense specified later, of the number of activated edge flaws in a specimen undergoing single fracture. Single fracture is defined as that resulting from a single instantaneous application of energy, such as the blow of a hammer whose motion is halted when fracture begins. Plural fracture consists of a short sequence of single fractures, as in a jaw or gyratory crusher; multiple fracture obtains when many applications of energy are made, as in a ball mill. The Griffith flaws associated with fracture of a specimen are flaws occurring within the volume of the body, on the new surfaces created within its interior and on the new edges (intersections of two surfaces) within the body. Fracture proceeds by activation of flaws in the volume of the specimen, in the fracture surfaces throughout its interior, and in the edges produced by intersection of fracture surfaces. If all the flaws and all groups of flaws are distributed at random, independently of other flaws or groups thereof, it has been shown by Gilvarry that the distribution function for fragment size has the form of a Poisson distribution. The fraction y, by weight or volume, passing a sieve of mesh size x is given by where k is an average spacing of activated edge flaws, j2 is an average amount of surface containing one activated surface flaw and i3 is an average volume corresponding to an activated volume flaw. The averages in question are referred to the mean dimension x of a fragment. On the basis of the theoretical

Jan 1, 1962

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 S. E. Bronisz, Dana L. Douglass

a Uranium was deformed by cold rolling, and the effects of this plustic deformation on the microstruc-ture of the metal were observed by the technique of transmission elecbon microscopy. The recrystalli-zation of selected cold-wmked samples was also studied by this method. The samples were prepared for examination by electrolytic thinning in a standard metallographic electrolyte which has not been previously reported as being used in producing thin foils of uranium. The microstructural changes observed in the a uranium were generally similar to those seen by optical microscopy. The changes in the microstructure of a uranium following cold work and recrystallization treatments have been the subjects of several studies employing optical and X-ray techniques.12-19 These studies indicated that twinning was the redominant mode of deformation at room temperature20 and that the regions of high energy associated with the deformation twins were preferred sites for the initiation of recrystallization.15,18, Recently several workers have prepared thin foils of uranium suitable for examination by transmission electron microscopy.'-lo The foils thus prepared have been used to study fission fragment damage; a particular dislocation reaction,' and gas precipitation." At this laboratory the microstructural changes occurring during the cold working and recrystallizing of a uranium samples have been studied by using transmission electron-microscopy techniques. The foils used in this work were electrolytically thinned in a polishing bath not reported previously in the literature as being used for this purpose. EXPERIMENTAL PROCEDURE The uranium used in this investigation was depleted a uranium sheet 4 mils thick. The results of a spectrochemical analysis of this material are shown in Table I. The as-received uranium was cut into strips, 1 in. wide by 5 in. long, and annealed in vacuo at 600°C for 2-1/2 hr. The resulting microstructure consisted of relatively uniform, equiaxed grains having diameters of 10 to 50 p. Cold Working. The cold-worked samples were prepared by rolling several of the annealed uranium strips to various thicknesses. The strips were reduced 3.5 ± 0.5, 6.0 ± 0.6, 10.0 ± 0.5, 19.5 ± 1.0, and 28.5 ± 2.0 pct in thickness. It was not possible to achieve these reductions with a single pass through the rolling mill, so the strips were turned end for end during the rolling procedure. The cold-rolled strips were stored in dry ice until samples from them could be thinned and examined (within 24 hr after rolling). Annealing. The samples used in studying the re-crystallization behavior of a uranium were obtained by annealing sections of the strips which had been reduced 10.0 and 28.5 pct. These samples were sealed in evacuated Pyrex capsules and were annealed in a tube furnace at temperatures from 200" to 600°C and for times from 1/2 to 2 hr, see Table l. Thinning: The samples from the rolling and heat treatments were prepared for examination in the electron microscope by electropolishing in a bath composed of 25 g chromic acid, 133 ml glacial acetic acid, and 25 ml water. This electrolyte has been used and recommended for preparing specimens for optical metallography It was necessary to mask the edges of the sample with a stop-off paint in order to prevent preferential attack in these areas. The paint used was made by dissolving polystyrene in benzene. The electrolytic cell consisted of a cylindrical stainless-steel cathode placed in a beaker containing the electrolyte, with the sample as the anode positioned at the cell axis. The temperature of the electrolyte was held below 15°C by cooling with liquid

Jan 1, 1963

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 Henry A. Wriedt, John Chipman

Equilibrium in the reaction of hydrogen gas with oxygen in liquid nickel, iron, and their alloys has been studied at temperatures of 1500° to 1700°C. The equilibrium con^stant, 0/p, [% O], is greater in nickel then in iron by a factor of 45 at 1594 C. In the alloy steel range, the activity coefficient of oxygen is slightly increased by the presence of nickel. NICKEL is more resistant to oxidation than is iron, in both the solid and the liquid state. The solubility of oxygen,' however, is greater in nickel than in iron. The two metals and their liquid alloys constitute an interesting series of solvents in which the chemical behavior of oxygen can be studied by methods similar to those used for liquid iron.' Experimental Method The furnace assembly, including the tube for preheating the entrant gas, was the same as that described by Gokcen and Chipman. " controlled mixture of water vapor, hydrogen, and argon preheated to approximately the melt temperature impinged upon the surface of the liquid metal. The charge, weighing about 35 g, was held in an alumina crucible at constant temperature in the gas stream for a prolonged period, then quenched in a stream of cold helium, sampled, and analyzed for oxygen by the vacuum fusion method. The temperature of the liquid metal was measured by a Leeds and Northrup optical pyrometer which had been checked against a similar instrument certified by the National Bureau of Standards. The freezing point of electrolytic iron in hydrogen, taken as 153j°C, was used for frequent checks. This point and the emissivity determined by Dastur and Gokcen" established the transmissivity of the optical system. The emissivities of Fe-Ni alloys at 1535" are known from the experiments of Smith and Chipman,h nd it was assumed that the variation with temperature is parallel with that of iron. Recorded temperatures are probably accurate within k5"C. Preparation of Gas Mixtures—Since the ratio H,O:H, had to be much higher for nickel than for iron, it was found convenient to construct two separate systems for gas preparation. In both cases, the ratio of argon to the sum of the other two gases was at least 4, the proportion of hydrogen and water vapor being varied to produce H,O:H, ratios from 0.2 to 70. The primary purpose of the argon was to diminish the errors caused by thermal separation' in the gas phase. In addition, its presence lowered the saturation temperature required for a given H,O:H, ratio and minimized the difficulty of preventing condensation. Both systems were modifications of the one used by Dastur and Chipman.' For H,O:HI ratios greater than 2, welding-grade argon purified by passage over magnesium chips at 620°C was mixed with hydrogen produced in the electrolytic cell shown in Fig. 1. The design was based on a cell used by Bodenstein and Pohl," the important features of which area nickel wire cathode and sodium hydroxide solution as the electrolyte. The rate of hydrogen evolution was controlled by a calibrated ammeter and a rheostat in the circuit, the current efficiency being taken as 100 pct. Precise control of the argon flow rate at about 300 ml per min was achieved by means of a capillary flowmeter having a controlled pressure drop. The molar flow rate was measured volumetrically before each run, and a constant flowmeter setting was maintained thereafter. The A-H, mixture passed through a water saturator of the type previously described," the saturation temperature being maintained automatically with a fluctuation of k0.04"C. The tube between saturator and furnace was heated by a resistance winding. The system was of all Pyrex construction with the exception of the stainless steel tube containing magnesium chips. This was connected with the glass tubing at the upstream end by

Jan 1, 1957

By R. E. Goodman, M. A. Mahtab

The, state of stress around a vertical wellbore in rock following nonlinear stress-strain laws is examined by means of finite element analysis. The wellbore, is considered an axisymmetrier body with axisymmetric loading. The initial vertcal and borizontal stresses are "locked" in the rock elements around the wellbore and a new state of stress is generated by the displacements which occur around the borehole. A point-uise variation of the elastic moduli is made on the basis of the new stress state and the triaxial data. The initial stresses tire now reintroduced along with tile changed moduli and original boundary constraints. This procedure is repeatrd until convergent stresses are reached. The effect of nonlinearity on streases is examined for a 6,000-ft wellbore in a scbistose gnezss and Berea sandstone using results of laborratory triaxial compression tests. The results sbou that the effect is restrited to one well radius from the bottom periphery of the bole. Beyond a distance of one-quarter radius, the elfect of nonlinearity on stresscs IS almost always less than percent for the cases conszdered, The consideration of a static pressure inside the well does not magnify the effect of monlinearity on borehole stresscs. INTRODUCTION The terms " wellbore" and "borehole" here designate cylindrical openings in the ground with vertical axis and a circular cross-section. A knowledge of the stress redistribution that occurs on excavating a wellbore is important in understanding (1) the behavior of the lined or unlined hole, (2) hydraulic fracture response, and (3) the effect of stress redistribution on drillability; also it is important in predicting initial stresses in the virgin ground, and in analyzing the response of measuring instruments placed in the borehole. Our knowledge of the state of stress around a wellbore has been restricted to homogeneous, isotropic, elastic material and derives chiefly from the analysis by Miles and Toppingl and the photoelastic work of Galle and Wilhoit2 and Word and Wilhoit.3 In this investigation the state of stress is examined for a nonlinear elastic material by means of finite element analysis. Many rocks possess stress-strain curves that depart notably from straight lines in their initial or final portions. While the literature contains abundant stress-strain data from triaxial tests (axisymmetric loading) on cylindrical rock specimens, there is little information on rock deformability under nonaxisymmetric loading conditions such as occur at each point around the bottom of a wellbore. Although there is some knowledge of the effect of intermediate principal stress on rock strength, there is virtually nothing known about its effect on rock deformability; therefore, we have assumed here that the effect of intermediate principal stress can be ignored. A schistose gneiss4 and Berea sandstone5 were selected as representative rocks for this analysis. The traditional graphs of deviator stress (ol -o3) vs axial strain were reworked to give the tangent modulus as a function of the deviator stress for varying values of the minor principal stress. The result is a nesting family of skewed, bell-shaped curves for tile gneiss (Fig. 1A) and the sandstone (Fig. 2A). A similar replotting of the lateral strain data defines the variation of Poisson's ratio (v) with the deviator stress and confining pressure. These curves, shown in Fig. 1B for the gneiss and in Fig. 2B for the sandstone, are not so well ordered as the tangent modulus curves. However, all of these display an increase of v with deviator stress application, but the rate of increase diminishes with confinement. The ET and v curves for the two rock types are tabulated in Tables 1A and 1A for use in a digital computer so that material properties corresponding to a given state of stress can be assigned by interpolation. STATEMENT OF THE PROBLEM Given a state of stress in the ground prior to drilling a hole (initial stresses), we assume that

Jan 1, 1969