Search Documents

By Joseph T. Singewald

AMONG the genetically interesting iron ores of the United States are those of the St. Francis Mountains near Ironton and Iron Mountain, Missouri. They are specular hematite in porphyry. The Iron Mountain and Shepherd Mountain deposits are generally spoken of as veins and those of Pilot Knob, Cedar Hill, and Russell Mountain as bedded deposits. Iron Mountain has produced nearly 4,000,000 tons of ore and Pilot Knob about 1,500,000 tons. The output of the other three localities has been insignificant. Pilot Knob is no longer worked as an iron mine, but the rock is now being crushed for use as road metal. Iron Mountain has continued as an active producer. In 1927, the property was acquired by the M. A. Hanna Co. and an extensive prospecting and exploration campaign was undertaken, so that the ultimate production of Iron Mountain will doubtless greatly exceed the figures mentioned. Previous Theories of Origin Descriptions of the deposits and theories of their origin have been concerned mainly with Iron Mountain and Pilot Knob, which represent the two contrasting types of occurrence. The Iron Mountain deposits have been considered of close magmatic origin by Whitney' in 18.54, Geijer2 in 1915, and Spurr3 in 1927. Whitney referred to them as of eruptive origin. Geijer considered them fissure fillings at magmatic temperatures. Spurr called them vein dikes formed from ore magma which was essentially fluid iron oxide. Schmidt4 in 1872 ascribed them to fissure filling by thermal chalybeate waters of meteoric origin. Nason5 in 1892 advanced essentially the same explanation. One of the best

Jan 1, 1929

By J. R. Long, C. E. Armantrout, A. H. Roberson, T. R. Graham, R. S. Dean

THE general program of the Federal Bureau of Mines on the study of alloys made with electrolytic manganese has been extended to copper-manganese-aluminum alloys. The initial results of the work are presented in this report and concern a determination of the mutual solubilities of these three metals to form the alpha solid solution. The principal studies of this system, reported in the literature, have been focused on the Heusler alloys containing 10 to 30 per cent Mn with 9 to 15 per cent Al. These investigations concern the ferromagnetic properties and the determination of the crystallization contours, with little emphasis on the phases present in the solid state outside of this area. Heusler1 recognized a breakdown of the beta solid solution of the magnetic alloy, which he says decomposes into the alpha, beta, and gamma phases of the copper-aluminum system, somewhat modified by the presence of manganese. Rosenhain and Lantsberry2 examined alloys ranging up to 10 per cent Mn and io per cent Al but were unable to reach definite conclusions concerning the phase boundaries in that area. Krings and Ostmann3 made magnetic studies and examined the system by thermal analysis. They present contours of the liquidus surface but do not place any of the phase boundaries. Corson4 [ ] has drawn a hypothetical boundary for the alpha field with few data to support such construction. The binary systems involved have received considerable attention by numerous investigators, and their general features are well established. The hopper-aluminum system (Fig, I) has been thoroughly examined by Stockdale,5 who places the limit of solubility of aluminum in copper at 7.4 per cent at the liquidus and 9.5 per cent at room temperature. Smiths established the eutectoid temperature at 570°C. by a study of isothermal decomposition of the beta phase. The copper-manganese system (Fig. 2) has been recently reviewed by Dean and co-workers.7 Using high-purity alloys, they place the solubility of alpha man-

Jan 1, 1946

By H. A. Wriedt

Zinc-vapor pressures in equilibrium with bcc solutions of zinc and iron at 703.5", 757", 793", and 900°C were measured by an isopiestic method. Thermody-na?nic properties of zinc in these solutions were derived directly and those of iron with a Gibbs-Duhem integration. Zirzc and iron exhibit large, positive, and moderately temperature-sensitive deviations from ideality. Excess enthalpies and Gibbs energies are positive in the ranges studied. It has been establzshed in relating published phase equilibria to the Gibbs energies of forming solid solutions that positive departures from ideality also exist in the fcc solutions. Tlze stability, as defined by Darken, of the bcc solutions is small but positive at 703.5"C. Its decrease with decreasing temperature indicates that a miscibility gap may exist in the stable or metastable bcc solutions at lower temperatures. The changes of thertnodynanlic properties relatable to changes in magnetic state are minor compared to other changes due to alloying. The accuracy with which the relation between lattice parameter and composition of the bcc solutions is known has been improved. Zinc increases the lattice parameter of a iron by O.ZOZA per unit atom action zinc to at least NZn = 0.25. TO our knowledge, only Ciganl has previously investigated thermodynamic properties of the bcc solutions of the Fe-Zn system. He measured zinc-vapor pressures in the temperature range from about 320" to 360°C by an effusion technique. We have measured the vapor pressures of zinc in equilibrium with bcc solutions of zinc and iron at 703.5", 757", 793", and 900°C by an isopiestic technique. This investigation and the pursuant thermodynamic computations are part of a program for studying several properties of these solid solutions. In the appendix we report a remeasurement of the relation at room temperature between the lattice parameter and zinc content. Related work not reported herein embraces measuring the diffusivity2 of zinc in ferrite at 700" to 900°C and the Curie point temperature3 as functions of composition. EXPERIMENTAL METHOD AND APPARATUS Purified iron at a series of constant temperatures was saturated with zinc at a series of constant zinc-vapor pressures, the range of conditions being so limited that the equilibrated solid solutions were ferri-tic. The equilibration was conducted in long capsules, Fig. 1, that had been evacuated, gettered with titanium, and sealed. The purified iron, part of a vacuum-melted ingot of electrolytic-iron stock (Glidden Co. 1044-grade), had the impurity contents listed in Table I. It was used in the form of millings about 0.001 in. thick. A zinc of 99.999 pct purity (United Mineral and Chemical Corp.) and later one of 99.9999 pct purity (Cominco Products, Inc.) were used, with no observed difference in results. Sufficient zinc was placed in the capsule (separate from the iron and kept at a different temperature) to supply the vapor phase and the iron alloy, while leaving a pool of pure liquid zinc throughout the run. The imposed temperature of the residual zinc pool fixed the vapor pressure of zinc throughout the capsule; condensation was avoided at other locations by making this the lowest temperature in the capsule. Two common complications of nonisothermal systems were, we believe, effectively absent from ours: first, molecular species other than monatomic have no appreciable concentration in zinc vapor at 700" to 900~~: second, the development of zinc-pressure differentials in the capsules was suppressed by use of relatively large-diameter tubing and the absence of diluents in the zinc vapor other than desorbed impurities. The specimens and zinc pool were located at two isothermal zones (each ranging within 1°C over several inches) of a wire-wound, 37-in.-long vertical-tube furnace. The required temperature distribution was achieved by external adjustment of the currents in three electrically isolated windings. The temperatures of the isothermal zones were controlled independently by ther-

