Search Documents

By L. A. Wilson, P. J. Root

The relative costs of heating a reservoir by steam injection and by combustion have been examined. The comparison was based on a model similar to that proposed by Chu.' The cost of boiler feed water, the price of fuel, pressure and plant capacity were parameters in determin-ing the costs of air compression and steam generation. The analyses compare the cost of heating to the same radius by the two methods. Results suggest that the two primary factors for comparison are the price of fuel and the amount of crude burned during underground combustion. The cost of fuel has a greater effect on the cost of heat from steam than it does on its cost by combustion. As a result, analyses indicate that when the price of fuel is low, steam may be unequivocally cheaper than air. The influence of heat loss is such, however, that as the heated radius increases combustion becomes relatively more competitive depending upon the amount of crude burned. This implies that steam may be cheaper for small stimulation jobs (huff and puff) but combustion may be more economically attractive for heating large areas (flooding). INTRODUCTION Use of thermal methods of recovery is an accepted fact today. After an induction period of several years, processes are being widely used that involve reservoir heating to augment recovery. Of the several techniques, steam injection and forward combustion appear to be destined to dominate the field. Although the objectives of both are the same, the basic differences between generating heat in situ and injecting heat after surface generation influence the cost in different ways. This study compares the cost of heating a reservoir either by steam injection or by forward combustion. There has been no consideration of recovery. Presumably, recovery from the swept region would be high in either case. The sole consideration was the cost of heating to the same radial distance by either process. PROCEDURE THE MODEL The basis for comparison was a mathematical model similar to that used by Chu' for combustion. The model simulates a radial heat wave in two-dimensional cylindricaI coordinates. It includes heat generation, conduction and convection within the reservoir and conduction in the bounding formations. Thus, heat losses from the formation are considered. Three significant modifications were made. 1. Equal logarithmic increments rather than equal increments were used for the mesh spacing in the r direction. By this technique large distances were simulated with relatively few mesh spaces. 2. A backward difference approximation to the convection term was used to avoid troublesome oscillations which result from a central difference approximation when the convection term is large. 3. The radial increments of the combustion zone motion were not necessarily uniquely related to the mesh configuration. The cumbersome step function introduced by the heat of vaporization of steam was circumvented by assuming the enthalpy of the steam to be a linear function of temperature between reservoir temperature and steam temperature. This is equivalent to assuming an average heat capacity numerically equal to the difference between the enthalpy of saturated steam and the enthalpy of water at reservoir temperature divided by the difference between the two temperatures. Heat losses obtained by this model are in essential agreement with those obtained by the analytical solution of Rubenshtein.' A detailed description of the model is presented in the Appendix. Using the model, the times required to heat to particular radial distances were obtained as a function of injection rate and other physical parameters. For the steam case, injected fluid was assumed to be saturated steam at pressures of either 500, 1,000 or 1,500 psia. The corresponding temperatures are 467, 544 and 596F, respectively. Thickness ranged from 10 to 50 ft and injection rate ranged from 100,000 to 1 million Ib/D. Reservoir and overburden temperatures at the injection well were assumed to be that of saturated steam at the injection pressure. The effect of maintaining the overburden temperature at the well at a different temperature (initial reservoir temperature) was examined with no significant change in behavior. The influence of wellbore heat losses for the steam case was determined in the following manner. The rates of heat loss as a function of time were estimated using an approach similar to that suggested by Ramey." he data were based on injection through 2%-in. tubing in 7-in. casing. Integration of these data over the entire iniection period yielded the total heat loss. Total heat losses were then corrected to their equivalents in steam (this number resulted from dividing the total heat loss by the latent heat). This was considered additional steam required to accomplish the reservoir heating and the total cost was increased accordingly.

Jan 1, 1967

By W. B. Gogarty, W. C. Tosch

