Search Documents

By B. Trafton, H. L. Hartman

An earlier paper1 introduced a rapid solution for mine ventilation network problems, employing the digital computer. In this paper, refinements and additions to the computer program which expand its scope of application and a solution to a ventilation network having a coal-mine configuration are presented. The revised program-written in FORTRAN for the IBM 7074 computer in Penn State's Computation Center-is applicable to any mine ventilation network employing natural splitting and having its fan located externally. It accomplishes the following: 1) Distributes an assumed total quantity of air flow among all airways of the network, performing iterations until an acceptable level of error is attained and the flow in each airway is known. 2) Computes the pressure drop (head loss) in all airways and the overall network and mine head losses. 3) Determines the operating point for the mine and a given fan, adjusting all air flows and head losses to the actual values. The basis of analysis used is the Hardy Cross method of successive approximation for solution of any fluid network, explained in detail in the earlier article. In brief, corrections to assumed air flows are computed for all the loops of a network; following distribution to the airways of each loop, the process is repeated until acceptable accuracy is attained. When a hand or desk-calculator solution is used, there is advantage to close approximation of the actual flows in the original assumption of air-flow distribution in the network. Since publication of the previous paper, in which assumption of initial quantities was also necessary in the computer solution, the program has been modified to make the first distribution itself, assigning the total quantity to one specified path in a split and then proceeding to correct the flows. The number of iterations required for a high-accuracy solution is little affected, at least in networks of few loops, 10 to 15 being sufficient to reduce the error to 1 pct or less. The principal limitation of the original program was that it contained no provision for handling a fan. Further, it did not determine the overall head loss in the network and was not designed for multiple parallel openings, such as occur in coal mines. The revised program was written to overcome these disadvantages and also to print out the adjusted flows for all airways, corrected for the system operating point. At the present time, the program can handle only an external fan. The procedure employed in solving a ventilation network is now as follows: 1) Input data consisting of airway characteristics and loop designation along with assumed direction of flow, for each branch of the network and any trunk airways, are fed into the computer. 2) Employing the resistance relationship, where R is resistance, K is friction factor, P is perimeter, L is length, A is area, and n is number of airways (where multiple intakes and returns are employed); values are computed for all airways. 3) An arbitrary, continuous route of flow through the circuit—from beginning to end—is selected and the path identified by airway numbers and read in. 4) Based on a total assumed quantity of flow Q (usually 100,000 cfm.) the computer employs the procedure described in the previous article in distributing quantities of flow among all the airways until a satisfactory level of accuracy has been attained (in the current program, the machine is turned off when the cumulative corrective error in cfm, is reduced to five times the number of airways). 5) Airway head losses H are calculated by the relation H = R Q2 and cumulated along the original designated flow route to find the overall head loss for the network. Trunk airways are then included to obtain the mine head.

Jan 1, 1964

By P. M. Kelly, E. P. Butler

The stability of a 70/30 Mn-Cu alloy aged to peak damping has been investigated using electron microscopy, X-ray diffraction, and torsional pendulum measurements. Storage at room temperature or at 100" C causes a reduction in damping capacity. This is shown to be due to the formation of subdomains which "lock" the main domain structure and so hinder their movement. After storage the damping capacity can be restored by a retransformation treatment, which generates structures identical to those found in newly aged specimens. Plastic deformation by rolling also causes a decrease in damping capacity. In this case the damping capacity can only be restored by retransformation if the amount of deformation is less than about 5 pct. THE general metallography and structure of these high damping capacity alloys has been discussed in Part I of this work.1 Part II deals with some of the variables that affect the damping capacity. In particular, the stability of the damping capacity of a 70/30 alloy aged for 2 hr at 425° C and quenched has been investigated. The temperature sensitivity of the damping capacity of a 70/30 alloy aged for 2 hr at 450° C has been reported by Birchon.2 The damping capacity begins to decrease with increasing temperature, and is very low above the M, temperature of 125°C. The cubic — tetragonal transformation is completely reversible and there is no apparent thermal hysteresis associated with the change in damping capacity. Another important effect was first reported by Dean and his coworkers.3 They noticed a decrease in the damping capacity with time for specimens of a quenched tetragonal alloy. More recently workers at the Admiralty Materials Laboratory 4 have studied the stability of an aged 70/30 alloy. They found that the damping capacity decreased to half its original value after about 1000 hr at room temperature. EXPERIMENTAL The techniques used for electron microscopy and X-ray diffraction have been described in Part I. The damping capacity was measured by means of a torsional pendulum, using rectangular specimens, 3 by 0.3 mm with a 6 cm gage length. The specific damping capacity, S.D.C., was calculated as a percent amplitude decay in free vibrations over 10 cycles, starting at an amplitude which corresponded to a 700 psi surface shear stress. Because of the stress dependence of the damping,' the values obtained were about a third those reported for a 5000 psi surface shear stress.' RESULTS a) Effect of Aging Temperature. The effect of aging temperature on the damping capacity is shown in Fig. 1. The damping capacity reaches a peak of 12 pct after aging at 400° to 425°C) then begins to decrease as the volume fraction of a Mn increases. At higher aging temperatures the structure becomes cubic again with a low damping capacity. b) Stability of the Aged Structure. The damping capacity of the aged alloy decreases with time at room temperature. This process could be accelerated by storing at 100°C, when about 100 hr were required for the damping to fall to half its original value. The results of room temperature and 100°C storage, for times up to 1000 hr, are shown in Fig. 2. The damping capacity starts to decrease after about 24 hr at 18° C and 1 to 2 hr at 100°C. It was found, however, that the damping could be almost completely restored by retransformation. The specimens were held at a temperature above the M, for a short time, 250° C for 5 min, and then quenched. The values of damping capacity for specimens that were aged, stored, and then retransformed are shown by the full points in Fig. 2. X-ray measurements of the c/a ratio and hardness readings on the specimens did

Jan 1, 1969

By J. Szekely, G. P. Martins

A formulation is given and a general analytical solution is presented for the problem of spherical phase growth in multicomponent systems. The results are valid for describing the rate of bubble growth in the absence of translational motion and, therefore, should be relevant to heterogeneous nucleation processes. Typical worked examples are given, which include bubble growth due to the simultaneous diffusion of hydrogen and nitrogen in molten steel, and the growth of CO bubbles in steel, due to the reaction: C + 0 = CO. The calculations predict that the presence of additional components would enhance the growth rate and also increase the absolute rate at which the individual components are transferred into the bubble. These findings should be relevant to degassing processes. THE problem of spherical bubble growth in the absence of translational motion is appropriate to the description of the heterogeneous nucleation and subsequent growth of gas bubbles while attached to a solid surface. Such problems are of considerable practical relevance to various vacuum degassing and steelmaking processes where bubble nucleation and subsequent growth play an important part in the overall reaction mechanism. A definitive formulation of spherical phase growth has been given in the now classical paper by Scriven1 and the analogous problem for solid phase growth has been discussed by Horvay,2 furthermore a thorough review of the relevant literature has also become available.3 All this previous work was concerned with one component or binary systems. The purpose of this paper is to show that the technique developed by Scriven may be generalized to alow the consideration of multicomponent systems. The results of this generalization would be particularly relevant, as gas bubble nucleation in the metallurgical systems previously mentioned frequently involves the simultaneous diffusion of several components. In the following a general formulation is given and typical worked examples are presented. FORMULATION Consider the growth of a spherical gas bubble due to diffusion in a quiescent liquid of infinite extent, as sketched in Fig. 1. The formulation of such problems requires the statement of the following conservation equations relating to the liquid phase: Equation of continuity Equation of motion Conservation of the diffusing species—Diffusion Equation. The following assumptions are made: a) The fluid is incompressible, quiescent, and of infinite extent b) The bubble is stationary, i.e., it undergoes no translational motion c) Spherical symmetry is observed throughout d) The mass diffusivity is constant e) Chemical equilibrium exists at the surface of the bubble, and any reactions occurring there are fast compared to the diffusion process f) Inertia and surface tension effects are neglected. Equations derived on the basis of these assumptions by previous investigators were found to provide a reasonable representation of bubble growth while attached to a surface in organic and aqueous systems4,5 For gas bubble growth in molten metal systems assumptions a) to d) chould hold reasonably well. The exclusion of chemical kinetics as a possible rate limiting step may not be universally valid, but should hold for many systems operating at high temperatures. The effect of hydrodynamic inertia may be neglected for moderate growth rates; finally, the assumption neglecting the effect of interfacial tension on the bubble growth was made for the sake of mathematical convenience; its inclusion would have required the use of quite cumbersome numerical techniques. While surface tension effects could be appreciable on the initial values of the overall growth rate, especially at low absolute pressures and in case of small bubbles, they are not thought to modify the relative effects present in multicomponent systems, the study of which is the objective of this paper. Within the framework of the above assumptions

