Search Documents

By S. Kado, G. Derge, L. Yang

Self-diffusion coefficients of silver in molten silver have been measured by means of the capillary-reservoir method in the temperature range 1002" to 1105°C. The results can be .represented by the equation D = (7.10 i 0.06) X 10-4 exp [(—8150 ± 1130)/RT]. The heat of activation agrees within the limit of experimental uncertainty with that for the viscous flow. The validity of the Eyring and the Stokes-Einstein equations are examined on the basis of the self-diffusion coefficients obtained and the known viscosity data of molten silver. MEASUREMENT of the self-diffusion coefficient of silver in molten silver was made in the temperature range 1275° to 1378°K, using the capillary-reservoir method and &"' as the radioactive tracer. The experimental arrangement is shown in Fig. 1. Silver of 99.99 pct purity (fine silver), was melted and degassed in vacuum in a graphite crucible situated inside a McDanel vacuum jacket. The crucible (4 in. tall and 1 in. ID) was protected both at the top and at the bottom with graphite shields and contained about 220 g of silver. The sample holder was a piece of graphite rod (1 in. long and 1/2 in. diam) with 4 capillaries drilled longitudinally into it. The holder was screwed onto a long graphite rod (12 in. long and 74 in. diam) which was connected to a steel disk hanging on a nylon string through a glass pulley. The steel disk kept the sample holder at the center of the McDanel tube and was heavy enough to overcome the buoyancy when the holder was lowered into the molten silver. By turning a handle connected to the pulley through a tapered glass joint, it was possible to raise or lower the sample holder without breaking the vacuum. Two such arrangements were constructed, sitting side by side in 2 Kanthal furnaces. One with non-radioactive silver in its crucible was used for the diffusion run and the other with radioactive silver in its crucible was used for filling the capillaries. To fill the capillaries with radioactive silver, the sample holder was degassed first and then lowered into the radioactive molten silver in vacuum. Purified argon gas was then admitted into the system to push the liquid silver into the capillaries. After cooling in argon, the filled capillaries were polished on emery paper to expose the top of the rod and used immediately as samples for the diffusion run which was carried out at the same temperature as that used for filling the capillaries. During the diffusion run, the samples were degassed in vacuum and lowered slowly into the nonradioactive molten silver before purified argon was admitted into the system. After a predetermined period of time, the samples were raised out of the melt and cooled in argon. The temperature of the melt was measured before and after the run with a chromel-alumel thermocouple inside a McDanel protection tube. Another chromel-alumel thermocouple attached to the outside wall of the McDanel vacuum jacket was used to indicate any fluctuation in the temperature of the furnace during the run, which was found to be less than ± 2°C for the longest run made (93 min). The temperature gradient inside the molten silver was found to be negligible. No significant difference was observed in the results obtained whether the melt was stagnant during the run or was stirred by raising and lowering the sample holder slowly a few times a minute for a distance of about 1/4 in. inside the melt.

Jan 1, 1959

By Julian Szekely

Through the simultaneous statement of Maxwell's equations. the turbulent Navier-Stokes equations and the differential thermal energy balance equation, a mathematical representation has been developed for heat transfer and fluid flow in electric arc furnaces. In broad outline the theoretical predictions appear to be consistent with experimental measurements. A discussion is presented of the general arc furnace control and optimization problem, together with the potential for the incorporation of the mathematical model given in the paper, into future control and optimization schemes. ELECTRIC ARC FURNACE STEELMAKING accounts for about 24% of the total steel produced in the U.S., consuming about 550kw-hr/ton of steel which amounts to about 1.4x1010 kw-hr annually. The percentage of steel produced by the electric arc furnace route has doubled since 1960 and is expected to increase further in the foreseeable future, particularly with the more extensive use of direct reduction technology. Table I shows, in terms of percentages, the principal items that contribute to the cost of running an electric arc furnace. It is seen that in addition to the fixed costs (labor, depreciation, taxes etc.) electrode costs and power consumption are the dominant inputs, and it follows that process changes that would ultimately reduce any of these cost components would have a significant impact on the efficiency of the operation. There are three main areas where progress in arc furnace design and operation must be directed in the next decade (Committee on Technology, 1981). 1. Increased furnace productivity 2. Reduction in processing cost 3. Improvement in the thermal efficiency of the furnace An estimate of the current operating efficiency of the electric arc furnaces may be obtained by comparing the theoretical minimum energy consumption per ton of steel with those achieved in operational practice. Theoretically 330-360 kw-hrs of electricity are required to melt one ton of cold scrap. For superheating, 8-12 kWh per 100°F are required. The actual energy requirements are given in Table II. It is seen that furnaces currently operate at 70-80% efficiency, with the major heat loss from the furnace resulting from the waste gases (Committee on Technology, 1981). While it is apparent that operational experience has produced a technology that is reasonably energy efficient this has been achieved in the almost complete absence of any detailed knowledge concerning the fundamental heat transfer mechanisms occuring in the furnace. The purpose of the present paper is to present a case for the usefulness of mathematical models in increasing the thermal efficiency of furnaces by investigating these heat transfer mechanisms and then combining this knowledge with existing electrode regulator mechanisms to produce truly efficient power cycles for melting and refining. In particular, existing mathematical models describing convection, radiation, and electro-magnetism in the furnace should be employed, and the development of optimization models which balance the cost of electrodes and energy against the thermal efficience of the furnace should be undertaken. Recent work in this laboratory (Szekely, McKelliget, 1981, 1982, 1983; Szekely, McKelliget, Choudhary, 1983) has made a useful start towards solving the problem of heat transfer within the electric arc furnace but the work still needs to be directly related to existing furnace power cycles and tested against plant scale measurements.

Jan 1, 1984

By I. Fatt

The chromalographic effect refers to the separation of constituents in a moving fluid phase which occurs when the phase is passed over a stationary phase, either solid or liquid, or large areal extent. In performing chemical analyses, the large surface contact areal is provided in the case of gas-solid chromatography by-passing the gas over a porous solid and for gas-liquid chronla tog rap hy by-passing the gas over a porous solid which has been coated with a nonvolatile liquid. The porous medium is usually placed in a cylindrical tube, although it may be a relatively flat sheet as in paper chroma-tography. According to one theory,' the separation results from the difference in time spent by various components in the stationary phahe. A gaseous component which is conlpletely insoluble in the stationary phase, or is not adsorbed on the solid, proceeds unimpeded through the column. The versatility and power of the various types of chromatogruphy for the analysis of complex mixtures are evident from their widespread use in performing chemical analyses and from the large number of symposia, books, and papers devoted to the subject in recent years. Associated developments are the application of chromatography to the determination of heats of solution, heats of adsorption, solubilities orK-values, etc. Gas-liquicl chromatography analyses may be characterized by the manner in which the samples are introduced and displaced from the column'. The three techniques used are known as (1) frontal analysis, (3) displacement analysis, and (3) elution analysis. The frontal analysis is performed by displacing an inert gas from the column by a continuous stream of the sample. As indicated in Fig. l(A), all but the first fraction leaving the column represent mixtures rather than pure components. The displacement analysis is conducted by filling the column with a sample, then displacing it with a more strongly adsorbed or absorbed vapor. The resulting chro-matogram is illustrated by Fig. 1(B). Finally, in the elution analysis a small sample of the mixture to be analyzed is introduced as a "plug" into a Rowing stream of carrier gas to give an elution-type chromato-gram. Fig. l(C). The three techniques may be used to provide models to which Row patterns in other physical situations may be compared. PURPOSE OF THE INVESTIGATION The purpose of this investigation was to determine whether natural reservoir rocks, when coated with crude oil, provided sufficient surface areas to give rise to a chromato-graphic separation of light hydrocarbon constituents. The tests were conducted with the light hydrocar-bons in the gaseous state. SANDS TESTED Two sands of widely different permeabilities have been studied: (1) Ventura core (16 md), and (2) Torpedo sandstone core (578 md).Table 1 gives a description of the cores that were tested. With the exception of one frontal-type flow arrangement which was performed on Core 3, the tests were of theelution-type in which a methane-propane mixture was eluted through the core by means of helium. The Ventura core, when impregnated with a weathered crude from the same field, became so impermeable to gas flow that it was necessary to conduct the tests on the granulated cores (Cores 2 and 3).

