Search Documents

By James E. Harding

AT just about the geographical center of the peninsula of Lower California is the El Arco gold mine. It is small and spotty, and three separate attempts to operate it in the past have failed. The only reason the present operation is possible is because of the higher price of gold and because the previous operators had spent about $300,000 developing the property with shafts, four levels, machinery, and a complete mill, all of which was acquired at a bargain price. All the known ore in the mine, to say nothing of the cost of recovering it, is

Jan 1, 1937

By George Krauss, M. Cohen

The reverse martensitic transfomzation (i.e., the conversion of martensite to austenite on heating) was investigated in Fe-Ni alloys containing 30.5 to 33.5 wt pct Ni. The reversed austenite was found to be sufficiently distorted to exhibit marked differences in mechanical properties and to undergo recovery and re crystallization on further heating. In particular, the yield strength of the reversed austenite was as much as 2.5 times higher than that of the virgin austenite: with five cycles of the direct and reverse martensitic transformation, the yield strength was 2.6 times higher. There was an accompanying decrease in strain-hardening capacity and in stable ductility or uniform elongation. The recrystallization was observed to follow sigmoidal kinetics, but was characterized by an abnormal1y large temperature dependence which was evidently due to carbon segregation at the migrating grain boundaries. Upon recrystallization, the strength Properties reverted to those of the virgin austenite, but the ductility did not fully recover until the annealing temperature was further raised to dissipate the carbon segregation at the grain boundaries. The transformation of reversed austenite to produce second-generation martensite resulted in an unusually imperfect structure, thereby leaditg to additional strengthening. Horlleuer, the martensitic phase was much less responsive in this respect than was the austenite. MOST martensitic transformations are capable of undergoing a reverse transformation back to the austenitic phase on heating.&apos; This reaction, like the direct transformation to martensite on cooling, is usually athermal. In analogy with the Ms andMf temperatures for the martensitic transformation (where M, > Mf), the reverse transformation is characterized by the As andAf temperatures (where As<Af). It is well-known that the finish (Af) of the martensite-to-austenite reaction on heating may occur at a temperature where austenite is not stable relative to other phases in the equilibrium diagram.&apos; Reverse martensitic transformations have been found to influence the properties of the parent austenitic phase in a number of ferrous alloys. The austenite thus produced after one or more reversals is significantly altered in strength, microstruc-ture, diffusivity, and annealing characteristics compared to well-annealed, virgin austenite. The first detailed observations of such phenomena were made by Wasserman3 on an Fe-30 wt pct alloy. He found that the reversed austenite (i.e., austenite formed by the reverse martensitic transformation) was much stronger than the original austenite, and the second-generation martensite (i.e., martensite produced from the reversed austenite) was somewhat stronger than the first-generation martensite. On heating the reversed austenite, recrystallization occurred in a manner similar to that of a cold-worked metal; the Laue spots sharpened and the strength decreased.4 In recent years, Russian workers have shown considerable interest in these reverse martensitic transformations. Most of such investigations5-10 have been concerned with the stabilization of reversed austenite relative to subsequent martensitic transformation, and because of its unusual nature, this phenomenon will be discussed in detail elsewhere." On the other hand, the annealing behavior of reversed austenite and the attendant variation in properties have received little attention beyond the experiments of Wasserman.3,4 Gridnev eta1.,8 in one of the few such studies. followed the hardness change of reversed austenite as a function of time at three annealing temperatures. The marked strengthening of austenite achieved by the reverse martensitic transformation has been referred to as "transformation strengthening" by Wasserman3,4 and as "phase hardening" by Malyshev et a1.7,9,10 to differentiate it from hardening produced by cold working. Both groups of investigators noted that further strengthening of the austenite could be attained by repeated or cyclic transformations, but the first cycle resulted in the most pronounced effect. Evidence of the altered microstructure of reversed austenite has been reported by Maximova and Nikonorova6 and by Edmonson and KO.12 The former found that reversed austenite was much more sensitive to etching than the original austenite, and the latter authors observed a (&apos;ghost" structure resembling martensite in certain areas of the reversed austenite. Enhancement of the self-diffusion rate in reversed

Jan 1, 1962

By O. W. Simmons, L. W. Eastwood, C. M. Craighead

Binary alloys of titanium with silver, lead, tin, nickel, copper, beryllium, boron, silicon, chromium, molybdenum, manganese, vanadium, iron, and cobalt were studied. One-half-pound ingots of the alloys were prepared in an arc furnace, employing a water-cooled copper crucible, an argon atmosphere, and a water-cooled tungsten electrode. The half-pound ingots were fabricated by forging at 1700°F in air to 1/4-in. slab, followed by hot rolling at 1450°F to 0.060-in. sheet. Tensile properties, minimum bend radii, hardnesses, response to heat treatment and aging treatment, and phase relationships were determined for these alloys. EARLY in 1947, as one phase of the evaluation of materials for Air Force Project RAND, Battelle successfully arc melted Bureau of Mines titanium and obtained a few basic properties of unalloyed titanium and several titanium-base alloys. The low density, excellent resistance to corrosion, and high tensile properties of these materials stimulated a great deal of interest among metallurgists and designers seeking better materials of construction. As a result of this early work, Battelle, under contract with the Air Materiel Command, Wright-Patterson Air Force Base, has continued extensive studies of titanium alloys. The result of the first year&apos;s work is described in three papers. The present one deals with binary alloys, the second with ternary alloys, and the third paper with quaternary alloys. Melting Method The arc furnace for melting titanium and other refractory metals was developed at Battelle Institute under Air Force Project RAND. During the present work, considerable improvement has been made in the design and operation of this furnace for making small ingots of titanium and titanium alloys. The improved furnace design is illustrated by fig. 1, which shows a cross-sectional view with the dimensions of various parts and their relationship to one another. The essential features of the furnace are the inert argon atmosphere, the water-cooled copper crucible, and the water-cooled tungsten elec- trode. The melting crucible is formed by spinning 0.064-in. annealed copper sheet. A groove for an O-ring seal is spun into the flange of the crucible. A fiber gasket between the crucible flange and the water-cooled brass cover provides electrical insulation. The brass cover is fitted with a sight glass, an opening for the water-cooled electrode, and a tube through which the material is charged. Direct current is used for melting, with the positive electrical terminal connected to the water jacket and the negative terminal to the electrode. A typical heat is made in the new furnace as follows: Approximately, 0.20 lb of titanium under 2 mesh per in. and the alloy addition, if any is to be used, are placed in the bottom of the crucible. The cover is clamped on the crucible so that the O-rings make a tight seal. The remainder of the charge is placed in a glass bottle, and this is connected by a rubber hose to the charging tube of the furnace. The crucible and the charge are evacuated with a mechanical pump to a pressure below 50 mm of mercury, or less if desired, and held at this reduced pressure for 5 min. Hot water is circulated through the jacket during the evacuation period to assist in outgassing the crucible. Following the evacuation, tank argon of 99.92+ pct purity is flushed through the crucible for 5 min and then adjusted to give a positive pressure of 1 to 2 psi. During subsequent operation, an argon pressure regulator introduces only enough argon to compensate for the small leakage which occurs through a mercury trap. Re-evacuation and reflushing the melting chamber with argon produced a negligible reduction in contamination. Two 550-amp generators, connected in parallel, constitute the power source. Normally, about 800 amp are used for melting alloys which are not extremely refractory. The generators are set for 90 v open circuit and 200 amp. The arc is struck and the electrode is withdrawn to produce an arc of 23 to 26 v. This voltage is maintained while the electrode tip is moved slowly around a circle 1 in. smaller than the diameter of the ingot. The current

