Search Documents

By M. Sheffield

This paper presents a technique lor predicting the flow of oil, gas and water through a petroleum reservoir. Gravitational, viscous arid capillary lorces are considered, and all fluids are considered to be slightly compressible. Some theoretical work concerning the fluid flow in one-, two- and three-space dimensions is given along with example performance predictions in one- and two-space dimensions. INTRODUCTION Since the introduction of high speed computing equipment, one of the goals of reservoir engineering research has been to develop more accurate methods of describing fluid movement through underground reservoirs. Various mathematical methods have been developed or used by reservoir engineers to predict reservoir performance, 193-5,7,8 The work reported in this paper extends previously published work on three-phase fluid flow (1) by including a rigorous treatment of capillary forces and (2) by showing certain theoretical mathematical results proving that these equations can be approximated by certain numerical techniques and that a unique solution exists. DISCUSSION The method of predicting three-phase compressible fluid flow. in a reservoir can be summarized briefly by the following steps. 1. The reservoir, or a section of a reservoir, is characterized by a series of mesh points with varying rock and fluid properties simulated at each mesh point. 2. Three partial differential equations are written to describe the movement at any point in the reservoir of each of the three compressible fluids. All forces influencing movement are considered in the equations. 3. At each mesh point, the partial differential equations are replaced by a system of analogous difference equations. 4. A numerical technique is used to solve the resulting system of difference equations. Capillary forces have been included in two-phase flow calculations.1,4,8 The literature, however, does not contain examples of prediction techniques for three-phase flow that include capillary forces. Where capillary forces are considered, each of the three partial differential equations previously discussed has a different dependent variable, namely pressure in one of the three fluid phases. Therefore, three difference equations must be solved at each point in the reservoir. Where large systems of equations must be solved simultaneously, an engineer might question whether a unique solution to this system of equations actually exists and, if so, what numerical techniques may be used to obtain a good approximation to the solution. It is shown in the Appendix that a unique solution to the three-phase flow problem, as formulated, always exists. It is also shown that several methods may be used to obtain a good approximation to the solution. The partial differential equations and differenceequations used are shown in the Appendix. Matrix notation has been used in developing the mathematical results. Two sample problems were solved on a CDC 1604. They illustrate the type of problems that can be solved using a three-phase prediction technique. SAMPLE ONE-DIMENSIONAL RESERVOIR PERFORMANCE PREDICTION A hypothetical reservoir was studied to provide an example of a one-dimensional problem that can be solved. Of the several techniques available, the direct method of solution as shown in the Appendix was used. The reservoir section studied was a truncated, wedge-shaped section, 2,400 ft long, with a 6' dip. (A schematic is shown in Fig. 1.) This section was represented by 49 mesh points, uniformly spaced at 50-ft intervals. The upper end of the wedge was 2 ft wide, and the lower end was 6 ft wide. The section

Jan 1, 1970

By E. E. Mason

The Canadian National Railway system is completing construction of the first major tunnel in North America to use the coarse-aggregate (+ ½ -in.) shotcrete technique of primary support and lining. The tunnel was driven 29 X 20 ft in cross section, 10,760 ft long, passing under an industrial and residential area of Vancouver. With the exception of 1400 ft in soil (mainly impervious tills), the tunnel passes through a very young series of flat-lying sedimentary beds, to a large extent so incompletely consolidated as to fall apart on exposure. For a conventional steel-support system, it was estimated 8 WF, 28-lb steel should be used at 5-ft centers. Cost of erection including lagging would amount to $120 per linear ft of tunnel at 304 per lb. With 8-in. concrete cover on the steel, lining would cost $280 per linear ft at $82 per cu yd. Thus, minimum cost of support and lining would have amounted to $400 per linear ft. It is now estimated that at least 90% of the rock tunnel would have required support and 100 % required lining. Specified shotcrete support was 6 in. in the arch and 4 in. in the walls. Cost of this work was roughly $100 per linear ft. In addition, another $10-$12 per ft has been incurred in drainage and waterproofing wet areas and in measures to correct experimental and operational errors, thus finalizing the support as a lining (Fig. 1). The tunnel was driven three shifts five days weekly, no blasting being allowed by city ordinance between midnight and 7.00 a.m. A six-drill jumbo was used to drill a 10-ft round, 110 holes blasted three times daily. A 2-in. layer of shotcrete was applied to the newly blasted arch commencing within 45 min of blasting. This was done from a flying deck that extended over the muck pile from the jumbo (Fig. 2). Also during the mucking cycle, the preceding incomplete arch support was brought up to 6-in. specification. The walls were sprayed 4 in. thick during the drilling cycle (Fig. 3). The shotcrete machines and aggregate bins were mounted on the jumbo.

Jan 1, 1970

By N. G. W. Cook

At present, tunnels can be driven in rock by either of two essentially different techniques, in one of which rock-breaking is accomplished by drilling and blasting and in the other by mechanical loading. The methods based on these techniques are comparable in terms of the rates at which tunnels may be advanced, although they differ widely in their range of applicability to different types of rock and lengths of tunnel. The method based on drilling and blasting has reached a high degree of development. It is applicable to even the hardest rocks and rates of advance in excess of 50 ft per day (fpd) are attainable. The main disadvantages of this method are that it involves a sequence of operations punctuated by blasting and that it is relatively uncontrolled, resulting in overbreaking the tunnel profile and damage to the surrounding rock. However, its flexibility, mobility, and low capital cost constitute real advantages in many situations, such as those involving short lengths of tunnel or low rates of advance. Tunneling machines based on mechanical breaking of the rock are a more recent development. They offer a continuous and controlled means of tunneling capable of rates of advance in excess of 100 fpd under favorable conditions. However, these machines weigh many tens of tons and consume several hundred horsepower, so that they are costly and lack mobility. At present, their applicability is limited to long tunnels in relatively soft rock where the high rates of advance and tunnel quality can offset their onerous capital cost. Many situations exist in which the superior quality of a machined tunnel and freedom from blasting would be real advantages, but where the rock is either too hard for existing machines or the length of the tunnel does not justify their capital cost: What are the possibilities of developing cheaper, more flexible, and mobile tunneling machines? Most existing tunneling machines fragment all the rock at the tunnel face to facilitate its disposal. If a machine were used to cut only the profile of a tunnel as one of a sequence of operations, it would be required to break