Jan 1, 1968

By F. G. Miller

Earlier studies of the mechanism of flotation showed that coal flotation products are contaminated by composite coal and mineral particles. The mineral matter of these particles has the dual effect of increasing the specific gravity and decreasing the coal surface of the particle available for bubble attachment. The present study showed that these particles are more loosely held in the froth column than the high grade coal particles, especially in the coarser sizes. Consequently, it is possible in the laboratory to produce a secondary concentration of the product by sprinkling the froth product with water to wash these loosely held particles from the froth during formation and recovery. To improve coal flotation, Bethlehem Steel Corp.'s research department has been working on a number of related problems based on both physical and chemical factors. The ideal is of course full understanding of the interrelated effects of both sets of factors. However, detailed studies1 of each physical factor have produced immediately usable results in commercial coal flotation and are providing the basis for our current work on the effect of chemical factors. Since these detailed studies provided the basis for the present investigation of froth sprinkling, their key results are highlighted here. For example, for a highly floatable bituminous coal, particle floatability was shown to be directly related to the percentage of coal surface available for bubble attachment and inversely related to particle mass. The data also underscored that in a given size range the addition of mineral matter to a particle makes flotation more difficult both by adding to the mass of the particle and by decreasing the amount of coal surface available for bubble attachment. Size also affects particle floatability in that as the size of particles of a given specific gravity increases, mass increases and the ratio between coal surface area and mass of the particle decreases. Flotation tests showed that any recovery of the coarsest particles in a coal feed is accompanied by progressively higher recovery of finer and finer particles. Consequently, fine high-specific-gravity particles are recovered with the coarse good coal. As a result, when the entire 14 mesh x 0 size range is processed by froth flotation, the good coarse coal is contaminated with high-specific-gravity fines. It is the small quantities of these fines that must be rejected from the froth product if flotation efficiency is to be increased above its present level. These principles and findings, together with work by Professor V. I. Klassen,' led to the froth sprinkling investigation reported in this paper. Working with Russian coals, Klassen had shown that froth sprinkling washes loosely held refuse particles from the froth and compensates for the water lost in the upper froth layers and increases its stability. As a result, the coal particles are retained in the froth and the loss of coarse coal from the froth is prevented. Decreases in ash of 0.1% to 1.5% were reported by Klassen. Working with a 28 mesh x 0 American coal, the author of the present paper found that single-stage flotation coupled with froth sprinkling produced a product comparable to that obtained by the more expensive method of roughing and cleaning.3 The investigation reported in the present paper expanded study of the effect of froth sprinkling to include each of the standard size fractions of a bituminous coal. The main findings are as follows: (1) For those size fractions between 48 mesh and 150 mesh, froth sprinkling produces a marked reduction in both sulfur and ash. In the case of those sizes coarser than 48 mesh, while this reduction is not as significant, the additional benefit of shorter retention time is still a plus for froth sprinkling. Furthermore, when froth sprinkling is coupled with normal fast flotation, the product purity approaches that obtained by starvation reagent feeding, i.e., gradual feeding of flotation reagent to float the recoverable coal from the pulp slowly in several increments. The specific results for these two size groups are: (a) 14 x 48 Mesh. By using froth sprinkling with fast flotation, product sulfur can be reduced to a level

Jan 1, 1970

By F. Keller, J. D. Edwards

The anodic treatment of aluminum presents problems of scientific as well as of commercial interest.l-3 Of particular interest is the fact that, during the anodic oxidation process, the oxide continues to form at the metal-oxide interface under any oxide previously formed. This has led to speculation as to the mechanism involved in the formation of the relatively thick oxide coatings that are used commercially for decorative or protective purposes. Furthermore, the ability of certain types of oxide coatings to adsorb dyes and other substances has stimulated research to determine the actual structure of these adsorptive coatings. It has been found that anodic oxide coatings on aluminum are composed essentially of aluminum oxide. They are formed by the action of oxygen ions penetrating to the metal surface during the electrolytic oxidation treatment. These coatings can be formed in a number of different electrolytes, as for example, those which contain sulphuric acid, chromic acid, oxalic acid, or boric acid. These are the commonest and most useful electrolytes employed commercially. Except for certain specific conditions of formation, the coatings in general have been found to be amorphous alumina, as far as can be determined by X-ray or electron-diffraction methods. Anodic oxidation processes can be arranged in three rather general classes if they are grouped in relation to the solvent action of the electrolyte on the coating.4 In the first class, the electrolyte has little or no solvent action on the coating that is formed. In general, coatings produced under such circumstances are nonporous and non-adsorptive. In the second class, the electrolyte exerts an appreciable solvent action on the coating. These coatings are porous and adsorptive. Finally, for the third class, the electrolyte tends to dissolve the coating about as rapidly as it is formed. This action produces electrolytic brightening or anodic polishing of the aluminum surface and at most leaves only a very thin film of oxide. The nonporous and nonadsorptive type of oxide film is represented by the coatings that are formed when solutions of boric acid are employed as the electrolyte. It is especially significant that these impervious films are formed in electrolytes that exert little or no solvent action on the coating. Where a boric acid electrolyte is used, the coating tends to form rapidly, with the result that the flow of current is soon reduced to substantially zero. This indicates that the growth of the coating has stopped. As a rule, the thickness of the coating is roughly proportional to the voltage employed for formation and the coatings of this type are exceedingly thin. In an electrolyte of this type, the coating has a high resistance when the aluminum is made anode; consequently, any current that may flow (leakage current) is very

Jan 1, 1944

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

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