Jan 1, 1951

By J. G. Savins

Certain types of macromolecules added to water and salt solutions flowing in turbulent motion can reduce the pressure gradient. Alternatively, the volumetric capacity of a pipe for these fluids is increased by the presence of these materials. Examples presented show that the drag reduction can become significant. Thus, the presence of 0.28 per cent of a gum derivative in a solution of sodium chloride flowing at 200 gal/min in a 1.89-in. pipe yields a pressure drop which is 0.44 of the single-phase drop measured under the same conditions of turbulent flow; the addition of 0.1 per cent of a vinyl derivative to a I-in. water line yields a throughput capacity which is 1.78 of the single-phase capacity at the same pressure drop. It is further shown that these phenomena are distinctly different from previous observations with other classes of non-Newtonian systems. There a simple lowering of friction factors below the levels predicted from the resistance laws for Newtonian fluids is associated with a suppression of turbulent motion. A rational physical explanation for drag reduction is advanced. Briefly, the proposed mechanism is a storage by the molecular elastic elements of the macromolecules in solution of the kinetic energy of the turbulent motion. INTRODUCTION This study was inspired by a recent review of some paradoxical drag reduction phenomena in turbulent pipe f1ow.I Under very moderate conditions of turbulent flow, the pressure gradient necessary to pump solutions containing certain specific kinds of polymers, fibers and metallic soaps may become appreciably lower than that required* to pump the solvent, i.e., water or a low-viscosity hydrocarbon, under identical flow rates in the same conduit. As shown by our review, this phenomenon of drag reduction in turbulent duct flow was first noted during the second world war, apparently arising in connection with the development of flame warfare weapons. Since that time several papers illustrating this phenomenon have appeared: Toms, Oldroyd,3 Agoston et al., Bundrant and Matthews, Robertson and Mason,6 Ousterhout and Hall,7 Daily and Bugliarello~ Lummus, Anderson, and Fox. 9 That there are practical applications for techniques which increase discharge or decrease the pressure necessary to transport a liquid through a pipeline is illustrated in the patents which have issued which take advantage of this peculiar phenomenon, e.g., Mysels,l0 Dever, Harbour, and Seifert.11 One also finds fragmentary evidence of this effect in the data pertaining to a few of the polymeric solutions studied by Shaver12 and Dodge.13 However, these investigators were concerned with the development of friction factor vs Reynolds number correlations for a variety of non-Newtonian solutions and suspensions, rather than in a study of drag reduction-. A similar kind of drag reduction effect has been observed in gases. Sproull, for example, reports that adding dust to air flowing in turbulent motion through a pipe results in a lowering of the pressure gradient at identical flow rates. There are also military applications for reducing the drag on hydrodynamic vehicles. For example, the possibility of injecting a rheologically complex fluid into the boundary layers of bodies to reduce the skin friction has been investigated by Fabulal5 and Granville.16 Along somewhat different lines are the drag reduction studies of Kramer.17118 He has shown that skin friction can be reduced by covering the surface of a vehicle with a flexible skin. The effect is apparently due to the boundary layer being stabilized by the presence of the skin. Drag reduction by means of coexisting gas and liquid boundary layers, e.g., film boiling and continuous gas injection, has been proposed by Bradfield, Barkdoll, and Byrne,19 Cess and Sparrow,20&apos;21 Sparrow, Jonsson, and Eckert.22 Here the skin friction occurs between a vapor and a surface rather than between a liquid and a surface. There are several references in the literature to friction-factor correlations for non-Newtonian solutions and suspensions: Shaver and Merri11,23 Dodge and Metzner, Clapp,25 and Thomas.26

Jan 1, 1965

By C. S. Finney, John Mitchell