By D. J. DeLazaro, W. Rostoker, R. E. Riley, M. Hansen

KNOWLEDGE of the isothermal transformation behavior and the TTT chart method of graphically summarizing such information has been of invaluable aid to the ferrous metallurgist in understanding and developing heat treatments intended to provide specific mechanical properties. It is to be expected that similar lines of study directed to titanium-base alloys will be as profitable. This paper summarizes, on conventional TTT charts, the isothermal transformation products and pertinent reaction rate data for binary alloys of titanium containing 1, 3, 5, 7, 9, and 11 pct Mo, respectively. The results so presented were derived from metal-lographic examination. Some X-ray diffraction work was done to clarify certain points of question. Ti-Mo System To the authors' k.nowledge, TTT charts have only been used to describe the reaction rates of eutectoidal decompositions. It would seem that the kinetics of decomposition of any unstable solid solution might be conveniently described in this manner. The equilibrium diagram of the Ti-Mo system' is shown in Fig. 1. The similarity with the Fe-Ni system will be recognized. Additions of molybdenum serve to stabilize the /3 phase to increasingly lower temperatures. The solubility of molybdenum in the a phase is restricted to less than 1 pct at 600°C. Under normal conditions. a homogeneous ,8 structure quenched to a, temperature in the a f & field will reject a of invariant composition. The rate at which this occurs depends on the degree of undercooling and the degree of supersaturation of /3. The structural character of the rejected phase may also be expected to be related to the undercooling. Materials and Experimental Methods A sponge titanium was used which has a nominal impurity content as reported in Table I. The nominal composition of the molybdenum used is also included. The alloys were prepared in 150 g melts by arc melting in an evacuated water-cooled copper crucible of the type reported in an earlier paper.' The pancake-type ingots were forged to 1/4x1/4 in. bars and homogenized for 24 hr at 1000°C in evacuated Vycor bulbs. Specimens Vi in. thick were cut from these bars for heat treatment. The heat treatment cycles followed an invariant procedure. The specimens were soaked at 1000°C for 20 min under a helium protective atmosphere prior to quenching into a lead bath held at a temperature corresponding to a chosen degree of undercooling. After a prescribed period of time at this temperature, the specimen was removed and water quenched. A series of specimens so treated for progressively longer periods served to establish the times for initiation and completion of the decomposition of the /3 phase at a particular temperature. Although heat treatment of unprotected specimens in a lead bath does permit surface contamination by lead diffusion, this method is necessary to insure sufficiently rapid quenching rates. As long as the bath treatments were of less than 1 hr duration. the contaminated surface was not unduly deep and could be removed without excessive grinding. Metal-lographically, the contaminated zone could be readily distinguished from the pure binary alloy. Transformation Characteristics of Ti-Mo Alloys Metallographic examination showed that a supersaturated p phase may decompose by one of two processes. Alloys containing 1 to 9 pct Mo, when quenched in water from the homogeneous /3 field,

Jan 1, 1953

By E. F. Sumner, C. K. Eilerts

A partial differential equation was developed and programmed to compute the cumulation and flow of the liquid phase of a gas-condensate fluid. It was found that such cumulation near the wellbore may attain saturations permitting liquid flow within periods of the order of 1 day. The radius at which this saturation obtains advances toward the external reservoir boundary with continued recovery, although it may move only 2 to 5 percent of that distance during the productive life of the well. This relatively high saturation can cause a serious dimunition in delivery capacity. The pressure dilferential required for flow of gas through a formation with such saturation may be as much as three times greater than the corresponding differential in a formation where saturation is less than 2 percent. Because of the radial geometry, total reservoir pressure differential required for flow with a zone of relatively high liquid saturation about the well was indicated to be two to three times greater than without the saturation. INTRODUCTION A computing program developed by the authors 1 to integrate the partial differential equation for transient radial flow of gas-condensate fluids required a prevalued function, W(Po) = kD(p0)/M(p0)z0(Po).......(1) related to cumulation of condensate in pores of the reservoir as pressure declines. Methods were available for relating mobility to saturation, saturation to liquid-vapor volume ratio, liquid-vapor volume ratio to pressure and, hence, mobility to pressure to obtain W(pD). However, certain assumptions had to be made in advance of program execution concerning movement of the liquid phase and change in radial distribution of liquid saturation as recovery progressed. Because reservoir fluid mobility and related pressure gradients are markedly dependent on liquid saturation, it was desirable to eliminate the necessity for these assumptions and in this respect improve the reliability of a result obtained by integration. Experimental research was begun to determine the relationship of variables that control the reservoir flow of low liquid-vapor ratio fluids. With the porous medium partly saturated, permeability to a two-phase flowing fluid could not be measured because the viscosity was unknown. Viscosity could not be measured without knowing permeability of the liquid-containing medium; but mobility, the ratio of permeability to viscosity, was measurable. This property, and liquid-vapor volume ratio of the flowing phases, have values determined by nature of the fluid and porous medium and, in addition, were found to be pressure, velocity and saturation dependent. A program developed by Kniazeff and Naville 2 for computing multicom-ponent, two-phase flow used polynomials to represent saturation dependence of permeabilities for gas and for liquid phases, and an effluent viscosity constant under reservoir conditions. The authors developed a program based on a partial differential equation of similar mass balance origin but designed to represent distribution of an unlimited number of components between gas and liquid phases utilizing the experimental information available as to variable dependence. Attention was given to liquid-vapor volume ratio of the flowing phases, a property which is essential as a liquid saturation boundary condition. A program for computing implicitly saturation increase at uniform intervals along the drainage radius was combined with an existing program for calculating pressures, and the two were used alternately in an iterative solution for both saturations and pressures. The function W(pD) needed for calculating pressures was replaced by the ratio of mobility to compressibility factor. BASIC EQUATIONS The partial differential equation for calculating

By W. D. Brown

