Search Documents
Search Again
Search Again
Refine Search
Refine Search
- Relevance
- Most Recent
- Alphabetically
Sort by
- Relevance
- Most Recent
- Alphabetically
-
Institute of Metals Division - New Method for Measuring Surface Energies and Torques of Solid SurfacesBy P. G. Shewmon
A novel technique for determining the surface energy (?) and its derivative with respect to orientation, (?') is described. Essentially it involves the 'floating" of a wedge on the substrate, said wedge being made of a material which is not wet or only slightly wet by the substrate, i. e., as a greased needle "floats" on water. A thermodynamic analysis of a system in which the wedge is supported entirely by surface energy is given. If the original suyface is not at a cusp orientation, the surface tension is directly measurable from the groove angle formed. If the original surface is at a cusp orientation, there may or may not be a groove depending on the relative value of ?' and the weight of the wedge. Experiments primarily on copper and silver showed that sapphire, quartz and refractory metal wedges were wet while graphite wedges were not. The technique was demonstrated to work using graphite wedges, but the results obtained were not as eccurate as those obtained by other workers using the wire-creep experiments. It is concluded that the technique might prove most useful with non-metals where ?' is large and filament creep experiments would be quite difficult. If an absolute value of the surface free energy (?) of a metal is to be determined, the most reliable methods used to date measure an average over the various orientations exposed on a polycrystalline sample. For example, ? for silver, gold, and copper have been measured by determining the force required to just keep a thin wire,' or foil,' specimen from contracting under the influence of ?. Herring 3 has predicted and experiment confirms, that the sensitivity of this method is inversely proportional to the grain size.' Thus it cannot be used to measure ? for a particular orientation by using a foil single crystal or a very coarse-grained specimen. An accurate value if ? for tungsten averaged over a range of orientations has been determined using a field emission technique. The same techniques cannot or have not been used to measure ? for non-metallic solids, and as a result the values available are much less accurate.4 This Paper resents a means of making an absolute determination of ? for a particular surface orientation on any solid, as long as the given surface orientation does not break up into other orientations during an anneal. Experimentally ? is found to vary with orientation and at a few low index orientations it is found to have a cusped minimum, i.e., the derivative of ? with respect to the orientation of the surface changes discontinuously at the low index orientation, see Fig. 1. The slope of a plot of ? vs orientation (herein designated ?') is called the torque on the surface, since it tends to rotate the exposed surface toward the low index orientation, or if the surface is at the cusp orientation it opposes any force tending to rotate the surface out of the low index orientation. The ratio ?'/? has been determined for a few metals, but in cases where this ratio is high there is presently no means of determining either ?'/? or the absolute value of ?' for the orientations present on an annealed surface. The technique discussed herein also provides a means of determining an absolute value of ?' for those orientations which deviate only infinitesimally from a cusp orientation. It should work best on surfaces where ?'/? is large; that is, for cases where no other technique is available for measuring ?'. Aside from trying to learn more about surfaces through measuring ? and ?', the primary reason for wanting values of ? or ?' is to study adsorption. From measurements of the variation of ? for a particular orientation with the concentration of an impurity, one can obtain the number of impurity atoms adsorbed per unit area (Ti) on that orientation using the Gibbs adsorption equation.' where µi is the chemical potential of the adsorbed impurity. Thus, if absolute values of ? could be obtained for the free surface of a given surface orientation as a function of µi, ri could be determined for the given orientation. Furthermore, by equilibrating a grain boundary with the given surface at various values of ki, one could also determine ri for the grain boundary. Similarly Robertson 6 has pointed out that if y is taken to be a continuous function of and µi, then a2 ?/a @a µ2 = a2 ?/a pi a +. Thus, at all orientations away from cusps the following equation holds From a measurement of ?' vs ki, it is thus possible
Jan 1, 1963
-
Part VII – July 1969 - Papers - Effect of Driving Force on the Migration of High-Angle Tilt Grain Boundaries in Aluminum BicrystaIsBy B. B. Rath, Hsun Hu
In wedge-shaped bicrystals of zone-refined aluminum it is observed that (111) pure tilt boundaries migrate under the driving force of their own inter-facial free energy. The boundary velocity is a power function of the driving force. The driving force exponent decreases with decreasing angle of misorien-tation. For example, at 64O°C, the exponent decreased from 4.0 for a 40 deg to 3.2 for a 16 deg tilt boundary. An evaluation of the driving force acting on the boundaries during their motion indicates that for low driv-forces, up to about 2 x l03 ergs per cu cm, the velocity is relatively independent of misorientation, whereas at higher driving forces a 40 deg tilt boundary exhibits the highest velocity. The measured activation energy for boundary migration approaches that for bulk self-diffusion at low driving forces, decreasing from 33 to 27 kcal per mole as the driving force is increased from 1 x l0 to 5 x l03 ergs per cu cm. These results are compared with current theories of grain-boundary migration. In previous experimental studies of grain boundary migration the driving force has been limited to a difference in stored energy across the boundary. This stored energy has been introduced into the crystal either by prior deformation1-3 or by grown-in lineage structure. A part of the energy stored in the deformed crystal is released by recovery either prior to or concurrently with grain boundary migration, thus introducing an uncertainty as to the magnitude of the driving force responsible for grain boundary migration. The grown-in lineage structure, though thermally stable during annealing, neither provides conditions under which different levels of energy may be stored in the imperfect crystal nor provides a control of orientation difference across the migrating boundary of a growing grain. Furthermore, because of variation in the lineage structure, it is difficult to determine accurately the energy stored in the imperfect crystal. Several investigations of grain boundary migration during normal grain growth have also suffered from difficulties in estimating the driving force because of uncertainties in the principal radii of curvature.~ In the present investigation the velocity of pure tilt boundaries in zone-refined aluminum bicrystals of selected orientation (40, 30, and 16 deg around the [Ill] tilt axis) has been measured in the absence of a dislocation density difference across the moving boundary, thus eliminating the previous experimental difficulties. The driving force for boundary migration is derived from a gradient of the total interfacial free energy of the migrating boundary in wedge-shaped bicrystals. A similar method was attempted by Bron and Machlin in a study of grain boundary migration in silver. However, they found that one of the crystals was deformed and consequently the motion of the boundary was partly due to a difference of stored energy across the boundary. The observed behavior of boundary velocities as affected by the driving force is examined in the light of the predictions of the current theories of grain boundary migration.7"10 The effect of boundary misorientation on velocity is compared with the theory of " which is based on a dislocation core model for high-angle boundaries. EXPERIMENTAL METHOD Seed-oriented bicrystals of zone-refined aluminum, 2.5 cm wide, 0.5 cm thick, and 12 cm long, containing tilt boundaries with a common (111) axis, were grown from the melt in the direction of this axis. Spectro-graphic analysis, reported earlier,'' indicated the purity of the crystals to be 99.999+pct. Three such bicrystals containing 16, 30, and 40 deg tilt boundaries were used. Wedge-shaped specimens were prepared from these bicrystals by spark cutting followed by electrolytic polishing. The angle of the wedge was usually 40 deg and the specimens were usually 0.25 cm thick. The intercrystalline boundary was located within 0.2 to 0.5 cm from the tip of the wedge. Fig. 1 shows a section of an oriented bicrystal containing an outline of a wedge-shaped specimen. The crystallographic directions shown in Fig. 1 represent the orientation of one of the crystals (the larger section of the bicrys-tal); the orientation of the other crystal differs only by rotation around the common [lil] axis. The parallel faces of the wedge always corresponded to the common (171) planes in both crystals, whereas the orientation of the side faces varied, depending on the misorientation angle. The bicrystal orientations were determined
Jan 1, 1970
-
Producing-Equipment, Methods and Materials - Emulsion Control Using Electrical Stability PotentialBy J. U. Messenger
A technique is described whereby the resistance of an emudian to breaking can be quantitatively determined. Produced ailfield emulsions are usually the water-in-oil type and, accordingly, do not conduct an electrical current. However, there is a threshold of A-C voltage pressure above which an emulsion will break and current will flow. The more stable an emulsion, the higher the required voltage. A Fann Emulsion Tester, modified so that low voltages (0 to 10 v) can be accurately measured, is suitable. This technique has application in evaluating the effect of a demuksifier on the stability of an emulsion. Emulsions can, in essence, be titrated with demulsifiers by adding a quuntity of demulsifier, stirring, and measuring the voltage required to cause current to flow. Any synergistic effect of two or more materials added simultaneously can be followed accurately. A demulsifier that significantly lowers the threshold voltage (from 100 to 400 v to 0 to 10 v for the emulsions in this study) is effective and can cause the enlulsion to break. A demulsifier that will bring about this drop in the threshold voltage at low concentration ir very desirable. The technique is also well adapted for rapidly screening demulsifiers. INTRODUCTION Stable emulsions in produced reservoir fluids resulting from certain well stimulation and completion procedures are common problems. The use of suitable demulsifiers can often mitigate these difficulties. At the present time, a rapid and efficient method for selecting satisfactory demulsifiers is not available. It is badly needed. Reliance is now placed primarily on trial-and-error procedures. A new test method has been developed which permits a more rapid and precise selection of demulsifiers. It involves measuring the electrical stability potential of an emulsion before and after a demulsifier has been added. This paper describes this method and shows where it should have application in field emulsion problems. NATURE OF OILFIELD EMULSIONS Two immiscible components must be present for an emuhion to form; we are concerned here with crude oil and water. An emulsifier must be present for tin emulsion to be stable. J Emulsifiers can be substances which are soluble in oil and /or mter and which lower interfacial tension. They can be colloidal solids such as bentonite, carbon, graphite, or asphalt which collect at the interface and are preferentially wet by one of these phases. Unrefined crude oils can contain both types of emulsifiers, A popular theory is that, of the two phases in an emulsion, the dispersed phase will be the one contributing most to the interfacial tension.' Usually this phase contains the least amount of emulsifier. The stability of a water-in-oil emulsion is affected by the fol1owing: l) viscosity; (2) particle or droplet size; (3) interfacial tension between the phases; (4) phase-volume ratios; and (5) the difference in density between the phases. A stable emulsion is usually characterized by high-viscosity, small droplets, low interfacial tensions, small differences in density between its phases, and slow separatian of the phases. It also has low conductivity (high electrical stability potential). Water-in-oil and oil-in-water emulsions"' are both common; however, oil field emulsions are predominantly water-in-oil emulsions. The emulsions which commonly occur during oompletion and stimulation operations contain a combination of several of the following: acids, fracturing fluids (oil, water, acid), and formation water and oil. Produced emulsions usually contain formation water and oil. Emulsions form in oil wells because oil and water are mixed together at a high rate of shear in the presence of a naturally occurring or unavoidably produced emulsifier. During the completion and stimulation of productive zones, and while formation fluids are being produced, oil and water are very often commingled. These mixtures are formed into emulsions by agitation which occurs when the fluids are pumped from the surface into the matrix of the formation or produced through the formation to the surface. Restrictions to flow (such as perforations, pumps, and chokes)".'" increase the level of agitation; tight emulsions are more likely to form under these conditions. Often an emulsified droplet is an emulsion itself.'" Therefore, emulsion-breaking problems can be quite complex. The complexity can be even greater if a third phase (gas) is included. Demulsifiers operate by tending to reverse the form of the emulsion. During this process, droplets of water become bigger, viscosity is lowered, color becomes darker, separation of the phases faster and electrical stability potential approaches zero. Any of these effects could be followed as a means of determining emulsion stability. However, electrical stability potential is the most reproducible and most easily measured parameter for following the stability of a water-in-oil emulsion.
Jan 1, 1966
-
Institute of Metals Division - Multistep Reactions in the Creep of CopperBy E. R. Gilbert, D. E. Munson
Creep of copper under 75 to 1.50 kg per sq cm stresses at temperatures near the melting point was found to he a complex reaction controlled by three mechanisms acting in parallel. In order of appearance with decreasing temperature, the auerage activation energies, Qc , are 168, 79, and 24 kcal per mole. Stress dependence of the minimum creep rate was found to he an exponential for the two high-Qc processes and a power law for the low-Qc, process. Transition of control occurs from one mechanism to another. The relative transition temperature depends upon the applied stress, and the range oiler which the transition occurs depends upon the difference in the activation energies of the mechanisms. The creep behavior at high temperature is explained by the climb of dislocations through thermally or mechanically formed jogs, CREEP in pure fcc metals at temperatures in excess of one half the absolute melting point is normally controlled by dislocation c1imb.1,2 A climb model which seems applicable, according to an extensive analysis of data,3 was derived by Weert-man.4 This model assumes jog-saturated dislocations and predicts an activation energy, Qc, nearly equal to that for self-diffusion, USd. Although the requirement of jog saturation is restrictive, agreement between theory and experiment seems adequate. Many other theoretical treatments, including an early model by Mott,5 include detailed consideration of jog formation as an initial requirement for climb. These models predict activation energies for creep which differ from those of self-diffusion. Seeger6,7 postulates an observed activation energy related to the stacking-fault energy. Thus, Usd <Qc<5Usd + Uj where Uj is the jog-formation energy. Seeger incorporated qualitatively the influence of the relative numbers of thermal and athermal jogs. Expanding this concept, Shoeck8 explicitly states a function based on formation mode and relative numbers of vacancies: e r ci exp {-uf/k T} exp {-Um/k T} [ 1 ] where Uf and Um are energies of vacancy formation and migration, respectively. The concentration of jogs, Cj, depends upon the manner of jog production. For intersection jogs, Cj is not sensitive E. R. GILBERT, Junior Member AIME, formerly with De- to temperature; for thermal jogs, Cj is proportional to exp {—Uj /k T}. Schoeck regards each mode as a distinct mechanism; therefore, the mechanisms may act together.299 The diversity predicted by theory, surprisingly, has not been substantiated by experimental results. A significant investigation must include the extremes in stacking-fault energy. Extensive creep studies of aluminum10 and nickel,11,12 high stacking-fault energy metals, have been made. Comparable studies on a low stacking-fault energy metal, such as copper, have not. It is the purpose of this paper to report the results of an investigation of the creep of copper under conditions which favor thermal jogs. EXPERIMENTAL Cylindrical compression creep specimens (0.240 in. in diameter by 0.400 in. long) were machined from cold-drawn rods of electrolytic tough-pitch copper containing 0.0007 Mg, 0.002 Fe, 0.001 Ni, 0.0005 Ag, Cd < 0.005, and Pb < 0.005 wt pct impurities. Undetected spectrographically was a nominal 0.04 wt pct 0, which occurs as a Cu2O constituent distributed discontinuously at grain boundaries. Vacuum annealing at 900°C for 15 min produced a stable 0.03-mm average grain diameter. Testing was carried out using apparatus similar to that described by sherby,13 modified by enclosing platens and a push rod in a vacuum cylinder. Normally this arrangement resulted in pressures less than 10 Only a slight surface tarnish, less than 0.0005 in. in thickness, occurred during the test. The applied stress, corrected for atmospheric pressure, was maintained within 2 pct of the desired true stress by the addition of lead shot at fixed strain increments. Creep strain was measured with dial gages as a relative displacement of the upper and lower platens; accuracy of measurement was 0.0001 in. Two creep-test methods were used, the differential or cyclic temperature14 and the isothermal, to obtain creep data at stress levels of 150, 100, and 75 kg per sq cm over the temperature range of 620° to 1032°C. Minimum creep rates were used from both test methods; this was considered proper because comparable temperature tests or cycling to the same temperature gave the same creep rate, within experimental error. The cold vacuum test chamber, with the unstressed specimen in place, was heated to temperature by placing a preheated furnace over the chamber. Temperature equilibrium was attained within 30 min. For the cyclic tests, the stress was removed during the 5 to 10 min necessary to effect the temperature change and re-
Jan 1, 1965
-
Part II – February 1968 - Papers - Influence of Work-Hardening Exponent on the Fracture Toughness of High-Strength MaterialsBy E. A. Steigerwald, G. L. Hanna
The influence of work-hardening exponent on the variation of fracture toughness with material thickness was studied for high-strength steel, aluminum, and titanium alloys. The results indicate that, when materials are compared at similar fracture toughness to yield strength ratios, the material with the lower work-hardening exponent undergoes the transition from flat to slant fracture at a larger thickness than material with a high work-hardening exponent. In the thickness range where complete slant fracture is obtained the reverse is true and a lower work-hardening exponent results in a lower fracture toughness. The influence of work-hardening exponent on fracture toughness is, therefore, dependent on the particular fracture mode. In the transition region a low work-hardening exponent is beneficial for fracture toughness while in the 100 pct slant region it is detrimental. THROUGH the use of fracture mechanics analyses, the influence of geometric variables on the crack propagation resistance of structures can be interpreted with reasonable consistency. However, in order to gain a more complete understanding of the fracture process, the influence of material parameters on crack propagation must be defined and coupled to the macroscopic fracture mechanics approach. The work-hardening exponent, which characterizes specific material behavior, may serve as an effective parameter to allow some degree of coupling to be accomplished. In the extension of a crack in a specimen of suitable dimensions the propagation process occurs in a stable manner when the crack extension force is balanced by the resistance to crack extension, which exists in the material at the crack tip. As the applied stress, and therefore the crack extension force, on the specimen increases, the resistance also increases primarily because the effective plastic zone at the crack tip, which is the main energy absorption process, becomes larger. Since the work-hardening rate of a material influences the stress-strain relationship, it will also influence the energy absorption process in the plastic enclave at the crack tip and hence should have an effect on crack propagation. A number of studies have been made correlating the strain-hardening exponent or the strain to tensile instability with the ability of a material to resist fracture. Gensamer1 concluded that a low-strain-hardening exponent would result in a steep strain gradient at the base of a notch. He reasoned that a large work-hardening coefficient would result in high-energy ab- sorption due to the increased area under the stress-strain curve. Larson and Nunes2 experimentally observed in Charpy tests on steels heat-treated to below 200,000 psi yield strength that the energy to failure in the fibrous mode, i.e., above the brittle-to-ductile transition temperature, was logarithmically related to the strain-hardening exponent. In order to avoid the complicating effects present in studying materials which undergo a brittle-to-ductile transition, Ripling evaluated the notch sensitivity of a variety of fcc metals with varying work-hardening exponents.3 The results indicated that the relative notch sensitivity, as determined from tests on specimens with a sharp notch, decreased with increasing values of strain hardening. Although the energy required for ductile or fibrous fracture increases with increasing work hardening, high-strength steels often exhibit improved crack propagation resistance when heat-treated to obtain low values of strain hardening.4,5 An analysis of whether low strain hardening is beneficial or detrimental to crack propagation resistance must depend on the particular fracture criterion involved. At temperatures where the material is relatively ductile and the development of a critical strain is required for fracture, high strain hardening increases the energy required to produce failure. In the transition region and below, however, a critical stress law appears to be valid6 and a low rate of work hardening may produce superior resistance to semibrittle crack propagation. The experimental program is aimed at studying these possibilities and determining the specific influence of strain hardening on the crack propagation resistance of several high-strength materials. MATERIALS AND PROCEDURE The alloys, chosen as representative of various classes of high-strength materials, are summarized in Table I. The heat treatments evaluated along with the smooth tensile properties are presented in Table 11. Pin-loaded sheet tensile specimens were employed to determine the smooth tensile properties. A strain gage extensometer (measuring range 0.200 in.) was used at a strain rate of 0.02 in. per in. per min. The work-hardening exponents were determined from the stress-strain curves generated in the smooth tensile tests and the assumption that the portion of the stress-strain curve beyond the yield point can be described by the power relationship: where a is the true stress, P is the true plastic strain,
Jan 1, 1969
-
Part III – March 1969 - Papers- Effects of Substrate Misorientation in Epitaxial GaAsBy A. E. Blakeslee
Morphological and electrical properties of GaAs epitaxial layers are influenced not only by changes in the nominal substrate orientation but also by small amounts of misorientation from the exact crystal planes. Deviations up to 5 deg from {11IA}, {11IB}, and (100) planes were investigated. Growth rates increase progressively with angle, approximately I u per hr per deg. Size and density of growth pyramids fall off with increasing angle, but other effects that are deleterious to the surface may occur which are heightened by increased misorientation. Carrier concentration decreases and electron mobility consequently increases as the angular offset increases, except in the case of strong compensation, where the mobility trend is reversed. It has long been known that changes in the crystallo-graphic orientation of the substrate may cause pronounced effects on the morphological properties of vapor grown semiconductor films. Reports of orienta-tion-dependent growth rates and surface characteristics are as old as the literature on epitaxy itself. shawl has recently published a comprehensive study of the dependence of growth rate on substrate temperature and orientation in epitaxial GaAs. It is also well-known that misorienting the substrate surface a few degrees away from the nominal low-index crystal-lographic plane often produces a much smoother epitaxial surface. This was reported by Tung2 for silicon, Reisman and Berkenblit3 for germanium, and by Kontrimas and Blakeslee4 for GaAs, and use is commonly made of this fact in the semiconductor industry to help guarantee smooth vapor deposits. The effects of substrate orientation on the carrier concentration and mobility of vapor grown GaAs were first documented by williams5 in 1964 and have been observed by several other authors since then,6,7 but no one has yet reported a careful study of how small changes influence these properties. We have made such a study and have found that sizable differences in growth rate, morphology, carrier concentration, and mobility can indeed be observed for epitaxial films grown on substrates that are oriented by progressive small increments away from the exact crystal plane. EXPERIMENTAL Early in the investigation an arsine synthesis system of conventional design8 was employed to produce growths on {111A}-oriented GaAs substrate crystals. In that early work, pronounced effects on carrier concentration and electron mobility were observed as a function of slight misorientation from this low index plane. That observation led to the more careful study that is reported here. An AsC13 system, differing in major aspect from those commonly in use today9 only in that the reactor is vertical rather than horizontal, was used for the detailed study. The gallium source was at 900°C and the substrates were at 750°C. The flow rate of pal-ladium-diffused H2 through the AsCl3 bubbler was 200 cu cm per min, and the flow rate of bypass H2 was also 200 cu cm per min. The substrates consisted of chro-mium-doped semiinsulating GaAs to facilitate elec-trical evaluation of the overgrowth by means of Hall and conductivity measurements on conventional eight-legged Hall bridges. They were misoriented by 0 to 5 deg from the {111A}, {111B}, and (100) planes, toward the (100) from the {111A} and {111B} and randomly toward the <111A> or <111B> from the {loo). The crystals were oriented for sawing by the Laue back-re-flection technique, which is good only to about ±1/2 deg; but after polishing or sometimes after epitaxial growth the wafers were checked by a diffractometer technique which is accurate to about * 0.1 deg. After lapping, the wafers were polished with NaOCl after the technique of Reisman and Rohr,10 and just before use they were cleaned in NaOC1, thoroughly rinsed with de-ionized water, and blown dry with nitrogen. Each run employed four wafers, each misoriented by differing amounts from one of the three major faces, and at least two runs were made for each orientation. The runs were continued long enough to provide at least a 15-µ or thicker layer. SURFACE MORPHOLOGY The appearance of all the films that were grown in a given run always changed from wafer to wafer as a function of increasing misorientation, but not always in the same regular fashion. At least three different trends were observed. These are more easily seen than described, and reference to the series of photo-
Jan 1, 1970
-
Coal - Mechanized Cutting and Face Stripping in the RuhrBy R. R. Estill
THE rank of the Ruhr coal ranges from a high volatile bituminous coal to an anthracite, depending to some extent on the original depth of the seam. The average Ruhr coal corresponds to a soft bituminous American coal of a coking quality. The average thicknesses of individual coal seams being mined are also comparable (59 in. against 65 in. in the United States). However, consideration of seam conditions and mining conditions other than those just mentioned emphasizes differences rather than similarities with United States soft coal. In general, the Ruhr seams now being mined are much more folded and inclined than American seams. Dips of 20' and 30" are common in seams now being worked, and 30 pct of the coal reserves in the district are in seams dipping more than 35". Only on the tops and bottoms of folds do we find rather flat coal seams. In addition to the folding there is extensive displacement by cross faulting plus a certain amount of strike faulting of an overthrust nature, which results locally in doubling or omission of seams. Because of the long history of mining in the Ruhr, nearly all coal lying near the surface has long since been mined out, and we find that the average depth of mining is at present about 2300 ft below the surface. Deep mining, folding, and faulting result in seam conditions requiring a great deal more roof support than one finds in American soft coal mines. In fact only in the anthracite district and the Rocky Mountain and Pacific coal fields do we find somewhat similar conditions. It is easy to say, therefore, that the problem of mechanization of coal cutting and loading in the German mines is quite different from that which we have so effectively met in America with our mobile cutters and loaders, duck bill loaders, and a room and pillar system of mining our drift and slope mines. Partly because of more limited coal reserves, the traditional German mining system is largely the longwall method, which gives an almost complete coal recovery. Backfilling must be extensively practiced to protect the longwall faces, the over and underlying seams and workings, and especially the surface industrialized areas and barge canals. The German engineers have accordingly concentrated their efforts on the design of cutters, loaders, and conveyors suitable to longwall methods rather than room and pillar methods. Undercutters with cutter bars like American models have been in use in the Ruhr since well before World War 11. In 1941 they accounted for 8.5 pct of the production. This percentage, of course, includes coal which was undercut but nevertheless had to be broken down with air hammers or with explosives. The most common of these cutters is the Eickhoff Standard cutter (see fig. 1). This machine does about 95 pct of the undercutting in the Ruhr today, and is available with either compressed air or electrical power and in at least four different sizes. A variation of the cutter is this one with two cutter bars (fig. 2). At the end of 1947 about 200 of these machines and similar cutters were accounting for 13.2 pct of the total production, a production which was, however, only 60 pct of the 1941 production rate, so that the actual cutter tonnage was only up to a small amount over 1941. In 1941 about 3 pct of the production was accounted for by shearing machines making their cut perpendicular to the longwall face. They were similar to those used in the States. These machines are today considered obsolete and now account for only 0.7 pct of the total production. They are located at only a few mines and at present do not seem to have much of a future in the Ruhr. For the future, the Ruhr miner is looking forward to rather extensive mechanization of face work, with two major types of equipment being developed almost simultaneously. On one hand there is the development of cutter loaders for use in relatively hard coal. They represent the further extension of ideas developed after relatively long experience with the Eickhoff cutter. On the other hand there has been since 1942 an intense interest in the Ruhr in the development of face-stripping methods, particularly by the Kohlenhobel (coal plow) and its modification. At the end of 1947 these cutter loaders, Kohlen-hobels and scrapers together were actually accounting for only about 1.4 pct of total production while air hammers still broke 77.1 pct and as much as 1.2 pct was actually broken by hand picks. However,
Jan 1, 1951
-
Coal - Mechanized Cutting and Face Stripping in the RuhrBy R. R. Estill
THE rank of the Ruhr coal ranges from a high volatile bituminous coal to an anthracite, depending to some extent on the original depth of the seam. The average Ruhr coal corresponds to a soft bituminous American coal of a coking quality. The average thicknesses of individual coal seams being mined are also comparable (59 in. against 65 in. in the United States). However, consideration of seam conditions and mining conditions other than those just mentioned emphasizes differences rather than similarities with United States soft coal. In general, the Ruhr seams now being mined are much more folded and inclined than American seams. Dips of 20' and 30" are common in seams now being worked, and 30 pct of the coal reserves in the district are in seams dipping more than 35". Only on the tops and bottoms of folds do we find rather flat coal seams. In addition to the folding there is extensive displacement by cross faulting plus a certain amount of strike faulting of an overthrust nature, which results locally in doubling or omission of seams. Because of the long history of mining in the Ruhr, nearly all coal lying near the surface has long since been mined out, and we find that the average depth of mining is at present about 2300 ft below the surface. Deep mining, folding, and faulting result in seam conditions requiring a great deal more roof support than one finds in American soft coal mines. In fact only in the anthracite district and the Rocky Mountain and Pacific coal fields do we find somewhat similar conditions. It is easy to say, therefore, that the problem of mechanization of coal cutting and loading in the German mines is quite different from that which we have so effectively met in America with our mobile cutters and loaders, duck bill loaders, and a room and pillar system of mining our drift and slope mines. Partly because of more limited coal reserves, the traditional German mining system is largely the longwall method, which gives an almost complete coal recovery. Backfilling must be extensively practiced to protect the longwall faces, the over and underlying seams and workings, and especially the surface industrialized areas and barge canals. The German engineers have accordingly concentrated their efforts on the design of cutters, loaders, and conveyors suitable to longwall methods rather than room and pillar methods. Undercutters with cutter bars like American models have been in use in the Ruhr since well before World War 11. In 1941 they accounted for 8.5 pct of the production. This percentage, of course, includes coal which was undercut but nevertheless had to be broken down with air hammers or with explosives. The most common of these cutters is the Eickhoff Standard cutter (see fig. 1). This machine does about 95 pct of the undercutting in the Ruhr today, and is available with either compressed air or electrical power and in at least four different sizes. A variation of the cutter is this one with two cutter bars (fig. 2). At the end of 1947 about 200 of these machines and similar cutters were accounting for 13.2 pct of the total production, a production which was, however, only 60 pct of the 1941 production rate, so that the actual cutter tonnage was only up to a small amount over 1941. In 1941 about 3 pct of the production was accounted for by shearing machines making their cut perpendicular to the longwall face. They were similar to those used in the States. These machines are today considered obsolete and now account for only 0.7 pct of the total production. They are located at only a few mines and at present do not seem to have much of a future in the Ruhr. For the future, the Ruhr miner is looking forward to rather extensive mechanization of face work, with two major types of equipment being developed almost simultaneously. On one hand there is the development of cutter loaders for use in relatively hard coal. They represent the further extension of ideas developed after relatively long experience with the Eickhoff cutter. On the other hand there has been since 1942 an intense interest in the Ruhr in the development of face-stripping methods, particularly by the Kohlenhobel (coal plow) and its modification. At the end of 1947 these cutter loaders, Kohlen-hobels and scrapers together were actually accounting for only about 1.4 pct of total production while air hammers still broke 77.1 pct and as much as 1.2 pct was actually broken by hand picks. However,
Jan 1, 1951
-
Institute of Metals Division - Plastic Deformation and Diffusionless Phase Changes in Metals-The Gold-Cadmium Beta PhaseBy L. C. Chang, T. A. Read
Diffusionless transformation in Au-Cd single crystals containing about 50 atomic pet Cd was investigated by means of X-ray analysis of the orientation relationships, electrical resistivity measurements, and motion picture studies of the movement of boundaries between the new and old phases during transformation. The nucleation of diffusionless transformation by imperfections and the generation of imperfections by diffusionless transformation were discussed. THAT connections exist between plastic deformation and diffusionless phase changes has long been recognized. Thus it is often possible to produce a diffusionless phase change in a temperature range, above that in which the change occurs spontaneously, by cold-working the initial phase. Certain aspects of the dislocation theory of the plastic deformation of crystalline solids also provide for a rather direct connection between the processes involved in plastic deformation and in diffusionless phase changes. Heidenreich and Shockleyl have pointed out that simple edge dislocations in f.c.c. metals are probably unstable, and that the more probable lattice imperfections, called extended edge dislocations, consist of two half dislocations separated by a distance of the order of magnitude of 100A. The region about two atomic planes thick between the half dislocations because of this stacking fault may be described as having the hexagonal close-packed structure. Presumably the stacking faults observed by Barrett" fter cold-working f.c.c. Cu-Si alloys resulted from the passage of such half dislocations through the lattice of the initial phase. It is now becoming clear that the development of a detailed theory of the atomic movements involved in diffusionless phase changes will require a consideration of the role played by lattice imperfections, just as such considerations are necessary to the understanding of plastic deformation mechanisms. This point of view has been recently set forth, for example, by Cohen, Machlin, and Paranjpe3 who pointed out the role which might be played by screw dislocations in nucleating diffusionless phase changes. The present paper reports on some aspects of the diffusionless phase change in single crystals of the beta phase alloy Au-Cd which serve to emphasize further the importance of lattice imperfections in diffusionless phase changes. The diffusionless phase change of Au-Cd possesses several remarkable features. One of these is that the interface between the high-temperature beta phase and the low-temperature orthorhombic phase typically moves with a low velocity, in contrast to the behavior observed in the transformation of austenite to martensite. Motion pictures of this slow interface motion have been prepared in the course of the work reported here. Another important feature of the Au-Cd transformation is the small amount of undercooling observed. The reverse transformation occurs on reheating to a temperature only 20" higher than the transformation temperature observed on cooling, and under some circumstances the hysteresis observed is substantially less than this. This narrow temperature range between transformation on heating and cooling is presumably in part a consequence of the fact that the transformation requires a lattice shear of only about 3". Finally, the orthorhombic product phase possesses unusual mechanical properties, as was first pointed out by olander' and Benedicks." After completion of the transformation on cooling the specimen can be severely deformed, yet on the release of load it springs back to its original shape in a rubber-like manner. Explanation of this phenomenon will require an understanding of the lattice imperfections in the orthorhombic structure and, correspondingly, of those in the initial body-centered cubic structure. Single crystals of Au-Cd alloy containing 47.5 and 49.0 atomic pct Cd were prepared from fine gold (99.95 pct purity) and chemically pure cadmium (99.99 pct purity) by melting the alloy in an evacuated and sealed fused quartz tubing and growing into single-crystal form by the Bridgman method. The Au-Cd alloy containing 47.5 atomic pct Cd undergoes a diffusionless transformation from an ordered body-centered cubic structure to an orthorhombic structure when it is cooled to about 60°C, while the reverse transformation takes place when the alloy is heated to about 80°C, according to electrical resistivity studies. The structures of these two phases have been studied by Blander,4 reinvestigated by Bystrom and Almin.e he lines of Debye photo-gram of powdered samples of Au-Cd alloy containing 47.5 atomic pct Cd prepared in this laboratory were identified and agreed fairly well with those of
Jan 1, 1952
-
Technical Notes - Effect of Feed Injection Position on Hydrocyclone PerformanceBy J. M. W. Mackenzie, C. J. Wood
In attempting to describe the size classification performance of a hydrocyclone, most workers have elected to use either an equilibrium orbit theory or an non-equilibrium orbit theory. The equilibrium orbit theory has been used by the majority of workers including Lilge,' Bradley; and Yoshioka and Hotta. In applying this theory, it is argued that particles in the body of a hydrocyclone attain an equilibrium radial position where the drag force on the particle resulting from the inward radial fluid velocity is balanced by the outward centrifugal force caused by the tangential component of fluid flow. When considered over the full height of the hydrocy-clone, attainment of this radial equilibrium orbit results in the particle following a conical equilibrium envelope. It is then argued that if this envelope lies outside the envelope of "zero vertical velocity," the particle will report to the underflow, while if the equilibrium envelope lies inside the envelope of "zero vertical velocity," the particle will report to the overflow or vortex finder product. The d50-sized particle which reports in equal quantities to the underflow and overflow is assumed to correspond to particles whose equilibrium envelope is coincident with the envelope of "zero vertical velocity." In considering the equilibrium orbit theory, it is apparent that the horizontal position of the particles in the feed inlet pipe should have no effect on their ultimate destination on the hydrocyclone. Each particle should attain an equilibrium position which depends on the density, size, and shape of the particle; the density and viscosity of the fluid; and the flow patterns within the hydrocy-clone. The nonequilibrium orbit or unsteady state theory has been largely developed by Rietema4 and Mizrahi.6 Mizrahi has listed four main objections to the equilibrium orbit theory. These objections center on the short residence time in the hydrocyclone, the fact that the experimental classification curve is much less sharp than is theoretically predicted, and the absence of negative efficiency conditions in hydro cyclones operating on a feed material which is much finer in size than d50. Proponents of the nonequilibrium orbit theory argue that for a particle to discharge with the underflow it must have sufficient outward radial velocity to reach the downward-flowing region close to the hydrocyclone wall in which the flow lines are parallel to the wall and the ratio of vertical to radial velocity is constant. It is then postulated that a d50 particle entering the cyclone at the center of the feed inlet will just reach this downward-flowing region as it reaches the apex. Thus for uniform distribution of particles across the feed inlet, half the d50 particles—that is, those injected in the half of the inlet area nearest the cyclone wall —will report to the underflow while those injected in the other half will not reach the downward-flowing region and will be carried inward to the center of the cyclone and thus report in the overflow. The exact thickness of the down-ward-flowing region of fluid adjacent to the outer wall of the hydrocyclone is uncertain but Mizrahi considers it to be equal to the apex radius minus the air core radius. Particles larger than d50 have a greater outward centrifugal force acting on them than the d50 particles and may reach the wall even if injected at a distance from the wall greater than Di/2 (Di is inlet diameter). Conversely, particles smaller than d50 may not reach the wall even if injected at a distance less than Di/2 from the cyclone wall. Since the equations put forward by the proponents of both theories yield approximately the same values of d50, it is not possible to decide between these theories by measurement of d50. It should be possible however to examine the theories by injecting a small stream of solids into the feed inlet of a hydrocyclone running on clear water. If the efficiency or classification curve is measured for various horizontal injection positions, then the curves should be coincident if the equilibrium orbit theory holds. If, however, the unsteady state theory describes the cyclone operation, then the classification curves should show finer d50 sizes for particles injected close to the cyclone wall. Experimental A 6-in.-diam hydrocyclone with geometry as in Figs. 1 and 2 was used. Quartz particles were injected as a 50% by wt pulp via an 1/8-in. steel probe. For each in-
Jan 1, 1971
-
Metal Mining - A New Incline in the Metaline DistrictBy Chas. A. R. Lambly
In the extreme northeast corner of the State of Washington, on the Canadian border, lies the Metaline mining district. This district is old in history, but young in production. Geology The Metaline district is a zinc-lead area of the replacement type in dolomite and limestone. The ore bodies of the Josephine horizon are in many ways similar to the ore bodies of the famous Tri-State zinc fields. The beds are faulted and folded and have varying low dips in varying directions, and underlie large areas of the district. History Production started in 1927 on a very limited basis. The property is now mining and milling 700 tons per day. The mine is opened by adit tunnels and a vertical shaft. As the ore horizons gained depth, it was necessary to sink inclines to follow the ore horizon (see Fig 1). From 1927 to date, approximately 600,000 ft of diamond drill was put down This work indicated that suficient tonnage existed to justify a redesigning of the whole operation, surface and underground. After four years of general study, the following program was planned: 1. A new mine entrance, which would be an incline, that could follow the ore body down at whatever pitch was necessary. The incline will be equipped with conveyors for the moving of ore and waste to the surface and with tractor-type locomotives for man and supply transportation. 2. The new incline also required a new type of mining which was developed and is now in use. It is called contour mining and will be described in a future paper. 3. The new incline exit would necessitate the moving of the mill and mine shops across the Pend Oreille River. This part of the program is now underway. The Incline The sinking of the incline was to start as soon as World War II ended and was as follows: The first leg of the incline was to be sunk from the surface 1600 ft on a 17" slope. The collar and first level at elevation 2180 ft, the second level at elevation 2000 ft, the third level at elevation 1875 ft, and the fourth level at elevation 1700 feet. From the 1700 ft elevation the incline was to flatten out to 12" for 400 ft to give the necessary depth for the ore pockets below the 1700 ft level and the necessary clearance for future sinking (see Fig 1 and 2). Due to lack of manpower in 1946, the program was changed and was as follows: A drift was driven from the old mine workings on the 1700 ft elevation in an easterly direction. At 1300 ft the drift was turned N 50" E and at this point a raise was driven 180 ft on a 50" slope. This raise intersected the Josephine horizon and commercial ore was encountered. At the 2000 ft mark, a main raise was driven, 245 ft on a 50" slope, and the 1875 ft elevation was cut. Exploration drifts were started on this level and production followed on a limited basis. The main drift at the 2500 ft point was turned N 35" E and ran parallel to and 10 ft east of and under the proposed incline line. At the proposed intersection of the drift and incline on the 1700 ft elevation, it was planned to raise the incline to intersect the 245 ft raise and to continue on to the surface, a distance of 1600 ft. When this proposed intersection point was reached, a heavy flow of water, approximately 800 gpm, was encountered and all work on the main drift face was stopped. This water flow flooded the main pump station in the old mine and the two lower levels with approximately 20,000,000 gal of water. The water was controlled and finally drained from the cave areas and lower levels after six months of pumping. After the heavy flow of water was encountered in the main heading, it was decided that the incline would have to be started from the surface, as originally planned, so that too much time would not be lost. The surface overburden had to be removed, a total of 6000 yards. A temporary dry house for 6 men was built. An 8 in. churn drill hole was intersected in the first raise driven from the 1700 foot elevation tunnel. Air and water lines were placed in this hole, and air and water were delivered to the collar of the incline from the mine working. The incline started down at 15 ft wide and 7 ft high through the Leadbetter slates. After sinking 4 sets, it was
Jan 1, 1950
-
Iron and Steel Division - Thermal Conductivity Method for Analysis of Hydrogen in Steel (Discussion page 1551)By J. Chipman, N. J. Grant, B. M. Shields
The vacuum tin-fusion method of analysis for hydrogen, developed by Carney, Chipman, and Grant, has been modified to permit the analysis of the evolved gases for hydrogen by means of a thermal conductivity cell. A properly prepared sample can be analyzed in 10 min with a probable error of ±0.12 ppm. A study of various methods for storage of hydrogen samples shows that samples can be safely held in a dry ice-acetone bath as long as six days. Storage in liquid nitrogen is necessary for samples to be held one week or more. HE vacuum tin-fusion method, as developed by I- Carney, Chipman and Grant,' is the only analytical procedure which has shown promise of being fast enough for use in the control of hydrogen during steelmaking. It was felt that further simplification and faster speed of operation could be effected by the use of thermal conductivity measurements for analysis of the gases evolved in the tin-fusion method. The application of conductivity measurements to the tin-fusion method is possible because: 1—the evolved gas is essentially a mixture of hydrogen, nitrogen and carbon monoxide with a hydrogen content usually over 50 pct, 2—the evolved gas is collected at a relatively low pressure, and 3— the thermal conductivities of CO and N2 are practically identical while that of hydrogen is very much greater. The major part of this research program was devoted to the construction and calibration of a vacuum tin-fusion apparatus which analyzes the evolved gases for hydrogen by means of a thermal conductivity cell. The second phase of the problem was associated with the development of a procedure for storage of samples prior to analysis. With the rapid quenching method for hydrogen sampling,' which seems to be the most practical for steel mill use, it is necessary that the samples be stored safely during the interval between sampling and analysis if the hydrogen content of the molten metal is to be maintained in the supersaturated solid samples. The thermal conductivity bridge has been used for a number of years in the analysis of certain gas mixtures. An elementary discussion of the theory and practice of gas analysis by thermal conductivity measurements is given by Minter.3 A more comprehensive discussion of the theory and of the various measuring circuits is presented by Daynes.' A complete knowledge of the theory and properties of the thermal conductivity of gases and gaseous mixtures can be gained by a study of the standard textbooks on the kinetic theory of gases."' The existing data on the thermal conductivity of single gases are reviewed by Hawkins: that for a number of binary gas mixtures by Daynes' and Lindsay." The thermal conductivity method may be applied to the determination of the composition of a binary mixture if: 1—the thermal conductivity of the mixture varies monotonically with composition, and 2— the two gases have measurably different thermal conductivities. The greater the difference between the two gases, the greater the sensitivity of the method.10 he method is applicable to the analysis of multicomponent mixtures when all of the gases in the mixture except one have nearly the same thermal conductivity. Fortunately, the mixture of hydrogen, nitrogen, and carbon monoxide evolved by the tin-fusion analysis' falls in this latter classification. The thermal conductivities of nitrogen and carbon monoxide are practically equal; and the thermal conductivity of hydrogen is approximately seven times that of the other two. Therefore, the thermal conductivity of a gaseous mixture of hydrogen, nitrogen, and carbon monoxide at known temperature and pressure can be related directly to the percentage of hydrogen in the mixture by suitable calibration. Usually the thermal conductivity of a mixture of gases is measured at atmospheric pressure where the thermal conductivity is independent of pressure over a wide pressure range. At very low pressures (below 1 mm Hg), the thermal conductivity of gases varies with the pressure. This phenomenon has been utilized in the Pirani vacuum gage for the measurement of pressures in the range of 10" to 10-0 mm of mercury.= Very little has been published concerning the variation of thermal conductivity with pressure at intermediate pressures between 1 mm Hg and 1 atm. However, preliminary measurements indicated that the thermal conductivities did vary with pressure over the range of pressures (up to 10 mm Hg) at which gases are delivered from the vacuum pump. Therefore, the calibration of the thermal conductivity cell had to be planned to include the effects of both gas composition and pressure. Such a calibration chart is shown in Fig. 4. Most industrial applications of the thermal conductivity method of gas analysis have used a compensated Wheatstone bridge circuit containing two
Jan 1, 1954
-
Part III – March 1968 - Papers - Formation of Phosphosilicate Glass Films on Silicon DioxideBy J. M. Eldridge, P. Balk
Phosphosilicate glass films were formed, by reacting gaseous P2O5 with SiO2, over a large range of temperature (800° to 1200°C) and gas phase composition (nearly two orders of magnitude of effective P2Ospressure). The film compositions generally corresponded with the liquidus curve, delineating the maximum solubility of the tridymite Phase of SiO 2 in phosphosilicate liquid solution at the temperature of film formation. It is shown that the P2O5 concentration of the phosphosilicate liquid film tends to decrease by reaction with the underlying SiO 2 layer until the liquidus curve is reached. The validity of the thermodynamic argument used in this explanation is supported by the results of a determination of the composition of borosili-cute films, prepared by reacting gaseous B2O3 with SiO2 at different temperatures. The kinetics of phosphosilicate film formation were described by a model predicated on a steady-state diffusion of P2O5 through the film. UNDERSTANDING of the processes leading to formation of phosphosilicate and borosilicate glasses is of great importance for producing passivating layers on FET devices. Passivating films with optimum characteristics are preferably formed in a separate step, independent of the doping of the semiconductor.' The results of an investigation carried out to gain improved insight into the mechanism of glass formation are presented in this paper. It would be expected that application of the known Pz05-Si02 and B 2 O 3-SiO2 phase diagrams should be useful in extending understanding of the glass-forming processes. However, the question of the propriety of treating thermally grown SiO2 in these binary oxide systems by the methods of equilibrium thermodynamics must be considered when this application is attempted. Although Sah et a1.' and Allen et al. 3 investigated the kinetics of formation of phosphosilicate glass (PSG), they failed to adequately relate their diffusion models to the occurrence of experimentally observed phases in the PSG/SiO 2/Si system. Horuichi and yamaguchi4 investigated the diffusion of boron through an oxide layer and described their results in terms of a model similar to that of Sah and coworkers. More recently, Kooi 5 and Snow and Deal6 reported the compositions of PSG films formed by depositing P2 O 5 onto SiO2. These compositions apparently coincide with those at the liquidus curve delineating the maximum solubility of crystalline SiO2 in phosphosilicate liquid solutions. These authors did not discuss the thermodynamic implications of their results on the structure of thermally grown SiO2 films. The structure of thermally grown Sio2 films and that of vitreous silica are generally thought to be quite similar. Since the solubility of a substance depends on its structure, it is relevant that the solubility of vitreous silica in water7 is highly reproducible, like the solubility of thermally grown SiOz in phosphosilicate liquid. Furthermore, the vitreous silica-water system appears to be in true thermodynamic equilibrium (viz., the same solubility value can be approached from both supersaturated and under-saturated solutions). Sosman7 suggested that a type of two-dimensional lattice may form at the silica/solution interface, resulting in the observed solubility behavior that is characteristic of a crystalline solid. An alternative explanation may be that vitreous silica has a microcrystalline grain structure. Other investigators have suggested that vitreous silica has essentially the structure of B cristobalite,' or is composed of microcrystals of p tridymite or cristobalite, or a mixture of both. Presumably the grain size would be sufficiently large to minimize any appreciable contribution of the grain boundaries to the solubility of the crystalline matrix. The present investigation was carried out to clarify the significance of the boundaries in the Pa,-SiO, and B2O3-SiO2 Systems in determining PSG and BSG (borosilicate) film compositions. Furthermore, the kinetic data for PSG film formation were extended, using a wider range of formation parameters than were previously reported. One model describing the kinetics of film formation will be presented that is compatible with the thermodynamics of the Pa5-Si02 system. EXPERIMENTAL PROCEDURE Glass Film Preparation. SiO2 films (1000 to 8000A thick) were obtained by oxidation of silicon substrates in dry O2 at 1100°C. PSG and BSG films were prepared by exposing these layers to gaseous oxides obtained by reacting high-purity POC13 and BBr3, respectively, with O2. A double-columned saturator was used to ensure complete saturation of the N 2 carrier
Jan 1, 1969
-
Reservoir Engineering-General - Effect of Vertical Fractures on Reservoir Behavior-Results on Oil and Gas FlowBy J. S. Levine, M. Prats
A homogeneous and uniform cylindrical reservoir containing oil and gas is fractured vertically on completion and is produced at a constant bottom-hole pressure. The fracture has an infinite flow capacity, is of limited lateral extent and is bounded above and below by the impermeable strata defining the vertical extent of the reservoir. Results show that such a fractured reservoir can be represented by a reservoir of circular symmetry having very nearly the same production history. The well radius of this circular reservoir is about 1/4 the fracture length and is essentially the same as that obtained previously for a single fluid of constant compressibility. At the same value of cumulative oil production, gas-oil ratios of fractured reservoirs producing at constant terzinal pressure are larger than those of reservoirs having no fractures. This leads to more inefficient use of the reservoir energy in fractured wells and results in lower reservoir pressures for the same cumulative oil production. The reduction in operating life due to fracturing a reservoir is not as great as that for a slightly compressible fluid. This diflerence can be accounted for by the lower reservoir pressure in the fractured reservoir and its adverse effect on the average mobility and compressibility of the oil. As anticipated, the reduction in operating life increases czs the reservoir permeability decreases. The type of results presented in this report can be used to determine the economic attractiveness of fracture treatments per se, to setect the initial spacing to be used in developing a field, and to compare the relative merits of fracturing available wells and infill drilling. INTRODUCTION The effect of vertical fractures on a reservoir producing either an incompressible or a compressible liquid has already been discussed in the 1iterature.l,2 Those results indicate that the production history of such a reservoir is essentially the same as that of a circular reservoir having an effective well radius of approximately one-fourth the fracture length. The present work reports on the effect of a vertical fracture on a reservoir producing two compressible fluids —oil and gas—by solution gas drive. Because of the empirical nature of the PVT and relative permeability data used to obtain the performance of such reservoirs, results can only be obtained numerically and with the aid of high-speed computers. Since reservoirs lose their radial symmetry when fractured vertically, pressure and saturation can no longer be given only in terms of distance from the well. Two coordinates (such as x and y) must now be used to describe the pressure and saturation within the reservoir, and, since we are dealing with compressible fluids, time is also a variable. Thus the solution of a vertically fractured reservoir requires finding two unknowns (pressure and saturation) in two space variables (say x and y) and in time (t). Since no means are readily and generally available for solving such problems at the present time, we have used the results of previous work1,2 to approximate the effect of a vertical fracture on a reservoir producing both oil and gas by depletion. The purpose of the present wmk, then, is to investigate the possibility of using available numerical techniques (limited at the moment to one space variable) to study the two-space-variable flow behavior resulting from a vertical fracture. Results obtained in the course of this investigation are also reported and discussed. Input and output data of the numerical methods used are given in practical units: BOPD, feet, psi, cp, and md. Results are discussed fist in terms of specific reservoir and crude properties and geometries. Later, dimensionless parameters are introduced in order to extend results to different values of some of the reservoir and fracture properties. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The reservoir has an impermeable, circular, cylindrical outer boundary of radius r,. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below' the producing layer (reservoir). It is assumed that there is no pressure drop in the fracture due to fluid flow. 1 indicates the general three-dimensional geometry of the fractured reservoir. Gravity effects and the effects of differential depletion resulting from variations in hydrostatic head (pressure) will be neglected. Thus, the flow behavior in the fractured reservoir is described by the
-
Part VIII – August 1969 – Papers - The Undercooling of Cu-20 Wt Pct Ag AlloyBy G. L. F. Powell
g samples of Cu-20 wt pct Ag alloy have been mdercooled to a maximum of 197°C by melting under a slag of commercial soda-lime glass in a vitreous silica crucible. No grain refinement of the primary copper was observed in samples undercooled to the maximum of 197°C. When the samples contained a small amount of oxygen, the copper dendrites were partially recrystallized at undercoolings greater than 97°C. In previous papers'-3 reporting the grain structure of undercooled silver and copper, it was observed that grain refinement was dependent on both undercooling and oxygen content. Grain refinement occurred in undercooled silver when the degree of undercooling exceeded the range 153" to 175"C, while in Ag-0 alloys (0.12 wt pct) fine equiaxed grains were exhibited when undercooling was greater than 50°C. Similarly, copper samples undercooled as much as 208°C displayed fan-shaped growth from a single nucleation site, while the grain structure of Cu-O alloys (0.08 wt pct) was fine and equiaxed at undercoolings larger than 150°C. Thus the presence of oxygen greatly reduced the undercooling at which grain refinement occurred. It was also observed that the change in grain size resulted from recrystallization and was not due to an enhanced nucleation rate in the liquid-solid transformation. It is possible that the influence of oxygen on recrystallization is due primarily to its presence as a solute element. walker4,' reported that, although a grain size change did not occur in pure nickel until the undercooling exceeded 150°C, small grains were observed in samples of Ni-Cu alloy solidified at small and large degrees of undercooling. Jackson et al.6 suggested that the fine grained structure of the Ni-Cu alloy resulted from the melting off of dendrite arms during recalescence. This remelting process may occur in alloys as a result of segregation during freezing which causes a variation in liquidus temperature from point to point within a dendrite. It was therefore decided to undercool copper with a metallic alloying element to ascertain whether the presence of a metallic solute would have a similar effect to oxygen in inducing grain refinement. A Cu-Ag alloy was chosen, since both metals had been shown to behave similarly on undercooling. The alloy Cu-20 wt pct Ag was selected since the eutectic constituent outlines the initial growth form of the primary copper, so that the as-frozen grain structure is not obscured if subsequent recrystallization occurs. This paper describes the results of undercooling experiments carried out with Cu-20 pct Ag samples undercooled to a maximum of 197°C and the effect of oxygen content on the grain structure of the undercooled samples. EXPERIMENTAL Melting was carried out in a small cylindrical resistance furnace using "fine" silver granulate and oxygen-free high conductivity copper. The procedure adopted was to melt the required quantity of silver in air in a clean vitreous silica crucible for approximately 15 min, freeze, and add granulated commercial soda-lime glass to form a complete surface slag cover, after which the sample was melted and frozen several times to reduce the oxygen content. The glass slag cover was approximately 3 in. thick. Pieces of copper (=50 g) were added to the crucible until the required quantity to make 350-g samples of alloy had been charged. Each piece was added quickly to the crucible which was held at a temperature slightly above the melting point of silver. The piece was quickly pushed beneath the glass to minimize oxidation and any oxide coating usually decomposed before the piece had settled down into the silver. After the full quantity of copper had been added, the melt was stirred with a silica rod to hasten homogenization and a Pt/Pt 13 pct Rh thermocouple enclosed in a vitreous silica sheath inserted for temperature measurement. Heating and cooling curves were recorded on a potentiometric chart recorder fitted with a zero suppression unit. The milli-voltage range of the recorder was adjusted so that temperatures could be read to 1°C. Heating and cooling curves were taken every hour until three consecutive readings gave the same solidus-liquidus range, consistent with the solidus-liquidus range for this alloy composition by reference to Hansen and Anderko.7 Metallographic examination of samples frozen at this stage, failed to show any variation in composition from bottom to top of the ingot. Consequently, it was considered that the melt was homogeneous at this stage and undercooling experiments were then car-
Jan 1, 1970
-
Part VII - Papers - Fatigue Crack Nucleation in a High-Strength Low-Alloy SteelBy Raymond C. Boettner
The present work had for its purpose: 1) the identification of crack nucleation sites in AISI 4340, quenched to martensite and tempered over a range of 'temperatures; and 2) the comparison of fatigue processes in AISI 4340 with processes observed previously in pure metals From constant def1ection-bending fatigue tests, martensite boundaries were identified as the favored crack nucleation sites in quenched and tempered AISI 4340. It, also, was concluded that the fatigue processes operating- in this lous-alloy steel were similar to Processes observed in pure tnetals. ALTHOUGH much engineering data has been accumulated on the fatigue properties of quenched and tempered martensitic steels,' fatigue as a process is not as well understood in martensite as it is in pure metals.' Important features of the fatigue process, such as the identity of the nucleation sites, have not been determined in the commercially important high-strength low-alloy structural steels. The present work had for its purpose: 1) the identification of crack nucleation sites in a low-alloy steel, i.e., AISI 4340, which had been quenched to martensite and tempered over a range of temperatures; and 2) the comparison of fatigue processes in the AISI 4340 with processes observed previously in pure metals. This comparison of the fatigue processes in the different tempers was restricted to the high-strain low-cycle part of the S-N curve. Under these test conditions, previous work on a number of metals has shown that a large number of cracks are nucleated in less than 30 pct of the fatigue life.3 Furthermore, crack nucleation sites are not restricted to inclusions but are also associated with intrinsic structural characteristics of the metal. MATERIAL A 20-lb ingot of vacuum-melted AISI 4340 (for composition see Table I) was hot-rolled to 1-in.-diam rod and then cold-rolled to a 1-in.-wide strip, 0.08 in. thick. Fatigue specimens, see insert of Fig. 1, were machined from the strip with the long dimension parallel to the rolling direction. m this orientation, the stringers of 1 to 2 p inclusions present in the sheet lay parallel to the stress axis in the specimens. The specimens were austenitited at 2050°F in order to obtain a large prior austenite grain size, i.e., 2 mm, which facilitated the subsequent identification of the prior austenite boundaries. A helium atmosphere was used to minimize decarburization. After austenitiza-tion at 2050°F, the specimens were transferred to a 1450°F furnace so that specimen distortion was held to a minimum in the subsequent oil quench. Previous work4 indicated that refrigeration in liquid nitrogen prior to tempering reduced the percentage of retained austenite in the quenched specimens to less than 5 pct. Tempering was carried out in air over the temperature interval of 200°to 800°F to produce a range of mechanical properties, Table I. The preparation of the fatigue specimen was completed by grinding about 0.005 in. from each surface and electropolishing in a chrome trioxide-acetic acid solution for 30 min. Examination of etched cross sections of specimens prepared in this fashion showed the foregoing specimen preparation to be adequate for the removal of the decarburized layer present after the heat treatment. Transmission electron microscopy showed that the as-quenched microstructure of this alloy consisted of a mixture of martensite plates containing either a high density of dislocations or microtwins. Previous work5'6 indicated that in the course of oil quenching autotem-pering resulted in the formation of E carbide on the martensite and microtwin boundaries. Tempering for 2 hr at temperatures up to about 400°F resulted in further precipitation of the E carbide. Finally, at about 400°F, cementite began to replace the E carbide on the martensite and microtwin boundaries in addition to forming a Widmanstatten structure within the plate matrix. EXPERIMENTAL S-N curves were obtained using electropolished specimens cycled at 1800 cpm as cantilever beams in fully reversed bending at selected constant deflections. The deflections were translated into surface strains by means of a calibration curve obtained through the use of strain gages. An argon atmosphere was used to minimize environmental effects. To investigate the development of fatigue slip bands, the specimens of the different tempers were unidirec-tionally bent to a surface strain of 0.005 to 0.007, photographed to record the location and appearance of slip bands so introduced, and then cycled to failure
Jan 1, 1968
-
Institute of Metals Division - Determination of Boundary Stresses during the Compression of Cylindrical Powder CompactBy M. E. Shank, J. Wulff
At the present time, the designer of dies for metal powder pressing is handicapped by relative ignorance of stress distribution and frictional effects at the interior surface of the die. Unckell was the first to develop a method for the study of wall friction. He used three Brinell balls on which the die rested during pressing. The total frictional wall force was determined by the size of impression these balls left on a soft metal plate. Since the method does not give radial pressures, or distribution of such pressures, coefficients of friction could not be determined. Although Unckel measured density distribution, he was not able to determine radial or shear stresses. Shaler2 has proposed theoretical expressions for the stress and density distribution within cylindrical compacts during pressing, in accordance with the experimental results of Kamm, Steinberg, and Wulff.3 By application of Siebel's method,4 Kamm et a13 plotted stress trajectories for two compacts. From the stress trajectories they calculated coefficients of friction from point to point along the die wall. As pointed out by Shaler in the discussion of Ref. 3, these values are based on progressive point-to-point calculations on finite size grid squares across the compact. In the region of the die wall such calculated values may therefore have considerable cumulative error. The purpose of the present paper is to develop an experimental method by which the nonhydrostatic pressures and shears acting on the interior wall of a cylindrical die can be measured. Such measurements can then he correlated with existing data to aid in the explanation of the pressing process. The method used is based on the elastic: properties of the thick-walled tube used as the die. The principle of super-position of force systems on an elastic body is assumed to hold. Electric strain gauges were mounted in adjacent positions on the exterior die wall in order to get an exact measurement of the variation of tangential strain over the length of the die during pressing. While in this paper, measurements in terms of only tangential strains are considered, it is well to note that similar calculations may be set up for axial strains. The latter are not preferred, since they tend to be smaller than the tangential strains and therefore permit less sensitive measurements. Discussion in this work is restricted to compacts pressed from both ends, since the elastic deformation of the die is then more amenable to analysis. Before choosing the electric strain gauge method, a more direct line of attack was considered and discarded. The discarded idea was the insertion of a pressure gauge through a hole in the die wall.* The gauge would have been in the form of a small piston. If pressure were exerted against such a gauge, it would move outward along a radius of the die. One disadvantage of the scheme is its inability to measure shears along the die wall. Another more serious disadvantage is the disturbance caused by the device itself. It would serve to change the forces it was designed to measure. No matter how small the movement of the gauge, when pressure is applied a discontinuity would exist in the wall surface at that point. Due to the stress concentration caused by the hole, abnormal deflections of the die wall would occur around the gauge. During pressing, powder would be forced into the resulting depression. The depression would then become larger with increasing compacting pressure. Powder, not being a fluid, is capable of supporting shear. The ease with which it would flow into the die wall depression to further move the piston is an indication, not of the radial pressure at that point, but of the state of shear retarding the movement. Thus the "pressure" gauge is really a criterion of flowability, and of the capability of the powder to support shear. For these reasons, it was decided that the electric strain method, herein employed, was more reliable, if more indirect. The gauges and lead wires, mounted on the external die wall do not in any way affect the behavior of the metal powder or the die during pressing. Theory of the Method THE EFFECT OF RADIAL PRESSURE ON THE DIE WALL Effect of a Single Small Band of Hydrostatic Pressure Consider a die which is a thick-walled cylinder of outer radius R. and inner radius Ri. If over a small finite length L there is a normal pressure P, a tangential strain distribution at the outer wall results. This is shown schematically in Fig 1. The exact shape of the curve may he predicted by an extension of the theory of a semi-infinite beam on an elastic foundation.6 This
Jan 1, 1950
-
Mining the San Juan Orebody El Mochito Mine, Honduras, Central AmericaBy Robert C. Paddock
INTRODUCTION A way of producing 3,000 tpd from the El Mochito Mine was needed. Of this production, 2,000 tpd must come from the San Juan orebody. The original sub-level stoping method did not give satisfactory results due to ground instability, and the highly irregular ore/waste contacts encountered . The experience gained from the initial system helped guide research into the ground instability problem. Results from this work, combined with knowledge gained about the orebcdy configuration, defined constraints that were previously not fully appreciated. These constraints, and others, combined with objectives, were considered together to develop a new mining method. No single technique was found to be suitable, so a hybrid mining system was developed. A combination of ramping, cut and fill, and vertical crater retreat, with an option to use top heading and benching was developed. To complement the mining system, the type of equipment needed was decided upoun. Also, to support the mining system at this expanded rate of product ion, major modifications of existing infrastructure were required. THE EL MOCHITO MINE The El Mochito Mine, of Rosario Resources Corporation, has been in continuous product ion since 198. The mine began operations in April of that yeas at a rate of 100 tpd. The reserves in 198 were 100,000 tons of silver ore assayed at 1,250 grams per tonne. As of the end of 1979, the El Mochito orebodies have produced over 5.6 million tonnes of ore averaging 516 grams per tonne silver, 6.8 lead, and 7.8% zinc. Present ore reserves are about 7.9 million tonnes, averaging 138 grams per tonne silver, 4.6% lead, and 8.7% zinc, with minor quantities of copper, cadmium and gold. An expansion plan to increase mill production two fold to 2,500 tonnes per day is underway. This expansion will require the mine to produce 3,000 tpd. The mine consists of numerous orebodies, all of which have been mined to a certain extent. Of all the orebodies, the San Juan contains 8% of known reserves. This amounts to about 6.7 million tonnes. The significance of the San Juan orebody to the future life of the El Mochito Mine is obvious. If the required mine production of 3,000 tpd is to be sustained, the San Juan must be the source of the majority of that production. Due to the mineability and overall logistics concerned with the other orebodies, the San Juan must be able to reach and maintain a production rate of 2,000 tpd by 1982. GEOLOGY OF THE SAN JUAN OREBODY The El Mochito Mine is a classic example of a chimney replacement deposit in limestone. Similar deposits axe found in Mexico, at the Naica, Providencia, and Santa Eulia Mines. The El Mochito Mine is located at the south- western end of the Sula Valley on the western edge of the Honduras Depression in the Central Cordillera and Central Highlands of Honduras in a setting of Mesozoic sediments. The orebodies occur in a structural basin developed between NNE trending normal faults and apparently hinged on the south end. Topographically, the Mochito Basin lies between the uplifted Santa Barbara mountain in the west and the Palmer Ridge on the east. The San Juan orebody occurs near the intersection of the NE trending San Juan fault and the ENE trending Porvenir fault. The downward continuation of the orebody is controlled by the westward rake of these NW and N dipping structures. The discovery of the San Juan orebody is attributed to analysis of structural evidence of known ore deposits by in-company geologists. The composition of the San Juan orebody is primarily garnet skarn, with local concentrations of hedenbergite and magnetite. The economically important sulfide mineralization consists of (in decreasing abundance), sphalerite , galena, pyrrhotite , and chalcopyrite. There is some indication that a Cu-Ag mineral such as tetrahedrite may also be present. The skarns were formed by replacement of the original limestone by hydrothermal water migrating upward roughly along the intersection between the Porvenir fault system and the San Juan fault system. Textural evidence suggests that the orebody is a composite of several pulses of hydrothermal activity which would explain, in pat, the great irregularity of the contacts and the large horizontal variation in mineralogy. A general pattern of skarn types can be seen in the orebody, partially accounting for the observed lateral variation in grades. This zonation is very generalized, and one or more zones may be missing in any given locality. The orebody is almost invaxiably surrounded by a 2 cm to 25 cm zone of bustamite skaxn with low values. The border skarn is usually
Jan 1, 1981
-
Part XI – November 1968 - Papers - The Effect of Dispersed Hard Particles on the High-Strain Fatigue Behavior of Nickel at Room TemperatureBy G. R. Leverant, C. P. Sullivan
To evaluate the effect of a dispersion of nondeform-able, incoherent, second-phase particles on high-strain cyclic deformation and fracture, recrystallized TD-nickel (Ni-2ThO2) and a commercially pure nickel, Ni-200, were fatigued under strain control at total strain ranges varying from 0.009 to 0.036. Relative to the Ni-200, the slip at the surface of the TD-nickel was more wavy and discontinuous due to the presence of the thoria particles. This made crevice formation (incipient cracking) within slip bands more difficult in TD-nickel than in Ni-200. Both materials cyclically hardened to a constant (saturation) flow stress which increased with increasing plastic strain amplitude. Cellular substructures were developed in both materials during cycling. The cell size in TD-nickel was controlled by the thoria particle distribution and was independent of plastic strain amplitude over the range investigated. The cell size in Ni-ZOO was larger than that in TD-nickel at similar plastic strain amplitudes and was a function of plastic strain amplitude. These results, together with the cyclic stress-strain curves for both materials, are discussed in terms of a model for fatigue strain accommodation at saturation recently proposed by Feltner and Laird. NUMEROUS fatigue investigations have considered the interrelation of slip character, dislocation substructure, and cracking in pure metals and solid-solution alloys. However, except for the studies of the low-strain fatigue of internally oxidized copper alloys1 and cast, dispersion-strengthened lead,' little is known about the effects which small, incoherent, nondeform-able, second-phase particles have on cyclic deformation and cracking processes. Effects due to the particles alone are often obscured by a dislocation substructure introduced during thermomechanical processing of dispersion-strengthened metals. In the present study, recrystallized TD-nickel and a commercially pure nickel, Ni-200, were employed to evaluate the effect of a thoria dispersion on high-strai fatigue deformation and cracking at room temperature. I) MATERIAL AND EXPERIMENTAL PROCEDURE The TD-nickel was supplied by DuPont as a 5/8-in.-thick stress-relieved plate which had been subjected to a proprietary schedule of thermomechanical treatments, and the Ni-200 as 3/4-in. bar which was subsequently annealed for 2 hr at 850°C in argon resulting in an average grain diameter of 0.05 mm. The compositions of these materials are given in Table I. The microstructure of the TD-nickel consisted of elongated grains parallel to the primary working direction with an average width of 0.16 mm, Fig. l(a). Many fine annealing twins were present indicating that the starting material was in a recrystallized condition; this supposition was confirmed by the absence of of any extensive dislocation substructure, Fig. l(b). Sheetlike stringers parallel to the rolling direction were occasionally seen both within grains and at grain boundaries. Some approximately spherical particles about 2 in diam, which may correspond to exceptionally large thoria particle aggregates, were also present. The average Young's modulus of the plate material in the rolling direction was 21.8 X 106 psi which is consistent with a {100}<001>recrystalliza-tion texture3'* being prominent. In transmission microscopy, the 2.3 vol pct of thoria particles generally appeared to be uniformly distributed although some clusters, 0.1 to 0.3 in diam, of larger particles were observed as previously reported for TD-nickel sheet,5 and stringering of particles was present in some areas as welt. The average diameter of the thoria particles was 450A with a calculated mean planar center-to-center spacing of 2100A, as determined by quantitative metallographic analysis.= The 0.2 pct offset yield stress was 36,000 psi which agrees with the value predicted by the modified Orowan relation7 for edge dislocations bowing between thoria particles of the size and spacing observed in the present investigation. Fig. 2 illustrates the specimen design employed for the axial high-strain fatigue testing. Adapters were screwed onto the threaded portions of each specimen so that testing could be performed in the same manner as that reported for buttonhead specimens.8 Stressing was coincident with the working direction for both materials. The gage section of each specimen was electropolished and lightly etched prior to testing. The total strain was controlled, being varied between zero and a maximum tensile strain ranging from 0.009 to 0.036. In addition to these tests, a circum-ferentially notched TD-nickel specimen was cycled over a total strain range of 0.0075. The same strain
Jan 1, 1969
-
Reservoir Engineering-General - Pressure Build-Up Analysis, Variable-Rate CaseBy F. Selig, A. S. Odeh
A second-order approximation to the exact solution of the diffusivity equation corresponding to the pressure build-up of a well producing at a variable rate is derived. This approximation is applicable when the well's shut-in time is larger than the total time elapsed since the well was first produced. The resulting equations are compact in form and easy to use. Thus, the need for Horner's' theoretically precise but rather laborious solution to the above problem is eliminated. In addition, these equations apply where the use of Horner's widely known approximate method is questionable. From a practical point of view, the reported method is best suited for analysis of drill-stem tests and short production tests conducted on new wells. INTRODUCTION The utility of drill-stem and short production tests in reservoir studies has long been recognized by the reservoir engineer. If interpreted correctly they could lead to a wealth of information upon which may depend the success or failure of reservoirs' analyses. Initial reservoir pressure and the average flow capacity are two quantities that are normally sought from a drill-stem and/or a short production test analysis. Pressures are the most valuable and useful data in reservoir engineering. Directly or indirectly, they enter into all phases of reservoir engineering calculations. Therefore, their accurate determination is of utmost importance. The flow capacity kh of the reservoir is indicative of its commercial capability. In addition, it can indicate the presence of a darnaged zone around the wellbore and, thus, the necessity for remedial measures. Of the several methods used to analyze drill-stem and short production tests, Horner's' method is by and large the most common. It applies to an infinite reservoir and or a limited reservoir where the effect of production has not been felt by the boundary. Horner's method makes use of the so-called "point-source" solution of the diffusivity equation. The point-source solution is approximated by a logarithmic function and the superposition theorem is utilized to give the familiar pressure build-up equation where is the shut-in time, q is in reservoir barrels per day and the rest of the symbols conform with AIME nomenclature. Eq. 1 was derived for a well which produced at a constant rate q from time zero to time t and was then shut in. In actuality, such a constant rate of production does not normally obtain. Therefore, a correction must be applied to Eq. 1 to account for the varying rates of production. Horner suggested two methods. The first, which results in a theoretically accurate solution, is rather lengthy and laborious and, thus, it is not suited for routine analysis. The second which has been termed a "good working approximation" is the one used by the majority of the reservoir engineers. In the second method, Eq. 1 is modified by simply introducing a corrected time t, and writing where q is the last established production rate prior to shut-in, and t, is obtained by dividing the total cumulative production by the last established rate. Horner's original paper does not give any indication that this method of correction is based on any theoretical justification. In addition, there is a question as to what constitutes the last established rate. In case of a drill-stem test some engineers use the average rate obtained by dividing the total fluid produced by the total flow time, while others calculate the average rate by dividing the total fluid produced by the last flow-period time. Obviously, different results obtain for the different flow rates used. Because of this, a simple method to the varying-rate case was developed which is theoretically sound and which defines clearly the flow rate and its associated time to be used in the calculations. The final equation arrived at is where q* and t* are a modified rate and time, respectively, and can be easily calculated. In addition, it is shown theoretically that Horner's approximate method, if used for a variable-rate case, gives the correct pressure but would not be expected to give the correct flow capacity. MATHEMATICAL ANALYSIS The general equation governing the flow of slightly compressible fluid in porous media may be written as The elementary solution to Eq. 4, representing an instantaneous withdrawal of Q units volume of fluid at the origin at t = 0, is known as the instantaneous sink