Search Documents

By T. B. Armstrong

Electrical energy conservation sees that service facilities are operated only when needed and that production equipment is operated under an optimum condition. Specific energy consumption (kw-hr per tonne) is not regularly recorded for comminution circuits because of variations in ore characteristics and the low grinding efficiency (0.1 to 2.0%) of the equipment. Improved data acquisition can produce meaningful results. Particle size monitor systems and density gauges permit comparison of energy (kw-hr) and production (tonne) for a "fixed" condition. Finding the operating condition producing the lowest unit energy consumption provides a target for an effective energy management program.

Jan 1, 1982

By W. E. Wickerman

The Magnetic Storm Monitor is an instrument that continuously records variations in the earth's total magnetic field at a fixed location. It is intended for use in conjunction with airborne magnetometer surveys. There are stormy periods when changes in the earth's magnetic field are so violent that data obtained with airborne equipment would be of little value. The Monitor provides this information so that flights may be cancelled or terminated, thus eliminating useless flying. Moreover, when some unpredictable change in the earth's field occurs during an airborne survey, the Monitor record, showing a similar event, may serve to salvage the data.

Oct 1, 1955

By Robert Uhley, Tsvi Meidav, L. Scharon

Before National Lead built an industrial plant on its Fredricktown property, some 100 miles south of St. Louis, a 750x500-ft area on the proposed site was investigated by electrical resistivity, seismic refraction, and churn drilling. Data were accumulated to prepare an isopach map of the overburden and to determine soundness of the dolomite bedrock and the depth to bedrock at fixed points. The site was on a gently sloping hillside. Elevations above sea level ranged from 760 ft in the northwest corner to 820 ft along the eastern edge, giving a total relief of 60 ft. There were no major drainage features, but in the northern part of the plant site there was an indicated depression area which later was interpreted as reflecting subsurface conditions in the bedrock. No rock was exposed in the entire plant site area; the surface was characterized by soil and overburden. GEOLOGY Prevailing geology in the area of investigation was relatively simple—Bonneterre dolomite overlain by residuum. The red overburden was that derived from the weathered Bonneterre. Some chert was present, interpreted in the general area as a weathering of both the Bonneterre dolomite and younger sedimentary rocks such as the Potosi dolomite. During the period of actual measurements the residuum varied from moist near surface on the western edge of the proposed site to dry near surface on the eastern edge, and for this reason overburden resistivities ranged from low on the west to high on the east. Seismic velocities were low, since the material was not well packed. The bedrock encountered was dolomite, specifically the Bonneterre of Upper Cambrian age which, according to McQueen,' consists of an upper dark gray member and a lower light gray member. These have been further divided into four zones, in ascending sequence a dark gray, sandy argillaceous dolomite; a dark brown or black, fine-grained dolomite: a light gray, finely crystalline, compact dolomite; and, at the top, a dark brown, finely crystalline, compact dolomite. In the area of investigation the upper two zones were encountered in the drilling tests as constituting the bedrock. The dolomite varies widely in porosity, which is low in the dense, fine-grained layers. In this area, therefore, both electrical resistivities and seismic velocities were high. Fracture zones in the Bonneterre are nearly Geology, and may extend to some depth, causing considerable deep weathering. It was realized that over deep fracturing the apparent resistivities

Jan 1, 1960

By Beverley S. Randolph

The coal-measures of George's Creek coal-region, in Maryland, lie in a comparatively flat synclinal basin, about 4 miles wide. The dip of the measures is, in places near the edge of the basin, as much as 20 feet in 100. The Pittsburgh seam, known locally as the " Big Vein," is the only seam worked at the present time. Ocean No. 3, one of the large mines of the Consolidation Coal Company, has its opening on the outcrop of the Pittsburgh seam, in a ravine which cuts into this seam to such an extent as to bring this opening to within a mile of the axis of the basin. After working all the coal to the rise of this point, the mine was further developed by extending a slope down the dip of the seam, making such an angle with the line of greatest dip as to give a grade varying from 1 to 10 per cent., down which the empty cars carry a rope from a single-drum winding-engine. The coal was worked from each side of this slope in successive Iifts during its progress. From the top of the slope, for a distance of 4000 feet to the fourth lift, the " make " of water is drained by gravity through a heading into the adjoining mine. Below this point the water was pumped into this same heading, at the fourth lift, from which point, as before mentioned, it flowed by gravity.

Jan 1, 1895

By David C. Jonson, W. Bruce MacKenzie, Arthur A. Bookstrom, Vaughn E. Surface, Neil K. Muncaster, Stewart R. Wallace