Jan 1, 1970

By E. E. Thurlow

The Geology Div. of Northern Pacific Railway Co. has adapted to its exploration efforts field procedures for geochemical sampling and semi-quantitative analysis which have been developed through research by the U. S. Geological Survey. The sampling is done in conjunction with reconnaissance geologic mapping on aerial photographs or during field checking of maps which have been prepared largely by photo-geologic interpretation. Stream sediment or soil samples are taken and readily soluble combined heavy metals are extracted with ammonium citrate followed by colorimetric titration with dithizone at the sample site. The effectiveness of this approach has been demonstrated both in positive and negative results, thereby assisting in the elimination of large areas of land from further immediate attention and in providing targets for followup exploration. An occa sional sample is submitted for spectrographic analysis as a check on the possible presence of other elements of economic significance. Field laboratory procedures are also being implemented for determining molybdenum, mercury, and other metals which may be valuable in themselves or useful as pathfinder elements in outlining mineralized areas. A lead-zinc occurrence in Oregon and a copper deposit in Montana were readily detected and outlined by stream sediment and soil sampling and a heavy metal anomaly that lacks obvious interpretation has been found in northwestern Montana. Detailed mapping and additional sampling are followed by other exploration techniques if warranted. Northern Pacific's program of reconnaissance mineral evaluation of company-owned lands and adjoining areas was outlined in a paper presented at the American Mining Congress meeting in Las Vegas four years ago. That report included policies and procedures with respect to mineral development on Railway Co. lands and may be referred to for information on these matters.1 In a sense this paper is a progress report on the use of geochemical prospecting in Northern Pacific's continuing mineral evaluation program. At this time there are positive results but as yet there is no real success story in the form of a minable ore body. GEOCHEMICAL PROSPECTING PROCEDURES The decision to use geochemistry as a reconnaissance exploration tool must be followed by research into sampling and analytical methods and into the interpretation of results. Various analytical methods were tried and evaluated by the Geology Div., including hot and cold acid extraction with and without fusion of the sample, panning and magnetic separation followed by dithizone testing of the magnetite, and cold aqueous extraction of metallic ions followed by indicator tests. Samples were first sent to the laboratory for analysis; later field laboratories were established and now virtually all testing is done at the sample site. It was concluded that, where applicable, the quickest dependable methods for reconnaissance use are the dithizone field tests for readily extractable heavy metals. This involves the simultaneous treatment of a sample, usually of soil or stream sediment, with two immiscible liquids; one a cold aqueous extractant (ammonium citrate) and the other a solution of dithizone in toluene or xylene.2 Dithizone is a colorimetric indicator that reacts quantitatively with copper, lead, zinc, and several other metals to form colored complexes in organic solvents. Results may be expressed in ml of 0.001 pct dithizone that are required to reach the blue end point, or calibration of the test results can be carried out with standard solutions of zinc; conversion factors would then be expressed in ppm, zinc equivalent. Calibration of dithizone solution with zinc standards indicates that 1 ml of 0.001 pct dithizone is equivalent to approximately 0.25 ppm of zinc. An authoritative report on dithizone field tests authored by H. E. Hawkes in the June-July 1963 issue of Economic Geology* states "The appeal of these techniques has been the simplicity and speed of the determinations as compared with more tedious laboratory determinations of trace metals in weathered rocks, soils, or plant materials. This obvious appeal, however, has often been more than matched by misuse of the method in such a manner that the results are

Jan 1, 1964

By William J. Potter, Roger N. Pyle

Introduction The primary ingredient for survival of mining companies during a recessionary period is to be innovative and resourceful in the structuring and financing of operations. Over the past decade, many mining companies have had to retrench and restructure their operations and to redeploy their assets just to survive. During the same period, traditional sources of financing for mining companies -- commercial banks, subsidized export credit loans and long term straight debt financing in the public markets -- became less available and more expensive. Meanwhile, new markets and new financing instruments provided new alternatives for mining companies. Those companies which exploited these alternatives greatly enhanced their capacity for survival and even for profitability. This article considers the strategies adopted by some companies to adapt to changes in this recessionary environment; the inadequacy of traditional sources of financing; certain of the new financing instruments, structures and markets which have become available; and how several mining companies have utilized these new instruments, structures and financial markets to improve their competitive positions. THE IMPACT OF RECESSION ON MINING COMPANIES AND TRADITIONAL FINANCING SOURCES In its 1983 year-end article "The New Economic Order" Mining Journal characterized many of the large private sector mining companies with operations in the industrialized countries as being "gripped between the jaws of high working costs and stagnating metal commodity prices." The article went on to characterize the Canadian mining industry in particular as being "in the difficult position of having to seek to sell the bulk of its output in a world market characterized by over-capacity, where production costs were determined by internal political and economic factors while prices were set by external competitive forces of supply and demand." Miners of a wide range of commodities in the industrialized countries have come to realize that their domestic mineral resources are in many cases no longer competitive with global sources of supply. To compete and survive they would have to rely on their technical and financial expertise in mineral exploration and development and in mining and processing. Strategies For Survival Global competition has compelled management in the extractive industries to improve their worldwide cost position. This is the only means whereby a mining company can differentiate its product and attain a substantial competitive advantage. During the past decade, for example, ASARCO Inc. and Inco Ltd. chopped their payrolls, postponed expansion plans and otherwise retrenched. ASARCO reduced its capital expenditures and cash requirements for operations in 1984 to a level almost equal to its internally generated cash. The company's finances were also strengthened by increasing total lines of credit from outside sources, by selling common stock, and by reducing annual pension funding through a change in actuarial assumptions. Other companies returned to their traditional organizational structure where they had a demonstrated competitive strength. Consolidated Gold Fields, for instance, resumed its operation as a mining finance house. Its capital investment programs are now devoted solely to resource development. In mid-1984, at a time of firming zinc prices, Union Miniere S.A. used the opportunity to integrate its smelting capacity and to direct investment to certain subsidiaries at the expense of less viable operations. Smaller corporations realized that it was not enough just to be the low-cost producer. For example, Drummond Coal Company in Alabama, which dominates its regional market due to its low-costs, has also developed a world-wide trading capability and has placed incremental tonnage from its own mines as well as brokered coal produced by others, thereby minimizing capital risk. In Canada, Campbell Resources has exploited difficult times by acquiring assets cheaply. In the process it improved its capacity to withstand cyclical economic swings by diversifying from solely a copper and gold producer to a company based upon asbestos, coal, copper, gold and natural gas.

Jan 1, 1985

By H. H. Hasler