Jan 1, 1970

By F. Keller, J. D. Edwards

The anodic treatment of aluminum presents problems of scientific as well as of commercial interest.l-3 Of particular interest is the fact that, during the anodic oxidation process, the oxide continues to form at the metal-oxide interface under any oxide previously formed. This has led to speculation as to the mechanism involved in the formation of the relatively thick oxide coatings that are used commercially for decorative or protective purposes. Furthermore, the ability of certain types of oxide coatings to adsorb dyes and other substances has stimulated research to determine the actual structure of these adsorptive coatings. It has been found that anodic oxide coatings on aluminum are composed essentially of aluminum oxide. They are formed by the action of oxygen ions penetrating to the metal surface during the electrolytic oxidation treatment. These coatings can be formed in a number of different electrolytes, as for example, those which contain sulphuric acid, chromic acid, oxalic acid, or boric acid. These are the commonest and most useful electrolytes employed commercially. Except for certain specific conditions of formation, the coatings in general have been found to be amorphous alumina, as far as can be determined by X-ray or electron-diffraction methods. Anodic oxidation processes can be arranged in three rather general classes if they are grouped in relation to the solvent action of the electrolyte on the coating.4 In the first class, the electrolyte has little or no solvent action on the coating that is formed. In general, coatings produced under such circumstances are nonporous and non-adsorptive. In the second class, the electrolyte exerts an appreciable solvent action on the coating. These coatings are porous and adsorptive. Finally, for the third class, the electrolyte tends to dissolve the coating about as rapidly as it is formed. This action produces electrolytic brightening or anodic polishing of the aluminum surface and at most leaves only a very thin film of oxide. The nonporous and nonadsorptive type of oxide film is represented by the coatings that are formed when solutions of boric acid are employed as the electrolyte. It is especially significant that these impervious films are formed in electrolytes that exert little or no solvent action on the coating. Where a boric acid electrolyte is used, the coating tends to form rapidly, with the result that the flow of current is soon reduced to substantially zero. This indicates that the growth of the coating has stopped. As a rule, the thickness of the coating is roughly proportional to the voltage employed for formation and the coatings of this type are exceedingly thin. In an electrolyte of this type, the coating has a high resistance when the aluminum is made anode; consequently, any current that may flow (leakage current) is very

Jan 1, 1944

By C. Bezemer, I. Havenaar

An investigation was carried out on the dynamic liltration behavior of drilling fluids. In a set-up consisting of a porous pipe through which the drilling mud was circulated, dynamic liltration rates and cake thicknesses were determined as a function of flow rate and filtration time. Eventually, a state of equilibrium was reached; i.e., both cake thickness and filtration rate attained constant values. Concurrently, dynamic filtration tests were carried out in a simple dynamic filter apparatus where mud flow was similar to that in a rotational vismmeter. In both arrangements the relation between equilibrium liltration data and rate of shear at the cake surlace was studied. Filter pipe and dynamic filter apparatus were lound to give identical relations. It was therelore concluded that the latter apparatus, which is suitable lor routine measurements on small amounts of mud, could be used to determine the dynamic liltration behavior of mud under borehole flow conditions. Some data are presented on the effect of filtration pressure and temperature on equilibrium filtration rate and cake thickness. INTRODUCTION In oil well drilling, the drilling fluid used generally forms a mud cake on the permeable sections of the borehole wall while filtrate is lost into the formation. As to the conditions under which this filtration process takes place, a distinction should be made between static mud filtration (mud is not circulating) and dynamic mud filtration (mud is circulating). In recent years several investigatorsl-5 have studied the process of dynamic mud filtration about which little had previously been known. This type of filtration is the subject of this paper which confirms several findings of previous investigators. It differs from already published works, however, in the following respects. 1. It contains a description of a simple laboratory apparatus developed for dynamic filtration measurements on small amounts of mud. 2. The effect of rate of shear at the cake surface on dynamic filtration has been studied. As a result of this study, dynamic filtration data pertaining to the borehole conditions of mud flow can be found from laboratory tests. 3. Some data are presented on the effect of filtration pressure and temperature on dynamic filtration. The work was limited to a small number of muds. No attempt was made to compare data on API filtrare loss with dynamic filtration rates. Previous workers have clearly shown that no correlation exists between static and dynamic filtration rates. EXPERIMENTAL ARRANGEMENT To establish whether the effect of mud flow conditions on dynamic filtrations is governed by the rate of shear at the cake surface, two different types of mud flow were used. In this connection the experiments were partially conducted by means of a porous pipe through which mud was circulated, partly in a dynamic filter apparatus which somewhat resembles a rotational viscometer. In both cases the flow conditions of the mud were well defined, thus making it possible to-calculate rates of shear at the cake surface. The filter pipe was 200 cm long with an internal diameter of 3 cm. Its perforated wall was covered on the inside with filter cloth on which the mud cake could form. Mud was pumped through the pipe at linear rates varying from 20 to 186 cm/sec. The required filtration pressure was obtained by means of a lack-pressure valve downstream of the filter pipe. In this arrangement the type of flow existing in the borehole is reasonably well simulated. The dynamic filter apparatus (Figs. I. and 2) consists mainly of two concentric cylinders with a clearance of about 0.5 cm. The outer cylinder (rotor) can be rotated at rates up to 1,000 rpm. The outer surface of the inner cylinder (radius about 2.5 cm) is covered with filter cloth or paper supported by copper gauze. The two cylinders are mounted inside a pressure

Jan 1, 1967

MINING; SOLUTION U.S. 4,066,297 - In the solution mining of uranium from an underground deposit of uranium ore using an alkaline leaching solution, precipitation of alkali metal salts is avoided by maintaining the Injected leach solution of an alkaline carbonate or bicarbonate and an oxidant undersaturated with respect to carbonates Undersaturation is effected by use of a water softening unit or by bleeding a portion of the recycled alkallne leach solution. A P. Spence, assigned to Atlantic Richfield Co Jan. 3, 1978 4 pp. (299-4). Filed: 6-1-76 (69 1,344)

Jan 1, 1979

By N. D. Shutler