Jan 1, 1970

By John W. Bader

Probing a mineralized area with a core drill to find out what's there is fundamentally a simple task, but a costly one. Quite commonly it costs $10, and sometimes more, to extract and assay each foot of drill core obtained. Yet, it is not unusual today to drill hundreds of holes on one mineral property. The initial objective of such probing is almost always simply to determine the potential worth of the ore- body, if any exists, and the question naturally arises: How little sampling is enough for this determination? The mathematics of statistics infer accuracy of estimates largely from the sheer number of samples taken. Mineral exploration usually involves accumulation of large numbers of samples, and thus it is usually not difficult to place narrow confidence limits on either side of the estimated average grade of a sampled orebody. The determination may, or may not, be as accurate as the derivations imply. If it is not, the fault is probably one of misapplication and not with the statistical processes themselves. Un- fortunately, indications are strong that the ordinary statistics of sampling usually cannot be applied to orebody exploration samples with actual attainment of the accuracy of estimate that they express mathematically. For this reason, a more deductive, less mathematical analysis of the sampling may well be worthwhile.

Jan 1, 1970

By William G. Fricke

MaCROSEGREGATION, as distinct from microseg-regation, is often inverse in nature. That is, the eu-tectic solute elements, which are expected to be con-

Jan 1, 1970

By R. J. Fruehan, F. D. Richardson

Electrochemical measurements have been made of the activity of oxygen in copper and its alloys with silver and tin at 1100" and 1200°C. The galvanic cell used was Pt, Ni + NiO/solid ellectrolyte/[O] in metal, cermet, Pt The results do not support any of the equations so far designed for predicting the activities of dilute solutes in ternary solutions from activities in the corresponding binaries. If, however, a quasichemical equation is used with the coordination number set to unity, agreement between observed and calculated activities shows that this empirical relationship can be useful over a fair range of conditions. SEVERAL solution models have been proposed for predicting the activity coefficients of dilute solutes in ternary alloys from a knowledge of the three binary systems involved. Alcock and Richardson1 have shown that a regular model, and a quasichemical model,' in which the dissolved oxygen is coordinated with eight or so metal atoms, can reasonably predict the behavior of both metal and nonmetal solutes when the heats of solution of the solute in the separate solvent metals are similar. But when this is not so, neither model gives useful predictions unless coordination numbers of one or two are assumed. Wada and Saito3 subsequently adopted a similar model to derive the interaction energies for two dilute solutes in a solvent metal. Belton and Tankins4 Rave proposed both regular and quasichemical type models in which the oxygen is bound into molecular species, such as NiO and CuO in mixtures of Cu + Ni + 0. However, their models have only been tested on systems in which the excess free energies of solution of the solute in the two separate metals differ by a few kilocalories. Ope of the reasons why more advanced models have not been proposed is because few complete sets of data exist for ternary systems in which the solute behaves very differently in the two separate metals. For this reason measurements have been made of the activities of oxygen dissolved in Cu + Ag and Cu + Sn. Measurements on both systems were made by means of the electrochemical cell, Pt, Ni + NiO/solid electrolyte/O(in alloy), cermet,Pt [1] The activity of oxygen was calculated from the electromotive force and the standard free energy of formation of NiO, which is accurately known.5 Before investigating the alloys, studies were made of oxygen in copper to test the reliability of the cell and to check the thermodynamics of the system. Of the previous studies those by Sano and Sakao,6 Tom-linson and Young,7 and Tankins et al.8,7 have been made with gas-metal equilibrium techniques; those by Diaz and Richardson,9 Osterwald,10 wilder," Plusch-kell and Engell,12 Rickert and wagner,13 and Fischer and Ackermann14 have been made by electrochemical methods. EXPERIMENTAL The apparatus employed was the same as described previously,9 apart from slight modification. The molten sample of approximately 40 g was held in a high grade alumina crucible 1.2 in. OD and 1.6 in. long. The solid electrolytes were ZrO2 + 7½ wt pct CaO and ZrO2 + 15 wt pct CaO; the tubes 4 in. OD and 6 in. long were supplied by the Zirconia Corp. of America. They were closed (flat) at one end. In one experiment with Cu + O, both electrolytes were used in the cell at the same time. The reference electrodes inside the electrolyte tubes consisted of a mixture of Ni + NiO. They were made by mixing the powdered materials and pressing them manually into the ends of the tubes, with a platinum lead embedded in them. The tubes were then sintered overnight in the electromotive force apparatus at 1100°C. By sintering the powders inside the tubes (instead of using a presintered pellet9) better contacts were obtained between the electrolyte, the powder, and the platinum lead. Troubles arising from polarization9 were thus much reduced. The electromotive force was measured by a Vibron Electrometer with an input impedence of 1017 ohm; the temperature was measured with a Pt:13 pct Rh + Pt thermocouple protected by an alumina sheath. The couple was calibrated against the melting point of copper. The cermet conducting lead of Cr + 28 pct Al2O3, previously found to be satisfactory9 for use with Cu + 0 was also found satisfactory with Cu + Ag + 0 and Cu + Sn + 0. Superficial oxidation was observed, but it did not interfere with the working of the cell. The reaction tube containing the cell was closed at each end with cooled brass heads and suspended in a platinum resistance furnace. The tube was electrically shielded by a Kanthal A-1 ribbon which was wound round it, and the ribbon was protected by a N2 atmosphere between the furnace and the reaction tube. The cell was protected by a stream of high purity argon which was dried and passed through copper gauze at 450°C and titanium chips at 900°C. All the metals used were of spectrographic standard. Procedure. In studies of the system Cu + 0, be-

Jan 1, 1970

By J. D. Hunt, R. J. McDonald