The introduction of ovens for the production of metallurgical coke is believed to be due to L. L. Norton who operated an iron foundry in the vicinity of Connellsville, Pa. Persuaded by his foreman, an English immigrant named Nickols from Durham, L. L. Norton put up a 12-ft square oven which produced coke in 1833. The coal used was taken from a local mine at Mounts Creek. The oven seems to have been used in conjunction with the customary method of coking in mounds. It was in the Connellsville district also, in 1841, that two carpenters, Provence McCormick and James Campbell, formed a partnership with John Taylor, a stone mason, for the manufacture and sale of oven coke. The task of the mason was to construct the ovens, while the carpenters were to build the arks by which the coke could be taken by water to the market at Cincinnati. The following account of the enterprise was given by McCormick: "James Campbell and myself heard in some way that I do not now recollect that the manufacturing of coke might be made a good business. Mr. John Taylor, a stone mason, who owned the farm on which the Fayette coke works now stand, and who was mining coal in a small way, was spoken to regarding our enterprise, and proposed a partnership-he to build the ovens and make the coke and Mr. Campbell and myself to build a boat and take the coke to Cincinnati, where we heard there was a good demand. This was in 1841. Mr. Taylor built two ovens. I think they were about 10 feet in diameter. My recollection is that the charge was 80 bushels. The ovens were built in the same style as those now used, but had no iron ring at the top to prevent the brick from falling in when filling the oven with coal, nor had we any iron frames at the mouth where the coke was drawn. The top and mouth had to be repaired when they fell in. In the spring of 1842 enough coke had been made to fill two boats 90 feet long-about 800 bushels each-and we took them to Cincinnati down the Youghiogheny, Monongahela, and Ohio, but when we got there we could not sell. Mr. Campbell, who went with the boats, lay at the landing some two or three weeks, retailing out one boatload and part of the other in small lots at about 8 cents a bushel. Miles Greenwood, a foundryman of that city, offered to take the balance if he would take a small patent flour mill at $125.00 hi pay, which Mr. Campbell did. He had it shipped here. We tried it, but it was no good, and we sold it to a man in the mountains for $30.00, and thus ended our coke business." So successful did the coke subsequently prove to be in use that the three partners were asked to deliver more. Evidently they had had enough of the coke business, however, for they refused to have anything more to do with it. Few ovens were built between 1841 and 1855, and it is reported that in the latter year, "there were only 26 coke ovens along the river above Pittsburgh". Successful coke makers of these years included Mordecai Cochran, Richard Brookius, and Colonel A. M. Hill. It was the use of coke in 1859 in the Clinton furnace erected by Graff, Bennett and Co. in a plant on West Carson Street, Pittsburgh, that brought the real beginning of the coke-iron era in America. Here the successful use of Connellsville coke as a blast-furnace fuel was demonstrated beyond all possible doubt, and from the year 1859 the coking industry expanded tremendously. The era of beehive coke ovens During the latter half of the nineteenth century and the early years of the twentieth, the major percentage of metallurgical coke produced in the United States came from beehive ovens. It was not until 1893 that the first battery of by-product ovens came

Jan 1, 1961

By D. M. James, J. C. Martin

The results are presented of a study of the application of analytical methods to the solution of two-phase flow into single wells. Approximate analytical expressions for the pressure distribution in two-phase flow are found ;;; a number of conditions. The results obtained from the analytical solutions me found to be in good agreement with results obtained by finite difference techniques using a high speed digital computer. Mathematical solutions for four sets of boundary conditions are presented. All of these solutions are composed of a short-term transient plus a steady or quasi-steady state. The rates of decay of the short-lived transients are analyzed. It is found that the durations of the short-term transients may be characterized by a parameter defined as the time constant which can be determined from simple relations. It is shown also that if the outer radius is much greater than the radius of the well. the short term transients decay at rates which are proportional to the square of the exterior radius, and the rates of decay are only slightly dependent upon the radius ratio. Numerical solutions based on finite difference techniques are presented for a number of conditions. The numerical solutions are in good agreement with the predictions based on the theoretical analysis for small and moderate drawdowns. Examples involving large drawdowns indicate that the nonlinearities in the equations of flow do not appreciably alter the longevity of the short-term transients. In all cases the time required for the short-tern transients to disappear is predicted satisfactorily. INTRODUCTION The mechanism by which oil and gas flow into a single well is of vital interest to the petroleum industry. The fundamental equations of two-phase flow which describe this mechanism are nonlinear partial differential equations. Numerical solutions of these equations describing pressure transients have been obtained with the aid of electronic computers. Although solutions obtained in this manner take into account a large number of effects, the reduction of this information to useful generalities is difficult. One method of obtaining generalities is the use of linearized approximations of the nonlinear equations. Since it is possible to obtain explicit solutions of the linearized equation, general properties of the role of pressure in the flow mechanism may be ascertained. Results obtained from this approach are limited to some extent by the linearizing assumptions. The severity of these limitations may be evaluated by comparing solutions of the linear equation with numerical solutions of the more exact nonlinear equations of two-phase flow. In the past considerable amount of work has been devoted to studying pressure build-up using the single-phase flow theory. Unfortunate1y, most prep sure build-up tests involve multiphase flow. A small amount of work has been done studying pressure build-up where the flow is two-phase. The encouraging results of these studies suggest that useful results may be found from additional studies of not only pressure build-up, but also the rapid transients associated with placing a well on production. This paper presents the basic theory of the pressure transients associated with placing a well on production and with closing it in. The paper is concerned chiefly with two-phase compressible flow; however, the results also apply to single-phase flow. The results are based on analytical solutions of the flow equations, and they are verified by numerical solutions using finite difference techniques. Much of the previous work on compressible flow into wells has been confined to single-phase flow. Some work has been done on compressible two-phase steady state flow, and solutions of the equations of flow have been found by finite difference techniques using high-speed computers. Muskat1 presents some rather general solutions to the equations of single-phase compressible flow into wells. Much work has been done on pressure build-up in wells (see, for example, Refs. 2-6). Almost all of the work on pressure build-up concerns single-phase flow with the exception of Ref. 5 and part of Ref. 2. Some work has been done on pressure fall-off in injection wells.718 Muskat9 presents the solution of the equations for radial steady state two-phase compressible flow. Handyl0 extends Muskat&apos;s solution to obtain

By John D. Price