In mid-Tertiary time a wet silici-alkalic magma penetrated the Precambrian rocks of what is now the Tenmile Range of Central Colorado and formed the Climax Stock. The stock is a composite one and was intruded in four separate irruptions, each giving rise to a major hydrothermal event. The first three of these produced large stockwork ore bodies that cap and flank the tops of their respective source intrusives; the fourth, and last, yielded almost no commercial mineralization. The ore bodies and their related zones of alteration are circular to ring-shaped in plan and arcuate in section. The different phases of the stock are so positioned with respect to each other that many of the circular and arcuate zones of mineralization and alteration overlap or are juxtaposed. Evidence of a time separation between the different ore bodies is found in certain dikes that were intruded between periods of mineralization and are well mineralized in one ore body but are barren in another. If it is assumed that the parent magma at Climax contained 10 per cent "water," that the "water" contained 0.1 grams per liter of molybdenum, and that the processes of extraction were 100 per cent efficient, then approximately 10 cubic miles of magma would be required for the formation of each ore body. The volume requirements, taken in conjunction with the cross-sectional areas of the different phases of the stock to explored depths, suggest to the authors that each magma chamber was columnar, extended to a depth of thousands of feet, and was connected with a much larger master reservoir. The repetition of similar intrusive, structural, and hydrothermal events, indicates periodic tapping of this reservoir. A regular chemical evolution of the parent magma is suggested by the systematic change in metal ratios and differences in the sites of deposition, relative to the source intrusive, in each of the four separate igneous hydrothermal stages. The preparation of ground was due to changes in magmatic and/ or hydrothermal pressures within the igneous-hydrothermal columns; small, repeated, up and down movements of the rock enclosing the tops of the columns formed the stockwork fractures into which the mineralizers penetrated. The forceful emplacement of each successive phase of the Climax Stock resulted in a slight doming and rotation to the westward of the earlier-formed geologic features. Postore normal displacement along the Mosquito fault is estimated at about 9000 feet. Part of one ore body was cut by the fault and lies at depth below the Paleozoic sedimentary strata on the hanging wall. On the footwall side of the fault, subaerial erosion and glaciation removed a large part of one ore body and uncovered the top of a second.

Jan 1, 1968

By Victor deWolfe

INTRODUCTION The Henderson Mine, located near Empire, Colorado, utilizes a continuous panel caving system to extract ore as one of the world's major producers of molybdenum. Any mine using a caving-by-gravity technique of mining must rely on closely controlled draw of the caved ore. This control is essential to insure proper caving action, to avoid damaging load concentrations of weight and to minimize the dilution of ore with waste material. Henderson is no exception. Draw control is a major factor in all production planning, from long- range plans to short-range and day-to-day ore scheduling. Draw control is reviewed constantly and administered daily in an effort to optimize production efficiency, ore recovery, and cave management. MINING METHOD The cave at Henderson is massive, moving slowly through large panels that are 244 m (800 ft.) wide by 610 m (2,000 ft.) long. Generally two cave areas are drawn at one time. The areas under active draw vary in size but can be as large as 244 m (800 ft.) by 244 m (800 ft. ) containing 400 draw points. Each draw point contains 45,360 mt (50,000 st) on the average and takes about two and one half years to exhaust. A complete panel is worked for seven to ten years. No pillar exists between panels, but rather a buffer zone of broken ore, or "static face," is left in each panel to be drawn with the adjacent, yet-to-be-caved panel in efforts of minimizing dilution of a working area from an exhausted one. (Figure 1) Production drifts are driven on 24.4 m (80 ft.) centers through the ore body. Between the production drifts are funnel-shaped draw bells on 12.2 m (40 ft.) x 24.4 m (80 ft.) centers to receive ore from the cave. Each bell is accessed by two draw points, one from the production drift on either side, thus forming a 12.2 m (40 ft.) x 12.2 m (40 ft.) draw pattern. Extraction of the ore is by rubber-tired, 3.8 m3 (5 yd3) load-haul-dump equipment. The LHDs then tram the ore a maximum of 49 m (160 ft.) to ore passes. Cave initiation and bell development are done from the undercut drifts which are parallel to and 17 m (55 ft.) directly above the production drifts. Longhole rings are drilled and blasted from the undercut drifts to define the bells and establish the undercut for caving. (Figure 2) DRAW CONTROL Since the cave line at Henderson is constantly advancing, it is necessary to be continually initiating new cave at one end while exhausting it at the opposite end. There must exist, therefore, an angle on the ore-waste contact in the broken rock from initiation to exhaustion. The basic concept of draw control is to keep this angle as smooth and even as possible, particularly at the time of exhaustion. If this is achieved, draw points are exhausted more or less in a line, avoiding pockets of remaining ore surrounded by exhausted areas. These pockets would cause spotty ore extraction at the time of exhaustion, increasing the amount of dilution occurring while introducing the potential for significant weight problems in the production area. To arrive at the desired angle on the ore- waste contact, maximum tonnage percentages are assigned to each row of draw points increasing at 10% or 15% increments (depending on cave size and velocity of draw) working away from the cave line. The available tonnage indicated by these percentages is the maximum allowable tonnage to be extracted from each draw point until the available tonnage percent- age is increased. As the cave moves, these percentages increase for each draw point regularly. However, in general the tonnage drawn from each draw point is kept at about 50% of this allowable maximum in order to maintain adequate available tonnage in the cave to sustain production for seven months if cave initiation were to cease. This available tonnage cushion is a safeguard built into the draw control program at Henderson to accommodate fluctuations in the rate of cave advance. When draw points move past the row of 100% tonnage availability, they are drawn past the desired 50% at the same increments per row until exhausted. (Figure 3) To achieve proper draw control, the number of LHD buckets to be taken from each draw point is assigned daily. The actual buckets taken, which may at times deviate from the

