Search Documents

By Rollin E. Hook, Attwell M. Adair

The yield and fracture behavior of recrystallized chromiztm have been studied in order to gain an understanding of the recrystallization embrittlement-rheotropic recovery phenomenon. The duc tile-brittle tensile transition temperature of recrystallized chromium was found to vary from 225 ° to -80°C as the ratio of specimen diameter/mean planar grain size decreased. The applicability of this observation to other bcc metals is briefly discussed. Brittle fracture in recrystallized chromium was found to be initiated predominately from intercrystalline rupture nuclei. Single crystals exhibit a high degree of ductility at 23 "C. The experimental data indicate that grain boundaries play a significant part in the yield and fracture behavim. Pres training significantly reduces the transition temperature, and high temperature aging does not result in re-embrittle-ment until complete recovery or re crystallization occurs. SEVERAL theories have been used to explain the recrystallization embrittlement-rheotropic recovery phenomenon in chromium. Wain, Henderson, and ohnstone' propose a Cottrell dislocation locking mechanism due to nitrogen, and Weaver and Gross and weaver 3 propose a strain-induced precipitation mechanism involving a chromium nitride precipitate. The researches of many investigators seem to leave little doubt that brittleness in chromium is promoted by interstitial impurities. Smith and seybolt4 have shown that the ductile-brittle transition temperature is lowered as the interstitial content is reduced. weaver3 provides evidence for a 150°C reduction in the tensile transition temperature as the nitrogen content is reduced from 20 to 5 ppm. Maykuth and affee' have suggested that the amount and orientation of grain boundary impurities may be responsible for the brittle nature of recrystallized chromium. This suggestion was based on their findings that wrought iodide chromium possessing a fiberous grain structure is ductile in tension at room temperature while a recrystallized grain structure is completely brittle. Previous investigations have not studied this phenomenon sufficiently to determine its cause. Perhaps this has been largely due to previous experimental findings that brittle fracture observed in chromium has been almost entirely transgranular.' In addition McNeil and imb,' who investigated coarse-grained as-cast chromium, describe the point of initiation of fracture as cleavage, close to the (100) plane, near or at grain bound-daries, or at included particles or voids. However, Hook, Adair, and ipsitt found that fracture in coarse-grained recrystallized chromium is initiated by intergranular rupture and then propagates nearly entirely by cleavage. This result, together with their finding that room temperature ductility in bending is observed for single crystals and coarse-grained polycrystalline specimens for which the number of grains through the specimen diameter is small, suggests that grain boundaries may play an important role in the brittle behavior of chromium. The recrystallization embrittlement problem can be resolved into two questions. These are, whether the controlling mechanism is predominately one of crack initiation or crack propagation, and whether brittle fracture is initiated predominately by slip or twin induced intergranular rupture or by cleavage. The following research was designed to study the role of grain boundaries on the brittle fracture of chromium in an effort to resolve these questions. EXPERIMENTAL PROCEDURES, RESULTS, AND DISCUSSION General Procedure. All tensile specimens (except 0.04 in. diam wires) were prepared from 1/4 in. diam swaged rod fabricated from one ingot. Details of fabrication are given elewhere. The results of impurity analysis for the swaged rod and for the

Jan 1, 1963

By E. L. Dougherty, J. W. Sheldon

Using the numerical techniques shown in this paper it is possible to compute the simultaneous dynamical behavior of multiple fluid-fluid interfaces in two dimensions. Hence, fluid-fluid inter face models of several oil recovery processes can be constructed which allow prediction of the effect of fluid flow in two dimensions on the behavior of these processes while taking into account various physical effects, such as saturation gradients, phase changes and thermal stimulation. Techniques for constructing fluid-fluid interface approximations for several oil recovery processes are reuiewed. The validity of these techniques is established by comparison to experimental results. The availability of a computer program for computing the behavior of multiple fluid-fluid interface problems makes the use of potentiometric models uneconomic. The results indicate the areal sweep efficiencies obtained from potentiometric model studies at high mobility ratios may be somewhat optimistic. The results also indicate that the areal sweep efficiency curves obtained at high viscosity ratios 2—4 using miscible fluids in porous media are misleading and that considerable care must be taken in defining the mobility ratio for a displacement process, especially for miscible fluids. INTRODUCTION Many of the methods used to estimate areal sweep efficiency assume that the oil recovery process can be represented by one or more fluid-fluid interfaces. When the mobilities of the displaced and displacing phases are assumed equal, this quantity can be computed using analytical mathematical methods.l.2 For non-unit mobility ratio, experimental methods based upon the fluid-fluid interface concept include potentiometric and electrolytic blotter models. Scaled models of porous media have also been used to study this problem. Fay and Prats3 and Aronofsky4 used numerical techniques to compute breakthrough sweep efficiency, but neither of these studies were extended beyond breakthrough. Numerical techniques which have been developed for solving a general fluid-fluid interface problem are presented in a companion paper.5 These techniques have been incorporated into a program for a large-scale digital computer which can treat simultaneously the dynamical behavior of as many as six fluid-fluid interfaces. Thus it is possible to simulate the behavior of several oil recovery processes taking into account such effects as saturation gradients. Techniques for representing several oil recovery processes by fluid-fluid interfaces are presented. Processes considered are conventional water flooding with and without a mobile gas phase, hot water flooding, miscible flooding and enriched gas drive. Other processes, whose effects could be studied by this technique, but are not considered in detail here, are underground combustion and steam injection processes and various slug processes such as gas-water injection and alcohol-solvent processes. Theoretical results obtained with single interfaces for mobility ratios of 0.25, 1, 2, 4, 8, 16, 32 and 64 in a repeated five-spot are compared to experimental results for this geometrical configuration obtained using miscible fluids and conductive analogs. Predicted recovery curves for water flooding in a repeated five-spot are compared to experimental results reported in the literature. Example calculations are shown for a conventional water flood in the presence of a free gas phase and for a hot-water flood. MATHEMATICAL PRELIMINARIES We assume that the rock matrix is homogeneous with constant thickness, absolute permeability and porosity. The fluids are incompressible. We also neglect the effects of gravity.

Jan 1, 1965

By E. Swartzman, G. C. Behnke