This paper describes a numerical mathematical model of the steamflood process that depends on fewer restrictive assumptions than models previously reported. The solution, however, is obtained economically. Example calculations are presented that, on comparison with experimental results, tend to validate the model. Results that expose certain process mechanics are discussed. The model describes the simultaneous flow of three phases — oil, water and gas — in one dimension. It includes the effects of three-phase relative permeabilities, capillary pressure, and temperature- and pressure - dependent fluid properties. Interphase mass transfer of water-steam is allowed, but the oil is assumed nonvolatile and the hydrocarbon gas insoluble in the liquid phases. The model allows heat convection in one dimension and two-dimensional heat conduction in a vertical cross-section spanning the oil sand and adjacent strata. The hydrocarbon-steam gas composition is tracked, but the effect of gas composition on water-steam phase behavior is neglected. The model is solved numerically in three separate stages. The three-phase mass balances are solved simultaneously using Newtonian iteration on nonlinearities occurring in the accumulation terms. The energy balance is solved separately by noniterative application of the alternating-direction implicit procedure. Separate solution of the composition balance is accomplished by straightforward solution of the finite difference equations. The method of effecting nonsimultaneous, stable solution of the mass and energy balances is the key to the success of the model. INTRODUCTION Mathematical tools as well as laboratory and field experiments are necessary to help us understand the complex steamflood process. A mathematical model can expose process mechanics and show the relative importance of process variables, but this ability is often limited by restrictive assumptions. Most known models of steam processes,l-5 with the exception of the model of Gottfried,6 are "simplified" in that they involve analytic approximations and require many restrictive assumptions. The primary utility of these methods lies in the routine use as an aid in engineering design. By contrast, the comprehensive model presented here is numerical and requires far fewer restrictive assumptions. It finds its primary utility as a research tool. It serves as an aid in understanding the nature of the process, in interpreting laboratory experiments and in evaluating and developing simpler mathematical models for engineering design. The major reason why previously presented models have been confined to the "simplified" class is evidenced by the one published exception. In 1965, Gottfried6 presented a numerical model for the combustion process of which the steam process is a subset. The result is a comprehensive tool (though it neglects capillarity and two-dimensional heat conduction) that is troubled with convergence problems and that requires 2 to 3 hours of IBM 7094 time to complete a calculation. Though our present model does not simulate combustion, it does consider capillarity and two-dimensional heat conduction and it overcomes the convergence and computer-time problems. MATHEMATICAL DESCRIPTION OF STEAMFLOODING FLUID FLOW The equations employed to describe three-phase fluid flow are of a familiar form. Darcy&apos;s. law provides expressions for the velocities of the three phases (oil, water and gas), which, when combined with oil, water and gas mass balances, give the partial differential equations governing flow of the

Jan 1, 1970

By B. C. Allen