Jan 1, 1981

By Lyall, Zickrick

As a result of the studies on recrystallization and crystal growth made in this laboratory, certain theories have been developed. These are expressed briefly in a paper by Dean and Hudson.' One of the principal points brought out in that paper is that in the deformation of a crystal, the energy supplied is distributed to the atoms of the lattice, probably by the forced formation of molecules. It is found that a Maxwellian distribution of the energy among these atoms will account for the grain growth at various temperatures. According to such a distribution the number of atoms with energy above a certain critical value, E, necessary for recrystallization at the absolute temperature T is given by n = Ne-E/kT where N is the total number of atoms. It is apparent that if this line of reasoning is correct, a metal which has been cold-worked and reheated to a temperature below that of complete annealing will have the atoms with the higher energy recrystallized, that is, reduced to normal energy and there will remain only atoms below a certain energy. If this metal is again cold-worked the peak in the number-energy distribution curve will, as compaied with its position after the first cold-working, be displaced toward the regions of lower energies, that is, the temperature of recrystallization will be higher. By a repetition of this process the position of this maximum and hence the recrystallization temperature of the metal could be raised considerably. In order to check this theory, experiments have been made on copper wire. The wire was drawn from 1/4-in. rod to 22 B & S gage and then sample lots were annealed at 400 and 220° C., respectively. The wire annealed at 400° C. gave an elongation of 34 per cent., whereas that annealed at 220° C. gave only 23 per cent. These wires were then drawn by identical practice to 34 B & S gage and the elongation determined after an anneal of 3 min. in an oil bath held at various temperatures. The figures obtained have been plotted in Fig. 1. It will be seen that the

Jan 1, 1927

By Robert H. Richards

The adaptation of European methods of concentrating ores to suit the conditions of this country has followed the lines that simplify machinery, diminish labor and increase capacity. Noteworthy instances are the substitution of hydraulic classifiers for the last one or two trommels in the series, and the substitution of a distributing-tank for the slimes in the place of the Spitzkasten. The question may be fairly asked, At what point do losses become serious in consequence of these changes? It is to study this question that the investigation herein described has been undertaken. After the losses are located and the methods of stopping them are pointed out, it will then lie with the individual mill to decide whether or not the plans proposed will save suficient value to warrant their adoption. In the discussion of sorting followed by sizing, as the terms are used in this paper, me have to deal with the preparation of slimes by the Spitzkasten for subsequent washing upon the slime-table. The Spitzkasten, or pointed box, is an inverted hollow pyramid, generally of wood, usually with four faces or sides coming together in a point below. The sides can have, according to Rittinger, a minimum angle of slope of 50" with the horizontal. The feed-water carrying fine slimes (less than 0.01 inch, or 0.25 mm. in diameter) for treatment is fed at one end and carefully distributed over the whole width. The outlet is over the opposite end. The water-current passes direct from inlet to outlet, and has usually been regarded as a surface-layer pass-

Jan 1, 1898

By William M. Shepard

The Tintic mining districts of central Utah com- prise one of the major silver-lead producing areas in the United States. Ore valued at nearly $450 million has been produced from these districts since the initial discoveries, in 1869, of veins cutting the Silver City monzonite stock. Subsequent production from the famous manto-type "ore runs", which characterize the Main Tintic district, was followed during the past 40 years by discovery of the blind orebodies of the East Tintic district, including the famous Tintic Standard mine and the most recent discovery, the Burgin mine, currently being brought to full production by Kennecott Copper. The occurrence of blind orebodies of exceptional size and grade in both the Main and East Tintic districts has made this area a field laboratory of major importance for the study of geochemical dispersion halos and alteration features associated with known ore, and has led to the discovery of additional concealed orebodies. The larger share of the basic data-gathering and interpretation has been done by members of the U.S. Geological Survey who have been studying the Tintic mining areas since before the turn of the century. These studies have been supplemented and expanded by the geological staffs of the local operating mining companies, by exploration geologists and by students.

Jan 4, 1966

By R. R. Goddard