THE free-swelling index, a numerical value for the free-swelling properties of coal, is being used to an increasing extent in specifying coal for burning in various types of equipment. The empirical standard ASTM method for determining this index uses, as a source of heat, a gas burner which requires frequent calibration because gas pressures are not constant, relatively small changes in gas pressure resulting in rather wide variations in the temperature attained in a crucible and thus in the resultant size of the coke button. Furthermore, the quartz crucibles specified for the test are not standard equipment, and relatively small variations in dimension and wall thickness from those specified often result in rather wide variations in the resultant size of the button. A method of analysis to be satisfactory not only for a scientific laboratory but also for use in industrial control laboratories under varying conditions should employ equipment which is rugged, easily calibrated, and readily replaceable by standard materials. Gas is not universally available as is electricity, and it is generally easier to calibrate and control a simple electric furnace or heater than a gas burner, the voltage not usually fluctuating to the degree that gas pressure does. Quartz crucibles are fragile, and the size specified in the test is not normally used for any other purpose in the laboratory. Although more expensive, the 15-cc platinum crucible, available in all laboratories where volatile matter in coal is determined by the standard ASTM method, is sturdy and standard and easily replaced. Because of these considerations, a modification of the standard method has been developed which employs a standard laboratory 350-w cone heater and standard 15-cc platinum crucibles. The results of numerous tests conducted with such apparatus show conclusively that the free-swelling index values as judged by the buttons produced with this equipment are exactly the same as those obtained by the standard gas burner method. Thus it is suggested that this electric cone heater method employing standard 15-cc platinum crucibles be considered as an alternative to the ASTM method, Designation: D720-46. The free-swelling index of coal is determined by means of the ASTM method, D720-46. This method, which is an adaptation of the crucible swelling test for coal of the British Standards Institution, consists of heating 1 g of coal in a special covered quartz crucible in a regulated gas flame for a period of not less than 2.5 min in which time the temperature of the inner surface of the bottom of the crucible reaches a temperature of 82025°C. The residual coke button is then compared in size and shape with a series of standard profiles shown in Fig. 1 and designated by number from 1 to 9 in increasing order of size. As the method is empirical, it requires rigid standardization of the apparatus and repeated calibration of the gas burner. The dimensions of the low form quartz crucible recommended must be rigidly adhered to as small variations in diameter or wall thickness result in variations in the rate of heat transmission and thus in the resultant size and shape of the coke button. The gas burner must be calibrated almost every time a new set of determina-

Jan 1, 1952

By J. Gurland

The effects of variations of carbon content on the constitution and properties of sintered WC-Co alloys were studied. The grain growth of tungsten carbide during sintering was measured and it was shown that it takes place mainly by solution and reprecipitation of the carbide in and from the binder at the sintering temperature. It was found that the decomposition of W3C03C by carburization during sintering is responsible for the formation of very large grains. THE properties of sintered WC-Co alloys are critically dependent upon their final composition and structure. Slight changes of cobalt and carbon contents and small variations of grain size result in marked changes of mechanical properties and performance of products. Previous investigators have studied the microstructure as a function of composition' and have dealt with the influence of cobalt upon properties.' The present work is concerned with the effect of carbon content upon constitution, microstructure, and grain growth during sintering, as well as upon the physical properties of sintered compacts. Phase Boundaries at the Sintering Temperature The WC-Co section through the ternary diagram W-Co-C has been studied by Wyman and Kelley,' Takeda,V andford and Trent," and Rautala and Norton." Two phases, namely WC and a cobalt-rich solid solution (P), appear in the sintered structure if the carbon composition corresponds very nearly to the theoretical carbon content of WC (6.12 pct C). But small deviations from the ideal carbon content will bring about the occurrence of either graphite or the double carbide W,Co3C (7) . The 7 phase was long believed to be metastable, this conclusion being based mainly on the work of Takeda. More recent studies have shown that it is a stable phase of the ternary system,"5 but the actual phase boundaries have not been established and there has been some discussion as to whether or not the 7 phase occurs at some stage of sintering in the fully carburized alloys." In view of this uncertainty it was thought advisable to undertake an investigation of the phases present during the sintering of WC-Co alloys. Commercially available tungsten, cobalt, and carbon powders were used in the present study. A typical analysis of impurities in a sintered compact is reported in Table I. The desired carbon compositions were attained by varying the amount of carbon added to the tungsten during carburization and controlling the composition during sintering by packing the compacts into beds of mixed alundum-carbon powders of suitably selected ratios. The alloys were sintered in a hydrogen tube furnace at 1400°C, then brought to equilibrium at the temperature under investigation and quenched in water. The high temperature phases were identified by metallographic and X-ray diffraction techniques. The samples were etched with a 10 pct solution of alkaline potassium ferricyanide. The phase boundaries, established by metallo-graphic examination, are shown in Fig. 1, which is a vertical section through the ternary diagram at 16 pct Co. The temperatures of the invariant planes, at 1298" and 1357°C were established by Rautala and Norton. The vertical section shows that the three-phase field, WC + 7 + L, extends into the two-phase field, WC + L, but does not include the stoichiometric composition of WC. In practice this means that the 7

Jan 1, 1955

By M. A. Cook

Sutherland used an ion-pair adsorption model to derive the author's hy-drolytic pee-acid) adsorption equation for the contact bubble curves of Wark and Cox. To do so it was necessary to postulate that the sum of coverages of surface anion sites by the collector and depressant anions is practically unity. It is here shown that the free energy of Sutherland's postulated collector ion (x-)-depressant ion (OH-) exchange reaction (for no other than the OH- to compete with X-), namely EX + (OH-) Z SOH- + X-, is in some cases positive and X- should not even adsorb. Also, in many cases the (OH-) adsorption must be quite negligible even though the collector covers only about a sixth of the mineral surface at the threshold of air bubble contact. This, however, contradicts the postulate that permits one to derive the correct bubble contact eqwztion for the ion-pair model. Furthermore, the requirement for charge neutrality leads in the ion-pair model to the (impossible) constancy of the ratio @+)/(OH-) along a contact bubble curve in which pH may change by more than seven units. HYDROLYTIC ADSORPTION EQUATION Hydrolytic (free-acid or free-base) and ion-pair (H', X-, or R', OH-) adsorption models of the collector in metal sulfide and similar flotation systems are superficially alike and somewhat difficult to differentiate. This is because they yield, under seemingly valid assumptions, equivalent forms of the equations describing the distribution of the collector between the aqueous solution and the mineral surface. This has been discussed by the author and associates1-4 and more recently by Sutherland,' both for simple collector-mineral and for more complicated collector-depressant-mineral systems. In the hydrolytic adsorption model the equations describing the contact bubble curves of Wark and coxeJ7 are as follows :2 Here X is the collector anion, D- the depressant anion, 9, and On are the fractions of the free-acid sites S on the surface of the mineral covered by collector and depressant free acids, respectively, and @ is the fraction of uncovered sites. Kx is the free-acid collector dissociation equilibrium con- stant, K1 that for the free-acid collector adsorption, KD that for the free-acid depressant dissociation, and K2 the equilibrium constant for the free-acid depressant adsorption. Eqs. 1 to 5 yield upon simultaneous solution the result is the critical free-acid concentration for bubble contact in the absence of depressant-free acid. Eq. 8 expresses the conditions along a simple collector-mineral contact bubble curve for the case HD = 0, or V = 0. Last and Cook showed that Eq. 6 applies generally and with good accuracy in the collector-depressant-metal sulfide mineral systems. Inci-dentally, it has been shown by infrared absorption studies that in some cases the free-acid adsorption process splits out water.' In this case, Eq. 2 might be written SHOH + HX Z SHX + H20, with a similar equation replacing Eq. 4 for depressant-free acid-water exchange. This possibility does not change, however, the basic equations or considerations of the hydrolytic adsorption model. Indeed, this type of free acid-water exchange reaction was predicted for fluorite and iron ore flotation by oleic acid in the dry grinding and boiling processes.410

