Search Documents

By Charles Malcolmson

AN elaborate series of experiments has been conducted during the past three years at Irvington, N. J., which has resulted in the perfection of a process for the manufacture of smokeless fuel from high-volatile coals, and for the recovery and refinement of the coal-tar products derived therefrom. These experiments have been financed by Messrs. Blair Co., of New York, and were conducted under the direction of Charles H. Smith, a member of this Institute, and the inventor of the process. The low-temperature distillation of coal has interested investigators for many years. Sporadic attempts have been made to solve the mechanical problems, but until the Smith process. was developed, they were not carried to conclusions of economic value. The present coal shortage and the increasing demand for smokeless fuels make this subject one of timely interest.

Jan 5, 1918

By S. W. Churchill

A digital computer was used to obtain an exact numerical solution for the transient behavior of the insulation and earth adjacent to an isothermal, submerged flat surface for a single set of parametric values. Comparison of the computed results with analytical solutions for limiting conditions revealed that a complete and general solution for all parametric values could be constructed from these limiting solutions. Complete and general solutions for insulated spheres and cylinders can also be constructed from limiting solutions. The procedure is illustrated in part for an insulated spherical tank. The heat flow was found to fall rapidly to a pseudo-steady state; after some time, it then decreased slowly to zero for a flat plate and cylinder, and to a low steady-state rate for a sphere. The accumulative heat flow during the initial falling-rate period may be a significant fraction of the heat flow during the entire first year. INTRODUCTION Underground storage of liquefied natural gas and liquefied petroleum gas has received considerable recent attention. The rate of heat flow from the earth to the storage cavity and the resulting temperature field in the earth are important factors in the technical and economic evaluation of such storage facilities. The objective of this paper is to indicate how complete solutions can be developed for the transient flux and temperature field for various geometrical configurations. The representative properties and dimensions, and the resulting parameters utilized in the illustrative results, are indicated in Table 1. The results presented herein are for dry earth. The effect of the latent heat of solidification of water in the soiI will be described in a subsequent paper. The temperature field in the insulation and earth is determined by energy balances, boundary conditions and initial conditions. The physical problem can be described mathematically with virtually no idealizations insofar as physical properties are known. It appears possible to solve the equations by numerical integration on a high-speed digital computer for any geometrical configuration and conditions. For complex situations, however, the computations are expensive and the results are highly specific. Analytical solutions have been developed for a few simple but important geometrical configurations and conditions, including one-dimensional heat transfer from (or to) earth at an initially uniform temperature to (or from) flat plates, spheres and circular cylinders. The solution for an insulated flat plate has also been derived but is in the form of a slowly converging series involving taklated functions. It was planned to use a computing machine to evaluate this series for a number of representative conditions. However, upon examination of the results of preliminary computations it was discovered that a complete, general and accurate solution could be developed by interpolation between several much simpler solutions for limiting conditions. This technique then was used to develop a a complete solution for an insulated sphere. The equations presented herein are derived in Carslaw and Jaegerl and other books on applied mathematics, or they are simple extensions of these previous results. Hence, all derivations are omitted. Emphasis is instead placed on the results and their

By R. Simon, L. Yarborough

This paper presents a correlation for predicting the critical pressures of gas-solvent-reservoir oil systems of the type encountered in gas and solvent injection processes. Experimental data from 14 systems were used to develop the correlation which expresses critical pressure as a function of composition, C7+ mol weight, and C7+ characterization factor. Critical pressures predicted by the correlation for the 14 gas-solvent-reservoir oil systems have an average deviation of 2.8 per cent of the experimental values. The correlation war used to study a wide variety of other hydrocarbon mixtures reported in the literature whose critical conditions are in the range of interest in oil recovery work (100 to 280F and 2,000 to 6,000 psia). Critical pressures predicted by the correlation had an average deviation, expressed as per cent of experimental value, of 5 per cent for 20 miscellaneous mixtures from the literature, and 6.4 per cent for I7 methane-butane-decane mixtures reported by Sage and Lucey. INTRODUCTION Predicting critical points of hydrocarbon mixtures is an important part of phase behavior calculations. If the critical point can be predicted for a specific mixture, the separation between the bubble-point and dew-point regions will be defined, and physical properties of the mixture can be calculated by using the principle of corresponding states. Numerous correlations for predicting critical points of hydrocarbon mixtures have been published. Six of these correlationsl-8 were for light hydrocarbons, and two correlations7,8 ere for high boiling-point petroleum distillates. Organick9 And Eilerts" presented correlations which combined the published data from both light and heavy mixtures, and these correlations became the only available methods for predicting critical pressures for wide boiling range mixtures. In the course of the authors' work on secondary recovery by gas and solvent injection, experimental data were obtained for 14 mixtures of dry gas, ethane-propane-butane solvents, and various reservoir oils. The critical pressures and critical temperatures of these gas-solvent-reservoir oil systems were in the range of 1,900 to 4,800 psia and 180 to 265F. The critical pressures predicted for the 14 mixtures by the correlations of Organick and Eilerts had an average deviation from the measured values of about 13 per cent. The authors felt that it was necessary to develop a better correlation for engineering calculations and the resultant correlation, which is presented on Figs. 1 and 2, predicts the critical pressures of the 14 mixtures with an average deviation of 2.8 per cent. The validity of the correlation was tested by applying the correlation to other hydrocarbon mixtures reported in the literature. This subject is discussed later in the paper. The authors feel that there is also a need for an improved critical temperature correlation for gas-solvent-reservoir oil systems. However, developing a critical temperature correlation is a significant study in itself and is not included in this paper. EXPERIMENTAL DATA The compositions and critical points obtained by the authors for 14 gas-solvent-reservoir oil mixtures are shown on Table 1. The methane concentrations in these mixtures varied from 20 to 63 mol per cent. The intermediate hydrocarbon group, consisting of ethane, propane, and butane, varied from 20 to 64 mol per cent, and the concentration of the C,, fraction from 10 to 20 mol per cent. The physical properties of the C,, fraction covered the range from light, paraffinic, nonasphaltic (Mixtures 5, 6, 7, 8, 9 and 10) to heavy, aromatic, and asphaltic (Mixtures 13 and 14).

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

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