By use of the frequency response method with a radioactive tracer, it was possible to determine fluid dispersion and distribution in a natural consolidated and an unconsolidated medium. Measurements were made in a linear flow system at oleic saturations of 69 per cent in the consolidated medium, and 100 per cent in both media. Dispersion and distribution were obtained by measuring the amplitude attenuation and the phase velocity of sinusoidal waves with a dual monitor apparatus. The gamma ray emissions permitted in situ measurements at any distance along the porous samples. One result of importance was that the effective diffusivity increased as the wave length increased. As a consequence, a dispersion coefficient appropriate for the injection of large slugs might exceed the value measured by use of small slugs. Since flow models based solely on fluid velocity and an effective diffusivity coefficient imply that the diffusivity should be independent of frequency, such representations were not adequate for the data of this study. A comparison was made with a capacitance model of porous media with deadend PVs, but even this model was not completely adequate. By using attenuation and phase velocity data, fluid dispersion can be predicted without postulating a differential equation satisfied by the tracer concentration, thereby eliminating the need of a complicated model to represent dispersion. INTRODUCTION The flow of similar miscible fluids through a porous medium can be fairly adequately described by two parameters: the average fluid velocity and the effective diffusivity.1-3 It has been pointed out recently, however, that significant discrepancies exist between this representation and the experimental data.4-7 An improved agreement can be obtained by introducing additional parameters based on the concept of dead-end pores. The purpose of the present investigation was to find out whether the frequency response method could be used to measure the relevant parameters. The method was used in the following form. A stream of fluid was flowed at a constant rate through a sample of porous material and the concentration of a radioactive tracer in the fluid was varied sinusoidally at a fixed frequency. The effects of flow through a porous medium are a decrease in the amplitude of the concentration wave and an increase in the velocity of the peaks of the waves above the average velocity. Attenuation and phase veIocity of the waves were measured as a function of frequency and fluid velocity. The simple two-parameter model implies that the diffusivity should be independent of frequency. Data reported in this paper show that the diffusivity decreases as the frequency increases. Hence, as shown also by many others, the two-parameter model is not completely adequate. Coats and Smith 5 used two additional parameters in their model: the volume of the dead-end pores and the rate of mass transfer between dead-end pores and the flowing stream. Their capacitance model of a porous medium containing some stagnant fluid, to which transfer occurs by molecular diffusion, did not explain the dispersion results of either the present study or of theirs. Instead, the capacitance effect can be better described as the result of extreme velocity variations within the pores of the medium, with transfer between the velocity zones by convection. EXPERIMENTAL PREPARATION OF POROUS MEDIA Data reported here were obtained with a natural Berea sandstone, and with two packings of 20/25 mesh sand. Pertinent information about these media is listed in Table 1. The sandstone core was first sealed with a non-

Jan 1, 1967

By M. H. Waxman, L. J. M. Smits

A simple physical model was used to develop an equation that relates the electrical conductivity of a water-saturatedshaly sand to the water conductivity and the cation-exchange capacity per unit pore volume of the rock. This equation fits both the experimental data of Hill and Milburn and data obtained recently on selected shaly sands with a wide range of cation-exchange capacities. This model was extended to cases where both oil and water are present in the shaly sand. This results in an additional expression, relating the resistivity ratio to water saturation, water conductivity and cation-exchange capacity per unit pore volume. The effect of shale content on the resistivity index - water saturation function is demonstrated by several numerical examples. INTRODUCTION A principal aim of well logging is to provide quantitative information concerning porosity and oil saturation of the permeable formations penetrated by the borehole. For clean sands, the relationships between measured physical quantities and porosity or saturation are well known. However, the presence of clay minerals greatly complicates log interpretation, particularly the electrical resistivity and SP logs, and considerably affects evaluation of hydrocarbon-bearing formations. The conductance and electrochemical behavior of shaly sands and their relation to log interpretation have been studied by many workers. wylliel and Lunch2 reviewed this work in some detail. Virtually all laboratory measurements of electrical resistivity and electrochemical potential of shaly sands published to date are the work of Hill and Milburn.3 Their measurements were made on about 300 cores covering a large variety of sedimentary rocks, and a wide range of equilibrating NaCl solution concentrations. Hill and Milburn described their conductivity data by an empirical equation in which the shaly sand conductivity Co was expressed as a function of the solution conductivity C and two parameters b and F01 . The quantity b was shown to be a measure of the effective clay content of the rock, being approximately proportional to the cation-exchange capacity of the rock divided by its pore volume. The latter ratio is designated as 9, in this paper and has the dimensions meq/ml or equiv/liter. Qv is identical with the term representing the concentration of fixed charges in the Meyer-Sieved and Teorel15 theory of permselective membrane behavior. F01 is a formation resistivity factor referred to a hypothetical equilibrating solution resistivity of 0.01 ohm m at 25C,* where clay effects presumably are minimized. F01 was correlated to porosity by an Archie-type equation.6 The Hill-Milburn equation describes their data with a standard deviation of approximately 1 percent and a maximum deviation of + 10 percent. Shaly sands behave as permselective cation-exchange membranes, their electrochemical efficiencies increasing with increasing clay content. The electrochemical potential data of Hill and Milburn were expressed graphically, and demonstrate that the membrane efficiency (or cation transport number) of these sands is a function only of the b value (i.e., Qv) and the respective salt concentrations of the two solutions forming the liquid junction in the sand. The diffusion potentials are not dependent on F01 or any parameter relating to the porosity or pore geometry of the rock. The Hill-Milburn resistivity equation correctly predicts a decreasing sand conductivity Co with decreasing solution conductivity C?. However, at some low value of C,, the calculated Co-C, function goes through a minimum; with further decrease in C?, the predicted sand conductivity increases sharply. As pointed out by Hill and Milburn, an increasing value of Co with decreasing C, is physically meaningless. Since this occurs below the range of practical values of C,, the usefulness of the empirical equation is not affected when