Jan 1, 1961

By J. F. Havard

Gypsum is a useful industrial mineral found abundantly on the earth's crust. It is inexpensive to mine and process, and its calcined products have a wide range of readily controllable properties such as strength, density and setting time. It has proven to be such a useful material indeed, that ten American companies, operating some 57 calcining plants in 25 states in 1958, produced or imported in that year a total of 13,649,000 tons of the crude mineral valued at $39,361,000, which they then converted into manufactured products with a value of $329,070,000. Gypsum (CaSO4.2H2O) is calcium sulfate with two molecules of combined water, which are responsible for most of its useful properties. It is closely associated in most deposits with anhydrite (CaSO4), which contains no combined water and has limited commercial value in the United States. Man has long recognized the useful qualities of gypsum. It was utilized by the Egyptians in the construction of the pyramids and, later, by the Romans for important structures. Many of the ancient peoples prized gypsum for artistic and ornamental purposes, carving statuary from alabaster and making casts with plaster. The tabular variety, selenite, was used for windows, and the light transmitted through these windows reminded people of light from the moon (selene), hence the name selenite. The use of gypsum continued through the Middle Ages and into modem times. Benjamin Franklin was a famed exponent of gypsum because he applied ground gypsum to clover fields in the form of large letters which later read, in the more luxuriant growth, "land plaster used here." Gypsum is still used as an effective conditioner for some soils, but in recent decades, particularly in North America and Europe, it has become much more important as a major material for the building industry. Description Gypsum and anhydrite are closely associated minerals. Gypsum is found in a number of forms: massive rock gypsum, alabaster, selenite, satin spar, and gypsite. Anhydrite is found as massive rock or as smaller bodies or isolated crystals included in gypsum. Massive Rock Gypsum comprises most of the commercial gypsum and occurs as a sedimentary rock, interbedded in most places with shale, limestone or dolomite. This rock is soft and crystalline. Its color is white or shades of gray, pink, yellow, brown or other colors. Alabaster is a pure, compact and fine-grained variety suitable for use in sculpture or turning into artistic forms; it may be white or selected for attractive coloration. Selenite is gypsum crystallized in flat, transparent, flexible foliated plates, found in clays or in surface vugs in gypsum. Satin Spar is a fine silky fibrous variety with the axes of the crystalline fibers perpendicular to the vein or joint in which it has formed. Gypsite is an earthy mixture of gypsum with clay and sand in secondary formations; it is sometimes useful as a soil conditioner and has been used in making base-coat plasters. Pure gypsum contains 20.9 pct combined water, 46.6 pct sulfur trioxide (SO3) and 32.5 pct lime (CaO). Standard specifications (American Society for Testing Materials Designation C 22-50) read that minimum purity for material to be called gypsum is 70 pct CaSO4.2H20. Most of the commercial gypsum is 85 to 95 pct pure as produced. Impurities in

Jan 1, 1960

By J. S. Aronofsky