Jan 1, 1970

By Edgar R. Kendall

A knowledge of the fundamental principles of locating, sampling, testing and evaluating materials is essential in prospecting for mineral aggregates. In this discussion, mineral aggregates will be considered sand, gravel, crushed stone and slag that are used as materials of construction. The quantity of mineral aggregates produced in the United States during the year 1935 for such purposes was roughly 165,000,000 tons, about 60 per cent of which was sand and gravel, 37 per cent crushed stone and 3 per cent crushed slag. The total value of this material was approximately $100,000,000. It is apparent therefore that the prospecting for suitable mineral aggregates is of considerable economic importance. In prospecting for mineral aggregates, the methods employed are the same whether the organization is large or small. In the large organization there may be available geologists for locating sources of supply, a land acquisition department to arrange for the necessary options or permits, engineers to do the sampling, a laboratory for testing, and a board of consulting materials engineers to make the final selection. Even small organizations must follow the same general procedure, though the responsibility for the selection of aggregates may rest with one man, who may or may not be an engineer, but who must depend entirely upon his own judgment. Locating the Supply The initial step in prospecting for mineral aggregates is the locating of suitable, adequate and economic sources of supply for the project in mind. Commercial plants in the vicinity should be investigated and a survey made of the existing structures that have been built of aggregates from these sources. If the structures are rather old and show no signs of deterioration, it is fairly good evidence that the aggregate is satisfactory. If materials of the proper quality can be obtained from commercial sources in sufficient quantity, and delivered to the job economically, it

Jan 1, 1942

By Edgar R. Kendall

A knowledge of the fundamental principles of locating, sampling, testing and evaluating materials is essential in prospecting for mineral aggregates. In this discussion, mineral aggregates will be considered sand, gravel, crushed stone and slag that are used as materials of construction. The quantity of mineral aggregates produced in the United States during the year 1935 for such purposes was roughly 165,000,000 tons, about 60 per cent of which was sand and gravel, 37 per cent crushed stone and 3 per cent crushed slag. The total value of this material was approximately $100,000,000. It is apparent therefore that the prospecting for suitable mineral aggregates is of considerable economic importance. In prospecting for mineral aggregates, the methods employed are the same whether the organization is large or small. In the large organization there may be available geologists for locating sources of supply, a land acquisition department to arrange for the necessary options or permits, engineers to do the sampling, a laboratory for testing, and a board of consulting materials engineers to make the final selection. Even small organizations must follow the same general procedure, though the responsibility for the selection of aggregates may rest with one man, who may or may not be an engineer, but who must depend entirely upon his own judgment. Locating the Supply The initial step in prospecting for mineral aggregates is the locating of suitable, adequate and economic sources of supply for the project in mind. Commercial plants in the vicinity should be investigated and a survey made of the existing structures that have been built of aggregates from these sources. If the structures are rather old and show no signs of deterioration, it is fairly good evidence that the aggregate is satisfactory. If materials of the proper quality can be obtained from commercial sources in sufficient quantity, and delivered to the job economically, it

Jan 1, 1942

By Stephen L. Bolivar, Susan H. Balog, Thomas A. Weaver