A new recovery process for producing oil under both secondary and tertiary conditions utilizes the unique properties of micellar solutions (also known as microemulsions, swollen micelles, and soluble oils). These solutions, which displace 100 percent of the oil in the reservoir contacted, can be driven through the reservoir with water and are stable in the presence of reservoir water and rock. Basic components of micellar solutions are surfactant, hydrocarbon and water. They may also contain small amounts of electrolytes and co surfactants such as a1cohol.r. The specific reservoir application dictates the type and concentration of each component. A salient feature of [he process is the capability for mobility control. Micellar solution slug mobility, by way of viscosity control, is made equal to or less than the combined oil and water mobility. Mobility control continues with a mobility buffer that prevents drive water from contacting the micellar solution. Laboratory and field flooding have proven that the process is technically feasible and that surfactant losses by adsorption on porous media are small. Introduction projects are under way to recover the maximum amount of oil under the most favorable economic conditions.' : New techniques are being developed to increase oil recovery,3" Polymer solutions are becoming an important means of controlling mobility in a waterflood. Thermal methods such as in-situ combustion and steam injection are being used in reservoirs containing highly viscous crudes. Surfactant flooding is receiving attention as a method of reducing interfacial tension to increase recovery.*'" Exotic recovery processes have been considered primarily for ' perations. Economics are unfavorable in most cases for tertiary recovery. studies at the Denver Research Center of the Marathon oil CO. have led to a new oil recovery method.* Micellar solutions (sometimes called microemulsions, swollen micelles, and soluble oils) are used to recover oil by miscible-type waterflooding. Basically, these solutions contain surfactant, hydrocarbon, and water. The method can be used in either secondary or tertiary operations. First, thc concept of thc process is considered in terms of the requirements for an effective miscible waterflood ing operation. Next, micellar solution properties are described including structure, composition, and phase behavior with reservoir fluids. Fluid characteristics are then considered as related to mobility control, and, finally, laboratory and field results are presented to illustrate the efficiency of the process. Concept of the Process Unit displacement efficiency and conformance determine the effectiveness of any oil recovery mechanism. In theory, a miscible waterflood should be capable of a 100-percent unit displacement efficiency with a correspondingly high conformance. Requirements for the slug of a miscible waterflood include (1) 100-percent displacement of oil in the reservoir contacted, (2) controllable mobility, (3) the capability of being driven through the reservoir with water, (4) a low unit cost to enhance economics, and (5) the ability to remain stable in the presence of reservoir water and rock. Micellar solutions satisfy requirements for the slug of a miscible waterflood process. Our discovery that these solutions acted as though they were miscible by displacing all fluids in the reservoir and by being displaced by water solved the miscibility problem. Adequate mobility control is possible by variations in solution viscosity through compositional changes. Economic requirements are met since micellar solution costs below $6/bbl appear possible, Mi cellar solutions stabilize surfactant in the presence of reservoir rock and water, thus reducing the importance of the problem of surfactant loss by adsorption. Fig. 1 illustrates schematically how these solutions are used. Operations start with injection of a micellar solution slug that serves as the oil displacing agent. Next, a mobility buffer of either a water-external emulsion or water solution containing polymer (thickened water) is injected to protect the slug from water invasion. Finally, drive water (water used in a regular waterflood) is injected to propel the slug and mobility buffer through the reservoir. Reservoir oil and water are displaced ahead of the slug, and a stabilized oil and water bank develops as shown in Fig. 1. Stabilized bank saturations are independent of original oil and water saturations. This means that, for a particular type of reservoir, the displacement mechanism is the same under secondary and tertiary recovery conditions. Oil is produced first in a secondary operation. For tertiary conditions, water is produced first. Movement of the slug through the reservoir is stabilized by the mobility buffer. An unfavorable mobility ratio usually exists at the interface between the buffer and drive

Jan 1, 1969

By James J. Bean

HE task of measuring the specific gravity of the -*- operating medium in a heavy-media separation system has never presented a particularly difficult problem because the medium is fairly stable and the overflow of the separatory vessel, as well as its underflow, can be sampled easily and accurately and the specific gravity of the suspension determined easily by weighing a known volume. However, while this method is simple and accurate it does require the operator to take the sample by hand and to weigh it and there is considerable temptation to avoid the periodic sampling if everything seems to be going well, or if something is occupying the attention of the operator. Furthermore all operators do not sample in exactly the same manner and considerable practice is required for two operators to be able to "check" each other to the last few hundredths, particularly if the sample is cut underneath the drainage screen where location of the point of sampling and load on the screens tends to influence the determination. While none of the above presents much of a problem, we have all recognized that some mechanical method of continuous measurement and recording would be advantageous since the operator would merely have to glance at the meter to check the gravity and to have an indication of the trend of any changes. Also if the instrument were of the recording type, a permanent record would be available for the guidance of the superintendent. The Eagle-Picher Mining and Smelting Co. was the first heavy-media user to actually install such a recording meter. In 1946 they installed in their Central Mill at Cardin, Okla., a specific gravity recorder manufactured by the Bristol Co. of Water-bury, Conn. R. A. Barnes, of the Bristol Co., working with E. H. Crabtree, Jr. and Elmer Isern, of Eagle-Picher, made the application and worked out the problems of sampling and measuring. Their attempts to measure the specific gravity of the medium in the cone itself were not entirely successful and they resorted to an outside sample tube for actually making the determination. Because of the particular flowsheet used, it was possible to tap off from the medium return pipeline a stream of medium and divert it into the sampling tube, which was provided with a constant level overflow and a spigot underflow, and into which the bubbler tubes dipped. The Eagle-Picher installation was successful and its possibilities were recognized by the Mineral Dressing Laboratory of the American Cyanamid Co. It was decided to install a similar unit in the heavy-media pilot plant to investigate further its possibilities. Chief among these was the continuous record which it was felt would be proof of the steadiness of the gravity in a heavy-media cone, something which is not always appreciated by POtential users. Because the heavy-media pilot plant is required to operate at a wide range of specific gravities, it was realized that the unit would have to record all gravities from 1.25 to 3.50, and do it to the nearest 0.01. It would not be necessary to record all of this wide range on a single chart and the method selected was to have 4 bands, each band range overlapping the other a small amount and calibrated so that with standard charts one division would represent 0.01 sp gr. A shift from one band to another could be arranged without alteration of the instrument itself, being accomplished by a simple change in the bubble-tube lengths, as described later. Accordingly, a recording type instrument was purchased and installed. Because there were some advantages in doing so, the first installation attempted to measure the gravity of the cone proper by placing the bubble tubes in the cone. This was not at all satisfactory and the second scheme utilized a fixed vertical screen at the surface of the cone, and an external sample-tube arrangement. We were particularly anxious to make this work as we felt it would be advantageous to measure the top level of medium where the separation was actually being made, but we were doomed to disappointment because it was impossible to keep the screen clean of float. Since the top gravity of the cone is the most convenient place to sample for control, a launder about 2 in. wide was installed longitudinally beneath the

Jan 1, 1951

By F. H. Buttner, H. Udin, J. Wulff

Using a modified Udin, Shaler, and Wulff technique, the surface tension of gold Udin, purified helium was found to be 1400 ± 65 dynes per cm for the temperature range 1017° to 1042°C. IN the original Udin, Shaler, and Wulff technique for measuring the surface tension of copper: variously weighted wires were allowed to extend or contract in a copper cell held at elevated temperatures in vacuum. By plotting stress vs. strain for a wire array in one test, the stress at zero strain is obtained. This is the point where the contractile forces resulting from surface tension are balanced by the applied load, according to the expression: y = e=o r [1] where y is the surface tension in dynes per cm; a,,,, the stress at zero strain in dynes per cm; and T, the radius of the wire in cm. The assumption that the wires deform viscously permits the drawing of a straight line through the points on the stress-strain plot. Justification of the assumption has received further experimental support recently.'-' The presence of grain boundaries in the wires requires a correction to the original expression used." Thus: y = d4=T [l- (dl) (ar)Y1 [2] where, n/l is the number of grain boundaries per unit length, and a, the ratio of grain boundary tension to free surface tension. Alexander, Kuczynski, and Dawson in studying the creep of gold wire in vacuum were unable to obtain reproducible values of the surface tension of gold. In plotting stress vs. strain for progressively longer times, they found that the stress at zero strain drifted with time from positive stress values to negative values. Similarly, for the surface tension of silver, reproducible values were obtained only when a purified helium atmosphere was substituted.' Evidently the evaporation rate of silver in vacuum is too high at the temperatures employed to obtain solid-gas equilibrium even in a similar metal enclosure. Thus reproducibility of results is lost. Experimental Procedure The experimental procedure was much the same as that originally developed by Udin, Shaler, and Wulff with a few modifications and improvements. For greater accuracy in strain measurements, knots gave way to cut gage marks as shown in Fig. 1. These were made with a hand-driven lathe in which razor blades serv'ed as cutting tools. Also a more precise cathetometer with a screw accurate to 0.00015 cm was used. The tests were conducted in an atmosphere of purified .helium rather than in vacuum in order to avoid possible evaporation difficulties. Five mil wire of high purity gold (99.98 pct) was used. After cutting in the gage marks, each wire of a series of about 12 was differently loaded by welding a gold ball to one end. This was done by dipping the end of the wire in a cooling gold droplet, previously melted on a charcoal block with a No. 2 acetylene torch. The other end of the wire was strung through a hole in a gold lid and twisted over the edge to hold the wires fixed and in suspension from the lid. The lid and mounted wires were then dipped in pure ethyl alcohol to dissolve any skin oils and dirt on the surface of the wires due to handling. Finally the lid was put in place on an alundum crucible lined with gold so that the wires hung freely within the gold-lined chamber. This whole assembly was next heated in a quartz nichrome wound tube furnace and heated for a few minutes at 600°C to soften the wires. After this anneal the wires were easily straightened with tweezers. The wire assembly was finally annealed 10" to 25°C above the subsequent test temperature for 2 hr. This treatment allowed the grains to grow to equilibrium size and shape. After the anneal, the lid was mounted in front of the cathetometer. The gage length was measured by sighting the 40 power microscope on the upper lip of the lower gage mark for the first reading, then traveling up to the lower lip of the upper gage mark for the final reading. This procedure was repeated four times to give an average gage length value. In this manner the annealed gage length and the final gage length could be measured to determine the strains. During all measurements, grain counts were made.

Jan 1, 1952

By James J. Bean

HE task of measuring the specific gravity of the -*- operating medium in a heavy-media separation system has never presented a particularly difficult problem because the medium is fairly stable and the overflow of the separatory vessel, as well as its underflow, can be sampled easily and accurately and the specific gravity of the suspension determined easily by weighing a known volume. However, while this method is simple and accurate it does require the operator to take the sample by hand and to weigh it and there is considerable temptation to avoid the periodic sampling if everything seems to be going well, or if something is occupying the attention of the operator. Furthermore all operators do not sample in exactly the same manner and considerable practice is required for two operators to be able to "check" each other to the last few hundredths, particularly if the sample is cut underneath the drainage screen where location of the point of sampling and load on the screens tends to influence the determination. While none of the above presents much of a problem, we have all recognized that some mechanical method of continuous measurement and recording would be advantageous since the operator would merely have to glance at the meter to check the gravity and to have an indication of the trend of any changes. Also if the instrument were of the recording type, a permanent record would be available for the guidance of the superintendent. The Eagle-Picher Mining and Smelting Co. was the first heavy-media user to actually install such a recording meter. In 1946 they installed in their Central Mill at Cardin, Okla., a specific gravity recorder manufactured by the Bristol Co. of Water-bury, Conn. R. A. Barnes, of the Bristol Co., working with E. H. Crabtree, Jr. and Elmer Isern, of Eagle-Picher, made the application and worked out the problems of sampling and measuring. Their attempts to measure the specific gravity of the medium in the cone itself were not entirely successful and they resorted to an outside sample tube for actually making the determination. Because of the particular flowsheet used, it was possible to tap off from the medium return pipeline a stream of medium and divert it into the sampling tube, which was provided with a constant level overflow and a spigot underflow, and into which the bubbler tubes dipped. The Eagle-Picher installation was successful and its possibilities were recognized by the Mineral Dressing Laboratory of the American Cyanamid Co. It was decided to install a similar unit in the heavy-media pilot plant to investigate further its possibilities. Chief among these was the continuous record which it was felt would be proof of the steadiness of the gravity in a heavy-media cone, something which is not always appreciated by POtential users. Because the heavy-media pilot plant is required to operate at a wide range of specific gravities, it was realized that the unit would have to record all gravities from 1.25 to 3.50, and do it to the nearest 0.01. It would not be necessary to record all of this wide range on a single chart and the method selected was to have 4 bands, each band range overlapping the other a small amount and calibrated so that with standard charts one division would represent 0.01 sp gr. A shift from one band to another could be arranged without alteration of the instrument itself, being accomplished by a simple change in the bubble-tube lengths, as described later. Accordingly, a recording type instrument was purchased and installed. Because there were some advantages in doing so, the first installation attempted to measure the gravity of the cone proper by placing the bubble tubes in the cone. This was not at all satisfactory and the second scheme utilized a fixed vertical screen at the surface of the cone, and an external sample-tube arrangement. We were particularly anxious to make this work as we felt it would be advantageous to measure the top level of medium where the separation was actually being made, but we were doomed to disappointment because it was impossible to keep the screen clean of float. Since the top gravity of the cone is the most convenient place to sample for control, a launder about 2 in. wide was installed longitudinally beneath the

Jan 1, 1951

By J. P. Zannaras

Referring to the article by R. J. Charles and P. L. de Bruyn, let us assume that W = weight of glass bar; P = weight of hammer; e = total deformation; K = unit of deformation; K = potential stress energy; E = modulus of elasticity; L = length of the bar; 7 = coefficient of inertia; h = height, ft; V = velocity, fps; and a = cross section area. The portion of kinetic energy which is effective in producing stress energy in a fixed bar struck horizontally is given by the formula" 1 + 1/3 W/P K =1+1/3w/p/(1+1/2w/p)2 P.V2/2g = Ph where 1 + 1/3 W/P ? (1 + 3W/P/(1=1/2w/p)2 [8] Putting e e = W/P =------------ From the above equation it can be seen that the maximum transfer of kinetic energy to stress energy is when e = 0 or W/P = 0 which indicates that the weight of the hammer must be very large as compared with the weight of the impacted rod. Eq. 8 diametrically opposes the conclusions reached by the authors of this article. In fact, if their suggestions were followed to the extreme when e = co when P = 0, there would be no transfer of kinetic energy to stress energy at all, as 7, becomes zero. Eq. 8 presumes that the velocity with which the stress is propagated through the bar is infinite, whereas the authors claim that the compression waves reflected are reaching the struck end of the bar prior to the complete transfer of the kinetic energy to cause such modification of the conditions there as to make them reach the reverse conclusions demonstrated by the above formula. That such interference exists is unquestionably demonstrated by the authors and others. However, if my observations are correct, such interference for this specific experiment and also for practical comminu- tion is insignificant, and the conclusions of the authors are in error and must be reversed to comply with Eq. 8. Eq. 5, w. = AE 2/2, given by the authors on page 51, is derived from the following equation (Eq. 9): K = 1/2Pe, where P = Sa, S = ?E, e = EL, and L = 1. The above formula, Eq. 9, cannot be applied in this case. This formula is applicable for static loads where the load increases from zero up to its final value, P, in such a way that the deformation at different instants is proportional to the loads acting at those instants and actually represents the area of a right triangle in the strain load diagram of base e and height P. The typical photographs shown in Figs. 3 and 4 represent the familiar strain load diagrams, and since the line of the wave marks the existence and intensity of the strain with the unquestionable conclusion that such strain has been caused by the action of a load acting continuously all along the wave until it reached the horizontal axis, the work stored at this point is represented by the area under the wave line and the horizontal axis and not by the area of the fictitious triangle given by the authors. Then if this is correct, even visual estimation of these areas at gage stations given in the typical photographs of Figs. 3 and 4 suffice to contradict the authors' calculation given in Figs. 6a and 6b and Figs. 7a and 7b. The typical photographs presented by Charles and de Bruyn show a considerable variation of the intensity of the strain at different stations but very small variation of areas which actually represent the stress energy at the corresponding stations. And, apparently, by squaring the small quantities, the authors magnified their error tenfold. J. M. Frankland's paperV iscusses the relative strain intensity and not the total energy for different types of impact loading. He states in his paper, "The reader is explicitly warned not to confuse the results in this report with those obtained when the load is applied by a blow as from a hammer. In this case the peak load rises to very large but mostly unknown values. The accompanying large deflections and stresses are the result of high values of P, not of the dynamic load factor n. According to Frankland "the dynamic load factor" is the numerical maximum of the response factor. It therefore appears that the authors followed the same procedure in obtaining the relative strain energy ab-

Jan 1, 1957

By Herman Dykstra, R. L. Parsons

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.' THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy's Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. 'Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If' .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951

By W. O. Philbrook, K. Natesan

The oxidation kinetics of spherical pellets of zinc sulfide made from Santander concentrates were studied using a thermogravimetric technique. The experiments covered a temperature range-. of 740" to 102O°C, 0-N mixtures varying from 20 to LOO pct O2, and pellet diameters between 0.4 and 1.6 cm. Mathematical models were formulated to Predict the reaction rate on the assumption that a single transport or interface reaction step was rate -controlling. Analysis of the data indicated that the process of oxidation was predominantly controlled by transport through the zinc oxide reaction-Product layer. ROASTING processes, which are reactions between solids and gases, are very important because they are employed in the production of a number of basic metals. These processes are highly complicated, and one needs to consider the transport phenomena of heat and mass between the solids and gases in addition to the kinetics of various chemical reactions involved. Because of such complications there is a lack of knowledge concerning the rate-limiting factors, which may strongly depend on temperature, particle size, gas composition, and solid structure. The oxidation of zinc sulfide, which is of commercial importance in zinc production, falls into this class of reactions. The major goal of this work was to elucidate the roles played by different process variables, such as reaction temperature, gas composition, pellet size, and pellet porosity, on the kinetics of oxidation of single pellets of zinc sulfide. Roasting of zinc sulfide single particles has been a subject of both experimental and theoretical investigations.&apos;-&apos; The reaction is exothermic and may be considered to be irreversible. Such a reaction has been found to proceed in a topochemical manner. In other words, as the reaction proceeds, a progressively thicker outer shell of zinc oxide is formed, while the inner core of unreacted sulfide decreases. It has been found experimentally, both in the present work and in the previous investigations,1-9 that the particle retains its original dimensions and the process requires transport of gaseous oxygen across the porous product layer for continued reaction. The reaction may be represented by ZnS(s) + 3/2 O2(g) = ZnO(s) + SO2(g) [1] The solid product considered here is only zinc oxide, since the diffraction patterns of zinc sulfide pellets oxidized partially at &apos;798" and 960°C showed K. NATESAN, Junior Member AIME, formerly St. Joseph Lead Fellow, Department of Metallurgy and Materials Science, Carnegie-Mellon University, Pittsburgh, Pa., is now at Argonne National Laboratory, Argonne, Ill. W. 0. PHILBROOK, Member AIME, is Professor of Metallurgy and Materials Science, Carnegie-Mellon University. This paper is based on a them submitted by K. NATESAN in partial fulfillment of the requirements for the Ph.D. degree in Metallurgy and Materials Science at Carnegie-Mellon University. Manuscript submitted December 2, 1968. EMD lines corresponding to original zinc sulfide and the newly formed zinc oxide. OXIDATION MODEL The generalized model for gaseous oxidation of zinc sulfide is illustrated in Fig. 1. This depicts a partially oxidized sphere of zinc sulfide in a gas stream surrounded by a laminar film of gas. The spherical sample of zinc sulfide of unchanging external radius r0 is suspended in a flowing gas stream of total pressure PT and composition specified by the partial pressures of the individual components. Partial pressures of the gaseous species in the bulk gas phase, at the exterior surface of the pellet, at the ZnS/ZnO interface, and at equilibrium for Reaction [I] are identified by the superscripts b, o, i, and eq, respectively. The overall reaction involves the following ~te~s:&apos;~&apos;~&apos; Step 1. Transfer of reactant gas (oxygen) from the bulk gas stream across the gas boundary layer to the exterior surface of the pellet and the reverse transfer of the product gas (sulfur dioxide). Step 2. Diffusion and bulk flow of oxygen from the pellet surface through the product shell (ZnO) onto the ZnS/ZnO interface and the reverse transfer of sulfur dioxide. Step 3. Chemical reaction at the interface, which results in consumption of oxygen gas and generation of sulfur dioxide gas and heat; at the same time the _________ PARTICLE SURFACE / x^^^^n^X / MOVING INTERFACE core-----/ sJSNxy " /T^T^02 \ V\V$\ ^NWX/ "*/— GAS BOUNDARY x. \\ZnO SHELLV/ / \^ . / &apos; ^ BULK GAS --------N_______ _____,------p£ w \ P« / (f> \ / a \ <i) / <----------------- _(o) ^^--------------(b) °- pso, pso2 ro ri 0 ri ro RADIAL POSITION Fig. 1—Generalized model for oxidation of a sphere of zinc sulfide.

Jan 1, 1970

By R. L. Parsons, Herman Dykstra

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.&apos; THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy&apos;s Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. &apos;Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If&apos; .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951

By V. M. Karve, K. K. Majundar, K. V. Viswanathan, J. Y. Somnay

The effect of calcium, magnesium, iron (both ferrous and ferric) and aluminum ions, which are commonly encountered in a typical beryl ore, was studied in the flotation of pure beryl, soda-feldspar and quartz. The vacuumatic flotation technique was employed. With sodium oleate as collector and in the absence of any activator, beryl floated in a pH range of 3 to 7.5, whereas feldspar and quartz did not float at any pH up to 11.5. The pH range of flotation increased in the presence of the ions studied. With calcium and magnesium ions beryl floated from 3 to 11.5 pH and beyond, soda-feldspar floated beyond pH 6 and quartz floated beyond pH 8. Ferrous ion activation was found to be similar to that of calcium and magnesium. Activation by ferric and aluminium ions was found to be complex and the lower and upper critical pH for all the three minerals was around 2 and 10 respectively. These studies indicated the possibility of separation of beryl from feldspar and quartz even in the presence of calcium, magnesium and ferrous ions between pH 4 and 6. Flotation tests on a mixed feed of pure minerals in a 10 g cell revealed that beryl can be selectively floated from feldspar and quartz if ferric ion is reduced to ferrous state or if it is complexed. Beryl occurs mostly in pegmatites, and hence is associated with feldspar, quartz and micas and small amounts of other minerals such as apatite and tourmaline. The separation of beryl from these minerals is difficult because all the silicates accompanying beryl have more or less the same physical properties. Specific gravities of beryl, feldspar and quartz are 2.70, 2.56 and 2.66 respectively. Electrostatic separation has been suggested but no work has been reported. &apos; The adsorption of sodium tri-decylate tagged with Cl4 on beryl, feldspar and quartz reveal similarity in surface properties. Much work has been reported on the flotation of beryl from ores, either directly or indirectly as a by-product, but little is known about the fundamental aspects of beryl flotation. Kennedy and O&apos;Meara3 laid emphasis on prior cleaning of the mineral surfaces with HF. Mica is removed first by flotation of beryl with oleic acid, around neutral pH. Runke4 introduced calcium hypochlorite conditioning in a final separation stage for activating beryl in a mixed beryl-feldspar concentrate, and after washing to remove the hypochlorite, floated beryl with petroleum sulphonate. The Snedden and Gibbs5 procedure is somewhat similar to that of Kennedy and O&apos;Meara. Emulsified oleic acid is used as collector. Recently Fuerstenau and Bhappu6 studied the flotation of beryl, feldspar and quartz with petroleum sulfonate in the presence of activators and stressed the importance of iron in the flotation of beryl. From the studies conducted in this laboratory, it was found that feldspar and quartz as such do not float with sodium oleate, but in practice selective flotation of beryl from feldspar and quartz in an ore is found to be impossible with sodium oleate as collector. A glance at the chemical analysis of typical beryl ore indicates the presence of several ions like Ca ++, Mg++, Al + + + and Fe+++ in abundance and Ti++++ and Mn++ in traces. Hence, in an attempt to explain the behaviour of feldspar in the beryl flotation, the effect of Ca++, Mg++, Al+++ and Fe+++, which are known as gangue mineral activators7&apos;8 has been investigated. Materials and Methods: Lumps of beryl ore (hand picked) were boiled with 10% sodium hydroxide and washed with distilled water. They were further boiled many times with 10% hydrochloric acid till no positive test for iron was obtained with ammonium thio cyanate. This was followed by thorough flushing with double distilled water. The lumps were crushed in a porcelain mortar and pestle under water. The minus 65 + 100 mesh fraction was used for testing and was always stored under distilled water. Pure feldspar and quartz were similarly prepared and the minus 65 + 100 mesh fractions collected. Inorganic ions tried as activators were ca++, Mg++ , Fe++, Fe ++ and A1 +++ . Calcium nitrate, magnesium chloride, ferrous ammonium sulfate, ferric ammonium sulfate and aluminum nitrate of G.R.E. Merck grade were used. B.D.H. technical grade sodium oleate was used as a collector. The vacuumatic flotation technique developed by Schuhmann and Prakash was used for studying the effect of pH on flotability. 7 The indications given by this work were confirmed by using 10 g miniature cell.&apos;

Jan 1, 1965

By D. Atkinson, C. Bodsworth, I. M. Davidson

A geometric model of interstitial solid solutions, which has been used previously as a basis for the prediction of carbon activities in Fe-C austenite, is shown to serve also for the calculation of nitrogen activities in Fe-N austenite. The model has been developed to enable predictions to be made of the activities of an interstitial element in the presence of two host atom species. The activities calculated via the model are shown to be in satisfactory agreement with the measured values in the austenite phase for carbon in Fe-C-Co, Fe-C-Cr, Fe-C-Ni, Fe-C-~n, Fe-C-Si, and Fe-C-V alloys and for nitrogen in Fe-N-Ni alloys. The effect of the second substitu-tional solute on the logarithm of the activity of the interstitial element is expressed as the product of a constant mad the atomic concentration of that solute. The constants so derived we related to the thermo-dynamic interaction coefficients which describe the effect on the activity coefficient of carbon of an added solute element. In recent years the thermodynamic activities of carbon and nitrogen in the single-phase austenite field have been determined for iron binary alloys and for several iron-base ternary alloys. In order to extend the use of these measurements, it is desirable to be able to predict with reasonable accuracy the activities of the interstitials at compositions and temperatures other than those which have been measured experimentally. In all the systems studied to date, the interstitial elements do not conform to ideal behavior. Hence, the available data cannot be extrapolated or interpolated using the simple thermodynamic concepts of solutions. Several models have, therefore, been formulated for the purpose of predicting the activity of an interstitial element in the presence of one species of host atom. These models can be divided into the geometric1"5 and energetic6-&apos; types. The former group is based on the assumption that at low concentrations the activity of the interstitial species is determined by a composition-dependent configurational entropy term and an excess free-energy term which is temperature-dependent but independent of composition. The purpose of this paper is to show that the treatment, based on a geometric model, can be extended to enable predictions to be made of interstitial activities in the presence of two substitutional host atom species. THE CONFIGURATIONAL ENTROPY OF MIXING ICaufman5 has shown that the configurational entropy, S,, for a binary solution comprising of a host atom species, A, and an interstitial species, I, can be expressed as: where NI is the atom fraction of the interstitial species, R is the gas constant, and (2 - 1) is the number of interstitial sites excluded from occupancy by the strain field around each added interstitial atom. The number of interstitial sites per host atom, p, is unityg for the fcc austenite solutions considered here. The configurational entropy of mixing for a ternary solution comprising two substitutional atom species, A and B, and one interstitial species, I, can be derived similarly. Let the number of atoms per mole of each of these species in the solution be represented by «a, ng, and nI. From geometric considerations, it is improbable that the addition of a few atom percent of a second host atom species will change the type of sites (i.e., octahedral) in which the interstitial atom can be accommodated in the austenite lattice. At higher concentrations (determined largely by the relative atomic radii of the atomic species present and any tendency to nonrandom occupancy of the host lattice sites) other types of interstitial sites may become energetically favorable. Restricting consideration to compositions below this limit, for 1 = 1 the number of suitable interstitial sites is given by (n + nB). However, if each interstitial atom excludes from occupancy (Z - 1) additional sites, the total number of sites available for occupation is reduced to (n + ng)/Z. The number of vacant interstitial sites is given by: The total number of recognizable permutations of the atoms must include the recognizable, different configurations of the A and B atoms on the host lattice. Assuming that these arrangements are purely random, and are not affected by the presence of the interstitial species, the total number of recognizable permutations in the ternary alloy is given by: The configurational entropy is obtained by expanding, using Stirling&apos;s approximation, and collecting like items, as:

Jan 1, 1969

By P. N. Richards, M. K. Ormay

The present Paper extends the previous work on cold reduced, low carbon steels to preferred orientations developed after various heat treatments. In recrystal-lized rimmed steel, cube-on-comer orientations increased with cold reductions up to 80 pct. Above that {111}<112> and a partial fiber texture with (1,6,11) in the rolling direction dominated. During grain growth, cube-on-corner orientations have been observed to grow at the expense of {210}<00l>. In re-crystallized Si-Fe (111) (112) and cube-on-edge type orientations are dominant near the surface and the (1,6,11) texture near the midplane for reductions up to 60 pct. With larger reductions {111)}<112> and the (1,6,11) texture are dominant. In cross rolled capped steel a relationship of 30 deg rotation was observed between the (100)[011] of the rolling texture and the main orientations after re crystallization. Most orientations present in recrystallized specimens can be related to components of the rolling texture by one of the following rotations: a) 25 to 35 deg about a (110) b) 55 deg about a (110) C) 30 deg about a (Ill) THE orientation texture of recrystallized steel is of interest where the product is to be deep drawn, because preferred orientation is related to anisotropy of mechanical properties such as the plastic strain ratio (r value);1,2 and in electrical steel applications where a high concentration of [loo] directions in the plane of the sheet improves the magnetic properties of the material. It is interesting to note that both these aims are to a large extent achieved commercially, even though the orientation texture of cold rolled steel does not show large variation3 and the recrystallized orientations are generally given as being related to the as rolled orientations mostly by 25 to 35 deg rotations about common (110) directions.4-6 There is, as yet, no single completely accepted theory on recrystallization. The three mechanisms that have been investigated and discussed are: a) Oriented growth b) Oriented nucleation c) Oriented nucleation, selective growth Largely from the observations of the recrystalliza-tion process by means of the electron microscope,7-11 there is now considerable evidence that the "nucleus" of the recrystallized grain is produced by the coalescence of a few subgrains to form a larger composite subgrain, which finally grows by high angle boundary migration into the deformed matrix. From the intensive work on the recrystallization of rolled single crystals of iron, Fe-A1 and Fe-Si al-loys4-" he following observations have been made: 1) The change in orientation during primary recrys-tallization can usually be described as a rotation of 25 to 36 deg about one of the (110) directions. 2) The (110) axes of rotation often coincide with poles of active (110) slip planes. 3) If several orientations are present in the cold rolled structure, the (110) axis of rotation will preferably be a (110) direction that is common to two or more of the orientations. 4) With larger amounts of cold reduction (70 pct or more) departure from these observations became more frequent. 5) After larger cold reductions, rotations on re-crystallization about (111) and (100) directions have been observed. K. Detert12 infers that a rotation relationship of 55 deg about (110) directions is also possible, by stating that the recrystallized orientation {111}<112> can form from the orientation {100}<011> of cold reduced partial fiber texture A.3 The observation by Michalak and schoone13 that (lll)[l10] formed during recrys-tallization in fully killed steel containing (111)[112],— as well as (001)[ 110] which is related to the {111}<011> by a 55 deg rotation about <110>-implies a possible 30 deg rotation relationship about the common [Ill]. Heyer, McCabe, and Elias14 have recrystallized rimmed steel after various amounts of cold reduction, by a rapid and by a slow heating cycle and found that the preferred orientations strengthened with increased cold reduction. The most pronounced orientation up to about 70 pct cold reduction was found to be {1 11}< 110>, after 80 pct cold reduction both {111}<110> and {111}<112>, after 85 and 90 pct cold reduction, {111}<112>, and after 97.5 pct cold reduction it was {111}<112> and (100)(012). In the present work, the orientation textures of the recrystallized specimens are examined under various conditions of steel composition, amount and method of cold reduction, and method of recrystallization. The relationships between the preferred orientations of the as rolled and recrystallized specimens, and the conditions for the formation of the various orientations during recrystallization are investigated.