Newell1 has shown that hydrogen is removed from steel in a vacuum at a temperature of 500° to 900° C. within 136 hr. Holm and Thompson2 also state that, especially when the hydrogen is high, the results obtained by extraction at 800°C. agree with those obtained by vacuum fusion. In this laboratory, an apparatus was designed (Oct. 1940) and used wherein the piece of steel was heated at 100°F. in a Torricellian vacuum (Fig. I). This vacuum was produced by raising and lowering a "leveling bottle or bulb" containing mercury. Extensive tests were made with this apparatus. The weakness of the procedure was the possibility of leakage of air into the tube around the rubber stopper at the bottom. The present plan is to raise and lower the mercury by applying vacuum from a mechanical pump, to the top of the tube and then to a bottle of mercury in which the bottom of the tube is immersed. The Apparatus The extraction is made in a Pyrex tube A, I in. i.d. It is drawn down at the top at N and an 8-mm. i.d. tube is welded on it. To this a Pyrex stopcock E, 4-mm. bore, is welded. The tube between N and E is preferably 8 in. long and is graduated in milliliters. The total length of A to the stopcock is 28 inches. Heat is supplied by a split-tube furnace. 0 in. long by 6 in. in diameter, hinged on the side. A short pyrometer couple P is inserted through the walls of the furnace. In use this is standardized by comparison with a standard couple in a I-in. tube, just opposite this short couple. In the I-in. tube A a glass rod R, flattened at one end, supports the sample S and the length of R is so adjusted that the sample is in proper place. (It may float in the mercury.) The tube A is inserted in a large rubber stopper to fit the bottle B and extends nearly to the bottom. It is so placed that the handle of the stopcock faces the operator. In this stopper also is a short tube D bent at right angles for the application of vacuum. Also in this stopper is an 8-mm. tube C reaching nearly to the bottom; this is bent so as to miss the heater and again bent to bring it in line with and as near as possible to the upper and reduced part of tube A. On this tube C is welded a 4-mm.-bore stopcock H. In order that the two tubes may be as near as possible, the handle of stopcock H faces the back, or is the reverse of stopcock E. The total height of II above the bottom is 28 in. As explained later, this is the pressure-measuring tube C. The bottle B is about half full of mercury. The mercury must cover the ends of A and C by about I in. when they are evacuated. By tees T and rubber tubing RT the tubes A, C and D, the latter through the three-way stopcock M, are connected to the vacuum pump. A sheet of asbestos at K prevents the stopcocks E and H from becoming so hot that they leak. A good grade of vacuum lubricant prevents leakage of E and H. In order to measure the pressure of the

Jan 1, 1945

By Walter J. Karplus

By means of a finite difference expansion a fluid flow field in cylindrical coordinates with axial symmetry, is simulated by a network of electrical resistors. A series of DC analog computing units, connected to the network bounrlary corresponding to the free surface effectively adjust the network parameterr to satisfy simultaneously Laplace's equation and the borrndary conditions. The system automatica1ly "re1axes" to the correct boundary shape in a fraction of a second. A detailed application of the method to the water coning problem in oil wells is presented. INTRODUCTION The water coning problem in petroleum production is a prime example of a class of problems known as free surface problems. A free surface arises in fluid flow regions when the shape of the boundary of a single phase flow region is affected in some manner by the pressure or potential gradients within the region. When there exists an interface of two fluids of different densities, the shape of this interface will be determined by the relative velocities of the two fluids and will therefore comprise a free surface. Other free surface problems of interest to petroleum engineers include the gas coning, water injection, and gravity drainage probIcms. In each case it is desired to predict the shape of the boundary of the oil zone as a function of the flow rate, the well penetration and the physical parameters and dimensions of the system. Such problems do not lend themselves to the classical analytical treatments of boundary value problems because in this case the boundary of the potential field and the potential distribution within the field are interdependent and must be calculated simultaneously. A further difficulty arises from geometric considerations. Since the fluid flow is of the radial type in most cases, such effective methods as conformal transformation and the hodograph are not applicable. It is the purpose of this paper to present a method whereby free surface problems of this type may be solved rapidly to the desired order of accuracy by means of a combination of a network of electrical resistors and conventional analog computer units. The water coning problem is treated in some detail. Then the procedure for extending the technique to other free surface problems is indicated. FORMULATION Most commercial oil reservoirs consist of an oil-saturated sand, bounded above by an impermeable strata and below by a water-saturated sand. Gravimetric separation results in a horizontal oil-water interface. When an oil well is drilled so that it partially penetrates the oil zone, and production is initiated at a constant rate, the oil-water interface will assume a new shape as a result of the pressure gradients in the oil zone. It should be emphasized that this is a steady phenomenon. After an initial transient period the water phase will be completely static and only the oil will be in motion. The new steady-state shape of the interface, as shown in Fig. 1, depends upon the depth of penetration, d, the thickness of the oil zone, I, the production rate, q, the hydraulic permeability zone of the sand, k, the viscosity of the oil, µ, and the densities of the oil and water, p.,. and pw.. The problem is to determine this relationship. Darcy's law + = - (P f g p z).....(1) it.

Jan 1, 1957

By Stanley L. Lopata, Eric B. Kula

The problem of preferred orientation has been considered in quantitative phase analysis by X-ray diffraction techniques. The average intensity of a diffraction peak can be obtained by integration over the whole pole figure, or by random rotation of the specimen in the X-ray beam, or by a combination of both. For rod or wire, where the pole figure shows rotational symmetry, the intensity is measured along a radius of the pole figure. Where no symmetry is present, rotational symmetry is imparted by rapid rotation of the sheet around its surface normal during measurement. These methods are applied to the determination of aus-tenite in an Fe-31 pct Ni rod and Fe-28 pct Ni sheet. X-RAY diffraction techniques are frequently used for quantitatively determining the relative amounts of two or more phases present in a structure.1,2 Briefly, the method requires a measurement of the intensity of a Bragg reflection for each phase. The method is based on the principle that the intensity of a Bragg reflection, IA, is directly proportional to the volume fraction of the given phase present, CA: T KRACA/2µ [1] where K is a constant independent of the diffracting phase, RA can be calculated and depends on the Bragg angle ? and the particular reflection measured, and µ is the linear absorption coefficient of the mixture.' For a mixture of two phases A and B, IA = RAcA [2] IA/IB=RAcA/Rbcb and CA+ CB = 1 [3] This allows the volume fraction of a phase to be calculated directly from the measured integrated intensities of the two phases present. This method is frequently used for quantitatively determining retained austenite in steels, but can be used for the quantitative determination of in titanium alloys, as well as phases in other alloy systems. The intensities are generally measured by diffracting only from planes whose normals are perpendicular to the specimen surface, with the implicit assumption being made that this intensity is typical of reflections from planes whose normals lie at other angles to the surface. If the specimen contains a texture or preferred orientation of the crys- tallographic planes by virtue of previous mechanical working, however, this measured intensity is correct only if the intensity level fortuitously is one times random. Examination of pole figures for actual rolled sheet shows that the intensity level at the center of the figure is seldom at the random level, and can very easily differ by a factor of five or more, leading to a commensurate error in the calculated amount of the phases. Textures have been commonly found in materials which are cold-rolled or recrystallized after cold rolling. This problem is becoming of greater importance since structural materials in use today, such as the high-strength stainless steels, frequently are further strengthened by cold rolling, and hence contain textures. Although high-strength carbon-containing alloy steels generally do not contain textures, some of the newer thermomechani-cal treatments will impart a texture even in these cases. In the case of hexagonal metals, such as titanium, probably every wrought sample will show some anisotropy. Since the standard X-ray techniques are inapplicable in all these cases, there is a clear need for some method of taking this anisotropy into account. The specific problem is to obtain the integrated intensity of a Bragg reflection for each phase, averaged over all orientations of the sample with respect to the X-ray beam. This can be done either by purely mechanical means or by a mathematical averaging. Mechanically, a small spherical speci-