It is proposed that extensive fluid motion occurs between the dendrites of a Partially solid casting and it is suggested the A segregates in steel castings are a direct result of this motion. The fluid flow arises because of density differences between liquid ahead of the dendrites and the interdendritic liquid. Experiments are described on a model system investigating this fluid motion. The importance of interdendritic fluid flow is discussed with reference to longitudinal segregation. SERIOUS macro segregation is often found in large killed steel ingots. Among the most common segregates are those known as A and V, Fig. 1. The A segregates are arranged over the surface of cones, the wide parts being at the base of the ingot. The A segregates are pencil like in shape and extend from the columnar into the equiaxed region. The V segregates are distributed over inverted cones but have a less well defined shape. They occur at the center of the ingot. Because of the commercial importance of macro inhomogenities these segregates have attracted wide-scale attention. The segregates contain alloying elements and impurities which depress the freezing temperature and are thus the last liquid to freeze. A number of explanations have been proposed to explain their presence. It has been suggested1'' that the segregates are the result of the equiaxed zone growing towards the columnar region together with the interaction of convection within the bulk liquid. It has also been proposed that the segregates are the result of nucleation and rapid growth ahead of the dendritic front.3 Neither of these explanations is completely satisfactory because A segregates are pencil shaped and are not found as continuous cones. In this paper an alternative mechanism is proposed, and experiments are described on a model system to determine whether the mechanism occurs in practice. It is thought that A and V segregates are not peculiar to steels but should occur in many other large castings. Fluid motion within the bulk liquid region of a casting has been extensively discussed. Little attention has been paid to fluid motion within the partially solid regions. The work415 that has been done has been mainly concerned with the fluid motion which results from the volume contraction on freezing. The volume contraction is used to explain various types of macro segregation. The present authors suggest that in many cases the fluid flow resulting from volume contraction will be small compared with that due to other causes. When the columnar region is growing, the liquid still remaining between the dendrites has a much lower temperature near the ingot walls than at the dendrite tips. The composition of the liquid in contact with solid will vary with temperature. The density of the liquid will thus vary both because of changes in temperature and composition. The layers of different density liquid will give rise to a slow but steady liquid flow through the interdendritic regions. If there is no density change with composition, the temperature gradient within the dendritic region should give rise to a downward flow between the dendrites, shown schematically in Fig. 2. When the liquid density changes appreciably with alloy concentration this motion could either be en-

Jan 1, 1970

By E. J. Rapperport, J. P. Pemsler

The thermodynamic activities of zinc in six solid solution Cu-Zn alloys ranging from 5 to 35 at. pct Zn were determined experimentally in the temperature range 400° to 600°C. This low temperature investigation was canducted in order to evaluate techniques developed to utilize the inherently high sensitivity of atomic absorption flocesses in the measurement of thermodynamic activities. Analytical expressions ,for the activity and actizlity coeflcient are given for the six alloys in the temperature ranges investigated. RELATIVELY few experimental methods are available for investigation of thermodynamic activities of alloys, especially in the solid state. The techniques most frequently used have been the electrochemical potential and the effusion methods, both of which have severe limitations in many instances. It is therefore desirable to expand the ability to perform such measurements in order to obtain new information as well as to provide an additional independent verification capability. In this work, we present a significant improvement in the spectrophotometric method for sensing small vapor pressures in static absorption cells. Similar techniques have been used previously;1"5 however, applications had been limited to relatively high pressures, often greater than 1 torr. Prior investigators have, for the most part, used broad spectral sources such as xenon or mercury lamps, and high intensity arcs. Hollow cathode sources were first suggested in 1956 6 and were used soon afterwards.4'5 These sources offer significant improvements in sensitivity and freedom from interfering spectral lines.'-' EXPERIMENTAL High purity zinc was obtained from Cominco Products, Inc., and copper from American Smelting and Refining Co. Both elements were of 99.999 pct purity. Copper-zinc alloys were vacuum melted in a high fired carbon crucible with each alloy pulled from the melt as a 4 -in. diam bar. The bars were swaged to -1/4 in. rods and vacuum annealed for 160 hr at 800° + 1°C. Samples for gross chemical analysis were taken at intervals along the length of the rods to ascertain the axial zinc gradient. Electron microprobe analysis of homogenized specimens indicated that the alloys had uniform compositions over their cross sections on a macro (200 p) and micro (1 u) scale to better than *1 pct (20) of the gross composition. This tolerance was determined by counting statistics, rather than assured composition fluctuations. All SiO 2 windows were high-ultraviolet-transmission grade to minimize intensity losses. Silica absorption cells were scrupulously cleaned consecutively in organic solvents, dilute HF, and distilled water before use. The empty cells were then flamed while under a dynamic vacuum, cooled, and removed to an argon-filled glove bag. Alloy pieces were cut and filed in the glove bag to produce fresh surfaces, and then loaded into the cells. The loaded cells were temporarily sealed, removed from the glove bag, reevacuated to 10-5 torr or better, and permanently sealed. The instrument used is schematically shown in Fig. 1. The spectral emission from a commercially made hollow cathode lamp (A) of a selected element is focused through an absorption cell (B) inside a well-controlled furnace (C). The intensity of the transmitted beam is measured using the spectrometer* (D) 'Techtron model AA4 atomic absorption spectrometer. which contains a grating (E) that disperses the light prior to impingement on the photomultiplier (F). The monochromator grating is adjusted so that only the wavelength of interest is measured. The power supply delivered an interrupted voltage to the lamp, causing a chopped radiation output to be transmitted. The detector read only the intermittent component of radiation incident upon it, so that all continuous noise signals (furnace radiation, and so forth) were eliminated. Three recording thermocouples contained in the muffle furnace were positioned along the length of the absorption cell: one at each end and one at the center. An effort was made to keep the ends of the cell several degrees hotter than the center to avoid window condensate. Appropriate thermal corrections were then necessary to relate cell pressure to radiation attenuation. Water-cooled heat shields, as shown in Fig. 1, were found to aid signal stability by protecting the hollow cathode and the photomultiplier from furnace radiation. The furnace had a 2-in. diam muffle, Kan-thal wound, with SiO 2 windows at its ends to minimize convective effects. The hollow cathode radiation was masked and focused to form a conic beam that was a maximum of { in. diam within the furnace. Thus, the 1.5 in. diam absorption cell easily contained the entire beam. The furnace was mounted on ball-bearing slides with positive positioning detents. This arrangement allowed the removal of the entire furnace assembly from the radiation path, position [I], Fig. 1, so that frequent sampling of the unattenuated beam intensity could be obtained. In all cases the beam intensity was kept constant to 0.1 pct as judged by readings taken immediately before and immediately after data collection. Only data for absorptions of less than 80 pct were utilized, as systematic deviations from linearity were found for greater absorptions.