A data-classification scheme was developed to detect potential uranium mineralization in the Montrose l º x 2º quadrangle, Colorado. The methodology developed is a rapid and efficient method of resource evaluation on a reconnaissance scale. The necessary techniques were developed and refined to digitize, integrate, and register various large geological, geochemical, and geophysical data sets for the Montrose quadrangle, Colorado, using a grid resolution of 1 km. All data sets for the Montrose quadrangle were registered to the Universal Trans- verse Mercator projection. The data sets include hydrogeochemical and stream sediment analyses for 23 elements, uranium-to-thorium ratios, airborne geo- physical survey data, the locations of 90 uranium occurrences, and a geologic map (scale 1:250 000). Geochemical samples were collected from 3965 locations in the 19 200 km2 quadrangle; aerial data were collected on flight lines flown with 3-5 km spacings. These data sets were smoothed by universal kriging and interpolated to a 179 x 119 rectangular grid (each grid block is 1 km2). A mylar transparency of the geologic map was prepared and digitized. All possible combinations of three, for all data sets, were examined for general geologic correlations by utilizing a color microfilm output. Subsets of data were further examined for selected test areas. A classification scheme for uranium mineralization, based on selected test areas in the Cochetopa uranium district, is presented. Areas favorable for uranium mineralization, based on this scheme, were identified and are discussed.

Jan 1, 1982

By C. H. Ault, D. D. Carr

Four counties in Indiana- LaPorte, Vigo, Vanderburgh, and Switzerland-have potential for deep underground limestone mining because they are in or near high population areas and have no surface-mineable limestone deposits. The principal markets are Gary-South Bend, Terre Haute, Evansville, and Cincinnati. Mineable seams suitable for aggregate are 120-150 m (400-500 ft) below the surface in LaPorte and Vigo Counties, about 90 m (300 ft) below the surface in Switzerland County, and about 600 m (2,000 ft) below the surface in Vanderburgh County. Some deposits are also suitable for high-purity uses or for cement raw materials, which enhances their potential for exploitation.

Jan 1, 1983

By E. O. Fromm, R. I. Jaffee, R. M. Evans

IN connection with a research on gallium alloys, a method for the thermal analysis of very low melting alloys was developed, which appeared to be worthy of more than passing interest. As used, the method was suitable for work in the temperature range, — 72° to 250°C. However, with suitable adjustment, the method could be adapted for use over broader temperature ranges. The Rosenhain temperature gradient furnace,&apos; adapted for low-temperature use, was the basis for the method. The gradient was set up in a stainless steel tube. The specimens were lowered through it at a constant rate while time-temperature or inverse rate data were taken. It was essential to maintain careful control of the top and bottom temperatures of the gradient. Also essential was good heat transfer from the walls to the specimen. This was accomplished by a liquid heat-transfer medium.

Jan 1, 1953

By Walter C. Labys

INTRODUCTION Interest in multicommodity stabilization schemes has increased recently for several reasons, a most important one being the experienced wide swings in primary commodity prices. These prices rose sharply between mid - 1972 and end - 1974 giving commodity exporting countries an opportunity to gain sizeable increases in export earnings. But since then, prices have tumbled and these countries are searching for a means to stabilize if not to increase their earnings. While this problem is particularly serious for developing countries, the developed countries are concerned about the impending upward price swing, which would increase inflation before recovery from the world recession is assured. As a direct consequence of the emergency meeting of the U.N. General Assembly devoted to commodity problems in April of 1974, interest has been renewed in the Keynesian (4) plan for an international multicommodity buffer stock program. In particular, UNCTAD (8) is designing an integrated commodity scheme to be discussed within that body. In terms of approaches that have appeared in designing such a program, none exist which relate to the interpretation of the independence and interdependence of prices of the commodities of interest. Typical criteria for selecting commodities for the scheme have included their importance in world trade, their particular origin from developing countries, and their feasibility for storage. But the financial management of a multicommodity buffer stock requires a consideration of commodity prices, yields, and risks as does as any portfolio of assets. The approach which we take in this study is that of application

Jan 1, 1977

By E. I. Organick, B. H. Golding