Jan 1, 1970

By J. E. Powell, F. H. Spedding

WHILE rare earths are reported to be widely distributed in nature and are not really rare," in practice, there are only a few minerals which are sufficiently rich in rare earths to serve as practical sources. Perhaps the best known of these is monazite which is a phosphate mineral containing rare earths and thorium. This mineral occurs as a dense brown sand in gravel beds and is particularly rich in the light rare earths of the cerium subgroup. This mineral is processed commercially for its thorium, cerium, and lanthanum content, and, consequently, furnishes rich concentrates from which neodymium, praseodymium, samarium, europium, and gadolinium may be obtained. Unfortunately, monazite is rather lean in rare earths heavier than gadolinium. A second mineral which is rich in the light rare earths is bastnasite, a fluoro-carbonate. Extensive deposits of this ore have been discovered in the western United States and have received considerable newspaper publicity in recent years. While bastnasite is very rich with respect to cerium, lanthanum, and neodymium, it contains even less heavy rare earths than does monazite. One of the better sources of heavy rare earths of the yttrium subgroup is gadolinite, a black silicate rock from which the rare-earth content can be extracted readily by acid leaching. It is obtained chiefly from Norway at the present time, although there are known deposits in the United States. Other sources of heavy rare earths include fergu-sonite, euxenite, and samarskite which are refractory tantalo-columbate ores. These minerals require caustic fusion or reduction to carbides with carbon before the rare-earth content can be extracted. All of the minerals which are rich in the heavy rare earths contain yttrium as a major constituent. After the rare earths have been extracted as a group from an ore by chemical means, it is generally convenient to precipitate them from acid media with oxalic acid in order to eliminate certain non-rare-earth impurities such as iron, beryllium, etc., which are usually present. The oxalate can then be readily ignited to R2O3. The oxide can be dissolved in acid and is the starting point for subsequent separation into the pure components. Perhaps the principal reason why the rare earths have not been studied as extensively as other elements of the periodic table, whose natural abundances are comparable, is that they are extremely difficult to separate from each other by the usual chemical means. Prior to 1945, the separation of one trivalent rare earth from another was a laborious process. All separations were based on repeated fractionation processes, i.e., fractional precipitation, fractional decomposition, fractional crystallization, etc. These processes were repeated from a few hundred to many thousands of times in order to obtain individual rare-earth salts of reasonable purity. Of course, mention should be made that, in the few cases where a rare earth could be oxidized or reduced to a valence state other than three, more conventional chemical means could be utilized to separate the oxidized or reduced ion from the other normally trivalent rare earths. The ionic states which deserve special mention are CeIV, SmII, Eu11, and Yb11. When it is possible to remove an element of the series efficiently, due to an optional valence state, its immediate neighbors also become easier to isolate. For example, binary mixtures of lanthanum and cerium, and praseodymium and cerium can be obtained by a relatively small number of fractional operations. The tetravalent state of cerium then allows the complete resolution of the binary mixtures by ordinary chemical means. Although the tetravalent state of cerium has been known for a long time, the divalent states of samarium, europium, and ytterbium were not used extensively in separations prior to 1930 because they are relatively unstable in aqueous media.&apos;-" No attempt will be made to give a comprehensive review of the extensive literature dealing with the separation of rare earths. Rather, this paper will be confined to a review of those methods which have been developed at Iowa State College during recent years, and which have proved extraordinarily successful for the isolation of highly pure rare earths in quantity. It was obvious that, if pure rare earths were to become generally available, methods would have to be developed wherein the thousands of fractional operations made necessary by the similarity of rare-earth properties could be performed automatically. The development of chromatographic techniques and ion-exchange resins appeared to offer a mechanism by which this objective could be accomplished. A number of early attempts were made to separate rare earths by these means; for example, Russell and Pearce12 passed a mixture of rare earths through a cation-exchange column and reported