Jan 1, 1970

By W. C. Winegard, C. J. Simpson

BNINGESSNER et al.1 have shown that the solidification substructure in dilute cadmium alloys persisted after repeated deformation, recrystallization, and annealing treatments. On the basis of this and similar evidence, Whelan has suggested in a recent publication' that the solidification structure might affect the mobility of grain boundaries even after primary re-crystallization. This suggestion is important, since the technique used by previous investigators3-5 to fabricate grain growth specimens consisted only of casting, deforming, and recrystallizing prior to the grain growth anneal. Usually no attempt was made to homogenize samples with prolonged annealing. Experiments. To test the Whelan hypothesis, grain growth experiments were performed with samples in which the degree of solute segregation was altered by controlling the solidification rate of the alloys. The initial experiments were performed on the Pb-Ag system. Alloys containing 0.01 at. pct Ag were produced from 99.9999 pct Pb and 99.999 pct Ag. Cylindrical samples, 40 cm long and 12 mm in diam, were solidified in graphite-coated Vycor tubes by lowering them through a vertical furnace at varying rates with the furnace power settings constant. Six-centimeter sections were discarded from each end, and the 28-cm section from the center was cut into sections 25 mm long for grain growth experiments. These samples were deformed with the aid of a drophammer to a uniform thickness of 5.7 mm. After almost instantaneous recrystallization, the samples were given a grain growth anneal for various times at 260°C in a salt bath. Experimental results are shown

Jan 1, 1970

By R. A. Greenkorn, R. C. Smith

Hele-Shaw cells are used to model creeping flow through porous media (where Darcy&apos;s law is valid). The effects of inertia on flow about obstructions in a Hele-Shaw cell can be calculated by a perturbation method if one can determine a solution to Laplace&apos;s equation. Results of a computer solution for flow about circular, square and elliptical obstructions are presented. These results show that or a modified Reynolds number of less than I, the inertia terms are small; and for values of less than 3, the average streamline predicts the ideal flow. Therefore, the analogy might be used for studying flow in porous media up to a modified Reynolds number of at least 3. INTRODUCTION The nature of fluid flow in porous media is of interest in the fields of soil mechanics, ground water flow, petroleum production, filtration and flow in packed beds. Because it is very difficult to study the phenomenological behavior of flow in porous media, homologs and analogs are used to study flow characteristics. A Hele-Shaw model, made of two closely spaced plates - usually glass — is often used as an analogy to two-dimensional flow in porous media. Hele-Shaw2 showed experimentally that the streamline configuration for creeping flow around an obstacle located between two closely spaced parallel plates is the same as for two-dimensional, ideal flow about the same obstacle. Stokes7 verified these observations mathematically. The usual equation of motion for flow in porous media is Darcy&apos;s law. The form of the mathematical statement of Darcy&apos;s law is identical, within a multiplicative constant, to the expression for the average velocity over the plate gap in the plane of a Hele-Shaw model. These models may be used to describe flow in both homogeneous and heterogeneous porous media. In the mathematical proof of the Hele-Shaw analogy it is assumed that the convective terms in the Navier-Stokes equations are negligible and that the equations of motion degenerate to Laplace&apos;s equation, with pressure the dependent variable. Whenever a Hele-Shaw model is used as an analogy to flow in porous media, .the validity of this assumption is in question. Riegels5 showed that if convection is not neglected, the velocity distribution around a cylindrical obstruction in the flow field depends on a Reynolds number, the plate spacing, and a dimension characteristic of the obstacle. Riegels&apos; solution, a perturbation solution, uses the boundary condition that the flow rate into the obstacle averaged over the plate gap at any point on the obstacle is zero. The method requires that a solution to Poisson&apos;s equation for the perturbation pressure be found. Riegels evaluated this solution for the case of the cylindrical obstruction. His method may be simplified by eliminating the need for solving Poisson&apos;s equation for the perturbation pressure. Instead, an analytic expression for the perturbation pressure gradient is obtained (valid for arbitrary shapes) and used to eliminate pressure from the equations for the perturbation velocities. The results show, for symmetrical shapes, that if NRe < 1, the convective acceleration terms are small, and that the average velocities represent ideal flow up to at least NRe = 3, where: L is a characteristic dimension of the obstacle perpendicular to flow, b is the plate spacing, u is viscosity, v, is velocity of approach and p is density. THEORY* Incompressible fluid flow through a Hele-Shaw model and a porous medium are analogous because the equations for the average velocity components in two-dimensional flow are identical if:

Jan 1, 1970

By D. M. Urich

New plant construction for pyrolysis and agglomeration during 1969 has been relatively quiet in the North American scene. U.S. pellet capacity approximates 52 million long tons and that of Canada 25.3 million tons. The foreign scene, however, has been extremely active during 1969, and this is expected to continue for the next several years. Changes in sinter practice have not been reported, other than a few changes towards increased bed depths and the increased use of screening facilities to size the product more closely. Interest in the field of briquetting in fields other than iron ore is high with a new fluorspar briquetting plant under construction at Brownsville, Texas. Considerable research appears to be pointed towards solving the problem of reclaiming steelmill waste products (flue dust, BOF dust, mill scale) by means other than sintering.

Jan 1, 1970

By B. S. Hickman