Jan 1, 1965

By J. S. Aronofsky, J. P. Heller

This paper presents a mathematical analysis of the fluid mixing which occurs during flow through porous media. The analysis is based on the well-known diffusion equation with mass transfer term. It is pointed out that the use of this equation is justified by a single general assumption which does not specify any particular mechanism for the mixing process. Formulas are given for tracer fluid concentrations for two different boundary and initial conditions. Calculated numerical values compare closely with some published results of experiments in which no viscosity or density difference existed between displaced and displacing fluids. INTRODUCTION A phenomenological theory of the mixing and diffusion process is presented for the flow in a porous medium of two miscible liquids of equal viscosity and density. This is based on the classical diffusion equation which has been used to explain such processes as Brownian motion and heat conduction. More recently, this same differential equation has been derived (in a manner analogous to the work in this paper) for the problem of heat transfer during fluid flow in porous media.' The purpose of this paper is to show that the available experimental evidence justifies the use of this simple diffusion equation. The recent data of Koch and Slobod' and of von Rosenberg' are used to test this diffusion model by comparing concentration curve shapes and calculating "effective diffusion coefficients". The same differential equation for miscible displacement, as well as some experimental confirmation, has been published in Japan by K. Yuhara,4 and recently in this country by Day.5 Yuhara arrives at his con-clusions through analogy with turbulent diffusion, de- scribing the microscopic velocities in terms of "eddy motion" — although unlike the ordinary eddies of hydraulics. Day proceeds directly from Scheidegger's6 concept of the dispersion of velocities. After a very clear description of Scheidegger's theory, he applies it to calculate the motion first of a drop and then of a thin layer of brine into a bed of sand. For the latter case, Day derives and solves the same differential equation as Eq. 1. In our presentation no detailed assumptions are made as to the mechanism of displacement. In fact, several microscopic processes may well be involved—the only general assumption made is that their cumulative microscopic effect may be described by the addition of a second-order term to the differentirl equation of continuity. DERIVATION OF EQUATIONS Let us consider the displacement in one spatial dimension of Liquid A from a sample of porous media by Liquid B. A and B are assumed completely miscible and of the same viscosity and density. Both fluids are also assumed to be incompressible and it is specified that the total flow is the same for any cross section. It is evident that the rate at which Liquid B would be carried across such a section would be given by the product of the fluid velocity, u, and the concentration, C, of B in the fluid mixture assuming the constancy of u and C over the cross section. Since the volume of Fluid B is conserved, this leads to the one-way wave equation - u ?C/?x = ? ?C/?l where x increases in the direction of u—the "Darcy velocity", Q/A (or velocity flux per unit area), and 9 is the fractional porosity. Consider the case where fluid B is displacing Fluid A and the concentration distribution is represented by a step funcLon (with a vertical discontinuity). The solution of Eq. 1 predicts no mixing

Jan 1, 1958

By J. F. Evers, S. Sadiq, F. W. Preston, G. W. Swift

Proper scaling of unsteady-state flow in laboratory gas reservoir models can be achieved if the permeability/fractional porosity ratio is less than 0.5 md. Porous media meeting this specification can be obtained from mixtures of sand and Portland cement provided (1) the water/cement ratio is in the range of 0.35 to 0.40, (2) the sand/cement ratio is in the range of 1/1 to 2/1, (3) the sand grain size is in the range of 48-100 to 18-28 mesh, (4) the mortar curing time is 7 days or more and (5) the cured mortar is dried under vacuum to reduce the temperature level employed during the drying process. The areal homogeneity of the mortars tested is reasonably good, although significant trends in vertical permeability were noted. The Klinkenberg and non-Darcy coefficients of mortar are such that proper scaling of these effects in the laboratory model could be expected. INTRODUCTION A long term deliverability prediction for a gas well is generally based on data determined from flow tests that are of short duration relative to the producing life of the well. These data, together with an estimate of gas reserves in place, form the basis for establishing the parameters of the mathematical model used to predict performance at future times. Results obtained from a mathematical model represent only an interpretation of the physical process being modeled. Experimental verification of the production process from discovery to abandonment is needed where production occurs under carefully controlled conditions. Ideally, one might consider such an experiment to be conducted in the field. Even if the cost of such an experiment could be justified, the time required would make it impractical to consider. However, if the time scaling problem could be handled properly, an experimental model of the gas reservoir production process could provide the necessary data to study the mechanisms of the production process in a reasonable time period. An experimental model has been designed and constructed using the data of this study as a basis. The model reservoir system consists of one producing well located at the center of a right-circular cylinder of porous material. Ten to 15 observation wells are located at strategic offset positions from the producing well. These wells, and the producing well, are fitted with pressure transducers so that the pressure-time history can be recorded at various radial and azimuthal positions in the drainage area during the entire production life of the system. Time is scaled in the model so that the time required to produce gas from discovery pressure to abandonment pressure is about 2 hours. This compares with a 20-year producing life in the model's field counterpart and gives a time scale factor of approximately 1:100,000. DESIGN CRITERIA To establish design criteria for the experimental model, the reservoir is considered to be a finite right-circular cylinder of radius rr, and height h with a producing well of radius rw located at the center and completed through the producing formation. This reservoir consists of a porous medium which has a uniform permeability k and gas porosity and the gas contained in the reservoir has density p, viscosity µ and molecular weight M. Initially, the pressure at all points within the reservoir is uniform at po. and the reservoir temperature is T. Since the reservoir production process has been demonstrated to occur in two distinct transient modes for either Darcy or non-Darcy flow,l it is

By M. E. Wadsworth, J. R. Lewis, J. M. Quets