The results of polentiometric model studies and numerical computations are described. The purpose of these studies was to determine the influence of the mobility ratio on flooding efficiencies during water encroachment in petroleum reservoirs, in systems of regular well geometry. The direct line drive well pattern having relative spacing distances of one and one-half to one was investigated for water to oil mobility ratios of 10. 1. and 0.1. The results of these studies indicate that the pattern Sweep efficiency i.; very dependent upon the mobility ratio. INTRODUCTION The almost universal occurrence of water in the immediate neighborhood of oil-hearing sands lends considerable emphasis to the desirability of having a rational analysis of water encroachment into petroleum reservoirs. This problem has 1)become increasingly important with the accelerated use of water injection for secondary recovery operations. The water encroachment problem has been studied from both the theoretical and experimental viewpoints by Muskat.' He has formulated, in a precise manner, differential equations to express flow in a water encroachment system, and rigorous analytical solutions have been developed for the linear, radial, and spherical cases.' However, it becomes exceedingly difficult to solve the encroachment problem rigorously for any general two-dimensional system in which the shape of the two-fluid interface is not immediately evident from the geometry of the system. The difficulty lies in the fart that a rigorous analytical solution requires that the shape of the two-fluid interface be known simultaneously with its instantaneous posititon and with the pressure distribution on both sides of the interface. Although this paper concerns the study of the two-fluid system of oil and water. the methods described in the paper are generally applicable to other systems of incompressible and immiscible fluids in which the displacing fluid sweep-the displaced fluid down to a residual value. In each region it is assumed there is only one mobile fluid and its mobility is constant throughout the region. The Same procedure. with minor alteration, is applicable for gas injection systems provided the assumption can be made that both the displaced and displacing gases can be described by steady state formulas. For the particular two-fluid problem considered in this paper, there is a relative permeability. K, and a viscosity. p. associated with the water zone, and a different permeability and viscosity associated with the oil zone. The K/p (mobility value associated with a particular volume element of the porous medium. at any instant of time. will depend upon whether the element is ahead of. or behind the advancing interface. It is assumed that all volume elements ahead of the interface are characterized by the K/p value for oil, while all the elements behind the interface are characterized by the K/p value for water. That is. the K/p value in the oil zone is constant throughout the oil zone. and the K/p value in the water zone is also constant but different from the value in the oil zone. The symbols are defined at the end of the paper. The well pattern chosen for initial study is the direct line drive; the swept-out area at breakthrough was determined by two different methods which are described in the section on theory. These are: 1. Potentiometric analyzer 2. Numerical computation The use of the potentiometric analyzer for studying water flooding problems is adequately described in the literature, particularly for cases where the mobility in the water zone is assumed equal to the mobility in the oil zone." ' 5

Jan 1, 1952

By L. G. Jones

An approximate method of calculation is developed in this paper which allows duplication of radial unsteady-state gas flow computer results where Darcy's law applies, such as those reported by Aronofsky and Jenkins and Bruce, Peaceman, Rachford, and Rice. Moreover, the new calculation method can be used to obtain results for radial unsteady-state gas flow obeying the quadratic flow law proposed by Duwez and Green. Means are discussed for predicting well behavior at single or superimposed flow rates in finite or infinite reservoirs, determining reservoir rock properties from well-test data, reproducing and interpreting back-pressure test data, and determining the radial extent and reserves of gas reservoirs from well-test data. Example calculations are presented for gas flow following both Darcy's law and the quadratic flow law. INTRODUCTION Since the publishing of U. S. Bureau of Mines Monograph 7,1 most gas-well testing methods have been based on the equation q = y (pl2 - pw2)a, where q = production rate, pl and pw, are formation pressure and sandface pressure, respectively, and y and a are constants to be obtained from test data. These methods, used for predicting both deliverability and "open-flow" capacity of gas wells, have been useful and accurate in many cases but unsatisfactory in others. Even at best, however, they do not supply information about the formation or lead to an understanding of nonsteady-state gas flow in porous media. Many theoretically based studies2-4 of gas flow obeying Darcy's law have been made.' Since the partial differential equations which result from combining Darcy's flow law with the continuity equation are nonlinear, most of the published research consists of either numerical solutions or analytical solutions for linear approximating equations. Such solutions have been of limited value in field work due to their unhandy form and their failure to correlate most field data. There is evidence5,6 which indicates that Darcy's law is inadequate to describe gas flow at some flow rates of practical interest. A quadratic flow law: which reduces to Darcy's law at low rates, is more successful in accounting for experimentally observed behavior. This flow law has been applied successfully to a few hypothetical reservoir cases in work which has recently been published.10 However, these numerical solutions of the equations involved have been successful only on a special analog computer. Routine applications to field cases would be awkward and have not been attempted. The present paper describes an approximate method for computing nonsteady-state gas flow solutions which has been completely successful in predicting the results for both Darcy flow and quadratic flow obtained by elaborate numerical methods. The new calculation method allows determination of the observable variables in gas-well testing at constant rates. It is similar to the scheme of using a succession of steady states suggested by Muskat7 in that it makes use of steady-state and material-balance equations. It also is similar to the work of Aronofsky and Jenkins2 in that the new method includes Aronofsky and Jenkins' stabilized flow equation as a special case. It improves upon both of these calculation schemes in that it accurately describes all portions of reservoir history and suggests means of determining reservoir rock properties from well-test data. This paper deals only with production from a reservoir, in which case the rate is defined as being negative. The reservoir studied here is a homogeneous, disk-shaped porous body of uniform thickness, with all boundaries sealed except the inner radial boundary of the well.

By E. S. Machlin, Morris Cohen

GRENINGER and Troiano' were the first to establish the fact that the habit planes of mar-tensitic products are usually planes of high indices. In steels containing 0.55 to 1.4 pct C, the habit plane indices are {225}A.* In steels with more than 1.5 pct C and in Fe-Ni alloys, the indices are {259In. The presence of both habits in certain alloy steels has been reported recently." Several investigators"-" have subsequently attempted to provide an understanding of these complex habit relations. Bowles'6 work is the most recent and the most successful. Like Greninger and Troiano, Bowles concluded that the transformation occurs in two steps, the first, a strain (or displacement) that is homogeneous throughout the volume of the martensitic plate, while the second is a strain that is homogeneous only within relatively small portions of the martensitic plate. The first strain yields the measured tilt of the specimen surface, but does not produce the final martensitic structure. This and the correct lattice orientations are achieved by the second strain. Bowles also reasoned that each strain could be described by the motion of a plane of no rotation and no distortion (hereinafter called the invariant plane) in a direction not necessarily contained in the plane.** In applying these considerations to steels with {225}a habits, Bowles found that the habit plane coincided with the invariant plane of the first strain, a result which signified that the atoms move along certain paths during the transformation. Whether this correlation results from the existence of some sort of a plate nucleus6 lying along the invariant plane of the first strain, or whether the invariant plane primarily defines the preferred directions of propagation of the first strain, will be discussed in a later section. There are two obstacles to the general acceptance of Bowles' treatment: 1—it has not been shown that elements of his double-strain analysis are unique, and 2—it is not known whether his treatment applies to transformations yielding the {2591A habit. Unless these two questions are resolved, the correlation between the habit and invariant planes in the {225}A systems might be regarded as coincidental. Experimental Procedure Using a 70 pct Fe-30 pct Ni alloy as an example of a {2591A habit system, preliminary attempts were made to obtain a double-strain solution based on a Nishiyamax** set of orientation relations and on a Greninger-Troianot set, adopting the stereographic technique described by Bowles" who used the Kurd-jumow-Sachstt set in connection with the {225}A habit. These attempts were unsuccessful because the method is essentially one of trial and error and requires either exceedingly good luck or a prior knowledge of the results being sought. Consequently it was decided to measure independently the strain associated with the formation of a martensitic plate as well as the habit-plane indices of the particular plate. By means of matrix algebra, the measured strain was applied to the austenite lattice, whose orientation with respect to the strain was determined, in order to ascertain whether the resultant structure would have the lattice of martensite. This was not found to be the case, and hence it was assumed that the measured strain was the first of two Bowles-type strains, characterized by the motion of an invariant plane. A solution of this displacement was obtained in terms of the invariant-plane indices and the direction of motion of the atoms. Then the second strain

Jan 1, 1952

By B. H. Caudle, A. R. Khan

The oil production performance of thin-oil-column type reservoirs producing entirely by natural water drive are discussed in this paper. The dimensionless production histories were obtained in a series of scaled porous-medium laboratory-model tests under varying effects of rate of production, mobility ratio across moving boundary, ratio of oil- and water-sand thickness and mobility ratio across fixed boundary including the effects of different residual oil saturations in the water-invaded regions. This scaled-model study of thin oil columns indicated that greater recovery of oil results for lower values of the following: rate of production, mobility ratio across the moving boundary, water-sand thickness relative to the oil-sand thickness, and mobility ratio across the fixed boundary or residual oil saturation in the water-invaded region of oil sand The results of the study are given as dimensionless production history plots which may be used to predict the performance of the field they represent. INTRODUCTION Oil fields producing under natural water drive can be classified into two major groups: (1) an "edge-water drive" reservoir group in which the water is encroaching from the reservoir periphery, and (2) the group of reservoirs in which the oil-sand is entirely underlain by an aquifer. Reservoirs in the second group are further subclassified as "bottom-water drive" reservoirs and "thin-oil-column" reservoirs. The latter reservoirs are characterized by their relatively equal oil- and water-sand thicknesses and by having both the oil and water zone relatively thin. When oil is produced from a water-drive reservoir by a partially penetrating well, water breaks through sooner or later before all of the displace-able oil is recovered. Many analyses of this general situation have been made with varying degrees of realism in their approaches, both theoretical and experimental5.6 in nature. Pertaining to the scaled laboratory model studies of natural water drives using porous media, detailed investigations have been made by Henley, et al.7 on bottom-water drive and by Caudle8 on edge-water drive reservoirs. No laboratory model studies, however, have been made to investigate the performance characteristics of thin-oil-column reservoirs. This paper presents the results of a study on scaled models of thin-oil-column reservoirs, operated under varying conditions corresponding to the actual field practices. SCALING PARAMETERS The object of scaling is to produce a one-to-one relationship between physical and dynamic parameters of the model -and prototype. The most comprehensive treatment of scaling techniques for oil reservoirs is given by Geertsma, et a1.,9 which covers two methods for arriving at dimensionless groups for scaling. Table 1 lists the scaling parameters suitable for creating kinematic and dynamic similarities between the model and the prototype in the variables considered for the thin-oil-column reservoirs. In order to achieve similarity it is essential that the dimensionless groups listed in Table 1 have equal values in the model and the protorype simulated. Since it has been established that the role of interfacial forces is neglible in gross fluid movement in a reservoir no attempt was made to scale capillary forces- Although there are reports describing Scaling criteria for these forces,10,ll consideration of the necessary but restrictive conditions for

Jan 1, 1970

By C. R. Clark, G. Swift, D. S. Roberts

The purpose of this investigation was to measure PVT behavior of various types and combinations of heavy hydrocarbon components from the paraffinic, naphthenic and aromatic classes where the relative proportions Of the various components were selected to approximate those of natural gas or gas condensate systems. Synthetic mixtures were used so that the compositions Of the various components could be measured accurately. n-Heptane, methylcyclohexane and methylbenzene were used to represent the paraffinic, naphthenic and aromatic components of the heavy fraction. In the mixtures studied, the heavy fraction composition was held constant at O. O5-mole lraction, the balance being methane except for one mixture where 0.10-mo1e fraction of an intermediate component, propane, was added. The PVT data for the mixtures were determined in a variable-volume, constant-mass apparatus at psuedo-reducedtemperatures from 1.84 to 2.00, over a pseudo-reduced pressure range from 2.3 to 12.0. The experimental results showed that, regardless of the type of heavy material (paraffinic, naphthenic, aromatic, or combinations thereof) mixed with methane or methane and propane, the compressibility factors at equal values of pseudo-reduced temperature and pressure varied by less than 2.2 percent. INTRODUCTION PVT data is used for natural gas and gas condensate fluids in determining reserves of reservoirs and performance of wells, in metering produced fluids, and in recombining samples for laboratory studies. While there are vast amounts of PVT data reported for natural gas and gas condensate systems from which useful correlations have been developed, ' the compositional analyses for these systems and the resultant correlations typically were made with components analyzed through some arbitrary carbon number, usually C6, with the residue reported as a lumped "heavy" fraction. The heavy fraction was stoichiometrically recombined on the basis of an apparent molecular weight to give the final compositional analysis. It is virtually impossible to make a systematic analysis of the effect that variation of the relative amounts of paraffinic, naphthenic and aromatic constituents of the heavy fraction of natural systems might have on PVT behavior. Because there is considerable latitude in the relative amounts of these constituents, one speaks of crudes or condensate liquids as being "paraffin base", etc. It is of interest to determine if changes in the relative amounts of these three types of hydrocarbons in the heavy fraction cause significant changes in PVT behavior. If so, steps should be taken to describe better the nature of the heavy fraction in PVT correlations. If not, more confidence can be placed in the correlations presently employed. DESIGN OF THE EXPERIMENTAL INVESTIGATION We attempted to determine if changes in the base of the heavy fraction would cause significant changes in the PVT behavior of gas or gas condensate fluids. Synthetic mixtures were used to enable us to systematically vary the base of the heavy fraction and to analyze accurately the fluids. We found that the heavy fraction of naturally occurring condensate systems seldom exceeds 0.05-mole fraction. Thus, the mole fraction of the heavy fraction in the synthetic systems studied was set arbitrarily at 0.05, to maximize whatever deviations in PVT behavior that might occur. Further, we judged that since the C7 hydrocarbons normally are present in greater quantities than higher carbon number hydrocarbons in naturally occurring systems, the use of n-heptane, methyl cyclohexane, and methylbenzene to represent the paraffinic, naphthenic, and aromatic species in the heavy fraction would be most appropriate. Finally, to enable us to compare our

Jan 1, 1970

By F. I. Stalkup

Vapor-liquid phase equilibrium experiments have been conducted in a static equilibrium cell on mixtures of a light, 45 API stock- tank gravity reservoir fluid and a rich hydrocarbon gas containing approximately 55 mole per cent of intermediate hydrocarbons. Both apressure-vs-composition study of the gas and a simulated reservoir fluid, and a multiple-batch contact simulation of the condensing-gas-drive oil recovery process were performed. In the latter experiments equilibrium gas and liquid compositions were analyzed. Also, conventional, "condensing-gas-drive", long-tube displacement experiments of the reservoir fluid and gases of various richness were performed. The results of these experiments could not be satisfactorily interpreted by the conventional pseudo-ternary-diagram representation of multicom-ponent phase behavior. The results seem to be explained better by considering a bubble-point surface and a dew - point surface joined in a plait-point locus. Portions of the plait-point locus cannot be "seen" directly by the rich hydrocarbon gas because of curvature of the dew-point surface. In such a system, continuous injection of the rich gas over stationary reservoir fluid might form a zone of contiguously miscible compositions from pure rich gas to pure reservoir fluid by: (1) saturating the reservoir fluid with injected gas to the bubble-point surface; (2) creating by mass transfer with fresh injected gas a path of contiguously miscible compositions along the bubble - point surface to the plait-point locus; and (3) creating by mass transfer with additional injected gas a path of gas compositions along the dew-point surface up to the point where direct miscibility results between dew-point fluid and the injected rich gas. INTRODUCTION The use of the pseudo-ternary-phase diagram to illustrate miscible displacement phase behavior has been discussed by several authors.5,7 Such a representation of phase behavior is not rigorous, but the ternary diagram nevertheless gives a qualitative picture of what actually occurs in a miscible displacement process. Fig. 1 is a typical illustration of miscible displacement phase behavior by a ternary diagram. The multicomponent hydrocarbon system is divided into three pseudo-components: a light fraction containing methane and nitrogen, an intermediate fraction containing ethane through hexanes plus carbon dioxide, and a heavy fraction containing heptane and heavier components. A two-phase region is bounded by a dew-point curve and a bubble-point curve, which are joined at the critical point. The concept deduced from such a representation for miscible displacement by a condensing-gas-drive process is as follows: a rich gas G, which lies to the right of the limiting tie line through the critical point C, is injected into the reservoir and contacts reservoir fluid L, saturating the reservoir fluid to give bubble-point fluid L1 and equilibrium dew-point gas GI. Continued injection of rich gas changes the composition of the saturated liquid L1 through a series of liquid compositions lying along the bubble-point curve, until the critical composition C is reached, at which point direct miscibility with the rich gas is achieved. Some equilibrium gas with compositions lying along the dew-point curve from G1 to C is also formed in this process. Thus, continuous injection of rich gas forms a

Jan 1, 1966

By P. K. Koh

Isothermal salt bath annealing of 0.014-in. thick 3 pct Si-Fe sheet was conducted at temperatures ranging from 927" to 1260°C in order to investigate the grain-growth behavior. Within the temperature range studied normal grain growth did not take place; instead nucleated grain growth prevailed. The nucleation frequency was found to increase with the annealing temperature, with corresponding decrease in grain diameter of ultimate secondaries after long-time annealing. Large secondaries grown by low-temperature isothermal annealing are essentially of (110)[001] orientation, and the (110)[001] texture is strong. The relatively small secondaries grown by high-temperature isothermal annealing also have orientations near (110)[001] as an average but widely scattered; thus the amount of (110)[001] texture is found to decxease with isothermal annealing temperature. Grain growth and texture development at various annealing temperatures can be described by families of sigmoidal curves with the time axis. ThE cube-on-edge (110)[001] texture and predominating large secondaries in fully annealed 3 pct Si-Fe are customarily produced by a 1100°C annealing of a primary recrystallized sheet. On annealing the primary recrystallized sheet above the temperature at which secondary recrystalliza-tion starts May and turnbull' have demonstrated that secondary grains decrease in size as the temperature is increased. The phenomenon was attributed to a decrease in the (G/N) ratio where G is the growth rate of secondary grains and N is their nucleation frequency. In addition, the growth of grains in orientations rather far from (110) [0011 occurred. The present investigation is devoted to detailed rate studies of the same phenomenon. Mi-crostructure, magnetic anisotropy, and texture data were obtained on strips of 3 pct Si-Fe isothermally annealed for various lengths of time in a BaC12 salt bath at temperatures ranging from 927' to 1260°C. The texture results are reported elsewhere.2 EXPERIMENTAL PROCEDURE Commercial 0.014-in. cold-rolled silicon-iron strip (3.16 pct Si), prepared by two stages of cold rolling with an intermediate short anneal, was given a decarburizing 3-min anneal at 800°C. Metallo- graphic studies indicated complete recrystallization. In order to provide a protective atmosphere and to obtain a rapid rise to an elevated temperature, a BaC12 and NaCl fused salt bath was used to heat samples isothermally at temperatures from 927" to 1260°C for various lengths of time. Torque magnetometer measurements were made on 1-in.-diam disc specimens before and after annealing, since magnetic torque curves supply information on textures and texture changes. Surface layers of the strip samples were polished and etched to reveal the grain structure. The diameters of grains were measured from micrographs at X100. The average value of diameters of matrix grains was established. Individual grains in the micrograph were cut out separately along grain boundaries. These grains were classified either as matrix grains or growth grains. The percent weight of growth grains in the cut-out patterns was reported as percent growth. Manv of the individual grains in the annealed samples were large enough to be X-rayed with a 5-mil beam for Lauegrams. Complete (110) pole figures were obtained for the primary recrystallized structure and short-time annealed structure with CoKa radiation at 30 kv in the back-reflection range, such as (220) for (110) poles, as described recently by Newkirk and ~ruce.~ Both the X-ray optics and the technique for pole figure differed slightly from those adapted previously.4

Jan 1, 1960

By Arthur B. Cummins

Diatomite is a siliceous rock of sedimentary origin, which may vary in degree of consolidation, but which consists mainly of the fossilized remains of diatoms. These are microscopic algae of the order Bacillaria. Diatomite is a unique form of silica, there being nothing else found in nature or produced by the manufacturing arts, that has its characteristic physical structure. There are numerous sediments which contain some diatom residues, admixed with a substantial amount of other materials-clay, carbonates or silica. While these are classified as diatomaceous silts, shales, or mudstones, they are not properly diatomite, a designation which is restricted to material of such qualities that it is suitable for commercial uses. The terms diatomaceous earth and kieselgur are synonymous with diatomite. The expressions infusorial earth and tripolite are obsolete. Diatomaceous silica is a definition best used to distinguish the basic component of diatomite, that is, the substance of the silica shell itself and is the major constituent of beneficiated diatomite or processed diatomaceous products. Some highly indurated materials, such as cherty shale, porcelanite, and opaline cherts which have resulted from diagenetic changes cannot be considered as diatomite although they may have been derived therefrom. On the other hand, diatomaceous oozes or loose sediments consisting mainly of diatom residues with or without much organic matter are source materials in process of becoming diatomite. Composition and Properties Pure diatomite is composed of opaline or hydrous silica, with small amounts of associated inorganic components-alumina, iron, alkaline earths, alkali metals, and minor constituents, all of which may occur in solid solution or as an essential part of the opal complex. Crude diatomite on the other hand may contain variable amounts of organic matter and soluble salts, and the usual contaminants of sedimentary formations-particles of rockforming minerals, sand, clay, or carbonates. As occurring in nature, diatomite contains an unusual amount of free water, which may vary from 10 to 60 pct or more. Diatomite, in the dry deposits which are characteristic of most of the important world deposits, occurs as a lightweight, comparatively light colored, "chalklike" material. The apparent density (dry basis) in situ varies from 20 to 40 lb per cu ft. The structure may be either massive or finely stratified. The earth is highly porous, usually soft and easily friable. It is nonreactive with acids and may be tentatively identified by the above characterizations. Confirmation is by microscopic examination which readily reveals the nature of the material. In deposits which contain organic matter the appearance of the material is totally different. The color may vary from light tan or gray to brown, greenish or nearly black. There is frequently difficulty in distinguishing it in the field from peatlike materials, but the difference is readily apparent after calcination. Powdered diatomite or diatomite products have an apparent density (loose weight) of 5 to 16 lb per cu ft, wet densities of 10 to 30 lb per cu ft, specific gravities 1.9 to 2.35, index of refraction 1.42 to 1.48. Diatomite is essentially micro-amorphous, but does show some tendency toward orientation in structure after a cristobalite pattern. Amorphous silica has been considered as a random association of cristobalite crystallites, the size of which are only slightly larger than a unit cell.20 The prop-

Jan 1, 1960

By H. M. Davis, T. S. Jones

A methodfor accelerating the reduction of a metal oxide by carbon monoxide is to incorporate in the reduction system a nonreducible oxide which has a strong affinity for carbon dioxide and will remove it from the product gases. The effectiveness of this method was investigated by veduction experiments on specially prepared powder mixtures of SrO-ZnO ("active'.) and SrCO3-ZnO (reference). Measurements at 750o and 677°C showed that the rate of reduction of ZnO by CO was as much as an order of magnitude greater in the presence of SrO. Under the experimental conditions at least as much CO2 as generated by veduction was absorbed when SrO was present. Enlzancernent of the veduction rate by SvO resulted from modification of' gaseous djfitsion processes. THE specific question posed in this investigation was whether the reduction reaction CO + ZnO - Zn(g) + COa [1] could be accelerated by incorporating strontium oxide into the oxide charge undergoing reduction. Provided the strontium oxide was active enough, it could remove carbon dioxide from the product gases of ( 11 and form strontium carbonate, thus modifying Reaction [ 1] to CO + ZnO + SrO = Zn(g) + SrCO3 [ 2] This study of the comparative rates of Reactions [ I] and [2] can be considered as a particular case of a general scheme for increasing the rates of reduction of metal oxides by carbon monoxide. The limits of applicability of this scheme are indicated in the simple thermodynamic analysis presented in the appendix. This analysis primarily dictated the choice of the system, CO-ZnO-SrO, employed in this investigation. In addition, it was thought that a study of this scheme by use of the CO-ZnO-SrO reaction system would be of special interest in view of indications in the literature1-3 that the gaseous reduction of zinc oxide is prone to rate control by diffusion in the gas phase. The distinction should be pointed out between this and other efforts to increase the rates of reduction reactions through the use of addition agents. Additions of a wide range of chemicals have been tried with de- grees of apparent success varying from slightly negative to highly positive. The reasons for making additions have included the use of salts to produce effects akin to catalysis,4, 5 the alteration of the microstruc-tural details of the reduction process (including solid-state diffusion),6-8 and the conversion of an en-dothermic reaction into one that is exothermic or autogenous.9-12 The means used to attempt alteration of the enthalpy change has commonly been to incorporate lime into the charge, so that heat is generated by reaction with either of the reaction products, carbon dioxide or water. The use of lime in this way makes endeavors such as that of Edstrom and Wiberg similar to the present one. Furthermore, the basic principle behind the present work has previously been set forth.13 However, a detailed investigation of the type of reaction described in this work does not appear to have been conducted. EXPERIMENTAL PROCEDURE From both conceptual and physical standpoints, Reaction [2] is to be considered as the sum of two reactions:

Jan 1, 1968

By Richard P. Hight

Abrasives include the substances, natural or artificial, that are used to grind, polish, abrade, scour, clean, or otherwise remove solid material, usually by rubbing action but also by impact (pressure blasting for example). They do not include abrasive tools, for instance, lathe tools and files--or polishing agents such as waxes, which act by filling pores. Detergents and cleaners whose action is chemical rather than physical are omitted al- though some chemical-action polishes and cleaners may also contain solid abrasives, for example, many automobile and metal polishes. General Considerations The most important physical properties of materials that qualify them for use as abrasives are hardness, toughness (or brittleness), grain shape and size, character of fracture or cleavage, purity or uniformity. For making bonded abrasive products such as grinding wheels, additional important factors are stability under high heat and bonding characteristics of grain surfaces. The economic factors of cost and availability are always important. No one single property is paramount for any use. For some uses extreme hardness and toughness are needed, as in diamonds for drill bits; for others, the factors of greatest importance are hardness and ability to break down slowly under use, to develop fresh cut- ting edges when grains become worn-for example: in garnet for sandpaper neither highly cleavable or friable grains nor extremely tough grains are wanted. For still other uses, great hardness is objectionable; for example, abrasives for dentifrices and for glass-cleaning soaps. For the most efficient use in the more critical applications, the different types of abrasives are rarely completely interchangeable; thus, while crushed quartz and garnet are both used in sandpaper, the papers are not at all interchangeable in their use applications. In the last analysis, the choice of a high grade abrasive depends upon the quality and quantity of work done by the abrasive per unit of cost. Initial cost of an artificial abrasive may be much greater than that of a natural abrasive but the artificial abrasive may do so much better work than the natural one, and do it so much faster that the ultimate cost will be less. It is for this reason that artificial abrasives have largely replaced natural abrasives. Abrasive Value Mineralogical hardness or "scratch" hardness as expressed in Mohs' scale is an important property in evaluating abrasive materials, but, as noted before, it is only one of several essential properties. The mineral hardness of pure crystal almandite garnet is about 7.5, but if the crystal is crossed by incipient fracture planes, or if it contains inclusions of other minerals, the apparent or useful hardness may be much lower. While the quartz grains in a sandstone have a hardness of 7, the bond holding the grains together may be so weak that the stone is valueless as a commercial abrasive. In artificially bonded wheels and stones, the hardness, strength, and character of the bond are fully as important as the hardness of the abrasive grains. Thus, in an overall consideration of abrasive hardness of loose abrasive grains, both "scratch" hardness and toughness must be considered. In naturally or artificially bonded abrasive stones, bond characteristics are a third factor, which is most important. The problem of abrasive hardness is further complicated by the inadequacies of methods of testing hardness and of expressing relative values. The Mohs' scale is inadequate both because the methods of testing are very crude and

Jan 1, 1975

By T. E. Scott, R. P. Zerwekh

The charactevistics of the deformatiotz of fcc thorium were examined between 78 and 473 K. Activation energy, actiuation volume, and the temperature dependence of the flow stress were determined and conlpared uith several theories for the rate-controlling mechnnistn of the tempemture-dependent component of the flou~ stress. It was established that the Cottrell-Stokes law does not hold, even approximately, for thorium. Through consideration of the above factors and other parameters, it was concluded that the vandom solute barrier theory proposed by Friedel demonstrates the best agreement with experimental observations of the deformation of thorium below ambient temperature. THORIUM, the first metal in the actinide series of rare-earth elements, is interesting and important from two points of view. In the more practical vein, it is a potential reactor fuel. As such, it is important to know its properties. Of specific interest here are the mechanical properties of thorium. Second, no detailed information concerning the deformation characteristics of the actinide elements has been gathered to date. In view of the interesting physical and chemical properties of these metals, their mechanical behavior might reveal characteristics which could improve our understanding of the mechanisms controlling deformation of the more common metals. For example, it has been shown1 that thorium exhibits a yield point which is greatly enhanced by interstitial solute elements such as carbon. This is not a normal occurrence in fcc metals. The goal of this investigation, therefore, was to characterize the deformation of polycrystalline thorium. MATERIAL AND PROCEDURES Thorium, made by the magnesium intermediate alloy process described by Peterson et al.,2 was rolled and swaged to 0.065-in.-diam wire. The metal analysis is given in Table I. After preparation, the oxygen and nitrogen contents had risen to 308 and 60 ppm, respectively. The 0.065-in.-diam wires were cut into 2-in. lengths and annealed in a vacuum of lo-' Torr at 750°C for 1 hr, producing a grain size of 2735 grains per sq mm. The ends of the 2-in. sections were masked off so that a 1-in. gage length could be electropolished. Specimens were polished at 35 v in a solution of 6 pct perchloric acid in methanol which was maintained at 0°C . The gage diameters after polishing were about 0.0625 in. The tests were run on a screw-driven tensile machine at a constant crosshead speed corresponding to a strain rate of 0.005 min-'. The strain rate was changed rapidly after each 0.01-in. extension between 0.005 min-' and one of three higher strain rates which were selected in order to achieve strain-rate ratios of 5, 10, and 20. Tests were run at several temperatures from 78" to 473°K; most were tested at 423", 300°, 215", and 78°K. EXPERIMENTAL RESULTS Flow Stress vs Strain and Temperature. Typical true stress-true strain curves for thorium are shown in Fig. 1 for several temperatures over the range 473" to 78°K. The curves for 78", 215", 300") and 423°K were obtained from strain-rate change tests using the lower (0.005 min-') or base strain-rate values of stress. Additional tests, using the same base strain rate at different strain-rate ratios, gave curves which virtually superposed on those shown, indicating that strain-rate changes did not materially influence the work-hardening characteristics at the base strain rate. The curves for 273") 373") and 473°K were standard tensile tests at 0.005 min-'. The behavior of the flow stress with temperature is