A simple correlation is presented for the prediction of saturation pressures in condensate-gas and volatile-oil mixtures. Saturation pressure is related directly to the composition of the mixture with the aid of two generalized composition characteristics, B, the molal average boiling point, and W, the modified weight average equivalent molecular weight. The need for hydrocarbon equilibrium constants is entirely eliminated. Only the values of the two conlposition parameters need be obtained, and these are calculated in a simple and reproducible manner from the stream analysis of the mixture and the ASTM distillation for the heavy ends. The correlation is given .in the form of a set of 14 working charts in which saturation pressures, either retrograde dew points or bubble points, are plotted cs temperature, forming partial phase envelopes for mixtures having discrete values of the composition parameters, B and Wul. A locus of critical states on each chart intersects the partial phase envelopes, distinguishing retrograde dew point pressures from bubble pressures. A comparison Between experimental and predicted saturation pressures reveals that: the probable error of the predicted values, when considering the data as a whole (214 experimental points), is five per cent; for high gasjoil ratio condensate fluids (approximately 40.000 or more cu ft/bbl) the probable error is about 8 per cent. In either case, however, the probable error is comparable with the overall uncertainty of the data from which the correlation is derived. INTRODUCTION There is a present need for a simplified, yet reasonably accurate method for predicting upper phase boundaries of reservoir fluids. In the case of high gas/oil ratio fluids, namely the so-called highly volatile oils and the condensate gases, no simple method, accurate within 5 to 10 per cent, is now available by which bubble point and retrograde dew point pressures may be predicted. In order to calculate the upper phase boundary of a volatile reservoir fluid, using correlations now available, it is necessary to find and employ the proper K-values (equilibrium volatility ratios of the hydrocarbons in that fluid). This process of finding the Proper K-values usually involves a double trial and error computation. Since the K-values are in themselves dependent upon the composition of the fluid, their use in the calculation of a saturation pressure actually represents an unnecessary complication. The correlation described in this paper relates saturation pressure directly to the stream composition, entirely eliminating the need for K-values. The only data required to calculate the Stream composition are obtained from the customary low temperature fractional analyses of separator gas and liquidt and from the ASTM distillation of the c,+ fraction of the separator liquid. A considerable simplification in the procedure for predicting saturation pressures is thereby effected without significant loss in accuracy. BASIS FOR A CORRELATION Empirical relationships between the saturation pressure and the composition of the fluid in question have already been partially explored by one of the authors.V n particular, a relationship between the critical pressure for a hydrocarbon mixture and its composition was developed. The critical pressure is in reality a saturation pressure, as it repregents a limiting case either of a dew point or of a buhhle point pressure for any given mixture. It was found that the critical pressure may be correlated as a function of two generalized composition characteristics or variables. One characteristic expressed a molal average prop-erty, the other a weight average property of the mixture. Specifically, the molal average boiling point (MABP, or B as it

Jan 1, 1952

By A. G. Weber, H. J. Welge

A published calculation method for predicting incompressible, multidimensional fluid displacement has been adapted to the problems of water and gas coning in oil wells. Since depth and radial distance from the wellbore are the two key dimensions affecting the shape of a gas or water cone, coning calculations are well suited to the use of two-dimensional methods. The alternating direction implicit procedure (ADIP) was used for relaxation calculations of two-phase potentials in a two-dimensional grid. From the potentials and the capillary pressure relationship, saturation and pressure distributions were calculated which trace cone growth with time. Predictions have been made for well producing histories both before and after core breakthrough. To check validity of the method, two-dimensional calculations have matched the coning behavior and produced water-oil ratio history of a laboratory sand-packed model. They have also matched the coning behavior of several producing wells, for which the calculations were compared with produced water or gas cuts and logs showing water or gas cone movements. INTRODUCTION The application of two-phase, two-dimensional calculations using ADIP to various reservoir flow problems has been described in the literature.1-3 When adapted for computer solution, this method has proved to be a powerful tool for simulating well and reservoir behavior. This paper discusses the method as applied to well coning calculations. Single-well and coning calculations comprise an especially difficult class of two-dimensional problems which require special techniques for computer calculation and determination of reservoir characteristics. Refs. 4 through 8 describe previous approaches to the coning problem. Several examples of water and gas coning calculations, including studies on both laboratory models and producing wells, are presented here. The two-dimensional method accounts realistically for the most critical parameters affecting coning behavior, including production rate, formation stratification, horizontal and vertical permeabilities, depth of well penetration, gravity and capillary forces. The method considers the different densities and viscosities of the two phases and the relative permeability and capillary pressure characteristics of the rock and fluids. In addition to tracing cone growth in the vicinity of the wellbore, the method calculates the overall movement of the fluid interface throughout the well&apos;s drainage volume. Incompressible fluid flow is assumed to occur between the producing interval and the well&apos;s limit of drainage. Calculations can be made for the producing history both before and after cone breakthrough. A typical two-dimensional grid or array of blocks used to solve a coning problem contains about 400 blocks. The well&apos;s cylindrical drainage volume can be represented by about 20 radial subdivisions and the formation thickness by 20 vertical subdivisions. The grid spacing is normally smaller near the withdrawal interval to define the cone shape accurately. For this work we used an IBM 7074 digital computer having a core memory of 10,000 ten-digit words. A typical study, covering 5 to 10 years of well producing history, required from three to six hours of computing time. MATHEMATICAL SIMULATION OF CONING BEHAVIOR BASIC METHOD In coning calculations, the reservoir volume drained by the producing well is represented by a two-dimensional system of blocks as shown in Fig. 1 for water coning studies. The horizontal dimensions of the blocks increase with radial distance from the well axis in geometric progression, i.e., the block size is small near the wellbore and large near the well&apos;s drainage radius (re). Vertically, the blocks are bounded by horizontal planes located at different depths through the

Jan 1, 1965

By R. J. Blackwell, A. L. Pozzi

Scaled laboratory-model studies provide a powerful method for evaluation of a proposed oil-recovery process. In recent years, models have been used extensively to evaluate processes in which solvents displace oil, both for general cases and for specific reservoir conditions. Since the performance of a miscible flood in a horizontal reservoir can be significantly affected by transverse mixing between solvent and oil, this displacement mechanism must be accurately simulated in the scaled model studies. Unfortunately, precise scaling of transverse dispersion coupled with the requirement of geometric similarity requires impractically large laboratory models and long times for experiments. If scaling requirements for miscible displacements could be relaxed while accurate simulation of essential displacement mechanisms is maintained, the utility of model studies would be greatly enhanced. The purpose of the work reported herein was to evaluate the relative importance of various mechanisms affecting miscible displacement and to ascer tain whether the essential features of the displacement process can be simulated even though some scaling groups are not satisfied. These studies were performed with completely miscible systems in linear, horizontal models packed with unconsolidated media. From the experimental results, a set of relaxed scaling criteria was formulated which allows the requirements of geometric similarity and equality of the ratio of viscous to gravity forces to be omitted for specified conditions. The relaxed criteria are valid whether transverse mixing is by molecular diffusion or by convective dispersion. Correlations which permit prediction of vertical sweep efficiencies in linear, horizontal reservoirs were developed from the experimental data when transverse mixing is by molecular diffusion. These same correlations may be used when transverse mixing is by convective dispersion if an empirically defined, effective, transverse dispersion coefficient is used in the description of the mixing process. The effective transverse dispersion coefficient correlation essentially duplicates the dispersion coefficient correlation for equal-viscosity, equal-density fluid systems. Experimental values for the effective transverse dispersion coefficient can be measured readily. INTRODUCTION One of the most effective methods for evaluation of miscible-displacement oil-recovery processes is that of displacements in laboratory models scaled to simulate reservoir conditions. For these laboratory studies to be meaningful, however, the essential displacement mechanisms affecting reservoir performance must be accurately simulated. Since the performance of a miscible flood in horizontal reservoirs, or in dipping reservoirs at high rates, can be significantly affected by transverse mixing of solvent and oil, this mechanism must be considered in the design of laboratory experiments. Unfortunately, precise scaling of transverse dispersion coupled with the requirement of geometric similarity requires impractically large laboratory models and long experiment times. This difficulty seriously limits the utility of laboratory model studies. Craig,et al,1 demonstrated that geometric similarity is not required when mixing is unimportant. Their experimental data indicate, for the cases studied, that the displacement is sufficiently characterized by scaling the ratio of viscous-to-gravitational forces. This work suggests that relaxation of the requirement of geometric similarity and, possibly, other criteria might also be permissible when mixing is important, provided suitable groups describing the mixing process are scaled. The purpose of the work reported here was to evaluate the relative importance of various mechanisms affecting miscible displacement and to ascertain whether the essential features of the displacement process can be simulated even though some scaling groups are not satisfied. Consequently, laboratory investigations were made of the dependence of recovery on dispersion and other variables,

By J. A. Danko, C. D. Pears, L. D. Schmidt, J. P. McGee

This paper describes a pilot plant built on the campus of West Virginia University by the Federal Bureau of Mines station at Morgantown, West Virginia, for gasifying pulverized coal in a mixture of oxygen and steam at pressures in the gasifier up to 450 psig. The work upon which this report is based was done under a cooperative agreement between the Bureau of Mines and the West Virginia University. The goal of Bureau of Mines work at Morgantown is development of a low-cost process for producing mixtures of carbon monoxide and hydrogen for the modified Fischer-Tropsch synthetic-fuel plants, which may be shifted to produce hydrogen for coal-hydrogenation plants. Gasification with oxygen offers the inducement that the gas may be generated at the pressure at which the gas will be utilized in a gas-synthesis plant. Since the optimum pressure for the synthesis appears to be about 400 psig, the gas should be produced at approximately 450 psig to take care of pressure drops in the purification train and control valves. If the gas is generated at the pressure of utilization, only the oxygen need be compressed, not the raw synthesis gas. This means a saving in power of 5 hp per hour per 1000 cu ft (CO + H2), which is approximately a saving of 4 cents per gallon of gasoline. For example, consider a commercial plant generating 18,000,000 standard cubic feet per hour of raw gas that would produce approximately 10,000bbl of liquid product per stream-day. Assume that this gas would be 85 pd CO + H2 and that the oxygen required would be 320 cu ft per 1000 cu ft of CO + H2. The following estimate is based on a commercial centrifugal compressor. If the 18,000,000 std cu ft per hour of synthesis gas is compressed to 450 psig with 74 pd efficiency, 107,000 hp is required for compression. If the synthesis gas is generated at 450 psig pressure, only

Jan 1, 1953

By Roy Lindgren

SINCE the year 1643, when the first blast furnace in America for treating iron ore was built at Saugus, Mass., out of mica schist quarried in the neighboring district, the procurement of a suitable refractory for furnace lining has been a problem of concern to the operators of furnaces. The stacks built of mica schist continued to smelt iron ore until about 1836, when, according to F. H. Norton, the first firebrick were produced1, at Queens Run, Pa. Other writers speak of brick having been molded and burned in Massachusetts about the year 1834. In 1841, Andrew Russell began to produce medium refractory plastic clay brick near East Liverpool, Ohio, that were used for lining blast furnaces1. The well-known Kentucky clay-producing district was not opened up until the year 1871, but since then it has produced a large percentage of the linings for iron blast furnaces. While some strides have been made by the refractories industry during the 100 years that have passed since the first firebrick were produced, it has been only during the last two decades that any real progress has been made towards bettering the product, even though the method of production had improved. Perhaps the fault lies with the user of the brick rather than with the producer, for not sooner demanding a supe-rior product. During the past 15 years the tonnage produced per lining has increased from 500,000 gross tons to 1,000,000 gross tons, and now some furnaces are producing 1,600,000 gross tons and better on a single lining. It is true that enlarged capacity of furnaces and improved practice have accounted for some of this increase in tonnage, nevertheless better quality in firebrick must be given credit for its share. However, we are not yet ready to say that we have reached a maximum life of furnace lining. We believe that a better product can be produced and that the refractories industries of America will, through their extensive research depart-

Jan 1, 1936

By John Dasher, Wilfred L. Freyberger, James H. Pannell, A. M. Gaudin

WHEN beryllium is bombarded by gamma rays under suitable conditions it evolves neutrons. This nuclear reaction, Bes + y + n + Be", [I] is utilized in the process that is described in this paper. The equation, which physicists often abbreviate to Bevy, n)Be" shows that a photon, y, is absorbed by the Be- sotope of beryllium which then emits a neutron, n, leaving Be8 as a residue." It has been determined that the threshold for this reaction is 1.63 mev; that is, gamma rays having an energy less than 1.63 million electron volts do not have enough energy to remove a neutron from the beryllium nucleus. The threshold value for the beryllium (y, n) reaction is lower than that for any other element. The next most sensitive element is heavy hydrogen (deuterium) whose gamma neutron reaction, D2 (-y, n) H1, has a threshold of 2.2 mev. Thus, gamma rays having an intensity between 1.63 and 2.2 mev will evolve neutrons from beryllium and from beryllium only. Nuclear reactions can be considered as occurring between a projectile and a target. The probability of the occurrence of any such encounter can be expressed in terms of the "cross section" of the target nucleus for the projectile and projectile energy. This cross section, then, is the area normal to the path of the projectile controlled by the target for the purpose of the reaction under consideration. Cross sections are usually expressed in barns, in which a barn is lo-" sq cm. The cross section for Eq 1 is very small, even for nuclear reactions, and of the order of 10"3 barns. It is therefore necessary that a high flux of gamma radiation at the appropriate energy be used to obtain an appreciable neutron yield. In view of the nuclear properties of beryllium it seemed that if pieces of beryllium ore on a moving belt were bombarded with an intense beam of gamma rays, those pieces which are rich in beryllium minerals would emit neutrons at a sufficiently high rate to make it possible to detect that emission, and to use it to actuate a mechanical device for removal of the active from the barren particles. The principal economic source of beryllium is the mineral, beryl. Beryl occurs mainly in pegmatites which are characteristically coarse grained. The suitability of mechanical concentration depends upon the size at which the desired mineral can be liber-

Jan 1, 1951

By James H. Pannell, John Dasher, Wilfred L. Freyberger, A. M. Gaudin

WHEN beryllium is bombarded by gamma rays under suitable conditions it evolves neutrons. This nuclear reaction, Bes + y + n + Be", [I] is utilized in the process that is described in this paper. The equation, which physicists often abbreviate to Bevy, n)Be" shows that a photon, y, is absorbed by the Be- sotope of beryllium which then emits a neutron, n, leaving Be8 as a residue." It has been determined that the threshold for this reaction is 1.63 mev; that is, gamma rays having an energy less than 1.63 million electron volts do not have enough energy to remove a neutron from the beryllium nucleus. The threshold value for the beryllium (y, n) reaction is lower than that for any other element. The next most sensitive element is heavy hydrogen (deuterium) whose gamma neutron reaction, D2 (-y, n) H1, has a threshold of 2.2 mev. Thus, gamma rays having an intensity between 1.63 and 2.2 mev will evolve neutrons from beryllium and from beryllium only. Nuclear reactions can be considered as occurring between a projectile and a target. The probability of the occurrence of any such encounter can be expressed in terms of the "cross section" of the target nucleus for the projectile and projectile energy. This cross section, then, is the area normal to the path of the projectile controlled by the target for the purpose of the reaction under consideration. Cross sections are usually expressed in barns, in which a barn is lo-" sq cm. The cross section for Eq 1 is very small, even for nuclear reactions, and of the order of 10"3 barns. It is therefore necessary that a high flux of gamma radiation at the appropriate energy be used to obtain an appreciable neutron yield. In view of the nuclear properties of beryllium it seemed that if pieces of beryllium ore on a moving belt were bombarded with an intense beam of gamma rays, those pieces which are rich in beryllium minerals would emit neutrons at a sufficiently high rate to make it possible to detect that emission, and to use it to actuate a mechanical device for removal of the active from the barren particles. The principal economic source of beryllium is the mineral, beryl. Beryl occurs mainly in pegmatites which are characteristically coarse grained. The suitability of mechanical concentration depends upon the size at which the desired mineral can be liber-

Jan 1, 1951

By Billy L. Crowder, John M. Fairfield

The implantation of B+ , P+, and As&apos; into silicon has been studied with the purpose of making shallow p-n junctions. The influence of such parameters as 1) ion energy, 2) target orientation and temperature, 3) total dose, and 4) annealing schedule was investigated. An energy range of 70 to 300 kev was used for boron and phosphorus implants and up to 500 kev for arsenic. It is found that the experimental projected range agrees well with theory and that shallow junction depths can be made reproducibly. ION implantation has received much attention recently as a technique for doping semiconductors. Specifically, it has the potential of supplementing or replacing the diffusion process as a method for making p-n junctions. In a few specific cases it has been used successfully to make semiconductor junction devices. Potential advantages of ion implantation doping over diffusion techniques are: 1) It affords greater control of shallow junction depths (< 0.2 µ) while maintaining high peak concentrations. This is particularly important for high-speed switching devices, since lower junction capacitances and resistances can be achieved. 2) More precise registration of small planar structures can be realized if proper masking procedures are employed. This advantage is especially useful in the design of high-density integrated circuits. It has been used to advantage in FET fabrication since the edge of the source or drain can be aligned precisely at the edge of the gate electrode.&apos; 3) Ion implanatation permits lower temperatures than diffusion techniques. This factor alleviates the problem of compatibility of diffusivities often encountered when designing multiple-junction structures. Also, the lower temperatures create fewer thermal defects and dislocations, which may account for the high efficiency of some ion-implanted solar cells.2 4) Impurity profiles can be more easily tailored to resemble ideal distributions. Successful exploitation of the potential advantages of ion implantation techniques will depend on increased knowledge and understanding of the subject. The factors likely to be influential in determining impurity distribution profiles in ion-implanted single-crystal targets have been reviewed by J. F. Gibbons.3 In addition to the mass and energy of the implanted ion, the total dose, target orientation, and target temperature are important parameters. The annealing temperature required for removing lattice damage and incorporating the implanted species on an electrically active site is very important. This paper describes an investigation of some of these factors. Implants of boron, phosphorus, and arsenic into silicon have been studied. Energy ranges of 50 to 300 kev were used for boron and phosphorus and up to 500 kev for arsenic. In addition to the implantation energy, the effects of total dose, target temperature, and post implant anneal have been investigated. EXPERIMENTAL PROCEDURE The implantation targets were silicon wafers cut from Czochralski-grown crystals, lapped, and chemically polished. The orientations were (111), (110). and (100) with misorientations of up to 7 deg from the principal axis. For this study, accurate target alignment (i.e., within 0.1 deg) was not available and quoted misorientation values should be regarded as approximate . The implantation equipment consisted of an ion source, a 300-kev linear accelerator tube, an electromagnetic separator, and the associated target supporting and beam focusing assemblies. The ion source was a simple oscillating electron type source,4 which has been described elsewhere.5 The gaseous compounds BF3, PF5, and AsH3 were used as ion sources for B+, P+, As+, and AS+&apos;. Analyzed current levels of up to 20 pamp could be obtained; however, for this investigation target current levels of 1-3 µ amp were usually employed. The analyzed ion beam was collimated through a double slit (1.4 x 0.4 cm) and swept perpendicularly to the long axis of the slit such that an area of about 2 sq cm on each target was covered. Dosages of around 1015 cm-2 were normally employed, but smaller amounts were also used for comparison. A uniform flux density over the bombarded area was assured by the continuous use of profile monitors similar to those described by Wegner and Feigenbaum.6 Post-implant annealing was accomplished in an argon atmosphere in a temperature range of 600" to 950°C. It was not part of the purpose of this investigation to study the annealing kinetics; however, some isochronal and isothermal anneal experiments were conducted to determine the time and temperature necessary to render a reasonably high portion of the implanted ions electrically active (i.e., higher than 50 pct). Post-implant anneal temperatures of around 900° and 600°C were required for boron, and arsenic and phosphorus implants, respectively. Arsenic and phosphorus implants increased in conductivity rather abruptly at the proper anneal temperature of the isochronal curve, but boron increased more gradually over a wider range. Isothermal anneal curves were reasonably flat after 10 min, so an anneal time of 1/2 hr was used for the experimental results described below. The profiling techniques were: 1) neutron activation analysis, 2) differential sheet resistance,7 and 3) junction staining.8 The differential sheet resistance technique is commonly employed in this type of study. Its principal disadvantage is the uncertainty of the ef-