R. W. Campbell (Jones and Laughlin Steel Carp., Pittsburgh)—As usual John Price has presented an excellent paper. I know of no one who has devoted more time and conscientious thought to this subject than he has, and his efforts have produced scientific data and conclusions of real value. Results of similar studies made by the Jones and Laughlin Steel Corp. on the blending of Pennsylvania coals may be of interest. During 1935 the late F. W. Wagner and the present writer developed an empirical physical fuel value formula for evaluating the combined physical tests of metallurgical coke, see Table IV. In the development of the formula, the tumbler test, shatter test, screen test, and porosity of coke have been taken into consideration. Each of these tests brings out some physical property not expressed by the others, and in estimation of the importance of these tests as a medium for evaluating coke for the blast furnace, their worth has been graded as follows: It was thought that some additional penalty should be applied to cokes that do not show a uniform screen size, assuming that coke larger than 4 in. was not desirable and that this size was an indication of irregular size. Therefore the screen test factor is further corrected by applying a penalty of 1 pct for each 1 pct material greater than 4-in. size as shown by the screen test. The physical fitness of coke for blast furnace use depends upon four fundamental characteristics: l—gen-era1 strength, as indicated by the shatter test; 2— abradability, or the resistance of cell walls to reduction in size by attrition, as indicated by the tumbler test; 3—porosity; and 4—size and uniformity. Each of these characteristics is not of equal importance, but surely, in an attempt to determine an accurate physical value, the individual value of each should be taken into consideration. The formula is used as a yardstick, so to speak, in determining variations in the normal production of metallurgical coke and in making comparisons between normal coke and experimental blends with other coals. Table IV gives results of coking tests on raw and washed high volatile coals and these same coals with blends of coke breeze, anthracite coal, char and low volatile coals, all from Pennsylvania. The table also gives the source of the coals and coke breeze, as well as their proximate analyses. It will be noted that the washing of the raw high volatile coal increased the physical fuel value from 26.17 to 39.16, showing the effect of removing shale and slate which produces a detrimental fracturing and weakening of the coke. It is a well known fact that many Pennsylvania high volatile coals contain an excess of resinous material or binder which has a tendency to weaken the coke, owing to its high evolution of gases during carbonization. Tests were therefore made with coke breeze, an inert material, in an attempt to alleviate this condition. Blends of washed coal with 5 pct and 7.5 pct coke breeze collected from the coke plant quenching sump (82.9 pct through 1/8-in. screen) resulted in cokes with a physical fuel value of 20.75 and 16.73 respectively, lower than results obtained from 100 pct raw high volatile coal. It will be noted that the percentage of + 4-in. coke increased but that the tumbler test was very low. Undoubtedly better results could have been obtained if the coke breeze could have been crushed extremely fine. During World War II the government asked the coking industry to consider the use of anthracite coal as a substitute for high quality low volatile coal, which was extremely scarce at the time. The next series of tests was made with blends of raw and washed high volatile coals and 5, 10, and 15 pct Pennsylvania No. 5 buckwheat anthracite (100 pct through l/8-in. screen). It was necessary to use some raw high volatile coal in these blends because of insufficient washer capacity at the time. It will be noted that the physical fuel values

Jan 1, 1954

Edward Dyer Peters, the only child of Henry Hunter Peters and Susan Barker Thaxter, was born in Dorchester, Mass., June 1, 1849. From his father he was a descendant of the Peters family of Ipswich and Andover, Mass., whose earliest ancestor settled in Essex County in 1659; while upon his mother&apos;s side he descended from several old Hingham families, as well as from the Quineys and from Governor Bradford. Between the ages of fourteen and sixteen he was a student at the Episcopal School for boys in Cheshire, Conn. Near the school was an old tin mine, and he used to spend his Saturday afternoons exploring it. When Edward was sixteen, his father sold his property in South-borough, and for the next three years lived in Europe. Sept. 26, 1865, Edward entered the Royal School of Mines at Freiberg, Saxony. In the summer of 1867 he joined his family on Lake Geneva, and spent a part of the season in that vicinity, going to Chamounix and taking other trips. Extract from a letter to his father: Elbingerode, Harz, Aug. 30, 1868: I left Freiberg with my friend Lilienthal, on Wednesday the 19th of August, reached Leipzig in the evening and spent the night in that city: the next morning we took the cars and in the afternoon reached Eisleben, in Prussia, which was to be our first stopping-place; there are immense copper-mines and smelting-works at this place, and the method of working the mines is very peculiar ... The next morning we got up at five and spent the whole day in the smelting-works where I obtained a number of drawings of the furnaces they use here and a good many valuable data. Saturday morning we got up at four o&apos;clock, and descended a mine about three miles from the town; the Director of this mine was one of the most intelligent and agreeable men that I ever met; he accompanied us in the mine and explained everything that was new to us. The system of mining at Eisleben is peculiar to this town, and is not used at any other place in the world; the copper ore comes in narrow, flat veins eighteen inches high, at the utmost, consequently the men work lying flat on their left side and using the pick with both hands: the ore is carried out of the side passages into the main ones, in very low waggons about twelve inches high, and four by two feet; the boys who pull these waggons have a board buckled to their left thigh and their left arm; they pass their bare foot through a strap in the waggon and wriggle along on their side dragging the waggon after them, with their foot, almost as fast as a man can walk above ground. In the afternoon we walked over to Hittstadt, about twelve miles, and remained there Sunday and Monday. Monday night, at ten o&apos;clock we went out to one of the smelting-works, to see the refining of the copper; this is only done at night, so we worked there from ten in the evening till eleven the next morning, and came home all blackened up and with blistered hands, for our skin was not quite so tough as that of the workmen: at noon the same day we took the stage to Ascherls-leben, and then the train to Stassfurt, where there are immense salt-mines, as well as factories for the manufacture of bromium, which latter process is, however, kept a profound secret, so we had very little hope of seeing anything but the salt mines; we were talking with our landlord about our disappointment, that evening, when he told us that the director of the Bromium works was sitting in the next room with no company but his glass of beer, so we went over to him, introduced ourselves as Americans, and actually succeeded in persuading him that it could do his business no earthly hurt if he showed us his works; at last he consented, after making us give our word of honor

Jan 1, 1919

By B. Yarar, J. Kaoma