Using primarily quantitative single crystal X-ray techniques studies have been made of the precipitation of the metastable w phase in alloys of titanium with Mo, Mn, Fe, Cr, and Nb. It is shown that, in agreement with earlier work on the Ti-V system, the composition of the w phase during aging approaches a constant; the value varies from one system to another in a systematic manner with the electron concentmtion. Results relevant to the mechanism of w precipitation during quenching and during aging are presented and discussed. Results are presented to confirm that the w phase is coherent and the nzorphology of the Necipitates is described and discussed. PRECIPITATION of the metastable w phase in alloys of titanium is of interest because of its severe embrittling effect and also because of the enhancement of the superconducting properties which accompanies precipitation.1 The general features of w phase precipitation have been discussed by cuillan and by Bagaratskii et d3 The work described in this paper forms part of a general study of w phase precipitation which was undertaken to elucidate various unknown or controversial features of the process, namely: i) What is the mechanism of formation of the phase ii) How does its composition, morphology, volume fraction, and particle size vary with the type and content of the alloying element and with the heat treatment conditions? iii) What is the relation between the w phase and the equilibrium a phase? In a previous paper4 the results of a detailed study of w phase in Ti-V alloys were presented. It was shown that: U During aging the composition of the w phase approached a saturation value of 13.5 to 14.0 at. pct V. ii) w precipitated as cube shaped particles and once precipitation was complete, particle coarsening did not occur. iii) a phase formed initially by direct conversion of w-phase precipitates. In this paper similar measurements are reported on Ti-Mo alloys and some more limited data given on the Ti-Cr, Ti-Fe, Ti-Mn, and Ti-Nb systems. The results i.1 the six systems are then compared. It has been shown previously by electron microscopy that the w precipitates as ellipsoids in both the Ti-Nb system,&apos; and in i-o.&apos; Blackburn and williamsz3 have reported the precipitates in Ti-Fe, Ti-Cr, and Ti-Mn are cubic, i.e., similar to the Ti-V system.476 EXPERIMENTAL METHODS 1) Specimen Preparation. Samples of the compositions shown in Table I were prepared by arc melting. High purity electron beam zone refined titanium (Materials Research Corporation Grade I—total impurity content =250 ppm) was used; the other source materials are given in the table. A commercial Ti-8 wt pct Mn alloy was used for studies of this system. After homogenization the materials were rolled or cut into strip approximately 1/4 in. by 2 in.; the strip was then heated at temperatures of 1200" to 1300°C for several hours by passage of a direct current through the strip in an ultra high vacuum system (<10-torr) and then quenched by admitting high purity helium and switching off the current. This treatment resulted in crystals up to 2 in. by 2 in. by A in., almost all of which had a [100] direction nearly perpendicular to the plane of the strip. The crystals were cut from the strip by spark machining and their composition checked by lattice parameter measurements supported by chemical analysis. 2) X-Ray Diffraction Methods. Lattice parameter measurements were made using a modification of the single crystal spectrometer described by ond. In this technique the crystal angles for reflection on either side of a highly collimated primary beam are measured. Generally the 400 CuK, or the 400 CuKp reflections of the bcc lattice were used. In the as-quenched crystals a precision of 1 in 105 was obtained, but due to line broadening the parameters in aged materials could only be measured to 1 in lo4. This accuracy was, however, quite adequate for the purposes of these experiments and enabled the composition of the p phase to be obtained to 10.1 pct in all systems except Ti-Nb where the lattice parameter change with composition is very small. The precipitation and growth of the w phase was examined using a single crystal goniometer attached to a Phillips diffractometer. The relative intensities and line shapes of b-phase and w-phase reflections

Jan 1, 1970

By W. R. Bull, J. R. Roos

A method is described whereby the true financial value of each constituent of each concentrate can be calculated. Each value is then compared with what this value would be in a perfect separation. The difference in each case represents the loss incurred by each aspect of an inefficient separation. The relative magnitudes of these losses indicate those avenues of plant investigation and research that offer the highest potential reward. The metallurgical efficiency of a concentration operation is difficult to express because there are two criteria of performance, namely, grade and recovery. Neither of these alone gives any indication of the effectiveness of an operation and, if they are combined in some formula, it must be empirical and unreliable because of many external factors. Since milling is an operation performed as one step in the most profitable overall method of producing a metal or other end product, the main objective of mill operation should be to make a maximum contribution to this maximum profit. When concentrates are sold to custom smelters, the income from the operation, and the costs, may be calculated. When concentrates are further processed by another division of the same company, and there is no formal financial relationship between, say, the milling and smelting divisions of the company, the mill superintendent&apos;s job is less well defined. But even so, if the mill superintendent works as though he were selling his concentrates, and uses reasonable, if fictitious, conditions of sale, then he can get a good indication of the effectiveness of his operation by calculating his "Economic Recovery". ECONOMIC RECOVERY This is a method used as a basis of operation by several mills and has been described in the literature. Very briefly, the "standard" of operation is taken as that of a mill separating the ore into perfectly pure concentrates with 100% recovery of each mineral in its respective concentrate. For example, a 4%Pb, 4% Zn ore (galena and sphalerite), would provide a lead concentrate weighing 4.62% of the feed and assaying 86.6% Pb and a zinc concentrate weighing 5.97% of the feed and assaying 67.1 % zinc. These perfect concentrates are (hypothetically) shipped to the appropriate smelters, freight and treatment charges and payments calculated, and the net return obtained. This is the optimum financial return for the operation of the mill and is used as the standard. Now the return for the actual, real-life, mill operation is calculated in the same way, and the ratio of the two is termed the Economic Recovery. Many additions and modifications may be made in any particular case. For example, if the zinc mineral is marmatite, the optimum grade of the zinc concentrate may be set lower than 67.1%. If silver is present, it will be recovered 100% in the lead concentrate in the perfect operation. Similarly, cadmium will be recovered 100% with the sphalerite. Mr. Weiss&apos; paper&apos; mentions a case where, because of extremely fine intergrowth between the valuable mineral, chalcocite, and pyrite, the optimum grade of copper concentrate was taken as 45% Cu instead of the theoretical 80%. Each operation will present its own possibilities in this way.

Jan 1, 1970

By R. K. Agarwal

Transforming from a mere novelty field of science to a well-recognized branch of engineering, the scientific discipline of rock mechanics continues to mature. The past decade has witnessed an unprecedented in- fusion of scientific research approach into mine design. In the preceding year alone, there were six national and international congresses which also indicated the increased interest, in rock mechanics and participation by scientific personnel with diversified backgrounds. It may be presumptuous to conclude that the time, effort and money spent thus far on procedures for mine structure design have necessarily resulted in the evolvement of design techniques; but certainly the prospects of evolving such techniques art: more promising at present than ever before.