Jan 1, 1970

By J. Federico Price, Joseph E. Aparicio

This paper describes the production of a specific iron oxide pellet that has given very good and consistent results as feedstock for HYC-Process Direct Reduction Plants. The sponge iron originated from these pellets has been used for the economic production of almost 15 Mt of low carbon steel via electric arc furnace meltshops. An overview of mining, beneficiation, pelletizing, direct reduction and steelmaking for this case is given with emphasis on the charactertitics of the self-fluxing oxide pellets.

Jan 1, 1982

By A. D. McMaster, M. L. Gimpl, N. Fuschillo

Low-power, high-speed, radiation-resistant, monolithic thin-film integrated circuits require thin-film resistors of high sheet resistance which are compatible with the processing requirements for monolithic silicon devices. For this purpose the structure and electrical resistance as a function of temperature were investigated for Nichrome-silicon monoxide cermet films on passivated single-crystal silicon substrates. A suitable material was found to be a mixture of 50 vol pct Nichrome-50 vol pct silicon monoxide. This mixture, when evaporated onto 300°C substrates. showed a low temperature coefficient up to 400°C, at both 270 and 1100 ohms per square. Resistors deposited on room-temperature, 300°, and 500°C substrates were heat-treated at 600°C in air for short periods of H- seed, low-power, radiation-re sis tant, monolithic thin-film integrated circuits&apos;-3 require thin-film resistors of high sheet resistance and low temperature coefficients of resistance which are compatible with the processing requirements for monolithic silicon active devices. The present investigation is directed towards the replacement of the Nichrome in the sandwich-type SO2-Nichrome-SiO resistors, currently used for DTL monolithic thin-film integrated circuits, with codeposited Nichrome-SiO cermet resistors. These Nichrome-SiO resistors are provided with an overcoat of SiO and are deposited onto thermally oxidized silicon wafers used in monolithic mi-crocircuit technology. The cermet resistors are shown to provide a lower temperature coefficient of time totaling up to 45 min. Electron nzicrogvaphs and diffraction patterns were taken of the as - deposited and heat-treated films. These results were then compared with electrical-resistance measurements for each condition. The room-temperature and 300°C substrate films showed a large decrease in resistance. The 500°C substrate film remained stable throughout the entire cycle. Electron-microscopy results show that the increased stability after heating to 600°C is largely due to crystallization of the films which begins in the vicinity of 400°C. The slight oxidation of the films when heated in air had a negligible effect on the resistance. The resistance of the film was a function of the substructure. resistance and a more versatile compatibility with silicon monolithic processing. Previous work3, on chromium-SiO coevaporated cermet resistors on glass and ceramic substrates without overcoats of SiO has been reported. For the present application, the thin films should possess the following ideal requirements: a) a temperature coefficient of resistance as near zero as possible up to 300°C; b) have a sheet resistance higher than 200 ohms per sq to conserve space on the chip; c) be easy to evaporate; d) be oxidiation- and radiation-resistant; e) show long-time stability at 100°C; } be stable in N2 atmospheres up to 600°C for short periods of time (5 to 45 min). Requirements a to e are well-known. The last requirement was optional and was due to a mask elimi-

Jan 1, 1967