[Introduction The critical surface tension of wetting of hydrophobic materials has been investigated extensively by Zisman et al. (1973) and relates the spreading of a liquid on a solid to the surface tension of this liquid. Thus, the Zisman ap¬proach shows that a plot of Cos 0 versus ?LV where, 0 is the contact angle and ?LV is the surface tension of that liquid in contact with its own vapor at equilibrium, produces a straight line which can be extrapolated to Cos 0 = 1. The intercept corresponds to a given ?LV, termed the critical surface tension of wetting. In instances, however, Cos 0 - ?LV plots show deviations from linearity. Alternative approaches have also included the plotting of adhesion tension (YLV COs 0) versus ?LV that can be interpreted in terms of critical surface tension data (Lucassen-Reynders, 1963). It has, on the other hand, been customary to correlate contact angle (0) with flotation data, although it is recognized that surface inhomogeneities and related contact angle hysteresis limit the interpretations pertaining to solid surface structures purely based on contact angle measurements. It is nonetheless true that 0 = 0, when the liquid is water, corresponds to a hydrophilic solid. Determination of ?c by contact angle measurements in conformity with the Zisman approach can encounter experimental difficulties, especially when solids are in a powder form. Following the reasoning that 0 = 0 corresponds to the hydrophilic state of a solid, froth flotation experiments can be expected to produce critical surface tension of wetting values for powders. The technique that consists of plotting recovery (% R) versus YLV data and extrapolating the linear part of the curve to % R = 0, has been described elsewhere (Yarar and Kaoma, 1984). This paper describes some of the results obtained by the extrapolation technique mentioned above, using various sulfide minerals. Experimental Materials Molybdenite (MoS2), sphalerite (ZnS), chalcopyrite (CuFeS2), chalcocite (Cu2S), and galena (PbS) were handsorted from massive samples and ground by mortar and pestle, and the -147 Fan fraction used in experiments; sulfur (S) was a C.P. Fisher Scientific Co. product. Method One gram powder samples were conditioned for three minutes in methanol-water solutions of predetermined surface tension values and the percents recovered by flotation using a Partridge-Smith Cell (Partridge and Smith, 1972), were plotted against each ?LV value. Results and Discussion The percent recovery versus ?LV for two naturally hydrophobic solids, i.e., natural molybdenite, MoS2 and sulfur are given in Fig. 1 where it is seen that ?c for sulfur is obtained as 31.5 mNm-1 whereas for MoS2, it is 26 mNm-1. Figure 1 also shows that the new flotation technique described, for estimating ?c is in good agreement with the Zisman method. Curves obtained by the same method for the minerals Cu2S, PbS, and ZnS are given in Fig. 2, where it is seen that for Cu2S, ZnS, and PbS, the extrapolation to % R = 0 produces ?c values of 36, 39, and 49 mNm-1, respectively. CuFeS2, on the other hand, gives a curve that does not permit such an extrapolation to obtain a ?c value with certainty. Nonetheless, treatment of the powder with 8.95 X 10-5 molar Na2S solution followed by washing with distilled water prior to flotation effects the curves as shown in the insert of Fig. 2. The type of curves given in Fig. 1 appear to be typical of hydrophobic solids as can also be observed in Fig. 3. The composition of solid and minute quantities of inclusions seem to play an important role as regards the surface properties of solids as reflected by the slope of the linear part of the % R - ?LV plot and thus the value of ?c obtained by this technique. For example, MoS2 of three different stoichiometric compositions produced the data given below. [MoS2MoS2.14Mo1.002S2 ?c mNm-126.029.031.0 dR/d?8.95.23.0]]

Jan 1, 1985

By R. F. Mehl, C. E. Birchenall

SINCE Maxwell1 first considered the self-diffusion process in 1872 its importance in the kinetic theory of matter has been recognized. Until the discovery of isotopes in 1913, a direct measurement of this quantity seemed impossible. The only information that could be gathered indirectly was for gaseous systems for which the kinetic theory was well developed. When methods of direct observation be-came available, investigations on self-diffusion were carried out in condensed phases. In metals these data are presumably of potential importance in the study of recovery, recrystallization, creep, sintering, and related phenomena. Following the pioneer work of Von Hevesy2 and his collaborators, who determined the rate of self-diffusion in lead with the naturally occurring thorium B isotope, several metals have been investigated. Pb, Ag³, Au.6,5, and Cu6,7,8,9 are examples of isotropic crystals existing in only one phase modification. Anisotropy of diffusion has been demonstrated in Bil6 and Znll. This paper is a study of a single metal in two allotropic forms and also of self-diffusion in a body-centered cubic lattice. Despite the fact that the measurements in each of the papers cited on self-diffusion in copper seem to be internally consistent to about 10 to 15 pet, the curves reported by different authors differ by factors of 2 to 4 at the same temperatures. This uncertainty may account, in part, for the failure to establish a successful correlation between the self-diffusion rates and other physical characteristics of the metals. No satisfactory theory of metallic self- diffusion has as yet been proposed to account for all the existing data and to permit estimation in other metals. Experimental Part: Two units of radioactive iron were used in these experiments, each a mixture of Fe53 and Fe59 (12). Fe53 decays by K electron capture and the emission of an X ray with a half life of about 4 years. Fe59 emits two beta spectra, one with a maximum energy of 0.26 Mev, the other 0.46 Mev, and gamma rays of 1.1 and 1.3 Mev, with a half life of about 44 days. They were present in nearly equal concentration initially. The composite half life started at about 50 days and increased as the Fe59 decayed, leaving Fe55 relatively more abundant. After aging a year, the half life was too long to measure significant decay in a month. The absorption properties also changed with time. Because of the importance of the absorption coefficient in the diffusion calculations, it was necessary to determine this quantity frequently over the period during which these experiments were carried out. This correction was unsuspected at the time of publication of notesl3 on this work and accounts for the discrepancy in alpha iron and for part of the discrepancy in gamma iron.* * The data in the present paper supersede those given in the notes entirely. In one note the scale of log D was inadvertently reversed. In addition to the correction for the drift in absorption coefficient, the D values for gamma iron at low temperatures were high owing to insuficient correction for diffusion occurring during heating and cooling through the high temperature part of the alpha range at a rate much slower than that employed in the runs reported here. The high temperature alpha rates are much higher than the low temperature gamma rates and correspond to large equivalent times at these temperatures. These experiments were discarded and new measurements taken. Independent experiments with the same activity units indicated a radioactive contaminant in the iron, but exhaustive attempts to isolate and identify it chemically failed. These experiments? indicate † These experiments will be discussed in a later publication from this laboratory. that the contaminant represented only a trace of little importance in the diffusion results. The active iron was plated from an iron chloride solution. Enough of the solution was dropped on a 1% in. square of filter paper to wet it thoroughly.

Jan 1, 1951

By L. S. Allen, R. L. Caldwell, W. R. Mills