Jan 1, 1955

By R. J. Reynik

A semiempirial small flunctation theory of diff- sion in liquids is presented, which employs a fluctuation energy assumed quadratic for a small atomic or molecular displacement and Einstein&apos;s random-iralh model. The resulting diffusion equation is given by In these equations. D is the diffusivity, is the average liquid shite coordination number (at interatomic distance d. cm. T is the absolute temperature, xu. em, is (the diffusive displacement. K, is the quadratic fluctuation energy force constant, and rg, cm, are the radii oj diffusing atoms A and B, respectively. The quantities Xn and K are calculated from the computer-filled values of the slope and intercept. respectively. The radius of self-diffusing atom or radii and of diffusing atoms A and B are eta United and compared with values reported in the literature.. The predicted linear variation of diffusivity with. It tempera lure htm been observed in approximately thirty-iire metallic liquid systems, and in over seventy-fiee other liquid systems, including the organic .alcohols, liquified inert gases, and the molten salts, ALTHOUGH the average density within a macroscopic volume element of liquid is constant for fixed total number of atoms. pressure. and temperature, there exist microscopic: density fluctuations within the respective volume element. As such the microscopic volume available to an atom and its Z first nearest neighbors at any instant of time fluctuates above and below the average volume available to these atoms. If one assumes that liquid state atoms vibrate as in a solid. and further postulates that the mean position of any atom in the liquid state is not stationary. but shifts during every .vibration a distance 0 5 j 5 xo. then every atom in the liquid state continuously undergoes diffusive displacements which vary in the range 0 5 j 5 ro. Mathematically. for a binary liquid system consisting of atcrms A and B. the maximum diffusive displacement. .YO, is defined by the equation: where d is the average liquid state interatomic distance at specified liquid state coordination number Z. and v~ \ and vg are the effective radii of diffusing atoms A and B: respectively. For self-diffusion. r^ equals rg , and Eq. [I.] reduces to: It is interesting to note that Eq. [l] or [2] can be used to compute the radii of the diffusing atoms, provided one had an experimental evaluation of xo. As such. the computed radii could be compared with metallic or crystallographic ionic radii to ascerlain the electronic character of the diffusing atoms. Thus it is proposed that in the liquid state the n~otion of an atom relative to its original equilibrium position of oscillation represents the thermal vibration of any atom and its Z first nearest neighbors. while the small and variable displacements. 0 5 1 5 xc,. of the centers of oscillation represent the complex diffusive motions of the atoms at constant temperature and pressure. This is consistent with data obtained from slow neutron scattering by liquids1 &apos; and resembles an itinerant oscillator model of the liquid state.&apos;" It is further postulated that the atomic displacements characterizing the liquid state diffusion process are essentially a random-walk process. As such. it nlay be described by Einstein&apos;s equation:&apos; where D is the diffusivity. sq cm sec-&apos;. j2 is the mean square value of the diffusive displacement. and i> is the frequency of density fluctuations giving rise to diffusion. FORMULATION OF DIFFUSION EQUATION The effective spherical volume occupied by an atom, as a consequence of a microscopic density fluctuation which enlarges the volume available to any atom, exceeds its average liquid state atomic volume by an amount: where AV is the enlarged spherical volume, v is the radius of the diffusing atom. and j is the elementary displacement distance from the original center of oscillation of the vibrating atom to a new center of oscillation position. For small atomic displacements. where c is a constant whose value depends upon the assumed geometry of the enlarged volume. For a spherical increase in volume, c equals 4nr2. Following the treatment of Furthl&apos; and ~walin." assuming the enlarged volu~nes AL7 for the diffusing atoms are distributed in a continuunl. the probability of finding a fluctuation in the size range 0 5 j 5 xo defined by Where c includes the geometric constant cl and Eij) is the fluctuation energy causing the volume change. But the proposed model assumes all the Z first nearest-neighbor atoms are centers of oscillation. and hence the probability that any of these atoms is adjacent to a fluctuation of magnitude 05j5xo is unity. Thus:

Jan 1, 1970

By D. Kunstelj, M. Stubicar, A. Tonejc, A. Bonefacic

A. Bonefafcic, D. Kunsfeli, M. Stubicar, and A. Tonejc The present communication is concerned with the variation of the texture in aluminum wires with die angle and temperature, at constant speed of extrusion. Experiments were carried out concurrently on refined samples, with a stated purity of 99.997 pct pure A1 and commercial samples of 99.5 pct Al. Ingots of refined aluminum samples were machined to 30-mm-diam by 150-mm-long billets. These billets were transformed to 5-mm-diam wires by drawing. The final form, suitable for examination, was obtained by extrusion through conical-face dies with a 1-mm hole diam and extrusion ratio 5:l in diam. The initial form of the commercial aluminum samples was drawn wire 5 mm in diam. These samples showed a poorly defined texture with (111) as a major and (001) as a minor component. A similar defined texture appeared in the refined aluminum samples after drawing to 5 mm diam. Conical-face dies with different angles (defined by the axis and the generating line of the cone) were used in our experiments. The values of the angles were 27, 35, 45, 57, 63, and 90 deg. The extrusion container was fitted with a heating element and controller permitting temperatures up to 600°C to be maintained within i5"C. Extrusion was performed at 250°, 300°, 350°, 400°, and 500°C at constant speed (approximately 1 mm per sec) and constant die reduction. The extrusion product was a wire 1 mm in diam and approximately 20 cm long. In order to remove the surface layer with the "conical" texture and to reduce the absorption by the X-ray examination of the samples, the extruded wire was etched to 0.22 mm in diam. Experiments were performed in the middle sections of the 20-cm-long wires. In addition to the die of l-mm hole diam, dies with a 1.5-, 0.7-, 0.6-, 0.5-, and 0.4-mm hole diam and 63-deg die angle were constructed. In our experiments we did not find in these ranges any important difference concerning the texture of the extruded wires and we continued our work solely on the 1.0-mm die. The diffracted X-rays (Cu K radiation) were recorded photographically. Diffracted intensities were measured on the (111) reflection with a microphotometer. The relative amounts of texture components were determined from the areas under the diffracted maxima. We found the texture of extruded aluminum wires to be strongly influenced not only by the temperature of extrusion and the purity of the sample but also by the form of the die. It is generally admitted that cold-drawn aluminum wires have mainly a (111) texture with a small amount of (001) component, Table I of Ref. 1. In our experiments with wires extruded in conditions represented by Fig. 1, in some cases a single (001) texture was obtained. If these wires were drawn repeatedly at room temperature, X-ray measurements revealed a duplex (001)-(111) fiber texture. Further drawings increased the (111) and decreased the (001) texture component. In Fig. 1 the percentage of material oriented with (001) parallel to the extrusion direction is represented as a function of the temperature and the die angle (a), for commercial and refined aluminum samples, respectively. From these diagrams we may draw the following conclusions. The slope of the die (a) influenced more strongly the texture at the lower rather than the higher temperatures. Again, a stronger influence was found in the case of the commercial in comparison with refined aluminum samples. In the case of the commercial aluminum samples the amount of material with (001) texture increases with increasing wire temperature in an approximately linear manner. This effect is less pronounced in the pure aluminum samples, with the exception of the die with a = 45 deg. In this case the (001) texture decreases with increasing temperature, as shown in Fig. 1. Component (001) is more pronounced in higher-purity aluminum samples. Our experiments led to the conclusion that both (001) and (111) components are essentially stable in extruded aluminum wires. As we obtained a single (001) texture starting with a sample of drawn wire in which the (001) component was very weak, our experiments revealed that the (001) component is not a remnant of the initial texture; this is in disagreement with the findings of Vandermeer and McHargue.1 We gratefully acknowledge discussions with Professor M. Paic. 1 R. A. Vandermeer and C. J. McHargue: Trans. 7MS-AME, 1964, vol. 230, p. 667.

Jan 1, 1969

By R. W. Ludt, C. C. DeWitt

The use of froth flotation for the separation of minerals has become one of the most important of ore dressing processes. Its particular adaptability to the enrichment of low grade ores has made the process an important factor in the national economy. The methods have been extended to the recovery of a great number of minerals. Among the few minerals which have resisted efforts toward industrial flotation is chrysocolla, a hydrated partly colloidal copper silicate. Chrysocolla, being a product of natural oxidation, has been found to occur in small quantities with many ores which are recovered by flotation methods. In present practice, these small quantities of copper silicate pass off with the tailings and are lost. The advantages to be gained by a satisfactory process for the recovery of chrysocolla is apparent. Any application of principles which points a way toward the satisfactory industrial flotation process for copper silicate would be of advantage. This paper presents an attack on this problem. Two methods for the recovery of chrysocolla have been developed by the United States Bureau of Mines.1,2 They have been successful on a laboratory scale but have been seriously restricted in industrial application by critical requirements in the procedure. In one of the Bureau of Mines methods,&apos; the ore is activated with sodium or hydrogen sulphide in an aqueous solution at a pH of 4. Amy1 xanthate is then used as a collector with pine oil as a frother in the flotation process. An excess of sulphide acts as a depressant and the state of optimum conditions is difficult to control industrially. In the second Bureau of Mines method,2 soap is used as the collector at a pH of between 8 and 9. The diffi- culties with this process are that soap is not a specific collector, that heavy metal or alkaline earth ions cause the formation of insoluble soaps, and that a more acid solution causes the formation of a free acid which does not act as a collector for chrysocolla. The problem of recovering chrysocolla by flotation involves the selection of a suitable collector. The collector molecule must be composed of an active polar group that has an attraction for chrysocolla, and of a hydrocarbon chain. Certain dyes have been shown to have an attraction for certain minerals. Suida3 found that hydrated silicates are colored by basic dyes. Dittler4 showed that chrysocolla, among other colloidal minerals of acid reaction, preferentially takes up such basic dyes as fuchsin B, methylene blue, and methyl green. Endell5 gave information to show that the colloidal material in clay may be determined by its selective adsorption of fuchsin. A simple experiment, likewise, illustrates the difference in the adsorptive power of chrysocolla and of silica for the basic triphenyl methane dyes. When a mixture of chrysocolla and silica is immersed in a very dilute dye solution, less than 5 ppm, the chryso-colla is rapidly dyed and the silica is dyed more slowly. The difference is substantial but one of degree. Dean2 showed that the dyes, crystal violet and toluidine blue, are taken up by quartz in an adsorption type process. The difference in the adsorptive power, however, offers the means by which a new collector may act. To form such a collector, a hydrocarbon chain must be attached to the dye molecule. This involves a process of organic synthesis. Butyl, hexyl, and octyl hydrocarbon chains were selected for substitution in the malachite green molecule. For the purpose of identification, the alkyl-substituted dyes formed are called: butyl-malachite green; hexyl- malachite green; and octyl-malachite green. An outline of the procedure for their synthesis is given in the appendix. It is generally recognized in the preparation of this type of dye that the chemical structure of some of the dye molecules varies. However, a uniform formula is attributed to the dye. Such a procedure has been followed in specifying the structure of these alkyl-substi-tuted malachite green dyes. The structure is given on the basis of their properties as an homologous series of dyes, on their method of preparation, and on the purity of intermediates used. Structure of substituted alkyl malachite green is: C6H4 N(CH3)2 p-R C6H4 CH C6H4 N(CH)2 Procedure The flotation cell is a Bureau of Mines 100-g, batch unit provided with an air inlet at the bottom above which is a variable speed agitator. The agi-

Jan 1, 1950

By John L. Miles

John L. Miles (Arthur D. Little, 1nc.)—Pollack and orris" have reported measurements on electron tunneling through A1-A12O3-A1 sandwiches in which the oxide was formed by gaseous oxidation in a glow discharge. From these measurements they deduced the asymmetry of the barrier and, since this is small, conclude that the mechanism suggested by Mott19 for the growth of oxide in thin A12O3 films is inapplicable. In earlier papers20 Pollack and Morris report similar work for oxide films grown thermally. In this case they find a greater asymmetry and conclude that the Mott mechanism is valid. I would like to point out that both these conclusions are quite unjustified. Mott suggests that the growth of the oxide film on aluminum results from the passage of ions through the already present film of oxide under the action of an electric field. This field results from a constant voltage which is in effect a contact potential between metal on one side of the barrier and adsorbed oxygen ions on the other side of the barrier. The theory does not require that the oxide grown is nonuniform either in stoichiometry or structure. It does however specifically assume that the partial layer of ionized oxygen on the surface remains adsorbed on the surface of the growing oxide. In other words, the so-called "built-in field" remains in the oxide only as long as the ionized oxygen is present. When a counter electrode of aluminum is deposited on the oxide, it will react with the adsorbed oxygen on the surface of the oxide, thus forming a small additional amount of oxide. It is clear, then, that there is no requirement in the Mott theory of oxide growth which would necessitate tunneling currents through an Al-A1203-A1 sample to be different when the polarity is reversed. Neither does the theory eliminate the possibility that some additional mechanism could cause the tunneling barrier to be asymmetric and hence tunneling currents to be a function of polarity in such a sandwich. Thus these tunneling-currents measurements are not germane to the question of whether the Mott mechanism is the true method of growth of aluminum oxide films. In fact, it is not surprising that there should be a difference between the oxide properties at the two interfaces (with resulting asymmetry in the tunneling barrier) since the growth conditions and growth rates must have been quite different at these two positions. S. R. Pollack and C. E. Morris (authors&apos; reply)— The point raised by Miles above is one has caused some confusion in the past. The following is an attempt to clarify this point. The built-in field which is responsible for the growth of the thermal oxide at low temperatures arises, according to Mott, because of the passage of electrons from the Fermi surface of the oxidizing metal to surface states introduced by the adsorbed oxygen. It is assumed that the energy of these surface states lies below the Fermi energy of the metal. Electrons therefore continue to flow from the metal to the surface until the built-in electric field raises the potential energy of the surface states to the value of the Fermi energy in the metal, at which time equilibrium is obtained between the surface states and the metal. That is in equilibrium as many excess electrons pass from the metal to the surface per unit time as vice versa. The surface of the oxide prior to deposition of a metallic counterelectrode can then be pictured as follows. The Fermi energy lies in the energy gap of the oxide and is essentially pinned at the energy of the oxygen surface states. The vacuum work function of the oxide is then given by the sum of the electron affinity of the oxide (i.e., the difference in energy between the vacuum and the conduction-band minimum) plus the energy difference between the conduction-band minimum and the Fermi energy. The deposition of a metal onto the surface of the oxide can result in a transfer of electrons across the extremely thin oxide only if there is a contact potential difference between the deposited metal and the parent metal or oxide. That is if the vacuum work function of the deposited metal differs from that of the parent metal, then charge can be redistributed across the oxide in order to equilibriate the Fermi energy across the structure. (It should be

Jan 1, 1965

By J. D. Verhoeven

In the past 15 years a considerable amount of experimental and theoretical work has been done concerning the onset of convection in liquids as a result of interm1 density gradients. This work, which has been doue in many different fields, is reviewed here and extended slightly to give a rrlore quantitative understanding to the probletrz of conzection in liquid metal dlffusion experinletzts. In liquid metal systems the capillary reservoir technique is currently used, almost exclusively, to measure diffusion coefficients. In this technique it is necessary that the liquid be stagnant in order to avoid mixing by means of convection currents. Convective mixing may result from: 1) convection produced as a result of the initial immersion of the capillary; 2) convection produced in the region of the capillary mouth as the result of the stirring frequency used to avoid solute buildup in the reservoir near the capillary mouth; 3) convection produced during solidification as a result of the volume change; and 4) convection produced as a result of local density differences within the liquid in the capillary. The first three types of convection have been discussed elsewhere1-a and are only mentioned for completeness here. This work is concerned only with the fourth type of convection. Local density differences will arise within the liquid as a result of either a temperature gradient or a concentration gradient. It is usually, but not always, recognized by those employing the capillary reservoir technique that the top of the capillary should be kept slightly hotter than the bottom and that the light element should be made to migrate downward in order to avoid convection. In the past 15 years a considerable amount of work, both theoretical and experimental, has been done in a number of different fields which bear on this problem. This work is reviewed here and extended slightly in an effort to give a more quantitative understanding of the convective motion produced in vertical capillaries by local density differences. The Stokes-Navier equations for an incompressible fluid of constant viscosity in a gravitational field may be written as: %L + (v?)v = - ?£ + Wv - g£ [1] where F is the velocity, t the time, P the pressure, p the density, v the kinematic viscosity, g the gravitational acceleration, and k a unit vector in the vertical direction. A successful diffusion experiment requires the liquid to be motionless, and under this condition Eq. [I] becomes: where a is the thermal expansion coefficient [a =-(l/po)(dp/d)], a&apos; is a solute expansion coefficient [a&apos; = -(l/po)(dp/d)], and the solute is taken as that component which makes a&apos; a positive number. Combining with Eq. [3] the following restriction is obtained: Since there is no fixed relation between VT and VC in a binary diffusion experiment, Eq. [5] shows that the condition of fluid motionlessness requires both the temperature gradient and the concentration gradient to be vertically directed. Given this condition of a density gradient in the vertical direction only, it is obvious that, as this vertical density gradient increases from negative to positive values, the motionless liquid will eventually become unstable and convective movement will begin. The classical treatment of this type of instability problem was given by aleih&apos; in 1916 for the case of a thin fluid film of infinite horizontal extent; and a very comprehensive text has recently been written on the subject by handrasekhar.&apos; It is found that convective motion does not begin until a dimensionless number involving the density gradient exceeds a certain critical value. This dimensionless number is generally referred to as the Rayleigh number, R, and it is equal to the product of the Prandtl and Grashof numbers. For the sake of clarity a distinction will be made between two types of free convection produced by internal density gradients. In the first case a density gradient is present in the vertical direction only, and, since the convection begins only after a critical gradient is attained, this case will be called threshold convection. In the second case a horizontal density gradient is present and in this case a finite convection velocity is produced by a finite density gradient so that it will be termed thresholdless convection. Some experimentalists have performed diffusion experiments using capillaries which were placed in a horizontal or inclined position in order to avoid convection. These positions do put the small capillary dimension in the vertical direction and, consequently, they would be less prone to threshold convection than the vertical position. However, if the diffusion process produced a density variation, as it usually does, it would not be theoretically possible to avoid thresh-

