Search Documents

By R. V. Higgins

This paper presents a computer method to obtain the shape factors and equal cell volumes of the channels for any well spacing pattern from a potentiometric model. By using this program the authors have processed the data for the seven-spot, direct line-drive and the staggered line-drive patterns. The data for the five-spot pattern had been previously processed by a noncomputer method and are included for completeness. The shape factors and volumes for the channels are presented in tables for those who want to use them to process data using their own permeability relationships and viscosities of their reservoir oils. The authors have used the data and sets of representative permeability curves to process sample calculationr of waterflood performances. The comparison of the calculated results shows that the influence of well spacing is small. The permeabilities of the reservoir rock to oil and water had a greater influence on oil recovery for a given pore-volume throughput of water than the well spacing pattern. The more water-wet the reservoir rock, the better the possibility of permeabilities which are conducive to good recovery. The viscosity of the reservoir oil also influences the recovery more than the well spacing pattern. The reduction in the percentage recovery of oil with increase in viscosity of the reservoir oils is small when oil viscosities are in the range of 0.1 to 3. Above this range the reductions in recoveries are extensive. Sample comparisons of the time required for different patterns to recover the oil are presented. Results of an example calculation are given to show the effect of the permeability profile on recovery. INTRODUCTION The effect of well spacing pattern on the recovery of oil when flooding with either gas or water has been studied by many investigators. Muskat et al.1 presented an analysis using conductivity, sweep efficiency and unit mobility to the time of breakthrough. Dyes et al.2 used experimental techniques (X-ray shadowgraphs) and different mobility ratios. They presented quantitatively the relationship between mobility and sweep efficiency at and after breakthrough. Hauber3 presented a method to predict waterflood performance for arbitrary well spacing patterns and mobility ratios. Craig et al.,4 using techniques similar to Dyes et al. to determine sweep efficiency for a five-spot pattern, added the use of relative permeability curves at breakthrough and thereafter. Douglas et al.5 sed rela- tive permeabilities and continuously changing saturations throughout the entire five-spot flood pattern. In obtaining their solutions they used finite-difierence equations. Hig-gins and Leighton6,7 also used relative permeabilities and continuously changing saturations throughout the pattern before and after breakthrough. They employed techniques that process a flood-pattern calculation on the computer in about one minute. The methods of Douglas et al. and Higgins and Leighton both checked closely the laboratory results for a wide range of mobility ratios. This paper presents some sample performances calculated by the Higgins and Leighton method that show the effect on recovery of different permeabilities and viscosities using the seven-spot, the line-drive and the staggered line-drive, as well as the fivespot flood pattern. No previous paper has presented these data using different permeability curves and continuously changing saturations throughout the flood patterns. The paper also presents (1) the results and analyses of the flood-pattern prediction, (2) the computer techniques for determining the shape factors and volumes from the potentiometric models for the foregoing flood patterns, and (3) the shape factors and volumes of the channels of the flood pattern in the event reservoir engineers may like to process waterflood calculations using their own permeability curves and reservoir oils. DESCRIPTION OF METHOD VOLUMES AND SHAPE FACTORS The use of channels taken from a potentiometric model (see Fig. 1) to aid in calculating the performances of water floods of nonlinear patterns has been thoroughly explained in the literature.0,7 Therefore, very little theory, discussion, or proof regarding this phase will be repeated in this paper. The computer method presented in this paper to calculate the volumes and the shape factors of the channels of potentiometric models employs the trapezoidal rule for the volumes and the Pythagorean theorem (the hypotenuse equals the square root of the sum of the squares of the two sides of a right triangle) for the shape factors. In calculating the volume of a channel, the area of each cell in the channel is determined and then multiplied by the thickness to obtain the volume. In determining the areas of the cells, trapezoids are constructed whose vertical sides are spaced Ax apart, as shown in Fig. 2. The length of the sides is the difference between an ordinate cut off by the top and bottom of the cell — usually equipo-tentials. The coordinates of the points along an equipo-tential or streamline are obtained by Lagrange's8 equation of interpolation for which the constants are coordinate points at the intersections. Three intersections for con-

Jan 1, 1965

By Alan T. Gorton, T. L. Joseph, Gust Bitsianes