Borehole models used in the study of nuclear well logging are reviewed and the merit of heterogeneous us homogeneous formation simulation is discussed. .A heterogeneous model for simulating sandstone formations of variable porosity is described. Simulation adequacy was checked by steady-state and pulsed-neutron measurements. INTRODUCTION In the course of developing new methods of nuclear well logging, the problem of testing new logging instruments under controlled conditions always arises. This testing is usually done in wells of known lithology, porosity and fluid saturation, or in laboratory test facilities which simulate the borehole environment. Although the former approach is possibly more satisfying, it is frequently difficult or inconvenient to find wells with enough variation in the cited parameters for a truly comprehensive study, On the other hand. laboratow- test facilities can simulate a wide variety of formation and borehole conditions if care is exercised to maintain adequate simulation under all conditions. Workers concerned with the development of nuclear well logging methods have reported several types of research facilities in the past. Tittle, et a1 used cubical tanks filled with Ottawa silica sand to study the steady-state behavior of thermal and epithemal neutrons in simulated sandstone logging environments.1 In addition to the effects of fresh vs salt water, a wide variety of casing sizes and compositions were investigated. Tittman also used a large tank of Ottawa sand to measure the slowing down length of sio2;2 for making the same measurement in CaCO3, however, an assembly of solid Vermont marble blocks was used. More recently a very useful and carefully planned test pit for nuclear logs was constructed under the auspices of the API.3 This facility consists of three massive limestone sections, each 6-ft thick and more than 5 ft in diameter. The blocks are buried, one on top of the other, in a concrete- lined pit having a total depth of 24 ft. Each section is made of a different limestone to achieve porosity variation, the selection being Carthage (2 per cent), Indiana (19 per cent) and Austin (26 per cent). The blocks are traversed by a 7 7/8-in. diameter borehole and are saturated with fresh water. Great care was taken to insure total saturation of the stone. Homogeneous mixtures and materials have not been used exclusively for simulating earth formations. Tittle, et a1 4 showed that a heterogeneous lattice of water and graphite could be used to simulate Liquid furfural with very little effecton measured thermal and epi thermal neutron distributions. Subsequently, a lattice-type model for simulating limestone earth formations was designed and built by Berry This limestone model was constructed of 7/8- x 11/4- x 36411. marble rods glued together with a plastic resin to form a 36-x 42- x 42in. cube. The cube was honeycombed vertically with 7/8- x 7/8- x 3641-1. holes for porosity and had a 111/2-in. diameter hole in the center to serve as the borehole. Lower porosities were simulated by placing marble rods in the square holes. The finished cube was used with the borehole vertical and set in a tank to hold liquid for filling the holes in the rock matrix. A second type of heterogeneous borehole model has been proposed by Kukharenko, et a1.6= In this concept the model is composed entirely of solid material. The hydrogen equivalent of formation fluids is simulated by making the model of alternating layers of dense rock and celluloid. Additional hydrogen is included by drilling vertical holes in the rock sections and filling them with solid forms of hydrogeneous material. The volume of the holes is equal to the volume of the celluloid layers, so an axial and radial dispersion of the hydrogeneous material is thus obtained. One such model was composed of alternate 50 mm-thick sand layers and 2 mm-thick celluloid sheets. Approximately 200 appropriately spaced plastic tubes were used in each sand section. The plastic tubes were filled with paraffin for additional hydrogen concentration. This arrangement yielded a model with a water-equivalent hydrogen concentration of 10 per cent. Experimental data taken in this model and in a model consisting of the same material but homogeneously dispersed showed that the heterogeneous system adequately

Jan 1, 1966

(Members are urged to send in for this column any notes of interest concerning themselves or their fellow-members) Members and guests who called at Institute headquarters during the period Jan. 10, 1917 to Feb. 10, 1917: Victor C. Alderson, Denver, Colo. G. Komar, New York, N. Y. Arthur W. Allen, New York, N. Y . H. H. Lauer, Allentown, Pa. George S. Blair, Pittsburgh, Pa. Frederick MacCoy, Antofagasta, Chile. F. J. Brulé, Toronto, Canada. E. P. Mathewson, Toronto, Canada. W. M. Corse, Buffalo, N. Y. Douglas Muir, San Antonio, Tex. John J Croston, Cullasaja, N.&apos; C. F. W. Osborn, Ajo, Ariz. W. A. Deane, Waxhaw, N. C. Mortimer A. Sears, Washington, D. C. P. H. Dudley, New York N. Y. W. H. Selden, Jr., Iron River, Mich. J. D. Gilchrist, Denver, Colo. W. H. Staver, Yonkers, N. Y. L. Hoffman, Pottersville, N. J. Ralph H. Sweetser, Easton, Pa. Guy R. Johnson, Florence, Wis. Benjamin F. Tibby, Salt Lake City, Utah. W. McA. Johnson, Hartford, Conn. Carl J. Trauerman, Butte, Mont. John S. Kennedy, Philadelphia, Pa. H. C. Wilmot, Arequipa, Peru. D. C. Bard has become a member of the geological department of the Chicago, Milwaukee & St. Paul R. R., with headquarters at Seattle. Jay A. Carpenter has resigned his position as metallurgist for the West End Consolidated at Tonopah and has been appointed superintendent of flip. Nevada. Packard mines it Rochester.

Jan 3, 1917

By R. L. Vaughn, A. G. Moon

In 1972 a major renovation of the flotation and regrinding circuits was completed at Kennecott Copper Corporation&apos;s Chino Mines Division concentrator at Hurley, New Mexico. At that time a new lime distribution and control system was placed in operation. Lime is used to depress pyrite, and flotation at a high pH or alkalinity is required for optimum recovery of principal copper minerals, chalcopyrite and chalcocite with good rejection of the pyrite. Regrinding is also required to liberate chalcocite, much of which occurs as intimately associated chalcocite-pyrite middling grains. Accurate control of the lime addition is necessary to obtain both optimum metallurgy and minimum lime usage. At low pH levels, flotation cell control is difficult due to excess froth, but excessively high pH levels can go unnoticed because froth characteristics do not change from those experienced at the desired pH. With manual control of the lime addition, the desired pH can be exceeded for periods of up to two hours. Test Work Plant test work on a lime control system using conventional pH probes gave disappointing results with high electrode maintenance and poor control characteristics. For example, the rougher flotation cells are operated at a pH between 11.4 and 11.6 and the probes required cleaning at least once a day, and occasionally every shift. These high maintenance requirements were considered unreasonable for a plant-wide system using six or more control loops.

Jan 1, 1979

By Hermsdorf, Richard P. E.