Samples of snythetic magnetite were reduced in hydrogen at various partial pressures and temperatures. The reaction mas found to be surface controlled and directly proportional to hydrogen partial pressure over the temperature range studied (400° to 1000°C). A transition in the mechanism occurred at approxirnately 590°C. This has been attributed to a change in the slow step of the reduction reaction resulting from the presence of stable wüstite at higher temperatures. Below 590°C the heat of activation for the reduction was found to be 14,700 cal per mole and it was 3,200 cal per mole above 590°C. THE reduction of the oxides of iron has been a subject of study and interest for many years. Edstrom1 has reviewed the work up to 1953 and his findings provide much basic information concerning details of the reduction of hematite (Fe2O3 ) and magnetite (Fe3O4) with both H2 and CO. He concluded that below 600°C the reduction of both hematite and magnetite was essentially identical, since the rate was controlled by a surface reaction at a magnetite gas interface and the iron formed was very porous. The decrease in the reduction rate for both hematite and magnetite immediately above 600°C was attributed to the formation of a stable wüstite phase across which diffusion was very slow. The wide difference in reduction kinetics for hematite and magnetite above 600°C was explained in terms of the character of the wüstite formed. The rate of magnetite reduction increased slowly with temperature beyond the minimum immediately above 600°C and rates equivalent to those observed at 550°C were not obtained until the temperature was increased beyond 850oC. Edstrom and Bitsianes2 studied the solid-state reduction of magnetite, the reaction of which can be written (4x-3)Fe + Fe3O4—4 FexO [1] where x is between 0.83 and 0.95. The rate was found to be controlled by the kinetics of transport of iron ions through the wüstite phase. Bogdandy and Riecke3 studied the kinetics of reduction of small grains of iron oxide in hydrogen and qualitatively attributed the slow step in the reduction process to the diffusion of hydrogen and water vapor through pores. McKewan4 has analysed the data of Joseph,5 Stalhane and Malmburg,6 and Lewis7 and has demonstrated that the kinetics of reduction of both hematite and magnetite are controlled by surface reactions and also that the variation in surface area during reduction must be accounted for. McKewan was able to correlate the reduction data by the equation rodo[1-(1-R)1/3] = Kt [2] where ro is the initial particle radius, do is the density of the iron oxide, R is the fraction reduced at time t, and K is a rate constant. The applicability of Eq. [2] for temperatures above 600°C where a stable wüstite phase is present, indicates the slow step in the over-all kinetics does not result from diffusion unless diffusion occurs laterally along the oxide metal interface or through a wüstite layer of constant thickness. Diffusional control through a layer of wüstite is unlikely, however, in view of the pressure effects demonstrated by Tenenbaum and Joseph,' Diepschlag,9 McKewan,10 and the results of the present study. McKewan has demonstrated a linear rate dependence on the hydrogen partial pressure over the temperature range 400" to 1050°C for the reduction of hematite.

Jan 1, 1961

By M. D. Hassialis