Jan 1, 1970

By M. J. Stewart, F. Weinberg

PREVIOUS investigations into convective flow in molten metals have examined systems in which the length-to-height ratio of the enclosure is large, usually using long graphite boats.&apos;&apos;2 In normal castings this ratio is generally small and it is not possible to extrapolate the results of previous work to casting configurations. The present investigation was undertaken to examine liquid flow in systems of low length to height ratio. The apparatus used for the experiments is shown in Fig. 1. The molten zone of tin is contained in an aluminum mold of inside dimensions 6.4 by 6.4 by 0.31 cm. A heating block containing a resistance coil is attached to one end of the mold and a cooling block is attached to the other end. These are used to produce a wide range of temperature differences across the mold. The mold is contained in a stainless steel tank which is heated from the outside by resistance heaters to give the required operating temperature during the experiment. Temperature measurements in the molten tin are made with an iron-constantan thermocouple PT equal to 5 atm. The phase diagram derived by us was the same as that found by Krivsky and Schuhmann, except for a similar discrepancy in one boundary line representing the gas composition in equilibrium with SO a. The vapor pressures of the pure components can be neglected at 1700°K, 1 atm and 2500°K, 5 atm. There-fore, in these cases PT is essentially the total pres-sure in the system. This is valid for all the diagrams presented by Krivsky and Schuhmann except for the one at 2500°K and 0.01 atm, where the vapor pressure of silicon is comparable to PT.6 This work was supported in part by the American Iron and Steel Institute. &apos;C. E. Wicks and F. E. Block: Thermodynamic Properties of 65 Elements, U S. Bur Mines, Bull. No. 605, 1963. R. P. Elliott: Constitution of Binary Alloys, 1st supplementary, McGraw-Hill Book Co., New York, 1965. 3W. A. Krivsky and R. Schuhmann, Jr.: Trans. TMS-AIME, 1961, vol. 221, pp. 898-904. 4M. F. Ancey-Moret, M. Olette, and F. D. Richardson: Mem. Sci. Rev. Met., 1964, vol. 61, pp. 169-76. 5J . F. Elliott and M. Gleiser: Thermochemistry for Steelmaking, vol. 1. Addison-Wesley, Reading, Mass. 1960. 6 J.Drowart, G. De Maria, and M. G. Ingraham: J. Chem. Phys., 1958, vol. 29, p. 1015. probe of 0.4 mm diam inserted into the open top of the mold. Radioactive sn113, in the form of a small sphere, was gently lowered into the molten tin at point A in Fig. 1 THERMOCOUPLE> WATER ARGON END HEATER QUENCH ARGON Fig. 1—Apparatus for detecting liquid flow. Point A indicates position the tracer was introduced.

Jan 1, 1970

By W. H. A. Timm

Lake Dufault Mines, Limited, is located ten miles north of Noranda, Quebec, Canada. Following discovery of a copper-zinc orebody late in 1961, metallurgical investigations of diamond drill core were conducted by Lakefield Research Limited during 1962 and 1963. These were followed by pilot-plant tests in May and June of 1964, utilizing underground development ore. Milling operations commenced on August 16, 1964, and reached rated capacity of 1300 tpd by October 1, 1964. Since that time, tonnage was gradually increased to a maximum of 1400 tpd and since January 1, 1968, gradually decreased to the present rate of 1100 tpd due to an increase in the hardness of the ore. The ore milled from start-up until January 1968, was mainly massive sulfides, consisting of a complex mixture of iron, Copper, zinc, and minor amounts of lead sulfides with important associated values in gold, silver and cadmium (average specific gravity of 3-8). From January 1968 until the present, the ore being processed is described geologically as a boxwork of sulfide stringers of chalcopyrite, pyrite, and pyrrhotite with some local dissemination in siliceous rhyolite (average specific gravity 2.9) supplemented by the massive copper-zinc ore. The limited size of the orebody and the mining method employed (open stoping by blasthole and slash benching methods) would have

Jan 1, 1970

By George S. Ansell, Thomas E. Davidson

It has been known for a considerable period of time that the ductility of even quite brittle materials can be enhanced if they are deformed under a superposed hydrostatic pressure of sufficient magnitude. The response of ductility to pressure, however, has been shown to vary considerably between materials. Prior work has shown that the effects of pressure upon the tensile ductility of Fe-C materials depend upon the amount, shape and distribution of the brittle cementite phase. In this current investigation, the effects of pressure upon the fracture mechanisms in a series of annealed and spheroidized Fe-C materials were examined. It was observed that the principal effect of pressure is to suppress void growth and coalescence, retard cleavage fracture and to enhance the ductility of cementite platelets in pearlite. Based upon the observed effects of pressure upon the fracture mechanisms, a proposed explanation for the enhancement in ductility by pressure and for the structure sensitivity of the phenomena is presented and discussed. THE effect of superposed pressure upon the tensile ductility of a variety of metals has been well documented.&apos;-&apos;&apos; Some of the results from several investigators are summarized in Fig. 1 where tensile ductility in terms of true strain to fracture (ef) is plotted as a function of the superposed pressure. As can be seen, a pressure of sufficient magnitude can significantly enhance the ductility of metals. However, Fig. 1 also demonstrates that the response of ductility to pressure and the form of the ductility-pressure relationship varies considerably between materials. Several explanations have been offered for the observed enhancement in ductility by a superposed pressure. Although no experimental evidence was provided, Bridgman13 and Bobrowsky10 proposed that the observed effect was due to the prevention or healing of microcracks or holes. Bulychev et a1.14 observed that cracks and voids in initially prestrained copper were healed in the necked region of a tensile specimen upon further straining while under a superposed pressure. Also, pugh5 observed that large cavities were suppressed in copper fractured in tension while under pressure. A second proposal has been forwarded by Beresnev et at al.6 This proposal is based upon the hypothesis that a material fails in a brittle manner because the normal tensile stress reaches a critical value before the shear stress is of sufficient magnitude to cause plastic flow. Since a superposed hydrostatic pressure will increase the ratio of shear to normal tensile stress, a sufficiently high hydrostatic pressure should favor plastic flow while retarding brittle fracture. Galli15 reported that a superposed pressure shifts the ductile-brittle transition temperature of molybdenum. This was explained based upon the reduction of the normal tensile stress by the superposed pressure. Pugh5 explained the occurrence of the observed pressure induced brittle-to-ductile transition in zinc in the same manner. Davidson et al.12 proposed an explanation for the enhancement of ductility by pressure based upon the effects of pressure upon the stress-state-sensitive stages of various fracture propagation mechanisms. Basically, they proposed that pressure will retard cleavage and intergranular fracture by counteracting the required normal tensile stress or will suppress void growth. They observed suppression of intergranular fracture and void growth in magnesium by pressure. Davidson and .Ansell16 reported ductility as a function of pressure for a series of annealed and spheroidized Fe-C alloys. Fig. 2, from this prior work, demonstrates that the effect of pressure upon ductility is structure sensitive in terms of the amount, shape and distribution of the brittle cementite phase. As shown in Fig. 2, in the absence of cementite or when the cementite is in isolated particle form (spheroidized), the ductility-pressure relationship is linear and the slope decreases with increasing carbon content. In the annealed carbon-bearing alloys wherein the cementite is in the form of closely spaced platelets (pearlite) or in the form of a continuous network along prior aus-tenite boundaries (1.1 pct C material), ductility as a function of pressure is nonlinear (polynomial relationship) in which the slope increases with increasing pressure. At the highest pressures studied (22.8 kbars), the slope of the curves for these materials tends to approach those for the spheroidized material of the same carbon content. In this current investigation the change in fracture mechanisms as a function of pressure for the materials shown in Fig. 2 has been examined. The possible connection between the observed effects of pressure upon the fracture mechanisms and the effect of pressure upon ductility is discussed.