Jan 1, 1968

By B. K. Larkin, H. R. Bailey

One of the most complicated and potentially one of the most promising secondary recovery methods is that of underground combustion. A number of field tests1,2,3 have been performed recently, apparently indicating that the underground combustion process is technically feasible. Relatively little has been published concerning laboratory or theoretical investigations of the underground combustion process. A few results of laboratory tube runs have been published4,5 and these data are valuable in understanding the process. The tube run data also emphasize the complexity of the thermal recovery process. Under various conditions, problems in heat transfer, multiphase fluid flow and combustion rates arise. For example, at low temperatures the combustion process proceeds so slowly that the front cannot be propagated. It is a problem in heat transfer to predict what injection rate of air is necessary to sustain combustion in a particular reservoir. Certain fluid flow aspects of underground combustion have been studied by Wilson, et a13. In the development of some of the flow equations, assumptions of adiabatic combustion were made. Vogel and Krueger1 have devised an analog model of the heat conduction process for a constant temperature cylindrical source moving radially outward in a homogeneous conducting medium. Jenkins and Ramey7, in a discussion of the Vogel and Krueger paper, present an analytical solution of a heat conduction problem related to problems considered by Vogel and Krueger. In the present paper an idealized model of the heat conduction problem is considered, which is more general than the model of Vogel and Krueger or Jenkins and Ramey. In particular, initial well heating, vertical heat losses and arbitrary frontal velocities are included. Equations for the temperature are obtained in terms of integrals. These integrals have been evaluated numerically for the case of frontal movement corresponding to a constant air injection rate and the case of constant frontal velocity. These results are presented graphically. From these graphs one may estimate conditions under which nearly adiabatic temperatures are obtained. The application of these results requires some knowledge of reservoir thermal properties, combustion rates and residual fuel content. Assuming no vertical losses it is possible to obtain explicit evaluations of these integrals in the linear con- stant velocity case and in the radial case corresponding to a constant air injection rate. PRESENTATION OF RESULTS The geometry of the radial system is illustrated in Fig. 1 where the reservoir is represented as a slab of thickness h. By making an energy balance for an infinitesimal volume, one obtains the following equation. T= x T+F(z.r).........(1) where F (z,r) =HWv1/pc (r-r1), if z h/2 = 0, if z >h/2 . In Eq. 1, V- represents the Laplacian operator and 8 is the Dirac delta function. Eq. 1 is the well-known diff-sion equation with a source term and may be applied to a linear system. For example, if the source is moving in the direction of the x axis, then the argument of the delta function is x - 1,. In deriving Eq. 1 the following approximations were made concerning the actual combustion process: (1) all heat transfer is by conduction in a homogeneous medium, (2) the formations bounding the reservoir have the same thermal properties as the reservoir, (3) the combustion fuel is consumed instantaneously as the front arrives, (4) there is a constant fuel concentration