[Surface Areas of Concentrates and Collector Coatings (T.P. 2002, by A. M. GAUDIN and G. S. PRELLER, Min. Tech., May. Discussion by M. D. HASSIALIS and the authors) . I Activation of Minerals and Adsorption of Collectors (T.P. 2082 by J. ROGERS and K. 1'. SUTHERLAND, Min. Tech., Nov. Discussion by A. F. TAGGART and the authors). . . . 3 Correlation between Mineral Behavior in Cataphoresis and in Flotation (T.P. 2005 by A. If. GAUDIN and S. C. SUN, Min. Tech., May. Discussion by M. D. HASSIALIS). . . . . . 8] Surface Areas of Concentrates and Collector Coatings M. D. HASSIALIS.*-In the current flotation theories there exists a common gap spanned by a common assumption pertaining to the arrangement and orientation of the collector ions at the solid surface. Also common is the belief that contact angle and flotation performance are dependent upon the arrangement and orientation of collector ions. This stimulating and informative paper sheds some light upon this general problem. One of the indices of arrangement is the area available per collector ion; in this connection, Gaudin and Preller have applied the B.E.T. method for surface determinations. In a surprising number of cases, their results show a greater specific surface for tailing than for concentrate. This cannot be accounted for by the differences in specific gravity and is attributed to the presence of clay minerals of "primary" sluices in the tailing. This is not a satisfactory explanation, for experience shows that when free slimed clay minerals are present in a flotation cell they will be carried mechanically into the concentrate. A part of the answer to the decreased specific surface of the concentrate would appear to be its usual flocculated condition, yet this should not entirely account for the large spreads reported. The possibility occurs that the nitrogen adsorption isotherms for clean mineral and collector-coated mineral are different. If that were true, and if the coated mineral surface exhibited decreased adsorptive powers, the discrepancy might vanish. In discussing the disposition of collector ions at mineral surfaces, the authors adopt the more consistent viewpoint that there exists some agreement between the ions of the substrate and those of the superstrate. In the case of xanthate ions on galena, a possible arrangement is shown in Fig. 2. This arrangement shows every xanthate ion being shared by two lead ions, whereas the chemical composition calls for two xanthate ions for each lead ion (p. g). This whole discussion is misleading. In the first place, it is known that xanthate ions are not abstracted by unoxidized galena, hence Fig. 2 should show as the substrate not Pbs but PbS„ Om. It would be equally erroneous to use the lattice arrangement of PGSO4, since the surface properties of xanthate-conditioned anglesite differ from those of xanthated galena. It is highly probable that the state of the galena surface exploited in its flotation is ephemeral, since weathered galena is only poorly, if at all, collected by xanthate. Just how a substrate of oxidized galena is to be pictured depends upon one's conception of the transformation from galena to anglesite. This writer is inclined to the following sequence of events. Surface sulphide ions undergo progressive oxidation. As more oxygen is taken up by the sulphide ions they become displaced from their equilibrium positions in the galena lattice, thus accommodating for the larger volume of the S„ Om ions. Concomittant with this change the forces acting on the lead ions

Jan 1, 1946

By H. R. Thresh, A. F. Crawley

THIS paper reports the results of viscosity measurements on cadmium and indium undertaken as part of a continuing program of liquid metals research at the Mines Branch, Department of Energy, Mines and Resources, Ottawa, Canada. For several years this work has continued with support from the Canadian Zinc and Lead Research Committee and the International Lead and Zinc Research Organization. In the literature, existing viscosity data for cadmium are conflicting and for indium are scarce. To achieve reliable viscosity values, further work was deemed desirable. Viscosities were measured using an absolute technique applied to an oscillating closed right cylinder of liquid metal. The apparatus, constructed in this laboratory, and the techniques involved in its use have been described in an earlier paper.' In calculating viscosity coefficients from the experimental data. Roscoe's2 mathematical analysis was used. Because accurate viscosity determinations are dependent on accurate density data, and since the available data in the literature are conflicting, the densities of liquid cadmium and indium were also determined in this laboratory by a pycnometric technique. This increased the control exercised over all the variables contributing to the accuracy of viscosity measurements. For both density and viscosity determinations, cadmium and indium of 99.999 pct purity were used. These materials were obtained from Cominco Ltd. In order to suppress excessive evaporation losses and help prevent oxidation, a technique developed in measuring the viscosity of zinc' was used for cadmium. This involved maintaining an atmosphere of purified hydrogen in the viscosimeter during the rela- tively long periods while the experimental temperature was changed and stabilized between individual determinations of the logarithmic decrement. Then, immediately prior to a decrement measurement. the viscosimeter was rapidly evacuated to a pressure of about 2 X 10-6 torr. On the other hand. with indium. no difficulties arise from vaporization or oxidation. In this case. therefore. the apparatus was held under a constant vacuum of about 2 x 10-6 torr. When all sources of random and systematic error have been evaluated, the technique is considered to yield data accurate within 10.5 pct. The viscosity of cadmium was measured to a maximum temperature of 451.8°C and that of indium to 346.5°C. These temperatures are. respectively. 130° and 190°C above the melting points. Two separate runs were made on indium. The results are given in Table I and, when plotted. yield the smooth curves evident in Figs. 1 and 2. Values of A and C to fit the well-known Andrade1 equation were determined by a regression analysis. The equation is: where n = coefficient of viscosity. i = specific volume. T = temperature in °K Despite its empirical nature. Andrade's equation has been shown by many workers to describe accurately the temperature dependence of viscosity. During the current research program, this equation has been found to fit the data precisely for all the metals and alloys that have been studied.5,6 The A and C' values for cadmium and indium are given in Table 11. Comparing the magnitude of the present data with those of previous workers,7-10 as has been done in Fig. 1, it is evident that the results of Fisher and Phillips10 deviate markedly from the rest. With the exception of one point from Gering and Sauerwald,7 all other data agree within ±2.5 pct. However, the agreement between the present results and the work of Menz and sauerwald9 is especially good over the entire experimental temperature range. Such agreement may be expected since Menz and Sauerwald's

Jan 1, 1970

By W. R. Foster, J. M. McMillen, A. S. Odeh

The complete equations of average linear momentum balance for a single-phase fluid in an incompreisible, homogeneous, porous medium are derived. The derivation begins with Euler's equation of motion for a continuum and uses an integral transform recently developed by Slattery. 7 For steady flow of a compressible, Newtonian fluid, the usual equations of motion result. For transient flow, the space-time description of the pressure is determined in the lowest approximation by the telegrapher's equation. From the analysis a new phenomenological coefficient results which connects the viscous traction to the deriv ative of the linear momentum density. The magnitude of this coefficient determines the velocity of sound through the pore structure in this approximation to the pressure field. INTRODUCTION The modification of Darcy's law of momentum balance for steady, single-phase flow through porous media has been discussed for many years. The first modification was suggested by Forcheimerl who added terms of higher order in the velocity. These can be expected to appear because the underlying microscopic equations of momentum balance are themselves nonlinear in the point velocity field. The Reynolds tensor pvv, which represents the convective flux of momentum density, appears in the momentum balance equation. Only in rectilinear flow (parallel stream lines) does the divergence of this tensor vanish. Since the steady flow stream lines in most porous media are not parallel, nonlinear dependence of the pressure gradient on the velocity should naturally appear. This nonlinearity has nothing to do with turbulence in the ordinary sense of random fluctuations in the pressure and velocity fields. It arises simply because the stream lines converge and diverge, even for steady flows. Klinkenberg2 demonstrated that the permeability coefficient in Darcy's law depends on the absolute pressure or, alternatively, on the density field. However, because he neglected inertial terms of the Forcheimer type, his correction coefficient could not be represented by a constant but tended toward a constant as the velocity decreased. Forcheimer's and Klinkenberg's modifications can be combined in a rigorous way to account for both inertia and slip during steady flow. This will be shown in a future paper. The transient change of pressure in porous media has been described by the diffusion equation.3 This form results from eliminating the velocity and density fields from a combination of the equations of motion in the form of Darcy's law, the continuity equation and an equation of state. Fatt4 suggested that the cause of deviations from the prediction of the diffusion equation for pressure transients lies not in the choice of Darcy's law as the equation of motion but on the existence of dead-end pores which might invalidate the averaged equation of continuity. On the other hand, Oroveanu and Pascal5 noted chat the time derivative of the momentum density must be included in the equations of motion since this measures the local rate of change of momentum density. Their differential equation for pressure is the telegrapher's equation (neglecting gravity). However, their form of this equation predicts that the speed of pressure propagation through the pore structure is the same as that through the bulk fluid. M. K. Hubbert6 attempted a derivation of Darcy's law by volume averaging the Navier-Stokes equations. Since these equations represent momentum balance at a point within an open set of points containing the fluid itself, Hubbert's volume averaging cannot lead to tenns involving transfer of momentum between the fluid and the walls of the pores. Once these viscous tractions are lost by choosing a control volume containing only fluid, they cannot

By R. A. Swalin, L. A. Neimark

The research work presented in this paper had initially a two-fold goal. First, further data concerning the low-temperature recovery of cold work was desired. Recovery phenomena have been extensively studied in face-centered-cubic metals but at the time of initiation of this research, no investigation concerning the low-temperature recovery in body-centered-cubic metals had been made. During the course of the investigation a paper by Martin on the recovery of point defects in molybdenum appeared.&apos; Comparison of our results on tungsten with Martin&apos;s results on molybdenum will be made. Secondly, knowledge concerning the influence of a particular type of impurity addition on the recovery kinetics was desired. It has been known for a long time that extremely small amounts of certain impurity additions exert a profound influence on the recrystallization behavior and the resulting microstructure of tungsten. Small amounts of potassium silicate and aluminum chloride added prior to the reduction of tungsten oxide to the metal have been found by Smithells to raise the recrystallization temperature markedly,2 and yield an ultimate microstructure in which grains several millimeters long are produced in fine wires of say 10 mil diam. If the impurities are added singly, ordinary secondary recrystallization results with none of the spectacular properties of the wire to which special additions were made. Swalin and Geisler3 and Rieck4 have studied the recrystallization behavior of such doped wire in some detail. The former workers found that small amounts of these impurities raised the recrystallization temperature over 400°C and electron microscopic examination revealed no obvious inclusions in the microstructures. In fact. chemical analysis showed that there was <0.001 pct K, <0.002 pct Al, trace of Si, and <0.0003 pct 0 in the wire studied. The possibility that the role of the special impurities is to keep the structure open by void formation and that submicroscopic porosity acts as the inhibiting agent has been discussed.2&apos;3 Meijering and Rieck5 have suggested that the impurities are strung out as second-phase particles in a direction parallel to the drawing direction. It was thought that a kinetic study of recovery could perhaps elucidate the role of the impurities to some extent. When a metal is deformed by cold work three basic types of defects are thought to be introduced into the crystal lattice, namely, interstitials, vacancies, and dislocations. Van Bueren6 has proposed that recovery of cold-worked metal is a four-stage process, each step being the result of the annealing of a primary defect or a combination of defects and having a unique activation energy. The steps in order of increasing activation energy are attributed to a) the diffusion of interstitials recombining with vacancies, b) diffusion of vacancy pairs, c) the diffusion of single vacancies, and d) climb of dislocations. In metals with relatively low- or medium-range melting points the temperatures for recovery are low enough so that the first three steps may occur during room-temperature deformation. For tungsten, because of the high melting point and subsequently high recovery temperature, it was calculated, using Nowick&apos;s approach7 that this recovery would occur in the temperature range 200" to 600°C. It is also known that recovery of X-ray line broadening, presumably due to dislocation climb, occurs at about 600°C.8 EXPERIMENTAL METHOD Two types of cold-rolled 10-mil-diam tungsten wire were used in the investigation; a) pure tungsten, (called undoped) and b) G. E. No. 218 wire (called