Jan 1, 1970

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

Heats of solution of gold and copper in dilute Au-Cu-Sn alloys have been determined using a liquid metal solution calorimeter. The self-interaction coefficient, Au - has been calculated at constant copper concentrations and n cu has likewise been determined at constant gold contents. Good experimental agreement is obtained between the interaction coefficients and nAu Cc thus demonsbating the reliability of the measured heat values. The measured data are compared with the Predictions of certain solution models. In previous publications,1,2 the results of calori-metric investigations of dilute Ag-Au-Sn and Ag-Cu-Sn alloys have been presented. The present work on the Au-Cu-Sn system concludes a program of studies of enthalpy interaction coefficients in dilute alloys of the Group IB metals with tin. Since the definition and derivation of an enthalpy interaction coefficient has been discussed previously,1,2 no restatement of this theory will be presented here. From the determination of the partial heat of solution of gold and copper in ternary alloys of various copper and gold contents, values of the interaction coefficients can be calculated. These coefficients give an insight into the various solute interactions that occur in the liquid solutions since changes in their magnitude and sign reflect bonding changes that are taking place in alloys of varying solute contents. EXPERIMENTAL Details of the design and operation of the liquid metal solution calorimeter used in this work may be found in a paper by Poo1.3 For the present studies copper of 99.999 pct purity was supplied by American Smelting and Refining Co., gold of 99.999 pct purity by A. D. Mackay, Inc., and tin of 99.99 pct purity by Baker Chemical Co. At the commencement of each series of experimental drops, a tin solvent bath consisting of between 70 and 90 g of the pure metal was inserted in the calorimeter. The weight of the bath was accurately determined and to it were added appropriate amounts of gold or copper to give alloys of the desired composition. For determinations of approximately 0.0015 g-atom samples of Cu were used and for measurements of ?HAu approximately 0.0025 g-atom additions of Au. The heat capacity of the bath was determined at regular intervals during a series of drops using tin calibration samples. Measurements were made of the heat of solution of copper in alloys containing a constant 0.01, 0.02, 0.03, and 0.04 mole fraction of Au, respectively, in order to determine ?HCu in each alloy, and the same mole fractions of copper were used to determine equivalent values for nAu at constant copper concentrations. The composition of the bath was maintained at the desired constant gold or copper content by making calculated additions of the appropriate solute throughout the experiments. The limiting values ?HAu in alloys of constant copper content and of %c, in alloys of constant gold content were studied as a function of the mole fraction of copper or gold respectively in order to determine and nCu. Heat content and heat capacity data used in calculating values of ?ºHAu and ?HCu at the experimental temperature of 720°K were obtained from Hultgren et a1.4 &apos; RESULTS AND DISCUSSION Determinations of ?HAu. The partial heat of solution of gold in pure tin as a function of gold concentration was determined in the previous study of dilute Ag-Au-Sn alloys1 and can be represented by the least-squares expression: ?HAu(l) =-8075 + 2413xAu [l] which is valid between XAu= 0.00 and xAu = 0.05. The standard error in the constant term, which represents the partial heat of solution of gold at infinite dilution in tin,?HºAu(l)is 35 cal per g-atom, while the standard deviation of the slope, which represents n Au is ± 619 cal per- agtom. Corresponding expressions for ?HAu(l) in alloys containing constant mole fractions of 0.01, 0.02, 0.03, and 0.04 copper were obtained from the data listed in Table I and are themselves given in Table II. Fig. 1 illustrates the partial heat of solution of gold as a function of its concentration in each of the alloys. For the four alloys of constant copper concentration, the values obtained for ?HºAU(l) (in order of increasing copper content) are -8141 i 36 cal per g-atom, -8210 ± 42 cal per g-atom, -8202 ± 46 cal per g-atom and -8268 ± 51 cal per g-atom. The corresponding values of the self-interaction coefficient, n Au, for these alloys are 3103 * 644 cal per g-atom, 2425 ± 676 cal per g-atom, 2574 * 717 cal per g-atom and 2523 ± 899 cal per g-atom. In Fig. 2 these values of n Au are plotted as a function of the copper content of the alloys and are seen to remain approximately constant within the experimental limits. The addition of increasing, small amounts of copper to dilute binary Au-Sn alloys thus has no apparent effect on Au-Au interactions in these dilute liquid solutions, although more exothermic values of ?HºAu(l) do result from an increase in the copper content of the alloys. Analogous behavior was observed with additions of silver to dilute Au-Sn alloys.&apos; By