Jan 1, 1969

By J. S. Levine, M. Prats

A homogeneous and uniform cylindrical reservoir containing oil and gas is fractured vertically on completion and is produced at a constant bottom-hole pressure. The fracture has an infinite flow capacity, is of limited lateral extent and is bounded above and below by the impermeable strata defining the vertical extent of the reservoir. Results show that such a fractured reservoir can be represented by a reservoir of circular symmetry having very nearly the same production history. The well radius of this circular reservoir is about 1/4 the fracture length and is essentially the same as that obtained previously for a single fluid of constant compressibility. At the same value of cumulative oil production, gas-oil ratios of fractured reservoirs producing at constant terzinal pressure are larger than those of reservoirs having no fractures. This leads to more inefficient use of the reservoir energy in fractured wells and results in lower reservoir pressures for the same cumulative oil production. The reduction in operating life due to fracturing a reservoir is not as great as that for a slightly compressible fluid. This diflerence can be accounted for by the lower reservoir pressure in the fractured reservoir and its adverse effect on the average mobility and compressibility of the oil. As anticipated, the reduction in operating life increases czs the reservoir permeability decreases. The type of results presented in this report can be used to determine the economic attractiveness of fracture treatments per se, to setect the initial spacing to be used in developing a field, and to compare the relative merits of fracturing available wells and infill drilling. INTRODUCTION The effect of vertical fractures on a reservoir producing either an incompressible or a compressible liquid has already been discussed in the 1iterature.l,2 Those results indicate that the production history of such a reservoir is essentially the same as that of a circular reservoir having an effective well radius of approximately one-fourth the fracture length. The present work reports on the effect of a vertical fracture on a reservoir producing two compressible fluids —oil and gas—by solution gas drive. Because of the empirical nature of the PVT and relative permeability data used to obtain the performance of such reservoirs, results can only be obtained numerically and with the aid of high-speed computers. Since reservoirs lose their radial symmetry when fractured vertically, pressure and saturation can no longer be given only in terms of distance from the well. Two coordinates (such as x and y) must now be used to describe the pressure and saturation within the reservoir, and, since we are dealing with compressible fluids, time is also a variable. Thus the solution of a vertically fractured reservoir requires finding two unknowns (pressure and saturation) in two space variables (say x and y) and in time (t). Since no means are readily and generally available for solving such problems at the present time, we have used the results of previous work1,2 to approximate the effect of a vertical fracture on a reservoir producing both oil and gas by depletion. The purpose of the present wmk, then, is to investigate the possibility of using available numerical techniques (limited at the moment to one space variable) to study the two-space-variable flow behavior resulting from a vertical fracture. Results obtained in the course of this investigation are also reported and discussed. Input and output data of the numerical methods used are given in practical units: BOPD, feet, psi, cp, and md. Results are discussed fist in terms of specific reservoir and crude properties and geometries. Later, dimensionless parameters are introduced in order to extend results to different values of some of the reservoir and fracture properties. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The reservoir has an impermeable, circular, cylindrical outer boundary of radius r,. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below&apos; the producing layer (reservoir). It is assumed that there is no pressure drop in the fracture due to fluid flow. 1 indicates the general three-dimensional geometry of the fractured reservoir. Gravity effects and the effects of differential depletion resulting from variations in hydrostatic head (pressure) will be neglected. Thus, the flow behavior in the fractured reservoir is described by the