Jan 1, 1969

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

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

Jan 1, 1952

By F. F. Skiba, J. E. Warren

An idealized miscible displacement process in a three-dimensional, heterogeneous porous medium has been studied via experimental computation based on a "Monte Carlo"' model. The macroscopic dispersion which results solely from variations in the permeability of the porous system is related to the scale of the heterogeneity as well as the distribution function for the permeabilities. The interpretation of laboratory data is often subjective and may bt? misleading since ambiguities are not always recognized. Furthermore, an experiment performed on a conventional oil-field core does not yield a valid measure of Macroscopic dispersion for reservoir engineering purposes since the scale of heterogeneity that is significant in the laboratory is not significant in the reservoir. The analysis of field-performance data is discussed and it is shown that a unique characterization of the reservoir can not be obtained from this information alone. A possible method based on a minimum of data, for estimating the spatial distribution of permeabilities in a reservoir is postulated INTRODUCTION Many studies concerned with diffusion and dispersion in porous media have been undertaken because of the potential importance of miscible displacement processes for increasing the recovery of oil. The results from a large number of these studies have been published; in a recent paper Perkins and Johnstonl presented a comprehensive survey of the pertinent literature and excellent appraisal of the "state of the art". In this investigation, dispersion resulting from macroscopic variations in the properties of the porous medium is of primary interest. Dispersion, in this sense, is caused by fluctuations in the velocities of the individual fluid elements as they move through the porous system. This phenomenon is the result of inhomogeneities in the permeability and/or porosity, but it can be amplified by differences in the characteristics of the fluids involved. The objectives of this study are the following: (1) to determine the qualitative manner in which macroscopic rock properties affect the observed dispersion coefficient; (2) to evaluate the possible influence of macroscopic dispersion on laboratory experiments; and (3) to appraise the significance of macroscopic dispersion with regard to reservoir performance. Investigations of this type are essential since many miscible displacement processes, while seemingly understood from a mechanistic point of view, do not behave in a predictable manner under field conditions. Although a plausible explanation is that the process is dominated by the environment, little serious effort has been devoted either to describing the reservoir itself or to determining the probable effect of the reservoir properties on the process. This study is intended to be exploratory in nature; however, it is hoped that it will provide additional insight into the physical problem and that it will suggest areas suitable for further investigation. EXPERIMENTAL COMPUTATION MODEL To limit the observed dispersion to that which results solely from variations in the properties of the porous medium, an idealized process in which one incompressible fluid displaces another identical fluid is simulated. This simulation is achieved by means of a modified Monte Carlo method — by tracking a number of mathematical particles (tracer) through the system according to the steady-state, single-phase velocity distribution. Then the system is characterized by the distribution of the residence times of the particles. Since any number of particles can be injected simultaneously (a true delta-function distribution) and only pressure conditions are imposed on the boundaries, the model is completely determinate. Furthermore molecular diffusion, microscopic dispersion, adsorption and/or reaction are automatically precluded from the model. The heterogeneous medium is represented by a rectangular, three-dimensional array of uniform, homogeneous porous blocks which are randomly arranged; it is assumed that any two opposing faces of the system are held at fixed potentials

Jan 1, 1965

By A. Gemperle