By J. H. Reynolds

The possibility of using a semiconductor, metal, semiconductor structure as a light detector is discussed. A brief theoretical argument is presented which predicts that this structure should have pho-toresponse characteristics similur to those of a metal-semiconductor diode, but with the magnitude of the current increased by a transistor gain factor. Measurements made at low frequencies on a particulur S-M-S device agree u~th this prediction aizd shoul that significant gain factors are attainable. The current gain is substantially higher than that predicted for the metal base transistor. It is proposed that the operation of these units is more closely approxinzated by the metal-gate model which predicts higher current gains. A qualitative discussion of its high-frequency performance indicates that detection of light modulated at microwave frequencies by a S-M-S deuice is possible. ADVANCES in the technology of chemical vapor deposition'-4 have made it possible to form a S-M-S structure consisting of a thin metallic layer sandwiched between layers of single-crystal semiconductor. It has been proposed5-7 that, if the metal layer is made sufficiently thin, this structure could be used in an amplifier having several advantages over presently available transistors including higher-frequency performance. Calculations8 predict that the maximum gain-bandwidth product is nearly an order of magnitude higher than that of a bipolar transistor having the same geometry. Maximum oscillating frequencies in the microwave range are predicted. Electrical measurements made on S-M-S devices have shown that, although the theoretical frequency characteristics are not presently attainable, current and power gain is obtained. Another possible application of the S-M-S structure is its use in a light detector. This structure is a three-layer relative of the metal-semiconductor diode. Light detection using M-S diodes is the subject of considerable present interest. They are known to have a short response time and the ability to detect light having a photon energy below the semiconductor bandgap. In this paper, the mechanism by which S-M-S photocurrent is generated is discussed. An expression is derived which shows that the response characteristics are similar to those of a M-S diode, but that the magnitude of the current is increased due to the transistorlike behavior of the device. The results of low-frequency response measurements on a particular S-M-S structure are presented. These agree with the theoretical predictions and show that significant gain factors are attainable. A qualitative discussion of the high-frequency characteristics is presented which predicts that detection of light at microwave frequencies by a S-M-S device is possible. THEORY In this discussion, the photocurrent, iT, measured across terminals attached to the semiconductor emitter and collector layers is related to the current pho-toemitted from the metal base layer. The latter will be referred to as the primary photocurrent. The behavior of the primary photocurrent is determined by the characteristics of the two metal-semiconductor diodes. The diode photoemission is illustrated in Fig. 1, which shows the emitter and collector potential barriers. The metal thickness has been greatly exaggerated. As will be discussed later, the actual barrier shape can be much more complex, but this diagram illustrates the essential features of the photoemission process. An absorbed photon having an energy greater than the barrier height can emit an electron into the emitter or collector. Since this process is discussed

Jan 1, 1968

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