High-temperature X-ray diffraction techniques were used to determine thermal expansion coefficients of iron and its oxides. Lattice parameters of a and iron, wiistite, magnetite, hematite, and goethite were determined and plotted as a function of temperature. From these results, the true and mean linear thermal expansion coefficients were calculated. In the case of noncubic crystals, the expansion coefficients were determined along each of the crystallographic axes. The effect of the magnetic transformation on lattice parameter was measured for a iron and magnetite at their respective curie temperatures. COEFFICIENTS of expansion for iron and its oxides were determined from the temperature dependence of their lattice parameters. From these data, the true linear coefficients were derived by means of the formula: at =[(l/ao )(dao/dt)] [1] where at is the true linear coefficient of expansion, a, is the lattice parameter in angstrom units, and all variables are referred to the temperature, t, in OC. Thus the true expansion coefficient represents the instantaneous change in lattice parameter on a unit basis and at a particular temperature. A second expansion coefficient is also often used. It is termed the mean linear coefficient and is defined in the following way: am = [ao(t) - ao(O)]/ao(0) l [21 This coefficient represents the average change in lattice parameter on a unit basis but at a specified parameter length, ao at O°C, and over a specified temperature range,1 to 0°C. Generally the true co-

Jan 1, 1965

By Luis Remy Valdivieso

Las minas de Marcona, principales productores de mineral de fierro del Perú, están ubicadas a 400 km, en línea recta al Sur-Este de Lima. Debido a la geología del terreno, los minerales eran virtualmente desconocidos hasta 1952. Ya, con la explotación geológica completa efectuada durante los úitimos diez años, se ha probado que la oxidación y el intemperismo han causado alteración en el tipo de minerales presentes. Asimismo, la gran mayoría de los minerales tiene un alto contenido de azufre que los hacen económicos. En ese sentido, el presente texto señala la comercialización de los minerales en Marcona, donde las pruebas metalúrgicas demostraron que tales minerales pueden ser sometidos a un tratamiento que los convierte en productos comerciales. Por ello, ante el requerimiento más exigente de los compradores referido a especificaciones físicas y químicas del mineral, la Marcona Mining Company decidió instalar una planta de beneficio capaz de impartir a los minerales provenientes de las minas que explota el más alto grado de pureza conocido actualmente. Así, la planta de San Nicolás produce cuatro tipos distintos de mineral de fierro para uso en altos hornos o bien en hornos de sinter, y como consecuencia, se aseguró la vida al distrito minero de Marcona.

Jan 1, 1965

By Harry C. Gatos, Frank J. Bachner, Mario D. Banus

The portion of the Nb-Sn phase diagram between 75 and 79 at. pct Nb at temperatures near the liquidus has been investigated by melting alloys of known composition and examining the microstmc-tzlres resulting when these melts were quenched. The Nb-Sn samples were contained in magnesia sleeves, and were mounted in a high-pressure apparatus. It was shown that Nb3Sn me1ts by a peri-tectic decomposition and that the composition of the compound Nb3Sn at the peritectic horizontal lies between 77 and 78 at. pct Nb. Analysis of the phases present in the samples was accomplished by anodic staining techniques and an electron-microbeam analyzer. DURING the past several years a great deal of work has been devoted to the Nb-Sn system, and the bulk of this work has been primarily concerned with the intermetallic phase Nb3Sn. However, the melting characteristics of Nb3Sn have not been investigated because of the difficulty in maintaining known compositions due to tin losses by vaporization, and the added difficulty of containing Nb-Sn alloys without sample-container reaction at the high temperatures involved. This paper reports on the development of a technique for containing Nb-Sn alloys under pressure high enough to prevent tin losses at temperatures in excess of the liquidus. Using this technique, melts of known composition were made and a phase diagram near the liquidus is proposed based on the microstructures resulting when these melts were quenched. EXPERIMENTAL The Nb-Sn alloys weighing 0.7 g each were prepared in the form of cylinders. Niobium powder, -200, +325 mesh, and tin shavings with a minute amount of tin powder for a binder were mixed in the desired proportions and cold-pressed at 50,000 psi in a 0.150-in. die. This cylinder was then enclosed as shown in Fig. 1, and the entire arrangement mounted in a 2.6-in. pyrophyllite tetrahedron.' The desired temperatures (-2500°C) presented several problems in high-pressure container materials. Boron nitride, the usual sample holder, was found to react with niobium to form the nitride. However, MgO cans and plugs machined from refractory rods did not melt nor react with the sample. Thick tantalum contact disks were used to carry the large currents (~450 amps) to the graphite heater, since tantalum did not melt and had enough strength to prevent extrusion of the sample. It was further observed that the graphite heater and tantalum disks could not be in direct contact with the pyrophyllite because it melts and decomposes at these temperatures to kyanite and coesite (SiO2),2 the former reacting with graphite and the latter with tantalum to form the silicide. The use of a MgO outer sleeve

Jan 1, 1965

By A. E. Dwight, A. F. Berndt

Binary and ternary alloys of uranium with transition metals were prepared with U2X stoichiometry. The compounds U2Tc, U2Rh, U2Os, and U2lr were formed by peritectic or peritectoid transformations, and were shown to be isostructural with U2Ru. The lattice constants are presented. Miscibility relationships obtained from the study of ternary alloys are graphically summarized in terms of an area-of-stability plot. The compound U2Re is not isostructural with either U2Ru or U2Mo. A series of compounds with stoichiometry U2X has been prepared in which X is technetium, rhodium, rhenium, osmium, or iridium. The only other U2X compounds, for which X is related to the above transition metals by proximity in the periodic table, are U2Mo and U2Ru. The structure of U2Mo is the tetragonal C11b type (MoSi2 type)1 and U2Ru is monoclinic, with a new structure type.' This investigation was undertaken to determine if any of the above compounds are isostructural with either U2Ru or U2Mo and to find the range of occurrence of these phases. EXPERIMENTAL Alloys were prepared from stoichiometric amounts of electrolytic uranium (99.99 pet pure) and transition metals of 99.5 pet purity (98+ pet pure Tc), by arc melting on a water-cooled copper hearth under an Ar-He atmosphere. Weight losses were negligible. The buttons were homogenized at temperatures between 700o and 800°C (900°C for U2Tc) in evacuated capsules and air-cooled. The X-ray powder specimens were annealed briefly at the temperature of homogenization. Except in the case of the technetium alloy, metallographic examination was used to determine the number of phases present. The radiation hazard involved in handling technetium precluded the use of metallography to examine U2Tc. The temperature and mode of formation of the U2X phases were determined by heating the metallographic specimens to successively higher temperatures, until examination indicated the presence of chilled liquid or a peritectoid transformation. RESULTS The U2X compounds investigated were found to be formed under the following conditions: UzTc unknown U2Ru peritectic 897o ± 3°C U2Rh peritectoid 755o ± 5°C U2Re peritectoid3 below 750°C Uas peritectic4 920°C U2Ir peritectoid 775o ± 5oC Debye-Scherrer patterns for the above U2X phases were made with CuKa radiation (nickel filter). The similarities in the positions and intensities of the strongest lines in the front-reflection region strongly suggested that U2Tc, U2Ru, U2Rh, U2Os, and U21r are isostructural with only minor differences in the unit-cell dimensions. The known lattice constants and intensities of the Debye-Scherrer rings for U2RU2 were used as a guide in order to index the low-angle X-ray reflections of the U2TC, U2Rh, U2Os, and U2Ir powder patterns. Approximate lattice constants for these compounds were thus obtained and were used along with the intensities observed for U2Ru to index the strongest lines in the back-reflection regions. A sin2 2? correction term was included and the lattice constants were refined and standard deviations estimated by least squares. Comparison of observed and calculated values of Q(= 4 sin2?/?2) used in these calculations are given in Table I. In each case when the positions of lines could not be resolved on the films they were not measured. The final lattice constants, with their standard deviations, are given in Table 11, along with the unit-cell volume and X-ray density. The values of the atomic volume5 of the transition ele-

Jan 1, 1965

By R. A. Sharma, F. D. Richardson

Manganese oxide acitivities and sulfide capacities have been measured at 1650°C in mixtures of MnO + A1203 and MnO + SiO2 - Al2O3. Sulfide capacities have also been measured at 165O0C for melts of MnO + SiO2, MgO + SiO2, MgO + SiO2 + Al2O3, MgO + MnO + SiO2, and SiO2 + Al2O3 and for crystalline CaO. The activity coefficients of the sulfides of manganese and magnesium have been derived where possible. It looks as though the activity coefficients of sulfides are greatest at compositions where the free energies of formation of the silicate mixtures have their maximum negative values. The amphoteric nature of alumina is reflected in the results for the binaries and ternaries in which it is present with silica. The limiting solubility of CaS in CaO at 1650°C is 6.3 x 10-4 mole fraction. THE activities of manganese oxide have been measured in a variety of silicate and aluminate melts by use of the equilibrium (MnO)slag = [Mn]pt .½O2 [l] and the experimental techniques described previously.' Sulfide capacities (CS = wt pct S (PO2/PS2)½) have also been determined in a series of melts, by use of an established method2 which exploits the equilibrium (O°) + ½S2 = (S") . ½O2 [2] If CS and the activity of a metal oxide are known in a melt, the activity coefficient of the metal sulfide can be calculated by use of the equilibrium constant K, for the pure phase reaction MO i ½S2 = MS + ½O2 [3] The activity coefficient of the sulfide is then given by the equation

Jan 1, 1965

By W. W. Mullins, F. A. Nichols

Solutions are developed, assuming surface diffusion and both internal and external volume diffusion, for the relaxation of bodies slightly perturbed from spherical and cylindrical geometries. Combined with those previously published for the nearly planar case, these results provide a means of gaging the relative contributions of the two diffusional proc-cesses in any given case. It is shown that in all sintering experiments to date, and probably in any attainable in practice, surface diffusion has played the dominant role, although most previous authors have assumed otherwise. It is also shown that surface diffusion predominates in normal field-emission tip blunting and also for the coalescence of gas bubbles introduced into metals by a bombardment. The surface-diffusion solutions for a perturbed sphere are combined with previous results for volume diffusion to show that the inclusion of interface diffusion permits considerably larger spheres to develop in diffusion-controlled precipitate growth before the onset of instability. A mechanism is also proposed for the spheroidization of precipitate platelets as well as rods. In a previous paper1 the relaxation of a nearly plane surface to flatness by the combined action of the transport processes of viscous flow, evaporation-condensation (in a closed system), volume diffusion, and surface diffusion has been analyzed under the assumption that all surface properties are independent of orientation. In this treatment, criteria were developed for deciding which process predominates, and solutions valid in the latter stages of the sintering of small wires and particles to a plane were obtained. A numerical solution, valid throughout the entire particle-sintering process for the case of surface diffusion, was subsequently obtained by the present authors.' It was found that the analytic solution (which assumed small slopes everywhere) is accurate to within -10 pct when the maximum slope of the profile is less than 0.3; the wire-sintering problem has also been solved nu- merically for the case of surface diffusion and here again the results converge to the analytic small-slope solution at late stages of the process, the two solutions agreeing in this case to within 10 pct when the maximum slope of the profile is less than -0.6. The purpose of this paper is to extend the perturbation solutions to nearly spherical and nearly cylindrical geometries. We treat only the two principal diffusional processes, i.e., surface and volume, but for these geometries we discuss volume diffusion both inside and outside of the solid. Our results, when coupled with Mullins' solutions1 for nearly planar surfaces, provide criteria for gaging the relative contributions of surface and volume diffusion to the over-all transport process in three basic geometries. A very interesting feature in the cylindrical case is the occurrence of instability for longitudinal perturbations with wavelengths greater than the cylindrical circumference, a classical result. This instability of the cylindrical surface is applied to give a quantitative explanation for the often-observed "erratic" pore closure in the late stages of the sintering of wire compacts; also, the theory previously presented for the spheroidization of rod-shaped precipitates by surface (interface) diffusion' is expanded now to include volume diffusion inside and outside of the particle. The results for circumferential perturbations on a long cylinder allow quantitative estimates for gaging the relative importance of surface or volume diffusion in the early stages of the sintering of spheres or wires. The results here demonstrate clearly that surface diffusion has played a very important, if not dominant, role in all sintering experiments discussed in the literature, although the surface-diffusion contribution to the kinetics has usually been ignored. The results for the sphere (surface-diffusion case) are added to the results obtained previously by Mullins and sekerka3 concerning instabilities of a growing spherical precipitate particle (with interface diffusion disallowed) to obtain a general solution to this problem including interface diffusion. The inclusion of interface diffusion is found to increase significantly the range of stability of a growing spherical precipitate for typical metallurgical cases. The following assumptions are made: (i) the initial surface lies everywhere near and has a slope differing only slightly from that of the reference

Jan 1, 1965

By J. L. Stokes, L. J. Demer

It has recently been shown' that epitaxial films of silicon carbide deposited on silicon single crystals possess desirable properties for semiconductor applications. There is one apparent disadvantage to such a system, however; since epitaxial deposition is a high-temperature process, it can result in plastic deformation of the silicon substrate. The accompanying Berg-Barrett X-ray diffraction photographs reveal plastic deformation in two silicon semiconductor crystals which occurred when the crystals were coated with an epitaxial layer of silicon carbide at 1150°C. Due to differences in thermal expansion between the silicon substrate and the epitaxial silicon carbide, the silicon was subjected to stresses while still in the plastic temperature range. Both crystals shown were 0.7 in. in diam and 0.0055 in. thick. The crystals were sliced from a larger crystal so that their polar axes were parallel to the [ill] direction. After cutting, the crystals were polished with diamond paste, and etched in HC1 vapor to remove surface damage re- sulting from the cutting and polishing operations. This provided clean, damage-free surfaces on which to deposit the silicon carbide. Figs. l(a) and l(b) are reflection and transmission photographs, respectively, of the same crystal. A smaller portion of the crystal is shown in the transmission geometry, because a vertical slit was

Jan 1, 1965

By D. R. Mathews

The solubility of carbon in electrorefined vanadium was determined to be approximately 0.27 wt pct at the eutectic temperature 1,650° C and 0.03 wt pct at 700° C. The tensile strength and hardness of vanadium increased with car-bon content to the solubility limit, but this effect was small, probably because of the narrow solubility range. Between the solubility limit and 0.3 weight-percent carbon, the second phase had essentially no strengthening effect in vanadium. The second phase, observed metallographically at the grain boundaries and within the grains, was shown by autoradiography to contain high-carbon concentrations. This second phase was identified as VgC by X-ray analysis. The VgC in the grain boundaries contributed largely to the poor cold-and-hot workability of vanadium-carbon alloys. Increased concentrations of V2C occurred at the grain boundaries of the high-carbon alloy ingots as a result of homogenization at 1,600° C for 1 hour. Alloy ingots containing more than 0.36 wt pet carbon (a calculated 2.9 wt pct VaC) following homogenization, cracked severely and failed by intergranular fracture upon swaging at 900° C.

Jan 1, 1965

Although the pattern of metamorphic zones south of the Thackaringa-Pinnacles Fault is much obscured by polymetamorphism, displacement of the orthopyroxene isograd by the fault agrees with the displacement of a horizon of "banded aplite" mapped by Andrews (1922). A dextral transcurrent component of 20,000 ft is indicated.INTRODUCTIONThe east-west Thackaringa-Pinnacles Fault (Andrews, 1922; King and Thomson, 1953), one of the major geological structures in the Broken Hill district, New South Wales, cuts across the general north-easterly strike of outcrops and fold structures of the Willyama Complex. The degree of shearing and recrystallization in its unusually wide crush zone and the manner in which adjacent rock units and structures have been dragged into conformity with its strike denote considerable movement on the fault. Information concerning the nature of this movement is of critical importance in the search for possible continuations or repetitions of lead-zinc mineralization south-west of Broken Hill.A transcurrent movement of 20,000 ft, with the southern block moving west, was estimated by Andrews (1922), who correlated prominent outcrops of "banded aplite" at The Pinnacles (north of the fault) and Big Hill (south of the fault). Displacement of this order, however, does not seem sufficient to explain the extent and character of the fault zone. Recent chemical and petrological studies of the banded aplites at The Pinnacles and Big Hill by Vernon (19B1) confirmed their similarity, but indicated an origin by metasomatic replacement of pre-existing folded and metamorphosed rocks. This raises the possibility that some...

Jan 1, 1965

By B P. Scott

Podsol soil horizons developed on Pleistocene glacial sediments overlying the Coronation Mine orebody were analyzed for copper and zinc spectro-photometrically, using biquinoline-isoamyl alcohol and dithizone carbon tetrachloride reagents. Copper and zinc are concentrated in the A. and C horizons, whereas the A1, A2 and B horizons are leached. Of economic significance is the fact that the C horizon contains a copper anomaly, in excess of 100 times back-ground, directly over the orebody. Due to the erratic distribution of high copper concentrations in the Ao horizon with respect to the orebody, that horizon is not considered as a reliable prospecting datum. U NDER certain conditions, metallic ore deposits are accompanied by primary halos of trace elements in unweathered rock surrounding the orebodies or by secondary halos in weathered rock, alluvium, soil, or vegetation above or near the mineralized zones. Where such halos occur, chemical analysis, as applied to the detection of trace elements, is a valuable prospecting aid, permitting the detection of mineralization of possible economic value.

Jan 1, 1965

By G. R. Speich, R. A. Oriani

The rate of coarsening (Ostwald ripening) of copper precipitate Particles in an a-iron matrix has been studied in three Fe-Cu alloys containing 2.3, 4.0, and 5.4 wt pct Cu in the temperature range 730° to 830°C. The kinetics of Ostwald ripening ohevs the law as deduced by Lifshitz and Wagner so that matrix diffusion controls the process. The equilibrium shape of the particles is that of rods with hemispherical caps, and is due to the BOTH Lifshitz and slyozov and wagner3 have recently solved the complex problem of determining the rate of coarsening of an assembly of precipitate particles. Both workers, although using somewhat different mathematical approaches, arrive at the same equation for the relation between the mean radius 7 of a distribution of spherical particles and the bulk solubility co, diffusivity D in the matrix, large anisotropy in the copper/a iron interfacial free energy (y,/y, = 4). A Lifshitz -Wagner equation, modified for both the nonsphericity of the diffusion field and the reformulation of the Gibbs-Thomson equation for the case of nonsphencal particles, is derived. Using the modified equation, the calculated values of the mean radius are smaller than those observed by a factor of about one -half. absolute temperature T, gas constant R, surface free energy y of the precipitate/matrix interface, and the time t of aging: where V, is the molar volume of the precipitate and v is a stoichiometric factor. Wagner's treatment is restricted to the case of spherical particles growing in a liquid matrix where strain energy is not a factor. Lifshitz and slyozov have discussed the effect of anisotropic particle shapes and internal strains on the rate of coarsening. Oriani4

Jan 1, 1965

By A. L. Jr. Pozzi, A. O. Garder, D. W. Peaceman

A new numerical method is proposed for the solution of multidimensional miscible displacement problems. Besides the usual stationary grid associated with numerical procedures, the method uses moving points for which positions and concentrations are computed each time step. Based on the method of characteristics for treating combined transport and dispersion, the method applies equally well to any number of dimensions and, in contrast with other numerical procedures, properly takes into account the physical dispersion, no matter how small. Extensive tests show that accurate solutions of one-dimensional problems can be obtained over a wide range of values of the dispersion coefficient, including zero. They also show that the moving points do not need to be uniformly spaced, and that increasing the number of moving points beyond 2/grid interval does not significantly improve the accuracy of the solution. Two-dimmsional calculations using the new method were carried out to simulate laboratory model displacements in which gravity caused overriding of solvent. Good agreement between observed and computed profiles and measured and computed concentrations of produced fluid was obtained for a viscosity ratio of 1.89; at higher viscosity ratios, fair agreement was obtained Computed results were found to be significantly affected by vertical dispersion, and negligibly affected by dispersion in the direction of flow, in agreement with conclusions based on experimental observations. INTRODUCTION Although miscible displacement processes are potentially of great economic significance, limited progress has been made in developing mathematical methods for predicting the outcome of a solvent flood. Because of the complexity of the partial differential equations which describe multidimensional miscible displacement, numerical methods are a natural approach to their solution. The task of finding suitable numerical approximations to one of the partial differential equations, which involves both dispersion and convective transport, has been particularly difficult. Dispersion in the absence of transport is described by the heat flow equation, a second-order equation of parabolic type. This equation has been very successfully treated by numerical methods. Transport in the absence of dispersion is described by a first-order equation of hyperbolic type. This equation has been treated numerically with some success in one dimension, both by Lagrangian and Eulerian techniques, 1.2 but extension to two or more dimensions has not been satisfactorily accomplished. Usually one of two things happens: either the numerical solution develops oscillations or it becomes smeared by artificial dispersion resulting from the numerical process. When transport and dispersion are considered simultaneously, this numerical dispersion may dominate the low physical dispersivicy that generally characterizes miscible displacement. Stone and Brian2 have developed a new difference equation for treating combined transport and dispersion in one dimension. By suitable choice of certain parameters, they succeeded in markedly reducing both oscillation and numerical dispersion. However, a satisfactory extension to higher dimensions has not yet been found. Peaceman and Rachford3 presented a centered-in-time and centered-in-distance difference equation combined with a 1,transfer of overshoot" procedure which was demonstrated to work well in one dimension. The technique was applied to the two-dimensional miscible displacement problem, but subsequent calculations with zero dispersivity showed little change in results compared with the calculations they presented. This indicates that their method involved a numerical dispersion of the same order as the physical dispersion. The present paper describes an alternate method for solving the combined transport-dispersion

Jan 1, 1965

By A R. Bullis

Metal Mines Limited (Bancroft Division) was incorporated m 1964 to take over the assets of Faraday Uranium Mines Ltd., including their uranium deposit near Bancroft. Ontario. This paper outlines the geology of the Bancroft area, and. then go on to describe the principle rock types m the mine workings, especially the pegmatite, which is the host rock for the uranium mineralization. The structure and mineralogy of the orebodies are thoroughly covered, and, after a summary of the production and ore re-serve figures, the paper concludes with suggestions for further exploration work. Location and History M ETAL Mines Limited, Bancroft Division, is located in Faraday township, Hastings county, Ontario. The main production shaft and concentrator are about: Y2 mile north of Highway 28 and about 5 miles west of the village of Bancroft (Figure 1). The presence of radioactive minerals in the Ban-croft-Wilberforce area has been known since the early 1920's. Uraninite was found at the Richardson deposit near Wilberforce in 1922, and some attempts were made to produce radium from the deposits of the area between 1929 and 1931.

Jan 1, 1965

By Edward W. Parsons

The explosive-anchored rockbolt was field tested in six mines where rock conditions varied from incompetent shale or mudstone through fairly soft sand-stone to firm sandy shales. The purpose of the testing was to document, under a wide variety of conditions and in comparison to other rockbolt anchors, the quality of the anchorage produced by the explosive-anchored bolt. In addition, the length, strength, and design of the explosive charge necessary to produce positive anchorage up to the yield strength of the weakest component of the bolt were determined. A l4-inch length beaded explosive charge was found to be suitable, in general, in relatively weak strata, and the length can be reduced in more competent strata. Of the 127 explosive?-anchored bolts tested, 118 bolts produced an evaluated anchorage in excess of 17,000 pounds. Only 4 of 43 expansion-shell bolts and 4 of 14 wedge bolts produced anchorage this good. In preliminary tests, properly coupled explosive-anchored bolts did not ignite an explosive gas mixture. Some minor phases included in the investigation were as follows: (1) effect of explosive-anchored rockbolt on surrounding rock; (2) effect of time on anchored bolts based on torque tests; (3) effect of submersion in water on rockbolt anchorage; (4) pull tests to higher loads; and (5) ignition tests.

Jan 1, 1965

By R W. Hutchinson

Cupreous pyrite deposits in Cyprus have many geologic similarities to Canadian massive sulphide copper-zinc ores as well as certain significant differences. A comparison of the Cyprus and Canadian deposits indicates a complex genetic history for the latter and suggests that their origin should be re-examined. The Cyprus ores occur in subaqueous extrusive rocks near and at the top of a Cretaceous volcanic pile or igneous complex. The rocks have not been deformed, nor have they been deeply buried or deeply eroded. Apparently the ores formed near or at surface in the top of the volcanic pile by a volcanic-exhalative process and from an ore-fluid of volcanic derivation. The Canadian deposits are of comparable composition and have similar geologic set-tings, but occur in the most ancient Precambrian rocks. Their enclosing Keewatin greenstones have, in most places, been deformed, regionally metamorphosed and invaded by granitic intrusions. Comparison suggests that the Canadian ores also formed initially near surface in ancient Keewatin volcanic complexes by a syngenetic, volcanic-exhalative process. Subsequently, during orogenesis, they were deformed, remobilized under conditions of elevated pressure-temperature and re-emplaced as epigenetic bodies into

Jan 1, 1965

The Shortage of Mining and Metallurgical Engineers I N many parts of the Western World attention has recently been drawn to the acute shortage of mining and metallurgical engineers produced by the universities. Naturally, the reasons for this shortage have been discussed, and various remedies for it have been proposed. In this paper, I want to examine some of them. It is now fairly generally ac-cepted among mining men that the easy attractions of modern urban life have left the career of mining in the metaphorical as well as in the literal backwoods. The industry finds it increasingly difficult to recruit adventurous and intelligent young men; for such men now seek their intellectual, if not physical, adventure in more glamorous occupations from a safe social base in suburbia. The salaries, opportunities and security offered by the mining industry do not complete effectively with the rewards offered by other industries and by commerce, taking into account the social disadvantages of a mining life. The excitement and other non-material attractions of mining

Jan 1, 1965

By W. F. Sullivan, J. R. Coleman

The use of radioisotopes in obtaining basic information on large scale industrial processing is illustrated with two examples. The first is concerned with the operating characteristics of a kiln scrubbing tower and the second with those of a Hardinge ore mill. The examples show the advantages associated with the industrial use of tracers as well as the limitations. The Atomic Energy Commission regulations governing the use of radioisotopes in this manner are also discussed. The nuclear properties of radioactive isotopes have provided a new technique for obtaining basic information on large scale industrial processing. A small portion of the process stream is made radioactive by mixing with a suitable radioisotope or by irradiation in a nuclear reactor. The "tagged" portion is then part of the normal process stream and can be followed by detection and measurement of the associated radioactivity. The advantages of radioactive tracing over other forms of tracers are well known, and include such factors as chemical similarity, capability of precise determination of very small quantities due to the high sensitivity of detection apparatus, and the absence of foreign element contamination of the final product. This report will illustrate two applications of the use of radioactive tracers in large scale industrial processing which have provided information that would have been hard to obtain by other means. The first is concerned with the operating characteristics of a kiln scrubbing tower and the second with those of a Hardinge ore mill. ATOMIC ENERGY COMMISSION REGULATIONS Atomic Energy Commission approval can be obtained by qualified personnel if the proposed quantity of radioactivity does not exceed AEC's standards for

Jan 1, 1965

By M A. Notte

Sand casting of steel is a process that offers the greatest design flexibility and highest metallurgical versatility. However, competition from other methods of fabrication is quite keen, constantly forcing the steel casting industry to the production of better quality castings at lower costs. New methods and materials have been put to use in the last ten years by foundries everywhere. This paper, in the form of a progress report, deals with new methods, processes and materials used by Canadian steel foundries for regular production, as well as with developments still in the initial but promising stages. T HE casting of steel to shape is the most direct method of producing steel parts in their final form. It is a basically economical process both as to cost and production time; it is also a process that offers design engineers the greatest flexibility and highest metallurgical versatility. In North America, the steel casting market has a normal yearly sales volume of 1,500,000 tons. In decreasing order, the largest customers are the railroads, construction machinery and earth-moving

Jan 1, 1965

During the past decade, the miner-al industries have been faced with an increasing shortage of engineers. The point has now been reached where the shortage may be termed critical. This shortage may be attributed directly to declining enrollment in the mineral science programs offered by the universities on this continent. Although the reasons for declining enrollment are undoubtedly many, two major hypotheses may be postulated. The first is that high-school graduates are either unaware of the mineral industries and the opportunities in them or that they are aware of them but have somehow obtained a false image of their true nature. The second is that they are aware of the true nature of these industries and the opportunities they offer but find them exceedingly unattractive. The problem of dwindling enrollments has, in recent years, been subjected to a concerted attack by the industries and professions concerned. However, their efforts to date have been based essentially on the first hypothesis. Sufficient results are now available to suggest that either the methods of attack are inappropriate

Jan 1, 1965