Electron diffraction and transmission electron microscopy on foils at room temperature were used to investigate the ordering of Fe-Si alloys containing 10 to 23 at. pct Si. A certain degree of DO3 order was found in all alloys. With the exception of the lowest silicon concentration for which the antiphase domains could not be clearly resolved, the alloys have a domain structure of two-domain type with boundaries having 1/4a01<111> displacement vectors for less than 12.3 at. pct Si and with boundaries having 1/2 a0<100> displacement vectors for more than 12.3 at. pct Si. The alloys with 12.3 at. pct Si have a domain structure consisting of fine domains with 1/2a&apos;o<100> boundaries within much larger domains with 1/4a&apos;o<111> boundaries. The development of these structures can be explained by transition of the alloy from the disordered state into the B2-type order and then into the D03-type order by the mechanism proposed previously for the FeSAl alloys. The existence of the B2 structure in the lower part of the investigated concentration range reported in some articles can be explained by fine domains with 1/2a&apos;o<100> boundaries formed by several disordered planes within large domains with 1/4a&apos;o<111> boundaries. The ordered structure predicted by the theory —with practically no domain boundaries —is found in the alloys having 12.3 at. pct Si where it develops in the B2 structure region. ORDERING in Fe-Si alloys was first studied by phragmenl who found that beginning with 13 at. pct Si the DO3 (Fe3Si) superlattice reflections appear in the diffraction patterns. The equilibrium diagrams constructed later by corson2 and Haughton3 from various measurements proposed the existence of a homogeneous solid solution (a phase) in the range from 0 to 25 at. pct Si. Jette and Greiner4 and Farquhar et al.5 measured the relation between lattice parameter and composition and they considered the break in the curve at 9 to 10 at. pct Si to be caused by the ordered solution a"(Fe3Si). Glaser and Ivanick6 Determined critical ordering temperatures of the alloys containing from 10.9 to 27.9 at. pct Si from the measurement of the electric resistivity of quenched samples. In all cases the critical temperature was lower than the melting point and it was highest for 25 at. pct Si. Lihl and Ebel7 measured the lattice parameter curves at various temperatures up to 1000°C. The region between two breaks on these curves, corresponding to 10 to 12.5 at. pct Si at room temperature, was considered by them to be two-phase (a + a"). They concluded by extrapolation of the measured values that a" in the alloy having 25 at. pct Si is stable up to the melting point. Davies8 studied superlattice reflections in the X-ray diffraction patterns of an alloy containing 8.7 at. pct Si. He found the B2 structure and short-range order in the slowly cooled samples and the DO3 A. GEMPERLE is Research Scientist, Institute of PhysicS, Czechoslovak Academy of Sciences, Prague, Czechoslovakia. __Manuscript submitted January 2, 1968. IMD structure in the quenched and annealed samples. This investigation first reports the presence of the B2 structure, phase a&apos;: in Fe-Si alloys. Meinhardt and krisement9,10 also found its existence in Fe-Si alloys by neutron diffraction. No order was detected by them in the alloy containing 8 at. pct Si. The onset of B2 order was observed at a composition of 9.2 at. pct Si. They found almost perfect B2-type order with partial DO3-type order at room temperature in the 10 to 12.5 at. pct Si range and almost perfect DO3-type order in the 12.5 to 25 at. pct Si range. They established the critical temperatures Tc of both the structures through measurement at higher temperatures. The critical temperature for the B2 structure was found to be always higher than the critical temperature for the DO3 structure of the same alloy. They extrapolated the curves of the critical temperatures and concluded that the alloys with more than 17 at. pct Si have the B2 structure up to the melting point and the alloys with more than 23 at. pct Si have the DO3 structure up to the melting point. The results of Meinhardt and krisement9,10 were confirmed by Dokken&apos;s measurement" of the temperature dependence of the electrical resistivity in the 10.8 to 15 at. pct Si range. On the other hand chessin12 detected by X-ray diffraction a considerably lower degree of order in the 12.7 at. pct Si alloy. ANTIPHASE DOMAIN STRUCTURE IN ALLOYS WITH B2- AND DO$-TYPE ORDER The ordered structures B2 and DO3 can be described in terms of a subdivision of the bcc lattice into four fcc sublattices with a parameter double that of the bcc lattice. Following Marcinkowski13 we will label them I. 11, 111, IV. The B2 structure in the AB alloy is formed by placing A atoms on sublattices I and II and B atoms on sublattices III and IV. In a non-stoichiometric perfectly ordered alloy having concentration (A) > (B), A atoms occupy sublattices I and 11. and A and B atoms distributed at random occupy sublattices III and IV. The DO3 structure in an A3B alloy is formed by placing A atoms on sublattices I, 11, and 111, and B atoms on sublattice IV. In a nonstoichiometric, perfectly ordered alloy having concentration (A)/3 > (B), A atoms occupy sublattices I, 11, and 111, and A and B atoms distributed at random occupy sublattice IV. As further shown in Ref. 13, two types of domains are possible in the B2 superlattice and the antiphase domain structure has associated with it boundaries with displacement vectors 1/4 a&apos;o<ll1> only. Four types of domains are possible in the DO3 superlattice and the antiphase domain structure has associated with it boundaries with displacement vectors 1/4a&apos;o<111> and 1/2a&apos;o<l00>. Bethe14 suggested on the basis of theoretical considerations that in a structure with two sublattices at low temperatures only one domain should be present in the whole crystal at equilibrium. Similarly Bragg15 concluded that at low temperatures the domain struc-

Jan 1, 1969

By J. E. Sandor, P. F. Schmidt