Jan 1, 1961

By N. T. Cotman, C. H. Pickering, P. B. Crawford

A beat-conduction model has been developed to study the flow of fluids in a stratified oil reservoir which is being subjected to unsteady-state depletion. To simulate stratification, plates of different metals were joined together so that cross flow could occur between strata. A description of the model, plate construction and necessary instruments is included. The data indicate that stratified oil reservoirs experiencing cross flow may be treated as a single, homogeneous, producing sand having properties intermediate between those of the layers making up the system. A technique for averaging core-analysis data to arrive at the proper average value of rock permeability and porosity is presented INTRODUCTION Although oil has been produced from reservoirs for 100 years and reservoir injection processes have been undertaken for 40 years, the industry is still inadequately informed on the problem of reservoir inhomogeneities. Stratification is perhaps the best-known type of reservoir heterogeneity. Urenl discussed the possible significance of this nonuniformity in his 1927 paper on the theoretical aspects of water flooding. Since this time the problem of vertical permeability variations has been the subject of much discussion and numerous papers. Levorsen2 simply describes the origin of layered beds in his Geology of Petroleum. Some beds are field-wide and may be traced from well to well, while others apparently pinch-out and become discontinuous. In addition, the physical characteristics of a layer may vary from place to place. For these reasons, a stratified system is not easily defined. In calculations involving oil reservoirs in which permeabi1ities vary, methods presently used contain many simplifying assumptions. The reservoir frequently is considered as consisting of permeable layers, each of uniform vertical and lateral permeability and each behaving as if separated from adjacent layers by thin impermeable layers. Other assumptions frequently made in reservoir calculations are steady-state flow conditions, unit mobility ratio, linear geometry, average or homogeneous porosity, and negligible capillary-pressure and gravity effects. Muskat3 has. developed solutions to the stratification problem for mobility ratios other than unity for both linear and exponential permeability distributions. He assumes lateral uniformity and continuity of all productive strata and does not consider porosity. This study will take into account transient effects, cross flow and varying porosity, as well as the variations in permeability, by means of a heat analogy. Muskat5 has shown the similarity of the equations which describe the flow of fluids in a porous medium to the flow of heat through a solid. Landrum, et al,6 made an analogy between the flow of fluid in an unsymmetrical reservoir and the conduction of heat in a similarly shaped metal plate. This analogy is summarized briefly as follows. If fluid is produced from an oil reservoir above the bubble point, the resulting pressure distribution will be identical to the temperature which would result in a metal or heat-conducting plate if heat were similarly removed. This means that a study of heat flow in stratified thermal models will behave like fluid flow in stratified oil reservoirs. For mathematical correspondencies, set- Ref. 6. The fluid mobility k/p is proportional to the thermal conductivity, and the porosity-times-compressibility product is proportional to the density times specific heat of the thermal plate. EXPERIMENTAL EQUIPMENT AND PROCEDURE The stratified reservoir model used in this study is illustrated schematically in Fig. 1. A Leeds and Northrup Speedomax Type G temperature recorder

By Jei Y. Choi, P. G. Shewmon

Using a scratch smoothing technique, the magnitude and the ratio of the surface diffusion coefficients (D,) in two orthogonal directions have been measured on copper for over 100 different surface orientations. Anisotropy was found for all surface orientations except those lying near the (100) and (111); however, these anisotropies were rarely greater than 2.0. The direction of maximum D, appears to coincide with the direction of step edges on any given surface, and if no such steps exist, no anisotropy is found. Thus the jump frequency along a step edge is shown to be faster than the frequency of a jump over the step. The values of Ds and AH for the (100) and the (111) orientations agree quite well with those obtained earlier by the grain boundary grooving technique. It is proven that the directions of maximum and minimum D, (the principal axes) for any surface orientation are known if the perfectly ordered plane corresponding to the given surface orientation possesses a plane of mirror symmetry. Thus the principal axes of D, for orientations lying along the sides of the unit stereo-graphic triangle are always known independent of the extent of thermal disorder or impurity adsorption. The merits and shortcomings of King and Mullins&apos; recently published equations relating Ds to the rate of scratch smoothing are discussed. 1N an earlier paper the authors have reported the variation of the surface diffusion coefficient (Ds) with surface orientation on copper.&apos; In that study Ds was determined from the rate of growth of grain boundary grooves on symmetric bicrystals. This gave values of Ds from 850° to 1070°C in a given direction on each of the six surfaces. It was found that the activation energy for surface diffusion was essentially the same for all six surfaces. The theory developed to explain these data assumed that the free energy of any given type of surface site depends only on the number of nearest neighbors the atom has. This same theory predicted that on any given surface Ds should be independent of direction. The primary aim of this study was to check the validity of this prediction for a wide range of surface orientations on copper. The growth of grain boundary grooves is ideal for studying Ds in a given direction on a given surface over a range of temperatures, because the same specimen can be repolished and used again. However, there is considerable work involved in growing a symmetric bicrystal, and for a survey of Ds in two or more directions on many surface orientations, the technique would be arduous. A better technique appeared to be the scratch method used primarily by Mykura. 2 If a scratch is placed in a surface whose surface tension is independent of orientation, the surface will smooth out in time under the influence of surface tension. To determine D, from this flattening, Blakely and Mykura2 used an approximate solution which was assumed to describe the rate of flattening. Recently King and Mullins3 (hereinafter K and M) have given a more exact analysis of the flattening of an isolated scratch. A second aim of this paper was to check the validity of this newer analysis and compare the values of Ds determined thereby with the values obtained from grain boundary grooving techniques. EXPERIMENTAL PROCEDURE In an effort to work with a particular predetermined surface orientation, light orthogonal scratches

Jan 1, 1963

By Thomas L. Kesler