THE mining and removal of coal from two or more beds, either simultaneously or consecutively, in vertically adjacent areas have always been matters of concern to mine operators from both operating and recovery points of view. However, with the progressive introduction and use of electrical and mechanical facilities, such problems assume even greater proportions. In this paper typical cases occuring in Cambia County, Pa., are given and procedures whereby such problems may be .greatly reduced, if not entirely eliminated, are suggested. Throughout a considerable portion of the Central Pennsylvania bituminous coal region, particularly in Cambria County (Fig. I), the portion of the region to which references in this paper are made, large areas are encountered where two or more beds of coal occur and have been mined extensively since the earliest days of mining in the county. Generally, operations were commenced in one bed of coal on a property, following which, as conditions appeared to warrant, operations were commenced in another bed or beds on the same property, in some instances by the same operator. Operating results by no means have been uniformly successful, as the following review of some typical cases indicates, partly because of lack of experience in earlier cases and in certain later cases to disregard for the experiences of earlier operators. It is suggested that repetitions of such failures may be reduced largely, if not eliminated entirely, by adherence to a few basic principles founded upon past experiences. Geology and Stratigraphy The geology, stratigraphy, and relief of almost all of Cambria County have been studied and reported at length by the U. S. Geological Survey in its Johnstown, Ebensburg and Barnesboro-Patton Folios and other publications, by the Pennsylvania Geological Survey and by various individuals. All exposed rocks occurring within the boundaries of the county are of sedimentary origin, consisting of alternate beds of sandstone, slate, shale, fireclay, coal, limestone, etc., and all have been classified by the U. S. Geological Survey as extending from the upper strata of the Catskill formation, Devonian system, to and including, in a few hilltop areas, the lower strata of the Monongahela formation, Carboniferous system. Included within that range is the Allegheny formation, Pennsylvania series, which embraces all of the commercially important beds of coal that occur in the county, viz., the A, or Brookville, to the E, or Upper Freeport, both inclusive, shown in Fig. 2. The measures are traversed by a series of anticlines and synclines, extending in a northeasterly and southwesterly direction, which form broad, gently sloping basins (flat to 15 pct pitch, exceptionally to 25 pct). However, the topography of the surface is such that long lines of outcrops of all beds of coal overlying bed B, and to a limited extent bed B also, occur over large areas, a factor which contributed to the early and extensive development of mining activities in the county. Earliest mining developments of commercial importance were made in B, C', D, and E beds of coal, and it is in these beds that the principal mining operations are being conducted. To a limited extent, mining operations also have been conducted in local areas in the A, B rider and C beds. Coal beds being mined today range in thickness from approximately 28 to 60 in. Intervals between the respective beds of coal are variable but conform somewhat to the following general averages: A to B, 60 to 100 ft.; B to C, 50 ft.; C to C', 70 ft.; C' to D, 45 ft.; and D to E, 45 ft. The following are typical mining cases obserced over a period of years in various parts of the county.

Jan 1, 1952

By C. R. McEwen

This paper presents a technique for calculating the original amount of hydrocarbon in place in a petroleum reservoir, and for determining the constants characterizing the aquifer performance, based on pressure-production data. A method for doing this based on a least-squares line-fitting computationwus proposed by vanEverdingen,Timmerman andMcMahon in 1953. We found that their method would not work when there is error in the reservoir pressure data — even moderate error. The technique presented here appears to give reasonable answers when pressure data are uncertain to the degree expected in reservoir pressure determinations. The major change introduced in the present analysis is to limit the least-squares line-fitting to yield only one constant— the amount of hydrocarbon in place. The water-influx constant is then taken as proportional to the oil (or gas) in place. The constant of proportionality can be computed from estimates of effective compressibility and reservoir water saturation. It is also pointed out that the commonly used least-squares analysis assumes all of the uncertainty to be in the dependent variable. The material balance should be arranged so that this condition is fulfilled if correct inferences are to be made from statistical calculations. Examples are shown of the application of the new technique to gas reservoirs—both hypothetical and real—and to the oil-reservoir example of van Ever-dingen, Timmerman and McMahon. INTRODUCTION The amount of hydrocarbon originally in place in a petroleum resenroir can be estimated by means of the material-balance calculation. Simultaneous observations of pressure and amounts of produced fluids are required, together with the PVT data applicable to the reservoir fluids. If water encroachment is occurring, it is desirable to try to infer the behavior of the aquifer, as well as the original hydrocarbon in place, from the pressure-production data. This imposes additional demands on the method of calculation, and uncertainty in the data can result in large uncertainty in the answer. In addition, if the size of a gas cap is to be established, the whole problem becomes indeterminate, as pointed out by Woods and Muskat.1 Brownscombe and Collins2 simulated a gas reservoir and its aquifer on a reservoir analyzer and derived quantitative information on the effect of uncertainty in pressure and aquifer permeability on computed gas in place. Among the various techniques of estimating the performance of an aquifer, the method of van Ever-dingen and Hurst? based on compressible flow theory, seems to have been the most generally successful (see Ref. 4, for example). In this paper we shall confine ourselves to their representation of the aquifer. In 1953,' van Everdingen, Timmerman and McMahon5 introduced a statistical technique for deriving the amount of oil originally in place and the parameters which describe the aquifer. (We shall refer to this technique as the "VTM method'', as suggested by Mueller.6) Their example reservoir had no gas cap. It has been our experience that the VTM method gives a reasonable answer when the data are very accurate, but that inaccuracy (particularly in pressure) can cause the method to break down. The effect was first observed in gas reservoirs, but has since been seen in oil reservoirs also. In this paper we present another statistical method which has been successful in achieving a reasonable answer where the VTM method has failed. In the new method, one less parameter is derived from material-balance computations. It is assumed that values can be established for effective compressibility in the aquifer and reservoir water saturation independently of the material-balance calculation. The water-influx constant can then be obtained from these data and the quantity of hydrocarbon in place.

By R. M. Asimow

Both smith1 and wells2 have shown that the dif-fusivity of carbon in austenite increases rapidly with increasing carbon concentration. Smith has shown, in addition, that part but not all of this increase can be attributed to the concentration dependence of the activity coefficient of carbon in austenite. From the work of Eshelby', it is known that, because of the image displacements, there is a uniform dilatation of the material in a body containing random centers of dilatation. It is thus tempting to explain the enhanced diffusivity as resulting from this apparent negative pressure caused by the interstitial carbon atoms.* That is, the activation free energy for diffusion is not directly a function of concentration if we consider that the alloy is subjected to a concentration-dependent pressure. To proceed with a quantitative solution, it is necessary to determine the activation volume AV for the diffusion process. An exact calculation of AV is quite a difficult problem and does not appear to have been solved as yet. A simpler procedure is to use the semiempirical, macroscopic theory of seyes' which relates the activation free energy AG to AV. Keyes has shown that for many types of diffusion processes the following relationship is valid Where p is the Gruneisen constant and x the compressibility. This relationship is based on the assumption that AG can be identified with strain-energy terms. In general there will also be contributions to AG due to a change in the energy of the valence electrons of the metal; thus the relation must be considered largely empirical. Surprisingly, for self-diffusion where it is well-known that the electronic contribution is quite important,' the Keyes relation works quite well for many systems. For interstitial diffusion where it might be expected that the relationship is on a firmer theoretical basis, the limited experimental data available is not in particularly good agreement as can be seen from Table I. However, it must be remembered that besides the uncertain basis of the Keyes formula there is considerable uncertainty in some of the experimental values of AV. The only real justification for using the Keyes formula in the present analysis is that it gives a result which satisfactorily solves the problem at hand and, further, there are no more reliable procedures for calculating AV available. The coefficient 3 in Eq. [I] is determined such that the room-temperature value for x is to be used. The activation entropy AS is estimated as 8.95 x 10"3 kcal per OK mole using the values in Table I1 and the Zener" relation AS = PAHIT, where p is a dimensionless thermal coefficient of Young's modulus which is to be evaluated at low values of T/T,, aH is the activation enthalpy, and T, is the melting temperature of the alloy. The resulting value of AG at 1000°C is 25.6 kcal per mole. The Gruneisen constant (p = 3cup/XC,), where p is the density, is evaluated from the room-temperature data in Table I1 after calculating C, from the relation. A value of p = 8.32 g per cu cm based on the value of a. in Table I is used. The result is C, = 5.0 x lo6 erg per g, H = 2.32, and AV = 3.41 cu cm per mole at 1000°C and infinite carbon dilution. It is assumed that the activation volume is not a function of carbon concentration. The theory could be refined to consider AV a function of concentration, but due to the uncertainties in the parameters there seems little justification for this. Carbon dissolved in austenite forms a nonideal solid solution, so that in writing an equation for. the diffusivity of carbon it is necessary to take into account the change in activity coefficient with concentration. The resulting equation is16