A general review of the present knowledge of the semiconductor/oxide interface is given, with emphasis on silicon and germanium. The electrical properties of clean surfaces and the changes that occur upon oxidation are described. A brief disczcssion of the methods and mechanisms of oxidation of silicon is given, and the interaction between the semiconductor substrate and the oxide phase is considered. Physical properties of the oxide film, charge storage effects due to the dielectric, and some device applications of oxide films on semiconductors are discussed. THIS paper will be heavily preoccupied with the elemental semiconductors germanium and silicon. Undoubtedly some injustice is done thereby to other semiconductor-oxide systems, but it is not readily possible to cover all aspects of so vast a field in the present review. In considering the semiconductor/oxide interface, we may start out by discussing briefly an idealized clean semiconductor surface, and may then follow the changes brought about by oxidation, as the oxide is allowed to grow progressively thicker. The references given are not a complete list, but are to be taken as guides to the literature. Tamml suggested in 1933, based upon the study of a one-dimensional asymmetrical potential model, that localized surface states may appear at the surface of a crystalline lattice; it was pointed out by Fowler that these states should occur in pairs. shockley3 examined the origin of such surface states on the basis of a more general one-dimensional model and found that surface states may appear in the forbidden region between the highest filled and the lowest vacant band. This happens when the two bands overlap, as, for instance, the "s" and (&apos;p" bands of the diamond structure. In a three-dimensional model these levels would produce a surface band which is half-filled for diamond, and the number of surface states should be approximately equal to the number of surface atoms. This picture of the surface has been refined by more recent theoretical work, like Handler&apos;s4 and Koutecky and Tomasek&apos;s~ without introducing essential modifications. When these concepts were applied to semiconduc-

Jan 1, 1965

By John W. Cahn

Tile avevage value of the mean curvature of surfaces in a specimen can be precisely delermined by sitrlple measurements performed on random sections or on 1 vojectiotzs of these surfaces. For surjaces enclosing particles, the quanlity is sllown to be velated to the probable tlurrrhev ot particles intevcepted by a random platne and, under certain conditions, to the perimeter of their Projected image and to the increuse in surjace area when the particles grow. Many of the velalionsl~ips are valid for any shape or size of the particles. THE recent growth of the field of stereology or quantitative microscopy is largely the result of the discovery of a limited number of rigorously defined geometrical parameters that can be precisely measured on random sections or projections of a structure without assumption about particle shape.&apos; Such parameters have been the volume fraction and surface area and line length in a unit volume. The present paper introduces the average of the mean curvature of surfaces in a volume, and shows that it can be rigorously measured by counting methods alone. AVERAGE CURVATURE ON A PLANE Consider a curve y = v (x) on the xy plane. The curvature Kp in the plane at any point is given by Kp = (dZy/dx2)/[1 + (dy/dx)2]3/2. The average curvature along the curve may be defined as Kp = (1/L) fKpdl |1] where dl = [(dx)&apos; + {dy)2]1/2 is the length of an element of the curve and .(dl = L is the total length of the curve. The integral on the right-hand side of Eq. [l] will be called the total curvature, It is dimensionless and readily evaluated by substituting d = arctan (dy/dx) to give ah which is the change in direction of the tangents to the ends of the curve, Fig. 1. Hence, Kp = (1/L)(?h) [2] It can also be seen that Eq. [2] is applicable in the limit of sharp corners. Even though the curvature is infinite at a corner, f Kpdl remains integrable and equal to the change of angle. The signs of the quantities Kp and dl are arbitrary and, if the areas adjoining the curve are equivalent, there seems to be no compelling reason to fix the signs. For a closed curve, however, it becomes worthwhile to use the signs to distinguish among different types of areas or phases enclosed. For such a closed curve AO = t2r and As a convention, we shall take the direction of the integration along the curve in such a way that the area or phase under consideration (henceforth called a) is kept to the left side of the curve, Fig. 1. Then Kp is positive where the curve is concave and negative where the curve is convex. For a closed curve Ad - +2n if the curve encloses a, and Ad = —2n if the curve is surrounded by a. PLANAR SECTIONS OF THREE-DIMENSIONAL STRUCTURE OF ISOLATED a PARTICLES Consider next a planar section of a microstructure containing particles of a in a matrix. We place no restriction on shape of connectedness of these particles, and begin by considering the curvature of only the boundary between a particles and matrix&apos; Clumps of contiguous a particles are then treated as a single a particle. The two-dimensional section will show in a unit area na simple closed curves enclosing a particles and nh closed curves of holes, in-

Jan 1, 1968

By H. H. Rachford

Numerical solutions of immiscible flow problems in which dispersive effects of capillarity are dominated by convection require excessively fine grid spacing with attendant high computing costs. The use of coarser spacing reduces cost but often produces oscillation or undue dispersion associated with displacement fronts, A numerical formulation is proposed here which should be applicable to two-dimensional flow problems. it is in part analogous to an approach previously tested for miscible systems. The convective transport is approximated using a change of variables to yield a coordinate system moving approximately with the local characteristic velocity. The capillarity-induced dispersive terms in the differential system describing the process are approximated with respect to a fixed coordinate system by the usual implicit formulation. One-dimensional tests of the procedure yielded results in which the saturation profiles tended smoothly to the zero-capillary pressure solution as the ratio of viscous to capillary forces was successively increased in a sequence of calculations. This contrasted favorably with solutions by other numerical procedures which would require attendant grid refinements to approach the zero capillary pressure results. INTRODUCTION Numerical solution of displacement problems has until recently relied on applying methods developed primarily for transient heat-flow problems. Such problems are classified as parabolic in type, and where the heat transport is purely by diffusion their solutions are characterized by a high degree of smoothness. It is not surprising, therefore, that for approximating these solutions available finite difference methods are quite adequate. In flow problems the transport is partly by diffusion, partly by convection or flow. Although the problem remains of parabolic type because the dispersive effects of capillary forces or diffusion play some role in every displacement, at high flow rates the problem is dominated by convection, and solutions tend toward those of equations of the hyperbolic type. Solutions of hyperbolic problems are characterized by the translation of fronts, or discontinuities, that may progressively increase in sharpness. Numerical methods for treating parabolic problems become less and less satisfactory as displacement rates increase and the role of dispersion due to concentration or capillary pressure gradients becomes small relative to transport due to flow. In computation the difficulty manifests itself as an error associated with the grid size chosen. 1-6 In summary, if the heat-flow type approximations are to include the terms arising due to convection, one of several choices may be made: (1) an upstream (to the direction of flow) approximation for the convection terms may be used; (2) a centered-in-distance (CID) approximation may be used; or (3) a recently developed approximation based on the theory of oscillation matrices may be chosen.6 The last appears to have significant promise for one-dimensional flow problems; its extendibility to two or three dimensions is an open question. In either of the first two approaches, a suitably small ratio of v&/D must be maintained, where v is the velocity, & is the grid spacing and D the effective dispersivity in the direction of flow. In the first choice, the approximation of the convective part is only first-order correct and errors introduced appear as a numerically induced dispersivity of magnitude proportional to v?x. In the CID choice, the approximation can be second-order correct, but the difference formulation fails to satisfy the maximum principle unless a condition on v?x/D is met. Practically, this means that for high flow rates oscillatory solutions may result in the neighborhood of a front unless exceedingly small grid intervals are taken. While the procedure proposed by Stone and Brian4 permits a less severe limitation to be placed on this ratio, ultimately the flow rates increase relative to the dispersivity the oscillation obtains. Further, extensions of their approach to higher dimensional systems may be attended by considerable

Jan 1, 1967

By V. N. Antaki

ELEMENTAL phosphorus (referred to as P,) is a true mineral product that has developed, in the relatively brief span of a century, from a laboratory curiosity to a heavy chemical. Annual production now exceeds 300,000 short tons. Continued rise in production has been predicted, consistent with new applications for P, derivatives and the expected population increase. In itself elemental P, is valuable chiefly as an &apos;intermediate chemical. The principal end uses at present are in inorganic phosphates, both simple and complex, of sodium and potassium, used in water treatment and detergent compounds, leavening agents, food supplements, and dentifrices. Organic P, compounds are important as plasticizers, oil additives, flotation agents, and more recently as gasoline additives (TCP being the outstanding example). The increasing use of high-analysis water-soluble and liquid solution fertilizers, as well as the direct application of phosphoric acid to the soil, will probably be the most significant factors in P, production growth rate. Coke oven plants now using sulfuric acid to fix byproduct ammonia values will gradually switch to phosphoric acid in order to market a fertilizer with completely available values. Recently sponsored experiments show that phosphoric acid can be directly admixed with cattle and poultry feed to produce rapid growth. Both P, compounds and modifiers for oil well drilling mud show increasingly promising significance.

Jan 3, 1957

(The following-list contains the names of those members of the Institute of whose connection with military service we have only recently. become acquainted.) ADAIR, ARTHUR C., Corporal, 4th Co., 1st Tr. Bn., 154th D. B., Camp Meade, Md. BATCHELOR, STILLMAN, 1st Lieut., Co. 6, E. 0. T. S., Camp Humphreys, Va. BATES, BENNETT R., 1st Lieut., Engineers, Co. 3, E. 0. T. S., Camp Humphreys, Va. BJORGE, GUY N., Co. 3, E. 0. T. S., Camp Humphreys, Va. BOAS, Ross, H., U. S. A. BRADLEY, . H., JR., Capt., 408th Engineers, Camp Cody, N. M. BRINTON, OWEN F., Capt., Engineers, U. S. A. BROWN, HAROLD C., 2d Lieut., A. S. A.P.

Jan 1, 1919

By C. Dinis da Gama

In order to reduce the complexity of mechanisms influencing rock. fragmentation by blasting a simulation approach is proposed, using the capabilities of microcomputer interactive graphics. Situations involving cratering end bench blasting are simulated, and both jointed and intact rocks can be subjected to explosive action, taking into consideration its breakage mechanisms. Size analysis curves relating rock distribution before and after the blasts are determined and compared with field data.

Jan 1, 1984