Present Raw Materials The lithium industry has had exceptional growth since publication of the last edition of this volume, and the present scale of mining and consumption of raw materials is a great many times that of 1949. The materials of present economic importance consist of four pegmatite minerals and lithium-bearing brine. Spodumene, a lithium pyroxene, occurs in greatest abundance, and is the chief source. It has the formula LiAlSi2O6 and theoretically may contain 8.03 pct Li2O, but in the more important deposits the content is 6.5 to 7.5 pct. The hardness is 6.5, the specific gravity 3.1 to 3.2, and the color mostly pale green to white although deep green, lilac, and yellow varieties exist. The crystals are prisms with strongly developed longitudinal cleavages, and crushing characteristically produces lath-shaped particles. Giant crystals up to 47 ft long are known, but in most occurrences the predominant size range is 0.5 to 12 in. long and 0.1 to 2 in. wide. Hypogene* alteration has formed lithia-bearing sericite and coarser muscovite and rarely eucryptite, while supergene* alteration yields kaolinite and hydrous mica. Lepidolite, of second importance, is a mica of varying composition approximating K2Li3-A14Si7O21(OH,F)3. It has lamellate cleavage, book structure, and a specific gravity of 2.8 to 3. The books range in size from microscopic to a diameter of several inches. The color is mostly violet, but other colors include pink, gray, yellow, and white. Most commercial lepidolite contains 3 to 4 pct Li20 and chemically is in the higher grade range of a series of lithium micas that includes polylithionite, zinnwaldite, and protolithionite. In this series, the content of Li2O is exceptionally as high as 7 pct. Petalite and amblygonite are also of economic importance, but are used in smaller amounts. Petalite has the formula LiA1Si4O10 and is white, gray, or pink. The hardness is 6, the specific gravity 2.4, the cleavage in two planes 114° apart, and the Li2O content commonly 3 to 4.7 pct with a maximum of 4.88. Amblygonite has the formula LiA1FPO4 in which OH may substitute for F to form montebrasite. It is white and has polysynthetic twinning in two directions at 90°. Amblygonite superficially resembles potash feldspar, but has a bluish or grayish tint unlike the cream or salmon tint of the feldspar.[t] The hardness is 6 and the specific gravity 3 to 3.1, and the content of Li2O is the highest among lithium minerals that occur in significant quantity, ranging mostly from 7.5 to 9 pct with a maximum of 10.2. The brine of Searles Lake in California30 has been a commercial source of lithium since 1938. See Chap. 7, p. 108. The "lake" is a saline evaporite deposit, and the exposed part has an area of 12 sq mi. To an average depth of 71 ft, this part of the deposit consists of crystalline salts with 50 pct interstitial voids. The voids are filled with brine that contains mainly sodium and potassium compounds, but the brine also contains 0.015 pct Li2O, which is recovered as dilithium sodium phosphate containing 21 pct Li20. This deposit overlies a

Jan 1, 1960

By R. C. Gifkins

R. C. Gifkins (CSIRO)—In this paper evidence is put forward to support the idea of grain boundary shearing in aluminum at 4.2°K and the phenomenon is explained in terms of a low-temperature "equicohesive temperature". It is stated that these findings show that grain boundary shear cannot be ruled out in any situation by temperature arguments alone. We believe, from our own recent study of "apparent grain boundary sliding&apos;&apos; at temperatures below 0.45 Tm, that it can be ruled out on temperature arguments alone.45 Our experiments were made principally with magnesium and a Ms-A1 alloy ("Magnox A12"). Careful examination of offsets of marker lines or steps at grain boundaries produced by deformation at room temperature (-0.3 Tm) showed many of these to result from curvature over narrow zones along the grain boundaries. A significant residuum did, however, appear as a sharp offset or step. It was found that the mean value of this step rose to a small value during the first 0.5 to 1 pct strain and then remained at this level with further strain. Moreover, exactly the same values were obtained (at the same stress) at 75°C instead of 21°C. That is, in this range of temperature the "apparent sliding" was not temperature-sensitive. At 125°C and above there was sliding which increased with strain and which was temperature-dependent. When this temperature-dependent sliding was used to calculate the sliding to be expected at room temperature, this underestimated the "apparent sliding" by a factor of 104. The conclusion drawn is that the "apparent sliding" consists of shear localized in a zone of a width measured perhaps in microns, whereas it seems that true sliding occurs on the boundary surface.46 A model for such localized shear has been developed in general terms, following a similar model of Urie and Wain.47 Other strong evidence for low-temper- ature sliding has been based on the offsets of marker lines consisting of a substructure seen using the electron microscope to resolve these.48 We have repeated this work and found that the effects can be accounted for in terms of cracks in an oxide film which is produced when the substructure is etched up; these cracks form over zones of shear rather than at actual interfaces. However, in aluminum at room temperature (-0.3 Tm) there is no "apparent sliding" which is optically detectable. We suggest that, when the temperature is much lower (and possibly the rate of strain higher, too), as in the experiments of Chin et al., then aluminum behaves like magnesium at room temperature, and a narrow zone of shear develops to accommodate differences in strain across the grain boundary. Such a zone could well give rise to triple-point cracks. It is hoped to investigate these suggestions further. B. Y. Chin, W. F. Hosford, Jr., and W. A. Backofen (authors&apos; reply)—It is interesting to have Dr. Gif-kin&apos;s comments, although it is somewhat difficult to reply in detail without the reference containing the work that he describes. Apparently though, he does find a low-temperature temperature-independent sliding that would seem to be related to the observations reported here. Perhaps some of the difficulty has semantic origins, growing out of the distinction between (&apos;true sliding" and "apparent sliding".

Jan 1, 1965

By S. S. Shen, M. J. Pool, P. J. Spencer, J. R. Guadagno

VERY few thermodynamic data are available for rare-earth alloys, partly because of the difficulties of obtaining the metals in sufficiently pure form and also because of their comparative scarcity up to the present time. The rare earths form a very interesting series in the periodic table in view of the progressive filling of their 4f subvalence electron shell and would be expected to provide useful information with regard to alloying characteristics of transition elements. It has been found in recent calorimetric investigations that the heats of solution of a number of the rare-earth metals in liquid tin are highly exothermic.1-4 Exceptionally strong bonding between rare-earth and tin atoms is clearly indicated in the liquid solutions. Continued thermodynamic studies of the behavior of rare earths in solution could thus contribute a great deal to present theoretical knowledge of alloying behavior and would yield an insight into bonding of a strength hitherto unencountered in liquid alloys. The determination of the heats of solution of the pure components in liquid tin is an essential preliminary to calorimetric work of this nature. Heats of solution of praseodymium, neodymium, and samarium in liquid tin at 750" &apos; and the heats of mixing of Pr-Nd solid solutions&apos; have already been determined in this laboratory. The present investigations with lanthanum continue the study of rare-earth elements. Details of the construction and operation of the liquid metal solution calorimeter used in this work have been published previously by Pool and Peluso and no further description will be given here. The lanthanum used in heat of solution determinations was supplied by Lunex who quoted the purity as 99.9 pct. Their detailed analysis of the material is given in Table I. Samples of 0.0003 to 0.0005 g-atom were used for each run and great care was taken to remove all surface oxide prior to weighing. The samples were stored at all other times in a beaker of oil to prevent their otherwise rapid oxidation. The compounds LaSn (Cudu, L1P type structure) and La&, (structure not classified) were prepared from the appropriate quantities of the two pure components by arc melting under a purified argon atmosphere and once again great precautions were taken to avoid any subsequent oxidation of the specimens. Although no structural examination of the compounds was made, the consistency of thermodynarnic measurements carried out on twenty or so individual samples of each indicates that the compounds were in fact homogeneous, single phase alloys. Samples of the compounds consisting of 0.0003 to 0.002 moles were used in the calorimetric measurements. A solvent bath consisting of approximately 80 g of 99.9 pct pure Sn was used for each series of runs and its heat capacity was determined at regular intervals throughout the experiments using pure tin or tungsten calibration samples. Measurements were made of the heat of solution of lanthanum in liquid tin at 725", 750C, and 775°K and at compositions up to 2.21, 4.76, and 5.25 at. pct La at the three temperatures respectively. In every case, the heat of solution of lanthanum remained constant at its limiting value, within experimental limits, up to the highest lanthanum concentration investigated. Since no change in was observed, it was concluded that the liquid/liquid + LaSn, phase boundary had not been

Jan 1, 1969