AMONG the newer lead smelting and refining plants of the country is that of the United States Metals Refining Co., at Carteret, N. J. Not only is the technical practice here modern and efficient, but exceptional care is taken to assure healthful operating conditions for the employees. One large building houses the entire department, which includes, as its chief units, one reverberatory smelting furnace, one blast furnace, two softening furnaces, one utility reverberatory furnace, two Faber du Faur retorts, five 120-ton kettles, three 60-ton antimouial lead kettles, and one 5-ton kettle for experimental purposes. There is a 10-ton overhead crane in the kettle aisle, which is also used for shipping lead. Other auxiliary equipment is referred to in the outline of operations. A general plan of the plant is shown in Fig. 2, and a general view in Fig. 1. Primary smelting operations are conducted in a 30 by 11-ft. reverberatory furnace with bottom, roof, and upper walls of clay brick, magnesite

Jan 1, 1934

By H. N. Treaftis, J. W. Cahn

THE previous determinations of the solvus of tin in solid lead disagree with one another by as much as 40°C or almost 10 at. pct. Even determinations that appear to be careful differ considerably in both the solubility and its temperature coefficient. Recent kinetic workl -3 on the rate of precipitation of tin from lead-tin alloys and thermodynamic work4 on the heat of solution of tin in lead have given some insight into the time necessary to reach equilibrium in this system. As this time is longer than that used in almost all previous determinations of the solvus, these determinations are questionable. Since the interpretation of recent kinetic results requires accurate values of the equilibrium solubility the present work was undertaken. DISCUSSION OF PREVIOUS RESULTS The methods used in previous studies may be placed into three classes: 1) Several alloys covering a range of compositions are allowed to reach "equilibrium" at a certain temperature.5-7 Some property, such as the X-ray lattice parameter of the lead-rich phase or the re- sistivity of the alloy, is then measured as a function of composition. The solubility limit for that temperature is taken to be that composition at which a discontinuity in the property occurs. However, the recent kinetic work has shown that tin precipitates in two stages. The first stage is rapid but leaves the lead still supersaturated. The second stage is very slow, and it seems from the equilibration times quoted by the investigators of the solvus that the end of the first stage of precipitation was usually mistaken for equilibrium. This is especially true at temperatures below 100°C where the first reaction is still quite rapid (from minutes to months), but where the second stage is so slow that very long holding times are required to reach equilibrium. Alternatively micrographic work8 on several alloys of different compositions seems to have given good results if done carefully. 2) The beginning of precipitation in a homogenized alloy is measured as it cools slowly8-10. Since reactions in solids undercool easily, such determinations will always give too low a temperature for the solvus. 3) Some property of an alloy is measured as a function of temperature while the alloy is above the solvus.5,9,11 The range of temperature can usually be extended slightly into the two phase region. Then the alloy is precipitated at a much lower temperature and slowly reheated. The solvus is taken to be the lowest temperature where the property of the reheated alloy is equal to that previously determined for the homogeneous alloy. Usually one also expects a discontinuity at the solvus in the temperature variation of the property for the reheated alloy. The difficulty in this method is the sluggishness of the dissolution process, and this method may give too high a solvus temperature. EXPERIMENTAL The phase boundary or solvus temperature was determined by a modification of method 3 for a series of chemically analyzed high-purity alloys, ranging from 5.1 to 26.4 at. pct Sn. The resistance, as a function of temperature, of a sample of homogenized alloy (which could easily be undercooled) was compared with the equilibrium resistance which that alloy sample attained after precipitation at room temperature followed by long isothermal anneals at various temperatures. Before each anneal the sample was rehomogenized and reprecipitated. The solvus temperature could be located to within 1°C as the intersection of the resistance vs temperature curves for the homogenized and equilibrium alloys. No attempt was made to locate the phase boundary more accurately since the limiting factor was the precision with which the alloys could be analyzed. This precision was ±0.1 wt pct Sn and corresponds to an error of 1 deg in temperature. This use of isothermal anneals differs from some of the previous investigations in which method 3 was used with baths of slowly rising temperatures. It was found early in the investigation that solvus temperatures tended to be high if the latter method was used even at extremely slow heating rates. For example a

Jan 1, 1961

By L. W. Rowe, S. S. Cole

A laboratory bench scale unit is described whereby finely divided chlorinatable residues are held for a short period by a restraining bed of a coarse-grained ore of comparable composition to permit &apos;flash&apos; chlorination of the high surface area powder. Results on chlorinating a finely divicled titania residue using rutile sand as the restraining bed showed that 92 pct utilization of chlorine and equivalent conversion of TiO2 in the fine residue to Ticl4 were attained with a dust loss of about 3 to 4 pct and an equivalent amount of ch1om&apos;nation of the titanium values in the restraining bed. THE chlorination of coarse grain titaniferous ores; i.e., 20 to 150 mesh sizes in fluidized beds is practiced commercially.&apos; One effective method found for the chlorination of a very fine material in a fluidized bed is the "flash" chlorination technique.&apos; By charging the finely divided charge material into the bottom of an established bed of optimum sized rutile, ilmenite, titanium beneficiate, or even an inert material such as sand, chlorination occurs before the dust can pass through the bed. The major requirement is that the bed be less chlorinatable than the charge. In this process, the charge consists of a finely divided (between one and 50 µ) titania residue obtained from a titaniferous ore beneficiation process. The charge is blended with a 20 pct wt basis addition of 20 to 60 .mesh size low volatile coke (to minimize hydrogen chloride formation) and fed to the bottom of the bed with the chlorine gas stream. A ball check device is used to permit entry of the charge and prevent the bed from discharging when the gas flow is cut off. The coarse pieces of coke serve to reduce compaction in the charge mixture and permit sat]sfactory feeding without plugging or jamming. Feed and gas rates are adjusted to provide for complete reaction of the charge material, with 90 to 95 pct chlorine efficiency, and maintain a gas flow in the range of 0.4 to 0.6 ft per sec. The equipment consisted of a 3-in. diam mullite tube enclosed in an electric furnace. A bow-shaped 3/4 in. Saran tube fed the charge from an air lock type of feed hopper (to prevent the escape of the chlorine gas into feed bin) to the bottom of the bed. The charge rate was set by providing a small awiliary hopper below the feed hopper where a weighed portion of the charge could be discharged in a specified time interval into the main chlorine stream. The chlorine rate was measured by a flowmeter. The reaction tube was filled with a bed of Australian rutile ore, this restraining bed serving two purposes; i.e., it rapidly heated the fresh charge to the reaction temperature and retained the particles in the bed for a sufficient time for selective chlorination to proceed. Due to the fine particle size and open, porous structure of the charge, this reaction was almost instantaneous, and, with proper adjustment of feed and gas rates, very little fine material passed through the bed without reacting and only a very small fraction of the restraining bed was chlorinated. A typical reaction procedure was to fill the bench scale fluidized bed unit with a 1 kg charge of Australian rutile ore analyzing 95.5 pct TiO2, 0.6 pct FeO, a beach sand deposit 40 to 150 mesh size. When fluidized, the bed was approximately 2 ft in depth. The ore was fluidized with a 1:l mixture of chlorine and nitrogen at a rate of 0.40 ft per sec. This gas mixture was necessitated by the comparatively short bed depth; a larger unit, with a deeper bed, could be fluidized with chlorine alone without significant free chlorine loss. Feed rate was 8.0 g per min of beneficiated ilmenite analyzing 85.0 pct TiO2, 2.6 pct FeO plus 1.6 g per min of 20 to 60 mesh petroleum coke. The