Jan 1, 1964

By F. G. Tahmoush, N. J. Grant, E. P. Abrahamson

The effect of grain size on the brittle -ductile transition temperature of pure iron and three dilute Fe-W alloys has been investigated by slow bend tests. The brittle-ductile transition temperature for pure iron increased with increasing grain sizes from 8000 to 500 grains per sq mm where the transition temperature then decreased with further increase in grain size. The dilute additions of tungsten lowered the transition temperature for fine grain sizes, suppressed a reversal of the curve in the coarse grained portion, and raised the transition temperature for this region. The frequency of twin occurrence increased with increasing grain size and decreasing temperature. The object of the present investigation was to re-examine the effect of grain size on the brittle -ductile transition of pure iron and several dilute Fe-W alloys. Current literature normally presents a simple relationship between grain size and the brittle-ductile transition temperature, namely, that the transition temperature increases with increasing grain size. Heslop and etch' have studied the transition temperature of mild steel by V-notch Charpy impact tests for the grain size range between 200 and 6000 grains per sq mm. The result was a linear plot of In 1-1/2 vs brittle-ductile transition temperature, where 1 is the grain diameter. Recent experiments by sipos2 revealed the possibility of a reversal in the slope of the grain size vs transition temperature plot for pure iron; however, his investigation was limited in scope, and the data obtained were insufficient to permit a definitive curve. The present investigation was undertaken in an effort to confirm or extend the grain size effect observed in iron by Sipos. PROCEDURE All alloys were vacuum induction melted, using 99.95 pct Fe containing 0.03 pct 0 and 0.004 pct C. Four ingots were made, ranging from pure iron to iron with 0.035 pct W. The chemical composition of each alloy, after processing, is given in Table I. The four alloys were prepared as 45 lb heats and cast into 4 in. diam by 16 in. molds. They were then extruded into 3/4 in. rods at 1800oF. Specimens for chemical analysis, taken from various positions along the rods, revealed no variation in tungsten content along the length. The extruded rods were reduced to 0.156 in. diam stock by a series of cold swaging operations involving intermediate vacuum annealing treatments. This procedure yielded a final cold reduction in area ranging from 10 to 90 pct. The swaged rods were chemically cleaned, cold pressed to a cross section of 0.10 by 0.18 by 1.0 in. long (2.54 by 4.73 by 25 mm). The specimens were wrapped in iron foil and vacuum sealed (10-4 mm Hg) in quartz. Various heat treatments at temperatures from 700" to 1300oC were utilized to produce grain sizes from 8000 to 1 grain per sq mm. Grain sizes were determined by measuring the average area of grains in a polished section. In specimens with grain size greater than 20 grains per sq mm, the total number of grains in the cross section was counted and then divided by the area to yield grains per sq mm. Grain size was measured both for longitudinal and transverse cross sections. The average of the two measurements was plotted vs transition temperature. An Instron Table Model testing machine was used to determine the transition temperature of specimens tested in three point horizontal bending. The specimen was supported by two 1/4 in. dowels one in. apart. Bending was accomplished with a push rod, the end of which had an 1/8 in. radius. The head speed was 50 in. per min. Temperatures were established through the use

Jan 1, 1963

By J. E. Hilliard