By O. C. Baptist

The Umiat oil field is in Naval Petroleum Reserve No. 4 between the Brooks Range and Arctic Ocean in far-northern Alaska. The Umiat anticline has been tested by 11 wells, six of which produced oil ; however, [lie productive capacity and recoverable reserves of the field are subject to considerable speculation because of unusual reservoir conditions and because several wells appear to have been .seriously damaged during drilling and completion. Oil is produced at depths of 275 to 1,100 ft; the depth to the bottom of the permanently frozen zone varies from about 800 to 1,100 ft, .so that most of the oil reserves are in the permafrost Reservoir pressures are estimated to range from 50 to 350 psi, increasing with depth, and the small amount of gas dissolved in the oil is the major source of energy for production. Laboratory tests were made on cores under simulated permafrost conditions to estimate oil recoverable by solution-gas expansion from low saturation pressures. The cores were also tested for clay content and susceptibility to productivity impaiment by swelling clays and increased water. content if exposed to fresh water. The results indicate that oil can be produced fronz reservoir rocks in the permafrost and that substantial amounts of oil can be produced from depletion-drive reservoirs by a pre.s.r~lrr drop of as little as 100 psi below the saturation pressure. Freezing of formation water reduces oil productivity much more than that due to increased oil viscosity: Failure of we1ls drilled with rtuter-base mud to produce is attributed to freezing of water in the urea immediately surrounding the wellbore. Swelling clays apparently contributed very little to the plugging of the wells. INTRODUCTION Naval Petroleum Reserve No. 4 lies between the Brooks Range and the Arctic Ocean in northern Alaska. The Umiat oil field is located in the southeastern part of the Reserve and is about 180 miles southeast of Point Barrow (the only permanent settlement in the Reserve and the primary supply point for drilling of the wells at Umiat). Eleven wells were drilled for the U. S. Department of the Navy, Office of Naval Petroleum and Oil Shale Reserves, between 1944 and 1953 to test the oil and gas possibilities of the Umiat anticline. Six of these wells produced oil in varying quantities and the best one pumped about 400 B/D.&apos; Estimates of recoverable oil range from 30 to 100 million bbl. The main oil-producing zones are two marine sandstone beds in the Grandstand formation of Cretaceous age: these are referred to as the upper and lower sands. Good oil shows were found throughout the sand settions in the first three wells drilled on the structure, but the highest rate of oil production obtained on any 01 the many tests was about 24 BOPD. These first wells were drilled with conventional rotary methods using water-base mud; later wells were drilled either with cablc tools using brine or rotary tools using oil or oil-base mud. These experiments were successful as is shown by comparing the oil production from Well No. 2 with that from No. 5. These two wells are only 200 ft apart and are located at about the same elevation on the structure. Well No. 2. drilled with a rotary rig using water-base mud, was abandoned as a dry hole after all formation tests were negative. Well No. 5. drilled with cable tools and reamed with a rotary using oil, pumped 400 BOPD which was the maximum capacity of the pump and less than the capacity of the well. These field results indicated that the producing sands were extremely "water sensitive" and it was assumed that the cause of this sensitivity was the presence of swelling clays in the sands. Because of the very unusual reservoir conditions and the difficulties encountered in completing oil wells in the permafrost. the Navy asked the U. S. Bureau of Mines to make laboratory studies under simulated permafrost conditions to assist them in estimating the production potential of the field and the recoverable reserves. These tests were designed to determine the cause of the plugging of wells in the permafrost and to test oil recovery from frozen sand by solution-gas expansion with the oil gas-saturated at very low pressures. EXPERIMENTAL METHODS AND PROCEDURES Samples Analyzed Core samples were analyzed that represent the lower sand in Umiat Well No. 7, the upper sand in No. 3. and both the upper and lower sands in No. 9. These sands should be productive in all of the wells because of their location on the structure. Core samples from