Jan 1, 1962

By R. M. S. B. Horta, W. T. Roberts, D. V. Wilson

Evaluation of results obtained by the Harris method for inverse pole figures is discussed. Two existing analyses and a new approach are compared. In the most frequently used analysis, different reflections are accorded equal weight. A second method involves area-weighting of reflections, mainly to take account of nonuniformity of Pole distribution. In the new method, reflections are weighted according to their multiplicity factors. Intensity ratios obtained from a sample of steel sheet by the three methods are compared. INVERSE pole figures offer a convenient method of depicting the proportions of material with various orientations referred to a unique specimen axis. The technique was first applied to fiber textures, but it is now frequently used as a description of texture in sheet metal, the unique axis in this case usually being the sheet normal. The Harris method&apos; for the determination of inverse pole figures has the great merit of simplicity. It suffers certain intrinsic limitations, but provided that these are borne in mind, the results allow useful conclusions to be drawn regarding the effects of texture components on certain mechanical properties. Much thought has been given to alternative methods of texture description, aimed at eliminating the ambiguity in the simpler methods.2"6 Refinements in the method of description are obtained at the price of increasing complexity, and it seems clear that a simple method, such as that of Harris, will continue to find application. The original derivation by Harris1 has been corrected by Morris7 and later by Mueller, Chernock, and Beck.2 Since the latter reference is more accessible, it is referred to more frequently than the work of Morris. Most users of the inverse pole figure method present their results on the basis of the treatment by Mueller et al., and it is the purpose of the present communication to comment on this treatment and to suggest a simple modification to the method of presentation which enhances the meaningfulness of the results. It will be convenient first to summarize the Mueller analysis as follows: The intensity Ihkl of a hkl reflection from a textured specimen is: Ihkl=CI0AL\Fhkl\2NhklPhkl [1] where C is a constant for a given sample I, is the intensity of the incident beam A is an absorption factor. With usual diffrac-tometer geometry and a flat specimen, A is inversely proportional to the density of the specimen. L is the Lorenz polarization factor Fhkl is the structure factor for the hkl reflection Nhkl is the multiplicity factor phkl is the fraction of crystals in the polycrystal-line aggregate with any particular (nkl) plane parallel to the surface. The corresponding equation for a specimen with randomly oriented grains is lR,hkl=CRl°ARL\fnkl\2Nnklk [2] the subscripts R indicating the factors that are different for the random sample, which will, in general, have a density different from that of the textured one. The intensity ratio is hkl , c A Phkl [3] IR,hkl CR AR PR The constants are eliminated by summing over all the measured reflections ^iRMi °R AR PR ljPkkl and from Eqs. [3] and [4] T IR,hkl The only unknown in Eq. [5] is Yj phkl and before this can be evaluated, it is necessary to adopt some form of normalizing procedure. Mueller et al. assume that when a large number, n, of reflections is used, the average value of p as defined by Eq. [6] below will be unity ? = £j**l3l [6] and the final result of their analysis follows: hkl Pkkl = i I&Jf~ [7] ±y* Ihkl ^ in,hki The normalizing procedure used by Mueller et al. involves an unrealistic assumption since it gives equal weight to each of the hkl reflections. Other normalizing procedures can now be examined with a view to selecting a more appropriate one. For an infinite number of reflections, the function p can be averaged over a continuous range of orientations covering the complete solid angle 4p by the integral relation: p =1 jp(a)dQ = l [8] This is the only exact procedure, but in practice a

Jan 1, 1970

By John B. Cook

Most of today&apos;s open pits take the form of conical-shaped excavations in the relatively flat or undulating terrain surrounding them. Ore is usually hauled uphill from the pit bottom by truck, rail, or skip to crushing, beneficiation, railroad loading, and stockpiling facilities near the edge of the pit. There are, however, a few orebodies whose locations in mountainous terrain causes special problems that generate different technical solutions. LAMCO&apos;s 10.5 million metric tons per year high-grade iron mine in the Nimba Range in tropical Liberia is probably one of the most extreme of these exceptions (Fig. 1). At Nimba there is a 2,000-ft difference in elevation between the upper- most levels of the orebodies and the start of the railroad below at the foot of the mountain. Haulage to the crusher is usually downhill at grades of up to 12% dictated by the topography of the rugged terrain. Ore is transported from the crusher down the mountainside by a conveyor sys- tem 10,096 ft long, whose average grade is 18%. The belt is braked by asynchronous generators, feeding power into the local network. The original crusher must be relocated at a lower elevation outside the final pit limit to allow that part of the mountain on which it now stands to be excavated. HISTORY The LAMCO Joint Venture is owned 75% by the Liberian-American- Swedish Minerals Co. (LAMCO) and 25% by Bethlehem Steel Corp. LAMCO is comprised of a 50% interest by the Liberian Government and 50% interest by Liberian Iron Ore Ltd (Fig. 2). This joint venture was launched to develop the high-grade iron ore deposits of the Nimba Range.

Jan 1, 1969

By Herman Dykstra, R. L. Parsons

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.&apos; THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy&apos;s Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. &apos;Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If&apos; .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951

By R. L. Parsons, Herman Dykstra

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.&apos; THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy&apos;s Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. &apos;Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If&apos; .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951