AS is well-known, the length per unit volume, Lv, of any lineal feature (such as grain edges in a polycrys-talline specimen) can be estimated from a count of the number of point intersections with a random plane of polish. Specifically:' Lv = 2NA in which NA is the expected number of intersections per unit area of the plane. As with all other measurements in quantitative metallography, the estimation of Lv is subject to a sampling error, the magnitude of which will depend on the number of observations. The analysis is usually performed by counting the number of intersections falling within a closed test figure of known area applied to a succession of different areas on the plane of polish. The test figure can be drawn directly on the screen of a projection microscope or else superimposed on the image of the structure by means of a reticule in a focusing eyepiece of an ordinary microscope. The use of micrographs is very inefficient unless, as in the present study, these are also needed for some other purpose. Consider an analysis in which n applications are made of a test figure within which the expected number of intersections is N. The expected total number of intersections counted is then given by NT = nN If (N) is the standard deviation in N for one application of the figure and if the latter is applied to non-overlapping regions so that the N's are independent in the statistical sense, then the standard deviation in Nt is given by: a(NT) = a(N)vn Eliminating n from Eqs. [2] and [3] we obtain: Thus, k will be a constant for a given specimen and test figure. From Eqs. [ I] and [4] we find As is to be expected, therefore, the fractional error in the estimation of Lv is inversely proportional to the square root of the total number of intersections counted. It is very desirable to have at least a rough estimate of k before starting an analysis, since it will then be possible to decide in advance by use of Eq. [5J how many counts should be made in order to achieve the desired accuracy in Lv. Without such information the investigator, unless lucky, will either do needless (and tedious) work or else obtain a result with too large an error. It is very easy to calculate k for the limiting case of N — 0, i.e., when the test figure is vanishingly small. In this case there is a negligible probability of two or more intersections falling within the figure. There are therefore only two possible outcomes: a count of one with a probability of, say, p and a zero count with a probability of (I-p). Thus and Therefore: We might expect k to depend on N and therefore on the size of the test figure. Such a dependence was observed2 in the analogous case involving the estimation of grain-boundary area per unit volume by a count of the number of intersections of a test circle with traces of the boundaries on the plane of polish. In this case the theoretical value of k at the limit N — 0 was 1.41. It was also predicted that k would display an initial sharp decrease with increasing N. Experimentally, k was observed to decrease to a value of 0.62 at N - 6 and thereafter remain approximately constant. Furthermore, the value of k for N > 6 was found to be insensitive to the form of the structure provided no pronounced anisotropy or grain-size variation was present. In order to determinewhether k would exhibit a similar dependence on N in the estimation of grain-

Jan 1, 1967

By L. H. Robinson

Triaxial compression tests on Indiana limestone indicate that chemical additives in the pore fluids can increase or decrease the strength in plastic failure but not in brittle failure. Tests using a series of the sodium salts of dicarboxylic acids show that odd-number carbon-atom chains decrease the yield strength as much as 3.0 x 10 3 psi, and even-number carbon-atom chains increase the yield strength as much as 9.0 x 102 psi. ,The most effective chemical (sodium citrate) decreased the yield strength of Indiana limestone from 16.0 x 103 to 12.0 x 103 psi. The chemical additives investigated have no influence on microbit drilling rate when drilling with zero differential pressure. However, when drilling at a 2.0 x 103-psi pressure differeniial and using an odd-number carbon chain compound as the circulating fluid, the roller bit drilling rates were slower than with water whereas the drag bit drilling rates were faster. One even-number carbon chain compound added to the circulating fluid produces the opposite effect: the roller bit drilling rate is slightly faster than with water and the drag bit slower than with water. INTRODUCTION A Russian book by Rehbinderl published in 1944 generated much interest in the possibility of aiding and accelerating the destruction of solids by the action of a liquid in which the solids are deformed. This book describes a process in which microcracks are generated during elastic deformation. To describe this process, the rock is divided into three zones: (1) the zone in which the destruction takes place, (2) the predestruc-tion zone (an affected area immediately below the destroyed zone) and (3) the undisturbed specimen. Consider a load applied normally to a surface of an infinite slab of rock. Immediately below the point of application of the load, the rock is destroyed by compression or shear forces. Below this area of destruction, because of a great number of discontinuities ranging from defects to weak points in crystals of discontinuities in crystal latices, very small cracks are created in or between the crystals (often referred to as microcracks). Upon removing the load, the broken rock in the destroyed zone could be removed. In the lower predestruction zone, however, the microcracks in and between the crystals could disappear by the process of rehealing, i.e., the molecular attraction from one side of the crack to the other side being sufficient to completely eliminate the crack. Thus, a large number of these microcracks possibly exist only for a relatively short period of time. The concept of small cracks being generated within the crystals of a rock is not unreasonable. If the rock matrix is deformed below the yield point, some of the defects present within the crystalline structure could start small cracks through the interior of the grains. Upon removal of the load, these microcracks reheal (or the newly generated surfaces come back together) so that all evidence of cracks completely disappears. Under various types of loading, and in various types of fluids, the rehealing of these cracks would be time dependent. If a liquid adhered to the surface of a rock, it could penetrate into the microcracks and thus delay or prevent the rehealing process. This suggests that the prevention of healing of the microcracks may be dependent upon the wetting characteristics of the solid. Thus, rocks which are water-wet should be more easily destroyed in water than in nonpolar hydrocarbon liquid. Adsorption of a liquid could also reduce the free surface energy which could explain the decrease in strengths of materials under nonrepetitive-type loading. Rehbinder's work prompted Boozer, Hiller and Serdengecti 2 to examine the effects of pore fluids on the deformation of rocks subjected to triaxial compression over a range of confining pressures, temperatures and strain rates. Fluids which strongly adsorbed on the rock grain surfaces decreased the strength of sandstone and limestone, which deformed plastically. For certain strain

By P. R. PaeIay, A. Slibar

Using three-dimensional, stress-deformation rate equations for a Bingham plastic, an approximate solution for the laminar flow of drilling mud between the drill pipe and casing is given for the case when the velocity gradients due to axial flow are large compared to velocity gradients due to rotation. INTRODUCTION The fact that drilling mud can be treated, to a first approximation, as a Bingham plastic has been established by several investigators12,3,4,5. A Bingham plastic is a material that is rigid until a certain combination of stresses exceeds a critical value. When the stresses are sufficiently high to initiate flow, the reduced stresses have a constant ratio to the corresponding deformation rates at each point. Prediction of pressure and torque requirements to produce the desired flow of the material have been made by analytical solutions due to Bingham6, Reiner7, and Buckingham8. One part of the present means of calculation that requires clarification is the flow between the drill pipe and casing. Present methods treat the flow resulting from the rotation of the drill pipe separate from the flow produced by the axial pressure gradient. Since the non-linear behavior of a Bingham plastic implies that solutions cannot be superposed, the present analytical investigation was undertaken to obtain a solution to the problem of combined axial and tangential flow. Equations have been proposed for the three-dimensional flow of a Bingham plastic9,10,11 and it has been shown that the proposed equations lead to the usual solutions from the one-dimensional theory". Applying these equations to the flow between two concentric cylinders the problem of axial flow only13 and tangential flow only12 have been solved. In addition to the separate problems of axial and tangential flow between concentric cylinders, the authors have determined the relationship between external torque on the inside cylinder and axial pressure gradient necessary to initiate flow". The solution given in Ref. 14 points out the complexity of the present problem. The problem considered here is the case of simultaneous axial and tangential flow. The axial flow, necessary to keep the drill bit free of removed material, results from the axial pressure gradient. Tangential flow is present since the drill pipe rotates in the hole. The solution given here is valid when the shear stresses in the drilling mud due to rotation are small compared to the shear stresses resulting from the axial flow. This case for laminar flow can be encountered in oilfield drilling. STRESS-DEFORMATION RATE EQUATIONS A detailed derivation of the equations governing the behavior of a Bingham plastic is given in Ref. 10. This section briefly outlines the physical reasoning leading to this three-dimensional formulation of the stress-deforma-tion rate law. Experiments on a certain class of materials have demonstrated the following properties: (1) when a certain combination of the stresses is below a critical value no deformation of the material occurs, (2) to a first approximation this combination of stresses is independent of the hydrostatic pressure (negative average normal stress), (3) if flow occurs the ratio of the reduced stress to the corresponding deformation rate is constant at each point of the flowing material, and (4) incom-pressibility. For the one dimensional case of shear deformation, Bingham6 formulated the governing equation as dxx, = 0 for sxx < to and µdxx = s xx = to for sxx to , . where dxx is the shear deformation rate, sxx the shear stress, to the critical shear stress, and p the slope of the sxx vs dxx curve for dxx > 0. Fig. 1 shows the behavior for pure shear between two parallel plates. Taking into account the physical properties given above and assuming the critical combination of stresses to be measured by the octahedral shear stress, the fol-lowing three-dimensional equations have been proposed.

Jan 1, 1958

By M. L. Hammond, G. M. Bowers

Silicon dioxide films have been prepared on silicon wafers by the oxidation of SiH, near 400°C. Results previously reported by Goldsmith and Kern7 have been qualitatively corrobovated using a system of different geometry, and significantly greater deposition rates have been achieved. The deposition rate is independent of thickness and nearly independent of the O2/SiH, ratio over a wide range; however, the rate is strongly temperature-dependent. The structure of the deposit depends upon the cleanliness of the substrate, and the electrical field strength of the deposited oxide depends upon the nature of the electrodes. MULTILEVEL interconnects in integrated circuits using aluminum metallization require an insulator that can be deposited at low temperature. Because the A1-Si eutectic is only 577°C, the maximum deposition temperature must be below about 500°C to prevent reaction of the aluminum with the silicon and with other materials in the integrated circuit. Many methods of depositing insulating films of SiO, and other dielectrics have been reported&apos;-6 but the technique described by Goldsmith and Kern7 offers the simplest method of preparing SiO, films. The technique involves the reaction of SiH, and 0, to form SiO, films in the 300" to 500°C temperature range. Previous work using this gas system has been reported8 and the initial development work at Philco Microelectronics was done at the East Coast R&D Laboratories. 9 This paper corroborates some of the results of Goldsmith and Kern,7 shows that much higher deposition rates can be obtained than previously reported, and presents some observations on the structure of deposited SiO, films. EXPERIMENTAL CONDITIONS Fig. 1 is a schematic drawing of the apparatus used in this study. It differs from that of Goldsmith and Kern in that 1) the gases are mixed before they enter the reaction chamber, 2) there is no fritted glass diffuser, and 3) the wafers are stationary during deposition. The control thermocouple was in a well located 0.06-in. from the substrate, and three gas supplies were used: 1 pct SiH4 in N 2, O 2, and N 2. Water vapor must be excluded from the plumbing to prevent the formation of colloidal SiO, and the SiH4 must be diluted with an inert gas because SiH4-0, mixtures are extremely flammable above about 4 pct SiH The wafer temperature was determined by attaching a thermocouple to a test wafer and the temperature was maintained to +5 C. The conditions used in this investigation are listed in Table I. The temperature range was chosen for compatibility with the A1-Si system. The volume ratio of 0, to SiH4 is the actual volume of 0, delivered to that of SiH4, taking into account the initial dilution of the SiH4. This ratio, which was used by Goldsmith and Kern,7 is a useful indicator of the experimental conditions. The active/inactive ratio is the ratio of the volume of active gases, SiH4 and O2, delivered to the volume of dilution gas. The active/inactive ratio is an important parameter to specify for deposition rate determinations since it determines the amounts of reactants that arrive at the sample per unit time. Fairly high dilution ratios were necessary to achieve uniform deposits of SiO,. Most of the process characterizations were performed on the polished (111) face of 3 to 5 ohm-cm, n-type silicon wafers which were 0.009 in. thick; but the data approximately apply to deposition of SiO, on thin films of aluminum. It is difficult to obtain accurate deposition rate data over aluminum films because of coarsening the aluminum surface during exposure to 300" to 400°C. The deposited SiO, film fol-

Jan 1, 1969

By J. D. Sudbury, J. E. Landers, D. A. Shock

The external corrosion of casing is one of the most important problems facing the present day production man. A 1953 NACE report&apos; estimated the annual cost of casing corrosion at $2,200,000 in 22 fields studied. The fields represented only 10 per cent of the wells in the United States. It was significant that 77 per cent of the failures reported were due to external corrosion. A variety of means for combatting the casing corrosion problem have been discussed in papers such as the NACE survey&apos; mentioned previously. No attempt to revicw these will be made in this paper. Thc scope of this paper has been limited to a discussion of the following major items: (I) an instrument for studying external electrolytic corrosion; (2) use of this instrument to determine the quantity of protective current required to prevent corrosion of casing; and (3) problem of accelerated external corrosion of the oil string at the base of the surface string. ELECTROCHEMICAL MECHANISM OF CORROSION It IS gcncrally belleved that one major cause of external corrosion is the existence of electrolytic cells operating between various formation waters. A theoretical discussion of the role of earth potentials in casing corrosion has been presented by L. de Witte and F. J. Radd.2 A brief picture of the electrolytic corrosion action between two differcnt formations (A and B) is presented in Fig. I. This picture is a review to many but is given to insure clarity. In the anode area, oxidation takes place as iron ions go into solution in formation A as a product of corrosion. In the cathode area. reduction takes place and hydrogen atoms are formed at the boundary between the pipe and formation B. The flow of electrons. as indicated. is into the pipe at the anode or corroding area and out of the pipe into the formation at the cathode. In order to conform with regulations of this journal, the flow of a conventional current (I! is also indicated in this sketch. All subsequent references will be to this so-called conventional current flow and will be indicated by the symbol (I). As this action takes place the layer of hydrogen developed on the cathode offers a resistance to the flow of current in the cell. It is said (hat the cathode has been polarized. This polarizing film can be rcmoved by addition of a so-called dcpolarizing agent which reacts with and removes this hydrogen film. Such an agent must be present in ordcr for this corrosion reaction to go at a rate sufficient to cause a severe corrosion problem. The two common depolarizing agents are: (I) oxy-sen. which readily reacts with hydro gen to form water, and (2) certain bacteria which use atomic hydrogen as part of their life process. It is not difficult to explain the presence of these agents. Oxygen may be present in some very shallow formations. Sulfate reducing bacteria are found to exist in nearly any ground where mud pits might be dug. THE CASING POTENTIAl-PROFILE Instrumentation A sketch of the equipment* uscd in determining casing potential profiles is shown in Fig. 2. The tool is run in the hole using a hoist truck with a three-conductor cable. This tool must be run with tubing and rods out of the hole and should be run with the flow line disconnected. By proper choice of electrode size it can be run in any size casing. The electrodes must be run in an open hole or a hole filled with oil. Measurements cannot be made while immersed in water. The tungsten carbide-tipped knives serve as the electrical robes and arc spaced 25 ft apart. The potential difference along a 25-ft section of casing is measured with a potentinmeter to an accuracy of about 5 mv. Current flowing in the pipe can then

Jan 1, 1956

By C. G. Dunn, F. W. Daniels, M. J. Bolton

IN recent investigations1. a data on relative grain boundary energies in silicon iron have been obtained. The present investigation is a continuation of this work along similar lines for the purpose of obtaining additional information on boundary energies. One aim has been to extend the data for a (110) series, which was formerly studied,&apos; to small differences in orientation, since a recent theory3,4 based on a dislocation model of grain boundaries predicts a E. ? (A-ln ? )* form of energy curve for small differences in orientation. Another aim has been to provide data for a (100) series—all grains having (100) planes parallel with the surface of sheet speci-mens—and to compare the results with the predicted form of energy curve. Finally, information on the nature of approach to equilibrium conditions was to be obtained through observations on the movements of grain boundaries. Recently published papers2,3,5 have treated the mathematical problem of expressing equilibrium conditions for grain boundaries for the general case when boundary energy depends upon boundary orientation. The equilibrium equations, which relate equilibrium angles and grain boundary -free energies, contain additional terms that express the variation of energy with boundary rotation. It has been necessary in the present investigation, however, to neglect these additional terms since no data for their evaluation were available. Dropping the additional terms leads to the following approximate equations which were used: E12/sin ?3 = E18/sin?2 =E28/sin?1 or those of the form: E28 + E13 cos ?3 + E12 cos ?2 = 0 Experimental Procedure Except for sample F1, all specimens were made from two lots of silicon iron called C and L respectively. The compositions of these are listed in table I. In the preparation of 12 three-grain specimens for the present investigation, a controlled grain growth technique, which has been described previously,&apos; was used. After preparation, the specimens were notched at the boundaries to shorten them as a means of providing more rapid approach to equilibrium in the anneals. Initial grain boundary angles were determined from micrographs taken at X500. Seven anneals, totaling about two days at 1300°C and two to four days at 1400°C, were run in an atmosphere of pure dry argon as described in a previous paper.&apos; Boundary changes generally could be followed without metallographic surface preparation, because thermal etching occurred during the anneals. These observations indicated angle changes of as much as 35" in some instances, with the major changes occurring during the first anneal. Annealing was discontinued when the angles reached stationary values. Fig. 1 gives a schematic diagram for purposes of defining crystallographic and boundary directions for three grains with common plane parallel with the sheet. Differences in orientations (expressed by A) are the angles between [00l] directions. They can be calculated with the aid of the ?&apos;s. The grain boundary angles can be calculated from the w&apos;S. Results Results are given in table 11. Types of crystal-lographic planes refer to those parallel with the sheet. Orientations are given in terms of [00l] direc-tions as shown in fig. 1. Except for grain 3 of specimen H6, which was unintentionally tilted 7" out of a (100) plane, all grains generally were within 1" or 2" of the plane specified. The directions of boundaries were measured generally to within 1" on both upper and lower surfaces, and the average direc- tions were determined for use in calculating the grain boundary angles. Variations in grain boundary directions refer to deviations from average directions. Boundaries near the junction point for all specimens generally were very close to perpendicu-

Jan 1, 1951

By R. F. Blakely, P. E. Potter

Any stratigraphic section or bedding sequence can be synthesized if there is a transition procedure from one lithology or bedding type to another, and if thickness distributions of the different lithologies are known. Stratigraphic sections of a fluvial sandstone body were synthesized with five bedding types: cross-bedding, massive beds, parting lineation, ripple mark and mudstone. The transition procedure from one bedding type to another used dependent, Markovian random processes which have a memory that extends one step backward in the depositional process. As observed in nature, median grain size and sand wave thickness (cross-bedding and ripple mark) decline upward in the synthesized sections as proportions of the different bedding types change. Grain size and permeability were also incorporated into the sections. By changing the transition procedures, bed thickness distributions, rate of upward decline or sand wave height and length, different types of sections can be synthesized, thus making it possible to model many different sedimentation problems. INTRODUCTION This paper describes a general method for synthesizing stratigraphic sections and bedding sequences of sedimentary, metamorphic or igneous origin. Synthetic generation is of interest for several reasons. Close correspondence between real and synthetic sections suggests that the factors used in the synthesizing model may indeed be the correct ones, thus giving the investigator a check on his assumptions. Rapid, inexpensive simulation of many stratigraphic sections permits one to synthesize a rock body (sandstone or carbonate reservoir) or, on a larger scale, the fill of a sedimentary basin. Harbaughl gives an example of mathematical simulation of a carbonate basin. He simulated the basin in the hope that improved prediction would follow better understanding of the depositional processes. From the petroleum engineer&apos;s viewpoint it seems reasonable to believe that the synthetic generation of rock properties and their distribution in a reservoir should be relevant in the study of reservoirs. Any stratigraphic section or bedding sequence can be generated provided there is a transition procedure from one lithology or bedding type to another, and provided the thickness distributions of the different units are known. The transition procedure involves random processes that are either independent or dependent. If the depositional process is independent, previous deposition will have no influence on present deposition. However, if it is dependent, past deposition will influence either present or future deposition. Such a dependent depositional process can be thought of as having a memory that extends backward in time through one or more pulses of deposition. A process with a memory can be described by a Markov process. Because the concept of memory or dependence appears to be in accord with our understanding of many depositional processes, Markov processes were used to synthesize the bedding sequences of this study (see Appendix). The above methods are perfectly general and are appropriate for any stratigraphic section or bedding sequence: bedding types in a beach deposit, an evolving carbonate bank or the changing lithologic fill of a thick geosyncline sequence. We chose to synthesize a vertical profile of a fluvial sandstone body because its characteristics were well documented, much was known about fluvial processes and fluvial-deltaic sandstone bodies constitute an important class of petroleum reservoirs. CHARACTERISTICS AND ORIGIN OF FLUVIAL CYCLE The fluvial cycle has been well documented in recent years by Bersier,2 Allen3-5and Visher.6,7 Deposits from fluvial cycles range from 10 to 150 ft or more in thickness and are characterized by a "fining upwards": coarse sandstones with occasional conglomerates grade upward into medium- to fine-grained sandstone, and hence into siltstone and mudstone. The dominant sedimentary structure is inclined bedding: thick cross-beds in the lower half

By J. B. Cheatham

Applications of the mathematical theory of plasticity promises to lead to the solution of many drilling and rock mechanics problems. ,Because of mathematical considerations, the inelastic behavior of rock has frequently been represented by a perfectly plastic model in conjunction with a yield criteria of the Coulomb or Mohr type. The totality of all stress states for which a solid ceases to behave elastically can be represented as a limit surface in stress space. ,Probing of such limit surfaces indicates details of strain hardening which are not provided by the standard triaxial testing procedure. Probing tests of the limit surfaces have been performed on Cardova Cream limestone to provide data for extending plasticity theory to cover situations in which consolidation and strain hardening are present. Test results indicate that this highly porous limestone undergoes a permanent volume decrease when it is subjected to hydrostatic pressures in excess of 3,500 psi. A virgin sample tested under a confining pressure of 1,500 psi has a yield strength of 1,700 psi; however, if the sample is subjected to a consolidation pressure of 5,000 psi, before testing at 1,500 psi, the yield strength is raised to 2,300 psi. Thus, both consolidation and strain hardening are important considerations in describing the mechanical behavior of this limestone. Tests conducted with the axis of the core having different orientations indicate that this rock is also anisotropic. Portions of the initial and subsequent limit surfaces are determined for samples loaded either perpendicular or parallel to the bedding planes. INTRODUCTION Previous experimental work in rock mechanics indicates that no mathematically tractable constitutive theory is inclusive enough to describe the mechanical behavior exhibited by all types of rocks under all conditions of stress and temperature. Indeed, the type of deformation encountered in a single type of rock is known to depend upon the stress and temperature conditions in the rock during deformation.l-5 Certain rocks and minerals, notably those minerals composed of ionic salts, have been shown to exhibit plastic deformation when tested under conditions of high confining pressure.2 Since the mathematical theory of plasticity provides simplifications over the theory of elasticity in certain types of problems, such as those in which limit analysis can be applied, it is of interest to know under what conditions plasticity theory may be applied to rock mechanics problems. The following factors determine the nature of the deformation a particular specimen will undergo: (1) the microscopic structure of the rock, i.e., the structure visible under an optical microscope, including number of phases, porosity, distribution of phases, (2) the mineralogical structure of the solid phases, (3) the conditions of stress and the rate of change of stress and (4) the temperature. Extensive experiments on Yule marble, Carthage marble, Solenhofen limestone and other calcerous rocks indicate that these relatively nonporous rocks deform plastically under certain conditions of loading, and creep under other conditions of loading.l-3 This study is concerned with the behavior of Cordova Cream limestone (Austin chalk) which is also composed almost entirely of calcite and thus has the same mineralogical composition but, because of a rather large porosity, it possesses a different microscopic structure. This investigation was undertaken to learn if Cordova Cream limestone deforms plastically despite the embrittling effect of the pore spaces, and to provide data which can be used to determine whether the mathematical theory of plasticity can describe the mechanical behavior of Cordova Cream limestone. A brief discussion is given of some of the underlying principles of plasticity theory concerning stress space and limit surfaces. Strain hardening is defined and a discussion is given of how the mechanical behavior of a material can depend on the previous history of the loading. This is followed by a description of the experimental procedures

By H. D. Outmans

In steady vertical flow, the interface of an immiscible liquid-liquid displacement is horizontal for any flow rate below the critical. In nonvertical flow, however, the shape of the interface in the steady state does depend on the flow rate, and the purpose of this paper is to calculate the unsteady interfaces during the transition of one steady state of flow to another. A knowledge of these transient interfaces is of considerable importance in reservoir engineering where the calculation of breakthrough recovery depends on the instant the interface reaches the producing wells and on the shape of the interface at that time. Although the emphasis is put on transient interlaces, which eventually approach stable equilibrium, it is shown that if the displacement exceeds a critical rate no equilibrium is possible. The interface is then unstable and viscous fingers are formed during the displacement. The critical rate and the shape of the transient and equilibrium interfaces are affected by the effective interfacial tension; but since this effective interfacial tension appears in the calculations only in combination with the inverse square of the thickness of the medium, its effect in the reservoir would appear to be negligible compared to its significance in model experiments. INTRODUCTION Stability criteria and the early growth of interfacial disturbances in a plane parallel to the boundaries of a dipping formation in which oil is displaced by an incompressible fluid were described in a previous paper.l This type of instability is significant in thin reservoirs. However, if the reservoir has appreciable thickness, then interfacial stability in vertical planes, normal to the upper and lower boundaries, also becomes important (the displacement is supposed to be parallel to these vertical planes). The difference between the two stability problems is that, in the first case, the intersections of the interface with planes parallel to the boundaries are normal to the direction of the displacement; in the second case, the intersections, this time with vertical planes, are not normal to the displacement. Instead, they are tilted at an angle which depends on the displacement rate. The tilt of steady interfaces was calculated by Dietz2 who also determined the critical rate of displacement for stability in the vertical plane by assuming that this rate would coincide with an interfacial tilt equal to the dip of the formation. The critical rate thus calculated is the same as has been found for thin reservoirs (see Eq. 1.1 of Ref. 1 and of the present paper). Dietz&apos;s calculation of the stable tilt was verified by laboratory experiments and the agreement was found to be fairly good.3 It is doubtful, however, that stable tilts actually exist in the reservoir because a change in production rate is not followed by an instantaneous adjustment of the interface to the new rate but, rather, by a transition period during which the interface changes from one equilibrium tilt to the other. The principal objective of this paper has been to describe these transient interfaces without putting any restrictions on the flow conditions or the shape of the interface, as had been done previously. The second objective was to compute the critical velocity, taking into account capillary effects, and the third was to evaluate, at least qualitatively, the shape of the front at rates above the critical, again without making the simplifying assumptions introduced by previous investigators.2,3 In the following sections two examples are given of the calculation of interfacial motion. The first describes this motion for an initially horizontal interface in a dipping layer, and the second for a vertical interface in a horizontal layer. The mathematical formulation of the problem is nonlinear in the boundary conditions, and this prohibits its solution in closed form. Instead, the solution is obtained in the form of higher-order approximations. 1 Before proceeding to a description of the mathematical model, however, we define two quantities

By R. M. Rose, D. P. Ferriss, J. Wulff

Single crystals of tantalum, tungsten, and the binary alloys thereof were grown by electron-beam zone melting and tested in tension between 77° and 373°K. The ductile-to-brittle transition temperature increased sharply with solute concentration, for both tungsten-rich and tantalum-rich alloys; apparently, the maximum transition temperature exists near the 50 pet concentration. Twinning was observed to be a consequence of cleavage-crack propagation. Cleavage occurred on(100)planes, and extensive tearing and "river markings" were present. The onset of brittleness could not be associated with any change in deformation mode, the yield point phenomenon, or the large rise in yield stress with decreasing temperature. Cleavage was observed to initiate at either the specimen surface, or subgrain boundaries. The observation was made that the grain boundaries of Polycrystalline materials are probably even better as nucleation sites for fracture, and are responsible for brittleness at high temperatures and lower solute concentrations than those of the corresponding single crystals. TUNGSTEN and tantalum, the two high-melting neighbors in the third transition period of the periodic table, present striking contrasts. Although both are bee, tungsten has far higher elastic moduli, and the highest ductile-to-brittle transition (DBT) of all the pure metals.&apos; Tantalum, however, is ductile at liquid helium temperatures. Large strain aging and yield point phenomena are observed in tantalum,2,3 but only to a lesser extent in tungsten.4 Tungsten work-hardens at a high rate,5 whereas this effect is negligible in tantalum.&apos; Since the mutual solid solubility of these two metals is complete,7 and no ordering has yet been reported in the binary system, the possibility of depressing tungsten&apos;s DBT by tantalum additions was considered. Also, tungsten additions to tantalum could produce hard, high-melting alloys which might be ductile at room temperature. Such alloys offered excellent bases for complex alloy development, particularly with regard to the effect of rhenium on Group VI metals.8"10 In addition, there was the opportunity to examine flow and fracture modes in the vicinity of the DBT, with particular emphasis on the correlation between brittleness and fundamental deformation mechanisms. Myers11 has reported that recrystallized alloys with up to 10 pet W in tantalum are ductile at room temperature, but alloys with 20 pct W could not be swaged at temperatures as high as 500°C. The whole spectrum of Ta-W alloys has been surveyed by Braun, Sedlatschek, and Kiefferl2 using sintered powder compacts. In addition to some physical properties they reported a maximum hardness at about 50 at. pet W, dropping to much lower hardness levels for the pure components. They were able to cold work tantalum-base alloys with as much as 30 pet W and to warm work tungsten-base alloys with up to 20 pet Ta by procedures similar to those used for pure tungsten. EXPERIMENTAL In order to study deformation modes and to keep the impurity concentrations at a minimum, electron-beam zone melting was used to prepare single crystals of high purity. The electron-beam zone melting apparatus is similar to that described by Calverley, Davis, and ever" and incorporates a grid control rectifier to maintain control of emission current. Pure metals were refined by making three zone passes at 4 mm per min up the rod. Alloys were prepared by melting together wire of one constituent on a previously melted rod of the major constituent. Subsequently, two zone passes were made in opposite directions to level the composition after the manner suggested by pfann.14 Purification of the starting materials was accomplished by evaporation of the several impurities, and the degree to which this was accomplished is indicated in the chemical analyses of Table I. The uniformity of composition along a zone levelled rod was within the i 1 pet accuracy of the fluorescent X-ray technique used to determine composition. Crystallographic orientation of the

Jan 1, 1962

By P. S. Kotval, H. Hahell

PRECIPITATION of M23C6 carbides at grain boundaries in austenitic steels and nickel base superalloys has received considerable attention. Hatwell and Berghezan,1 in a study using carbon extraction replicas, observed that dendritic growth of the M23C6 carbide phase occurred at grain boundaries in austenitic steels. Merrick and Nicholson2 have observed that in thin foils of "Nimonic 75" alloy the precipitation of M23C6 carbide was discontinuous at grain boundaries. Lewis and Hattersley3 have suggested that the type of M23C6 precipitation observed in Ni-Cr austenites represented the early stages of discontinuous precipitation. In a more recent study, Singhal and Martin4 have concluded that grain boundary M23C6 precipitation involves the migration of an austenite grain boundary by a mechanism very similar to that proposed by Tu and Turnbull.5 Tu and Turnbull&apos;s suggestion essentially involves grain boundary migration under the driving force arising from the difference between the surface free energies of the opposite sides of the plate. Thus, Singhal and Martin&apos;s observation4 implies that the first stage of discontinuous precipitation has taken place and that the grain boundary has provided a diffusion short circuit for the sub-stitutional solute atoms (principally chromium) which assist the preferential growth of the M23C6 plate in contact with the grain boundary. However, no evidence of the characteristic cellular morphology typical of discontinuous precipitation was found in Singhal and Martin&apos;s studies. Presumably, this is due to the low volume fraction (i.e., low supersaturation of solute) of the precipitating carbide reported by Singhal and Martin.4 Recently, Raymond6 has studied precipitation of the M23C6 carbide in Incoloy alloy 825 and shown that the growth of this phase is perpendicular to grain boundaries at 1400°F and in the plane of the grain boundary at 1700°F. The former aging condition resulting in discontinuous or "cellular" precipitation. Raymond has suggested an explanation for this behavior based on the relative rates of chromium diffusion by volume and grain boundary diffusion at the two different aging temperatures. In the course of a recent investigation on precipitation behavior in several &apos; strengthened superalloys, precipitation of M23C6 carbide at grain boundaries was studied. Under conditions of a high degree of supersaturation, the M23C6 carbide precipitates in the discontinuous or cellular morphology in Waspalloy (Ni-18.51 Cr-2.0 Fe-3.53 Mo-1.16 Al-2.78 Ti-0.09 C). This note presents electron microscopy observations of this mode of precipitation and its effect on mechanical properties. Specimens of Waspalloy after "solution" heat treatment at 1950°F show essentially single phase matrix with a few residual particles of "primary" MC carbide. In normal commercial heat treatments, precipitation of Ni3(Al,Ti) ordered fcc r&apos; phase occurs throughout the matrix and M23C6 precipitates as globular particles at the grain boundaries. In these commercial heat treatments the degree of supersaturation during precipitation is low. However, if specimens of Waspalloy are quenched from the solution heat treatment temperature of 1950°F and then aged at an aging temperature where the degree of supersaturation is high, discontinuous precipitation of M23C6 occurs. Figs. 1 and 2 show thin foil micrographs of Waspalloy, solution heat treated at 1950°F, water quenched, and aged at 1292° and 1598°F, respectively, for 500 Fig. 1—Thin foil of Waspalloy, solution heat treated at 1950°F for 1 hr, water quenched, and aged for 500 hr at 1292°F. Magnification 15,500 times. Fig. 2—Thin foil of Waspalloy, solution heat treated at 1950°F for 1 hr, water quenched, and aged for 114 hr at 1598° F. Magnification 13,000 times.

Jan 1, 1970