The interface energies of chronzium, molybdenunz. hugsten, and their solid-solution alloys Cv-35Re, MO-33Re, and UJ-25Re were studied at 0.6 to 1.0 of the absolllte liquidus ter)zpe,vature using fiz&apos;e )izethods. Liquid surface tension, yv , was deter mined clsing the pendant-drop and drop-weight methods. Results are, respectizlely, 1700, 2370, and 2480 +100 dynes per ct for the rhernium -containing alloys and essentially the same as tlwse reported for liquid chro)riln, trolybdenum, and tungsten. Average solid slrjace energy, rsv< xias ))zeasured using tlre fiber-extetlsion method. The ratio of ysS, the acerage high-angle grain-boundary energy, to ySV cclas jolnd fronz grain-bolzdary grooue angles fort)zed at the surface in an inert atrfizosphere. Absolllte iute?:face energies were deterawined using ?nultip/rase equilibria involzing suitable liquids of known surface tension (tin, silver). Interpretation of the experimented results in view of pvobable tenzperatzcre, orientation, and purity effects giz,e the follouling approximations in ergs per sq ctn: ysv (i2lo. Mo-33Re) - 2100, ySS (Mo, Mo-33Re) - 800, rr (defornzation twins in MO-33Re at 1200"C) - 800. ysV (Cr. Cr-35Re) - 2400. YSS (CY, Cr-35Re) - 1000. Probably Ylv- YSV- 2500 for tungsten and W-25Re, giving yss (It&apos;, W-25Re) - 900. The interface energies of solid and liqid ch?&apos;omiu?.z, nolybdenu?rr, and tungsten are not geatly aff&apos;ected by rhenium and therefore are not a ttlajor factor in the ductili zing rhenium effect in Croup VI-A metals. THE interface energies of the refractory Group VI-A metals, chromium, molybdenum, and tungsten, are not well-established. The objective of this investigation was to study the liquid surface tension, solid surface energy, and grain-boundary energy of these metals and compare them to those found for the bcc solid-solution alloys, Cr-35e,&apos; 0-33e,&apos; and -25e. Five techniques were used to measure interface energies in high-purity polycrystal rod, wire, and sheet at 0.6 to 1.0 of the absolute liquidus temperature. The alloys were chosen to see if there was any connection between interface-energy behavior and the ductilizing rhenium effecL4j5 EXPERIMENTAL WORK Materials. A description of the materials used is presented in Table I. Chromium rod was prepared by arc melting iodide process crystals supplied by Chromallo Cor., hot extruding, and warm swaging to 0.63-cm-diam rod.6 The sheet was prepared by rolling as-extruded rod to 95 pct reduction in area from a hydrogen furnace at 800" to 900°C and surface grinding off 0.02 cm from each face. Cr-35Re rod was prepared by arc melting sintered rhenium powder and iodide chromium crystals, warm rod rolling to 50 pct reduction in area in cans, and swaging to 60 pct reduction in area at 1100" to 1200°C. Some of the rod was warm-rolled to sheet and then surface-ground. Portions of swaged chromium and Cr-35Re were further reduced by swaging and drawing to 0.013-cm-diam wire by the General Electric Co. Mo-33Re and W-25Re rod, sheet, and wire were provided by Chase Brass and Copper Co. The molybdenum sheet consisted of two lots, both essentially the same except for the carbon content. Liquid Surface Tension. The liquid surface tension of Cr-35Re, Mo-33Re, and W-25Re was measured by a combination of pendant-drop and drop-weight methods using techniques already decribed." Following out-gassing, molten drops were formed on the ends of centerless-ground Mo-33Re and W-25Re rods by electron bombardment at 5 x 106 mm. Similar drops were formed on outgassed Cr-35Re rods by induction heating under 1 atm of 99.995 pct Ar. Solid Surface Energy. Solid surface energy was measured by conducting microcreep experiments on molybdenum, Mo-33Re, chromium, and Cr-35Re wires at 2350°, 2306, 1550°, and 180O°C, respectively. In preparation, gage marks -2.5 cm apart and -0.001 cm deep were circumferentially scribed on the wire with a razor blade. Weights of the wire material were then attached. Five to seven reasonably straight wires were hung in a container made out of the wire material. The free end was placed through a small hole in the removable top and secured by bending a small portion 90 deg. The containers not only tended to provide vapor-solid equilibrium for the wires but also protected them from gaseous impurities. They were nominally 2.5 cm in diam by 5 cm high and were made from extruded chromium rod, Cr-35Re arc casting, molybdenum bar stock, or welded Mo-33Re sheet. After deg re as ing, the assembly was outgassed at a relatively low temperature to 2 x 10"5 mm and then recrystallized 2 to 8 hr at the creep temperature in a rhenium-element resistance furnace. The static argon atmosphere was gettered by tantalum radiation shielding. Specimen temperature was measured optically to 25"C using calibration with known melting points and blackbody conditions. The wires generally developed a stable bamboo-type structure according to Fig, l(b), (c), and (d) and retained their gage marks [upper portion of Fig. l(d)]. One or two of the weights were clipped off to provide a low load for the creep anneal. To minimize the possibility of bending or breakage, the wires remained attached to the top of the annealing container which was held to keep the wires vertical. The distance between gage marks was

Jan 1, 1967

By R. N. Thurston, E. A. Nesbitt, G. Y. Chin

A general relalionship between the amount of glide shear (due to slip) and the macroscopic shape change has been developed. Since the deformation can be large, finite strain analysis is employed. In this treatnzent, the shape change is expressed by a deformation gradient matrix, F - [axi/aXj1, where X and x refer to the initial and final positions of a particle. This matrix is readily evaluated in terms of the amounl of glide shear, a , and the direction cosines of the slip-plane normal and the slip direction for a single active slip system. In the case of deformation from slip on several systems, the product F(a) of the several deformation gradient matrices is first calculated. Then, by assuming that the final configuration is reached by a long series of small shears of magnitude a, occurring more or less alternately in the several slip sys- IN problems of crystal plasticity, it is often necessary to relate the amount of glide in the operating slip systems to the macroscopic strain components. This relationship is relatively simple when only small strains are considered. In this case the separate strain contributions from several slip systems are additive. Such a procedure is incorrect in the case of large plastic deformation. To the authors&apos; knowledge, no general treatment of the latter problem has appeared in the literature. In specific instances, Mark, Polanyi, and schmid3 have derived the relationship between glide and axial elongation during tensile pulling of a single crystal when only one slip system is active. A similar relationship for duplex slip was worked out by v. Goler and Sachs.4 Taylor and Elam5 have also studied in detail problems of large plastic deformation. Due to uncertainties as to the slip systems at the time, terns, the final deformation gradient matrix can be obtained. Mathematically, the problem is reduced to finding the limit of F(dN as N — m while a — 0 in such a way that the Product Na = 0, a ftnite conslant designatzng the accunzulated amount of shear. It turns out that this limit is simply , where F, is the matrix whose elements are the coefficients of a in the F(a) matix. Application is made of the present treatment to a fcc crystal of Perlrzalloy compressed on the (110) plane and contrained to elongate zn the [112] directzon. In addition, it is shown that one rnay readily obtain from the general analyszs the well-known formulas relating elongation, amount of glide shear, and amount of lattice rotation for crystals deforming by single and double slip under tension. however, they were mainly concerned with proving that slip occurs on {111}(110) systems in fcc metals. The present general treatment arose out of recent investigations of magnetic anisotropy in cold-worked Fe-Ni alloys6 as well as some related problems in strength anisotropy of single crystals.&apos; Detailed application of the general analysis will be provided in the case of a Permalloy single crystal compressed on the (110) plane and constrained to elongate only in the [112] direction. It will be seen that, as expected, the small strain approximation leads to significant errors after moderate straining. It will also be shown that the present general treatment yields the formulas of Mark et al. and of v. Goler and Sachs in those specific cases. GENERAL CONSIDERATIONS Consider a homogeneous deformation in which a material point initially at moves to both positions being referred to the same set of Cartesian axes. Such a deformation can be

Jan 1, 1967

By Ernest F. Burchard

The Pao deposits of iron ore are in the State of Bolivar, Venezuela, near Upata, about 30 miles south of the Orinoco River and 90 miles east of Ciudad Bolivar. They were discovered 4 or 5 years ago and have not yet been developed to the stage of mining. In March-April, 1929, the writer devoted a week to the examination of these deposits. The maps, drill records and chemical analyses presented herewith were derived from original sources, and the field photographs were made by the writer. The petrographic and metallographic data are from notes by C. S. Ross and M. N. Short, of the U. S. Geological Survey, who have kindly collaborated in these special phases of the study. The Iron Obebodies Topography and Accessibility The deposits are in a region of low mountains representing northern outliers of the Guayanan Highlands. These mountains reach altitudes of 1800 to more than 2000 ft. above sea level, as determined by aneroid barometer. The iron ore is found on mountain crests that reach altitudes of 1915 to 1930 ft., and may be followed down precipitous slopes for 200 • to 300 ft. or more. The greatest vertical range in outcrop is nearly 400 ft., and drilling has indicated that the ore extends to this depth below the surface in one place, at least. The superior hardness of the iron ore accounts for the preservation of the local mountain crests with the intermediate valleys eroded in softer granitic rock. The deposits of iron ore are reached by a private road for automobiles branching northeast from the Ciudad Bolivar-Upata road, a distance of about 12 miles. This road was used in 1928-29 by trucks for transportation of food and camp supplies, gasoline, drilling equipment, etc., but was not suitable for the transportation of iron ore. Between Ciudad Bolivar and the iron-ore area the western two-thirds of the automobile road passes through nearly level, open country known as the savanna, largely used for grazing cattle. About 96 km. (60 miles) east of Ciudad Bolivar at

Jan 1, 1931

By Paul Merica

THOSE of you who, like myself, have had the privilege of hearing previous Howe Memorial Lectures will note, and I hope without too much disaffection, that for the first time in this series of lectures iron has been chosen as the subject. I feel sure that the predominant choice by previous lecturers of subjects relating to steel has been principally in deference to the well-known interest this subject held for Professor Howe, and that their choice has not reflected any lack of respect for the material, iron. Nevertheless, in the face of this record, I feel some diffidence in introducing this "Cinderella" of materials into the circle of distinguished subjects hitherto treated by Howe lecturers, which is relieved only by the fact that I hope to be able to present some evidence that this Cinderella has at last undergone some of those magical trans-formations which render her acceptable even in the society of her dis-tinguished sister, steel. Fortunately I can cite the old master himself in support of my choice of subject and recall to you that Professor Howe by no means overlooked or neglected cast iron in the long course of his metallurgical meditations. Indeed his attention was early arrested by a metallurgical aspect of cast iron which he believed and we still recognize as of the most funda-mental importance-namely, the genesis of graphite in gray iron and the means for the control of its quantity and distribution. He made some comments over 20 years ago regarding the nature of J. E. Johnson&apos;s oxygenated iron which indicate that even then his thoughts were taking directions which today we are finding it useful to explore in connection with the theory of high-strength iron. Had we had in the subsequent years the continuing benefit of Professor Howe&apos;s illuminating consider-ation of this important subject, our progress in the improvement of cast iron would undoubtedly have been even more rapid.

Jan 1, 1937

By J. M. W. Mackenzie

The interactions between oil drops and mineral surfaces have been examined for the cetyltrimethylammonium bromide (CTAB)-quartz and sodium dodecyl sulfate (SDS)-hematite systems. The results have been interpreted using a electrical double layer interaction model, a "sweeping-penetration" model based on steric interactions and a free energy of adhesion model. The adsorption of long-chain ions at the mineral interface was also studied and found to follow a Fuerstenau-Modi model of combined Coulombic and chain-chain interaction forces. The basis of the emulsion flotation process and the feature which makes it differ from conventional flotation is: encounter, followed by adhesion, between oil drops and mineral particles which are suspended in an aqueous environment. The manner in which the adhesion between these two phases contributes to the flotation process may be through the formation of min-eral-oil drop-mineral linkages, mineral-oil drop-air bubble linkages, or a combination of both. Possible flotation mechanisms derived from these adhesion processes are shown schematically in Fig. 2 while Fig. 1 shows an example of the first mechanism for quartz particles which have been floated using a cetyltrimethyl-ammonium bromide (CTAB)-Nujol emulsion as collector. In Fig. 2 the three phases, mineral, oil, and air, for convenience have been shown as having approximately equal size. This would not always be the case in an actual emulsion flotation operation. For the flotation of sand-sized particles of approximately 0.5 mm diam, the mineral particles would be at least ten times the size of the largest emulsion drops present in the system, while in the flotation of —10u slime particles, the mineral phase may be of the same order of size as the oil phase. The air bubbles in a flotation cell are normally about 1 to 2 mm in diam and therefore, in general, would be larger than the oil drops. Although the relative size of the three phases does not alter the basic adhesion processes between oil drops, mineral particles, and air bubbles, which depend on the physical and chemical properties of the interfaces involved rather than on the extent of these interfaces, it does alter the degree of deformation of the oil drops and air bubbles which will occur on adhesion. For mechanisms 1 and 2 to take place, the oil-water and mineral-water interfaces must be partially replaced by a mineral-oil interface while in mechanism 2 there is the additional requirement that the oil-water and air-water interfaces must be replaced by an oil-air interface. In view of the importance of the oil drop-mineral particle encounter and subsequent adhesion stages to both of the emulsion flotation mechanisms discussed, a study has been made of some physicochemical factors which control these processes for Nujol-cetytri-methylammonium bromide (CTAB)-quartz and Nujol-sodium dodecyl sulfate (SDS) -hematite systems. Experimental Electrokinetic Measurements.—The zeta potential of mineral particles was measured using a streaming potential technique similar to that described by Fuerste-nau.1 The streaming potential was measured using a multirange electrometer while the plug conductivity was measured using a 50 cycle sec1 bridge. Measurements were made on 28 x 65 mesh particles. Contact Angle Measurements.—A captive bubble technique was used to measure the contact angles for air bubbles and oil drops. The angles were measured on fracture faces of unmounted minerals which proved easier to clean than polished mounted specimens and should more closely represent a flotation system. Reagents: These have been described in a previous publication.0 Minerals: Rock crystal quartz from the New England district, N.S.W., was used for contact angle and streaming potential measurements. The quartz was hand

Jan 1, 1971

By D. B. Robinson, C. A. Macrygeorgos, G. W. Govier

Experimental data have been obtained on the volurrletric behavior of ternary mixtures of methane, hydrogen sulfide and carbon dioxide at temperalures of 40°, 100" and 160°F up to pressures of 3,000 psia. The results indicate that the compressibility factors for this system do not agree with compressibility factors for sweet natural gases at the same pseudo-reduced conditions. The deviation increases as the temperature and methane content decrease. Discrepancies of up to 35 per cent were observed. A careful analysis has been made of the existing pUrblished data on compressibility factors for binary systems containing light hydrocnrbons and hydrogen sulfide or carbon dioxide. It has been found that the deviation of actual from predicted compressibility factors for methane-acid gas mixtures is a function of the methane content and the pseudo-critical properties,.v of the mixture. The ratio between actual compressibility factors for methane-acid gas mixtures and compressibility factors for sweet natrlral gases at the same pseudo-reduced conditions has been currelated over a range of pP,, from 0 to at least 7 arid a range of pT, from about 1.15 to at 1east 2 0 with an error not exceeding 3 per cent and over most of the range within I per cent. The validity of the correlation for mixtures containing appreciable hearvier hydrocorbons has not been fully established, but it is shown to be preferable than the use of a corretation based only on hydrocarbons. INTRODUCTION Although a relatively accurate method for predicting compressibility factors of pure materials is provided by charts based on reduced properties and the assumption that the compressibility factor is a unique function of T P and z the determination of the correct values of compressibility factors for gas mixtures is somewhat difficult. Two general methods of dealing with gaseous mixtures have been proposed. The first assumes a direct or modified additivity of certain properties of the mixture in terms of the properties of the individual components. Examples of this method are based on the familiar laws of Dalton and Amagat. The second method averages the constants of an equation of state applicable to the pure components. Both of these methods are of limited value in engineering calculations because the first usually provides reliable answers only over narrow ranges of pressure and temperature and the second is cumbersome to handle. In petroleum engineering practice accurate estimations of the volumetric behavior of natural gases arc frequently required. To fulfill this need, several generalized compressibility charts have been developed.&apos; &apos; Of these, the one prepared by Standing, el al is widely used at present. In the construction of charts of this type a third method for dealing with mixtures has been followed. It is based on correlation of pseudo-critical properties as outlined by Kay and calculated from the critical properties of the individual components in a mixture. Although these charts provide relatively accurate information on the compressibility of dry or wet sweet natural gases, they are less reliable when used for gases containing high concentrations of hydrogen sulfide or carbon dioxide or both. Thus, an experimental program, although time consuming, is the best means now available for the determination of the volumetric behavior of sour or acid gas mixtures. An increased interest in the behavior of these gas mixtures, particularly in connection with some of the fields in Western Canada where the acid gas concentration of the reservoirs may be as high as 55 per cent and where hydrogen sulfide alone may be as high as 36 per cent, provided the incentive for this study. It was the purpose of the investigation to determine the volumetric behavior of selected mixtures of methane, hydrogen sulfide and carbon dioxide over a range of temperature from 40" to 160°F and at pressures up to 3,000 psi. EXPERIMENTAL METHOD The apparatus used in this investigation was basically the same as that described by Lorenzo.&apos;" The amount of each pure component used in preparing the gas mixtures was measured over mercury in a glass-windowed pressure vessel. The pure components were then transferred individually in the desired amounts to a second glass-windowed pressure vessel where the volumetric behavior of the mixture was determined. Volume was varied by mercury injection or withdrawal. The capacity of the cell was about 125 cc. Temperatures in the cells were measured with copper-constantan thermocouples and a Leeds Northrup semi-

By Henry M. Howe, Joseph Winlock, Francis B. Foley

This investigation was undertaken for the purpose of determining, in a systematic way, the effect of the rate of cooling of steels, heated to above the transformation range, on their various mechanical properties. This paper gives the results obtained in the investigation of the heat treatment of carbon steels carried out under the auspices of the Committee on The Heat Treatment of Carbon Steels of the Engineering Division, National Research Council. The objects of this committee, of which Dr. Henry M. Howe was chairman, were to increase knowledge of the influence of heat treatment on the mechanical properties of carbon steel. The mcthod pursued was to subject specimens of steel of the carbon contents most used for engineering work, 0.34 per cent., 0.52 per cent., and 0.75 per cent., to various heat treatments, and to test the mechanical properties thus induced. The steels were the same as were used by the committee on fatigue phenomena. This project was financed jointly by the U. S. Bureaus of Mines and Standards; industries and university laboratories have cooperated. A large part of the work was done at Doctor Howe&apos;s private laboratory. The committee consisted of F. B. Foley, H. AT. Boylston, H. C. Boynton, G. K. Burgess, H. J. French, J. H. Hall, Zay Jeffries, R. S. Johnston, F. C. Langenberg, H. F. Moore, E. H. Yeirce, H. S. Rawdon, W. E. Ruder and 11. J. Stagg. The work was undertaken at the suggestion of Doctor Howe, who thought, at the outset, that steel in the sorbitic state might have properties which, for most engineering purposes, would be better than those of the usual or pearlitic state. By a systematic variation in the heating

Jan 1, 1923

By Joseph Winlock, Francis B. Foley, Henry M. Howe

This investigation was undertaken for the purpose of determining, in a systematic way, the effect of the rate of cooling of steels, heated to above the transformation range, on their various mechanical properties. This paper gives the results obtained in the investigation of the heat treatment of carbon steels carried out under the auspices of the Committee on The Heat Treatment of Carbon Steels of the Engineering Division, National Research Council. The objects of this committee, of which Dr. Henry M. Howe was chairman, were to increase knowledge of the influence of heat treatment on the mechanical properties of carbon steel. The mcthod pursued was to subject specimens of steel of the carbon contents most used for engineering work, 0.34 per cent., 0.52 per cent., and 0.75 per cent., to various heat treatments, and to test the mechanical properties thus induced. The steels were the same as were used by the committee on fatigue phenomena. This project was financed jointly by the U. S. Bureaus of Mines and Standards; industries and university laboratories have cooperated. A large part of the work was done at Doctor Howe&apos;s private laboratory. The committee consisted of F. B. Foley, H. AT. Boylston, H. C. Boynton, G. K. Burgess, H. J. French, J. H. Hall, Zay Jeffries, R. S. Johnston, F. C. Langenberg, H. F. Moore, E. H. Yeirce, H. S. Rawdon, W. E. Ruder and 11. J. Stagg. The work was undertaken at the suggestion of Doctor Howe, who thought, at the outset, that steel in the sorbitic state might have properties which, for most engineering purposes, would be better than those of the usual or pearlitic state. By a systematic variation in the heating

Jan 1, 1923

By John C. Martin

The discovery of a new gas reservoir demands that the planning of a sottnd well-spacing program be initiated early in the development stage. It is the purpose of this discussion to illustrate by actual field examples the application of basic well-spacing principles, previously developed for oil reservoirs, to the problem of well spacing in natural gas fields. These studies are presented for the field use of geologists and engineers who are concerned with the initial planning of the proper development of the newly discovered gas reservoir. INTRODUCTION The phenomenal growth of a vigorous natural gas industry emphasizes the increasing importance of natural gas as a source of energy, fuel, and raw materials to our nation&apos;s economy. Since 1945 marketed production of natural gas for the U. S. has increased 21/2 times to a record high of 10.6 trillion cu ft during 1957. As major participants in the gas industry, we share an added interest to develop and produce our natural gas reserves with constantly improved efficiency. The subject of well spacing is vitally important to the gas industry, for the well itself plays a significant role in the development of the natural gas reservoir and in control of the recovery process. Maximum utilization of wells is an integral part of sound conservation practices. The discovery of a new gas reservoir demands that careful choice of well location and well spacing be made and that the planning of a sound well spacing program be initiated early an the development stage. With the drilling of the first development well, efforts of the geologist and engineer must be directed toward acquisition of adequate technical evidence upon which a firm recommendation for a spacing program may be based. With this technical appraisal as a foundation, operators and state regulatory agencies jointly can go far in providing a framework for sound development of gas fields to achieve a program of conservation that avoids the unnecessary well. WELL-SPACING CONCEPTS Through laboratory and field investigations of the mechanism of the recovery of oil and gas, of fluid behavior, and of effective control of reservoir and well, a crystallization of ideas regarding reservoir behavior has emerged as a well-developed technology. Associated with a better understanding of the fundamental principles underlying reservoir and well behavior has been the growth of concepts concerning the role of wells and their spacing in the development and operation of an oil or gas reservoir. In addition to serving as outlets for the withdrawal of fluids from the reservoir, wells are recognized as having two other important functions: (1) providing access to the reservoir to obtain information concerning the characteristics of the reservoir and its fluids, and (2) serving as a means by which the natural or induced recovery mechanism may be effectively controlled. Beyond a minimum number of wells required to fulfill these two functions, additional wells will not increase recovery. With particular emphasis upon well spacing in oil reservoirs, many studies of the well spacing-recovery relationship have evolved the concept that the ultimate oil recovery is essentially independent of the well spacing.&apos; These fundamental concepts are no different when regarding the role of wells and their spacing in the natural gas reservoir. They are equally applicable to the consideration of well spacing in gas reservoirs. For the gas reservoir, the problem of well spacing then revolves around the question of drainage and the degree or extent to which a well may drain gas from its surrounding reservoir environment. Theoretical and Experimental Work During the past 30 years, theoretical and experimental work carried on to study the physical principles involved in the flow of fluids through porous media has shed light upon the matter of drainage. Fundamental mathematical equations have been derived to describe the mechanism of flow of oil and gas through porous rocks. With the recent advent of high-speed digital computers, attempts have been made with mounting success to develop solutions, employing numerical techniques, to mathematical expressions that describe more rigorously the physical behavior and mechanism involved in the unsteady-state flow of compressible fluids, such as a gas, through porous rock. In 1953, Bruce, Peace-man, Rachford and Ricez published a stable numerical procedure for solving the equation for production of gas at constant rate. The results of these calculations are significant with respect to this matter of drainage, for they indicated (1) that depletion of the gas reservoir resulted in a drop in pressure at the extremity of the

By H. H. Rachford

Numerical solutions of immiscible flow problems in which dispersive effects of capillarity are dominated by convection require excessively fine grid spacing with attendant high computing costs. The use of coarser spacing reduces cost but often produces oscillation or undue dispersion associated with displacement fronts, A numerical formulation is proposed here which should be applicable to two-dimensional flow problems. it is in part analogous to an approach previously tested for miscible systems. The convective transport is approximated using a change of variables to yield a coordinate system moving approximately with the local characteristic velocity. The capillarity-induced dispersive terms in the differential system describing the process are approximated with respect to a fixed coordinate system by the usual implicit formulation. One-dimensional tests of the procedure yielded results in which the saturation profiles tended smoothly to the zero-capillary pressure solution as the ratio of viscous to capillary forces was successively increased in a sequence of calculations. This contrasted favorably with solutions by other numerical procedures which would require attendant grid refinements to approach the zero capillary pressure results. INTRODUCTION Numerical solution of displacement problems has until recently relied on applying methods developed primarily for transient heat-flow problems. Such problems are classified as parabolic in type, and where the heat transport is purely by diffusion their solutions are characterized by a high degree of smoothness. It is not surprising, therefore, that for approximating these solutions available finite difference methods are quite adequate. In flow problems the transport is partly by diffusion, partly by convection or flow. Although the problem remains of parabolic type because the dispersive effects of capillary forces or diffusion play some role in every displacement, at high flow rates the problem is dominated by convection, and solutions tend toward those of equations of the hyperbolic type. Solutions of hyperbolic problems are characterized by the translation of fronts, or discontinuities, that may progressively increase in sharpness. Numerical methods for treating parabolic problems become less and less satisfactory as displacement rates increase and the role of dispersion due to concentration or capillary pressure gradients becomes small relative to transport due to flow. In computation the difficulty manifests itself as an error associated with the grid size chosen. 1-6 In summary, if the heat-flow type approximations are to include the terms arising due to convection, one of several choices may be made: (1) an upstream (to the direction of flow) approximation for the convection terms may be used; (2) a centered-in-distance (CID) approximation may be used; or (3) a recently developed approximation based on the theory of oscillation matrices may be chosen.6 The last appears to have significant promise for one-dimensional flow problems; its extendibility to two or three dimensions is an open question. In either of the first two approaches, a suitably small ratio of v&/D must be maintained, where v is the velocity, & is the grid spacing and D the effective dispersivity in the direction of flow. In the first choice, the approximation of the convective part is only first-order correct and errors introduced appear as a numerically induced dispersivity of magnitude proportional to v?x. In the CID choice, the approximation can be second-order correct, but the difference formulation fails to satisfy the maximum principle unless a condition on v?x/D is met. Practically, this means that for high flow rates oscillatory solutions may result in the neighborhood of a front unless exceedingly small grid intervals are taken. While the procedure proposed by Stone and Brian4 permits a less severe limitation to be placed on this ratio, ultimately the flow rates increase relative to the dispersivity the oscillation obtains. Further, extensions of their approach to higher dimensional systems may be attended by considerable

Jan 1, 1967