Jan 1, 1970

By J. B. Cohen, S. L. Sass

The nucleation process for y&apos;, Ni3Ti, is shou&apos;n to change from heterogeneous to uniform as the undercooling within the phase boundary increases. As unifornz nucleation beconres copious, (100) alignment is observed in the early stages of aging; howecer. whether or not the alignment is due to spinodal decomposition could not be determined. The presence oi a large-scale reversible segregation of titanium and not the intermediate precipitate y&apos; may be the cause of the initial strengthening in this alloy system and the discrete regions detected magnetically by Ben-Israel and Fine. Evidence is presented for the role of dislocations in the formation at long aging tines oJ the equilibrium precipitate . ThE nickel-rich portion of the Ni-Ti phase diagram consists of a terminal fcc solid solution of titanium in nickel called y.&apos; For temperatures below 1290°C if the solubility limit is exceeded and the temperature is high enough, the excess titanium will combine with nickel to form q, Ni3Ti, a four-layer hcp ordered structure with stacking sequence ABAC."~ Since a change in structure is involved, formation of q is slow and a metastable precipitate referred to as y&apos; forms. This metastable precipitate is an ordered cubic phase in the unconstrained state having a structure similar to Cu3Au (LIZ) with the composition Ni,Ti."&apos;5 From X-ray studies,8, 7 the following sequence is indicated for the precipitation reaction at low temperatures: 1) The appearance of satellites around low-index diffraction spots and splitting of high-index spots, indicating the formation of a slightly tetragonal phase in an imperfectly periodic arrangement. The tetragonal phase is y&apos; or the depleted matrix, depending on the volume fractions of matrix and precipitate (and hence, the alloy composition); tetragonality apparently arises from coherency strains. 2) The appearance of diffraction spots from the equilibrium second phase, q. Ben-Israel and Fine&apos;" used magnetic methods to measure the changes in matrix composition during aging, and to estimate the precipitate composition in a Ni-10.1 at. pct Ti alloy. They observed that the alloy, solution-treated at 1270°C and quenched, contained heterogeneities with discrete compositions which gradually vanished upon aging. It was suggested that y&apos; may initially form as a defect structure with several possible compositions deficient in titanium as compared to Ni3Ti; the composition Ni3Ti develops with prolonged aging. The initial composition was Ni6Ti which would require a very large deviation in stoichiometry for y&apos;. It was noted that at 700°C the precipitation reaction was 80 pct complete after 2 hr and that the volume fraction of second phase was 0.2 after 1 hr. Mihalisin and Decker&apos; suggested that the equilibrium precipitate, besides being nucleated at grain boundaries, could also form at stacking faults. ~errick" examined the formation of in Ni-Cr-Ti, and suggested that a collapsed vacancy cluster within y&apos; acted as a nucleus for intragranular q at high temperatures, 820" and 900°C. After long aging times at lower temperatures, 650" and 750°C, diffraction patterns taken from regions containing ribbons of stacking faults showed spots associated with q. Aging greatly enhances the mechanical properties of these alloys, even at high temperatures, and the metastable precipitate y&apos; is thought to be responsible.&apos; At 600°C the hardness of a Ni-10.1 at. pct Ti alloy increases rapidly during the first 20 hr of aging and continues to increase slowly thereafter.&apos;" At 700°C the hardness reaches a maximum at 1 to 2 hr and then decreases.l1 In this paper the kinetics of formation of y&apos; and the mechanism of transformation to in a Ni-Ti alloy are described. Evidence is presented that suggests that large-scale regions of titanium segregation are present initially and these may have an important influence on the initial strengthening of the alloy: these regions may be the heterogeneities detected by Ben-Israel and Fine. I) EXPERIMENTAL TECHNIQUES The Ni-10.3 at. pct Ti-0.6 at. pct A1 alloy was prepared by vacuum melting. Its chemical analysis is given in Table I. A 0.007-in.-thick strip was sandwiched between foils of the alloy and placed in an open Vycor tube with titanium chips (to get oxygen) and solution-treated for ; hr at 1150°C in a Globar Furnace under a flowing argon atmosphere. Following a quench into an ice-brine mixture at -5" to -2"C, a sample to be aged at temperatures up to and including 700°C was sealed in an evacuated Vycor capsule before being put into a fused salt bath. Aging above 700°C was carried out on bare samples in a salt bath. When removed from this bath and quenched into an ice-water mixture, the sample was normally coated with a thin layer of oxide easily removed with emery

Jan 1, 1970

By Warren E. Holland, Leon H. Robinson

The pore pressure must be included in any discussion of the state of stress within a rock. Rock cylinders at elevated stress levels, loaded uniaxially to a stress less than the yield strength, can be fractured by holding the constraints constant and increasing the pore fluid pressure. This can be described in terms of the Mohr envelopes using effective stresses. The effective stress concept can simplify data with various confining and pore pressures. When the pore and confining pressures are equal, the pore pressure appears to be totally effective in reducing the confining pressure effects. The area over which the pore pressure acts has been the subject of some prior discussions. Terzaghi 1 defined boundary porosity as ". . . the ratio of that part of the area of the potential surface which is in contact with the interstitial liquid and the total area of this surface." His method of calculation included a discussion of splitting failure which was objected to by Mr. Kessler, National Bureau of Standards. The values calculated by Terzaghi were close to unity. Leliavsky,2 using unjacketed cement cores, reported that the interstitial fluid pressure appeared to act upon about 92% of the surface of failure. He concluded: "The fraction of the area of cross-section over which the internal pressure acts can therefore have nothing to do with ordinary porosity, and the results of the experiments tend to support Professor Terzaghi&apos;s views." In 1947 Balmer3 determined the boundary porosity of concrete when the confining and pore pressures were equal. He also found values close to unity. In 1948 McHenry4 published some of the results from the same equipment and included the first data in which the pore pressure

Jan 1, 1970