Search Documents

By L. J. Kashar

By means of a continuum-mechanics approach, the following equation for calculuting the pressure for axisymmetric cold forging has been derived: where P is the extrusion pressure; K and n are the strength coefficient and the strain-hardening exponent of the steel, respectively; a and b are strain factors, primarily dependent on the die design; and R is the extrusion ratio. Data from compression and backward-extrusion tests performed at a crosshead speed of 0.5 in. per min on a nonwork-hardening material pure lead) and on a work-hardening material (1018 steel) indicate that the predicted extrusion pressures obtained from the derived equation agree very closely with the actual extrusion pressures. Therefore, the equation can be used to predict backward-extrusion pressures from standard mechanical-test data. The decision to cold-forge a part is based on several considerations, most of them related to economics. However, the production feasibility is usually determined by the pressure needed to extrude the part and the availability of tooling that can repeatedly sustain this pressure for over several hundred thousand cycles. A major contribution to the cold forging of steel would result from the development of steels that could be forged at reduced pressures and from the development of techniques that permit the accurate prediction of the pressure necessary to cold-forge the part. The ability to predict extrusion pressures is necessary to properly design the tooling required for the cold-forging process. Attempts to theoretically predict the cold-forging pressure, especially for work-hardening materials, have been unsuccessful in previous investigations of the cold-forging process.1-3 Many empirical relations have, therefore, been proposed relating the extrusion pressure to the chemical composition,' the extrusion ratio,2 the tensile strength,273 and hard- ness, 2 and the yield strength3 of the material being extruded. However, these empirical relations are valid only for the specific conditions used in their development, and, in some cases, the empirical relations had a scatter band greater than * 10 pct. A study of the extrusion process was undertaken with the objective of developing predictor equations for extrusion pressure by using a continuum-mechanics approach to the problem of metal flow. The derivation of such an equation is summarized in this paper, and data are presented that confirm the validity of the derived equation. DERIVATION OF THE EXTRUSION EQUATION If the extrusion of a material through a die is considered in terms of a work balance, it has been shown that it is possible to derive an equation relating the work necessary to move the material through the die to the deformation work undergone by the material. The work done, W, in moving the material is given by where F is the extrusion force, I, is the distance the material was moved, A, is the original cross-sectional area of the material, P is the extrusion pressure acting on A,, Ap is the cross-sectional area of the extrusion punch, and Pp is the punch pressure acting on Ap. The deformation work, WD, is given by where V is the volume of metal deformed, ay is the flow stress of the metal, and E is the total strain undergone in the extrusion. If there are no frictional losses, the work done, W, is equal to the deformation work, WD; therefore, Since the volume of a metal does not change during deformation,

Jan 1, 1968

By R. H. McLean

Velocity, kinetic energy and shear in crossflow beneath three-cone jet bits may influence cleaning of the bottom of the borehole and the teeth of the bit. Laboratory investigation shows that each of these parameters is a function of the diameter of the borehole and the product of the volume rate of flow and velocity through the nozzles (QVn). Increasing QV, or decreasing the diameter of the borehole increases each parameter. These functions provide means of predicting the magnitude of each parameter and of scaling the cleaning forces. INTRODUCTION In drilling operations using conventional jet-type rock bits, the impinging jels create an important flow mechanism. Called crossflow, this flow mechanism originates in the impact area of the jets, spreads across the bottom of the hole and supplies the principal source of energy to clean the teeth of the bit and most of the bottom.' Besides providing the means of cleaning cuttings from the bit and the hole bottom, crossflow may also have other, less direct, effects on the rate of penetration. The shear stress generated on the bottom by crossflow will inIluence the thickness and permeability of any filter cake of mud solids or crushed material which forms on the bottom.' These factors may affect the rate of penetration. A previous pulblication introduced some fundamental concepts of crossflow.' Crossllow was shown to occur in a thin layer adjacent to the bottom, and to cover the bottom completely. The maximum velocity in the crossflow above any position on the bottom was found to be directly proporlional to the square root of the product of the volume rate of flow and velocity through the nozzles—the jct QV,,—and inversely proportional to the diameter of the borehole. This information indicated that the effectiveness of crossllow in scavenging the bottom can be improved by maximizing the jet QV. The investigation reported herein amplifies the definition of crossflow. Complete velocity profile data above a representative position on the bottom are analyzed. These data better illustrate control of the capacity of crossflow to scavenge the bottom, and also relate shear stress on the bottom to known, controllable parameters. The conclusion reached in the previous publication that maximization of the jet QV. produces the maximunl cleaning beneath current jet bits is unchanged by these new data; rather, it is strcngthened. Data presented here show that the kinetic energy flux above a representative position on the bottom is maximized by maximizing the jet QV. The shear stress on the bottom will also be shown to be maximized in the same manner. Since the functions relating the jet QV. to velocity, shear stress and kinetic energy also involve the diameter of the borehole, means of equating, or scaling, these quantities in different sizes of boreholes will be illustrated. EXPERIMENTAL EQUIPMENT AND TECHNIQUE JET BIT MODEL Data were recorded from the same laboratory model as used in the aforementioned investigation of the flow around a jet bit.' The model consisted of a 43/4-in. three-cone, jet-type rock bit in a smooth, flat-bottomed borehole constructed of lucite. The bit had a shape and nozzle placement closely resembling larger three-cone jet bits commonly used in field operations. Fig. 1 illustrates the impact area on the bottom of a jet from this bit. Details of the orientation of the jets may be found in the previous publication. TECHNIOU'E OF MEASUREMENT Measurements of the crossflow were made by inserting a very small Pitot tube through the bottom of the simulatetl boreholc. Extreme thinness of the layer of crossflow necessitated accurate measurements of the height of the Pitot tube above the bottom to achieve close definition of the velocity profile. A cathetometer, which could be read to the nearest 0.005 cm, was used to make this measurement.

Jan 1, 1966

By Margaret Dienes, Karl L. Fetters

It is recognized that ingot structure is important in determining the quality of finished steel. Such elements of ingot structure as the size, shape and distribution of primary crystals; the size, shape and distribution of shrinkage cavities; the degree and kind of segregation (both macro and micro); the distribution of nonmetallic impurities and gas cavities are each the result of solidification phenomena. 'These structural features vary with variations in ingot practice. Their control requires a knowledge of their correlation with the variables of ingot practice. The structures of full-sized commercial steel ingots have been studied extensively by the Steel Ingots Committee of the Iron and Steel Institute.l This work has contributed much to our knowledge of ingot structure and the factors that control it. These investigations of necessity have been limited in scope because of their costliness and the great amount of time required. Studies of the structure of small steel ingots have not been very satisfactory because of the rapid rate of freezing of the ingots. Their structures are thus difficult to resolve and to interpret. Medium sizes of steel ingots in the ranges of 150 to 500 lb. do offer opportunity for study, but their cost also becomes greater and their use less expeditious. Experimental models for the study of steel-ingot structure have been made of nonferrous metals. Gathmann2 used lead-tin alloys and reported that they were satisfactory for determination of the location of the major shrinkage cavities and the Liplanes of weaknessv resulting from the pattern of crystallization. Certain features of mold design have been illustrated by these results. Bezdenezhnihk3 has used copper and aluminum ingots similarly and has shown the effect of mold-wall thickness on the depth of columnar grains. These studies of nonferrous models have been useful, but their results have been limited because the analogy with steel ingots obviously is not a close one. Model wax ingots have also been useful in pointing out important features of mold design. Gathmann4 watched the solidification of wax ingots, through a glass-walled mold, and has reported the effects of mold type on the location of the shrinkage cavities. Segregation was simulated by using a mixture of paraffin and stearine and adding a red dye that on freezing segregated preferentially to the stearine. Brearley5 also used wax models to determine the effects of some of the principal features of mold design. Carlsson and Hultgren6 dsed stearic acid to show that the columnar crystals of ingots tend to slope upward progressively from the mold wall toward the center of the ingot. As with the nonferrous ingots, the wax ingots have

Jan 1, 1943

By Margaret Dienes, Karl L. Fetters

It is recognized that ingot structure is important in determining the quality of finished steel. Such elements of ingot structure as the size, shape and distribution of primary crystals; the size, shape and distribution of shrinkage cavities; the degree and kind of segregation (both macro and micro); the distribution of nonmetallic impurities and gas cavities are each the result of solidification phenomena. 'These structural features vary with variations in ingot practice. Their control requires a knowledge of their correlation with the variables of ingot practice. The structures of full-sized commercial steel ingots have been studied extensively by the Steel Ingots Committee of the Iron and Steel Institute.l This work has contributed much to our knowledge of ingot structure and the factors that control it. These investigations of necessity have been limited in scope because of their costliness and the great amount of time required. Studies of the structure of small steel ingots have not been very satisfactory because of the rapid rate of freezing of the ingots. Their structures are thus difficult to resolve and to interpret. Medium sizes of steel ingots in the ranges of 150 to 500 lb. do offer opportunity for study, but their cost also becomes greater and their use less expeditious. Experimental models for the study of steel-ingot structure have been made of nonferrous metals. Gathmann2 used lead-tin alloys and reported that they were satisfactory for determination of the location of the major shrinkage cavities and the Liplanes of weaknessv resulting from the pattern of crystallization. Certain features of mold design have been illustrated by these results. Bezdenezhnihk3 has used copper and aluminum ingots similarly and has shown the effect of mold-wall thickness on the depth of columnar grains. These studies of nonferrous models have been useful, but their results have been limited because the analogy with steel ingots obviously is not a close one. Model wax ingots have also been useful in pointing out important features of mold design. Gathmann4 watched the solidification of wax ingots, through a glass-walled mold, and has reported the effects of mold type on the location of the shrinkage cavities. Segregation was simulated by using a mixture of paraffin and stearine and adding a red dye that on freezing segregated preferentially to the stearine. Brearley5 also used wax models to determine the effects of some of the principal features of mold design. Carlsson and Hultgren6 dsed stearic acid to show that the columnar crystals of ingots tend to slope upward progressively from the mold wall toward the center of the ingot. As with the nonferrous ingots, the wax ingots have

Jan 1, 1943

By R. Schuhmann, G. W. Powell

The concept of local equilibrium is examined and, in particular, the applicability of the concept to two-phase binary diffusion couples is discussed. It is concluded that, if binary solid solutions are treated thermodynamically as ternary solutions, vacancies being the third component, then recent experimental observations of nonequilibrium interface compositions (i.e., compositions different from the phase diagram values) are consistent with a state of complete dynamic equilibrium at the interface. CLASSICAL approaches to diffusion and to other irreversible processes often utilize a postulate of "local equilibrium". That is, one assumes even in the presence of gradients, fluxes, and time changes of properties, a) that the Second Law can be formulated locally, and b) that thermodynamic equations of state among intensive properties are identical with those for equilibrium phases. A direct consequence of this assumption, for example, is that one would calculate the total entropy of a nonequilibrium phase simply by integrating the entropy density (eu per unit volume) over the volume of the phase, assuming that entropy density is the same function of energy and particle density as for equilibrium phases. In the study of diffusion and of diffusion-controlled phase transformations, metallurgists take this concept of local equilibrium for granted and utilize it automatically. In fact, experimental evidence of all kinds is overwhelming that "local equilibrium" is a valid limiting law, approached more and more closely as gradients and rates of change approach zero. On the other hand, we have neither experimental nor statistical thermodynamic proof that this principle of "local equilibrium" should remain rigorously valid at high gradients and for rapid processes. Theoretically, in fact, one would expect deviations because interactions between adjacent volume increments of material affect the evaluation of entropy. Also one observes experimentally structural and small energetic differences between equilibrium and nonequilibrium phases. The purpose of this paper is to present a thermodynamic framework for the study and correlation of deviations from strict local equilibrium in isothermal diffusion and related processes. This procedure, in effect, preserves the principle of local equilibrium by thermodynamically assigning the deviations from static equilibrium properties to a defect component. Although this theory is entirely phenomenological, the possibility does exist of correlating the theory with experimental observations of defects under favorable circumstances. The most immediate application appears to be to the equilibria at the interfaces in two-phase binary diffusion couples. A recent pertinent work by Eifert et al.1 has shown that the concentrations of aluminum in the a(fcc) and B(bcc) phases at the 0-8 interface in Cu:Cu-12.5 wt pct A1 diffusion couples are greater than those given by the phase diagram. It is postulated by these authors that the supersaturation of the o phase at the interface provides the driving force for the a-to-B lattice transformation; the excess free energy of the a phase, ?Gaxs, arising from the greater-than-equilibrium aluminum concentration is about 3 cal per g-atom. On the basis of the schematic free energy-composition diagram shown in Fig. 1, Eifert et al. concluded that the thermodynamic situation at the interface is one of the partial equilibrium; that is, = for the more mobile aluminum whereas will be reexamined in the latter part of this paper. THERMODYNAMIC PROPERTIES OF DEFECTS A rigorous analysis of the thermodynamics of defects requires an added degree of freedom in the Phase Rule sense, because the standard formulations of solution thermodynamics (partial molal properties of components, free energy minimization, equality of chemical potentials, and so forth) involve the implicit

Jan 1, 1970

By M. G. Corson

INFORMATION regarding the use of manganese alloys has hitherto been incomplete and available only from widely scattered sources. This paper attempts a systematic description of properties .and uses of alloys other than manganese bronze, manganin and duralumin, which are probably the only non-ferrous alloys associated with manganese in the minds of most metallurgists.1 PROPERTIES OF MANGANESE Manganese can react either acid or basic; it is hard to obtain in the pure state; it has two different crystalline lattice structures, both of the cubic system, one stable above 850° C. and the other below 650° C., both coexisting in the intermediate range. Although the commercial product is hard and brittle, the pure metal seems to be ductile. In fact Zhemchuzhny, in Russia, was able to obtain ductile manganese and draw it to fine wires by the simple expedient of adding about 3 per cent. copper to the alumino-thermic product. Copper most probably brings the silicides, aluminides, nitrides and carbides present in the commercial product into a much less harmful form of distribution. One industrial use that can be foreseen for wrought manganese articles is as anodes for electrolytic processes in sulfate and nitrate solutions, particularly with high voltages. A practically insoluble skin of manganese peroxide is likely to form in such cases and an insoluble and strong anode material will be available.

Jan 1, 1927

By W. L. Mamme, G. Y. Chin

The problem of selection of the active slip systems for a crystal undergoing an arbitrary strain was analyzed by Taylor and by Bishop and Hill in terms of a minimum (internal) and a maximum (external) work criterion, respectively. These two criteria have now been generalized to include crystallographic slip on several sets of slip systems, twinning mixed with slip, and slip by (noncrystallographic) pencil glide. The generalized treatment also takes into account the possibility of a Bauschinger effect and of unequal hardening among the shear systems, which were considered in the Bishop and Hill work. Optimization techniques of linear and nonlinear programming are shown to be applicable for the numerical calculation of the minimum or maximum work. In the case of crystallographic shear, the constraint functions are linear and hence the optimal work is obtained as the saddle value of the lagrangian function Wi(y) e minimum and W,(u) + (a) for the maximum, where Wi is the (internal) work, We is the (external) work, Y is the crystallographic shear strain, u is the applied stress, and and are constraints. It is shown that the Lagrangians are functionally the same and the saddle value of one problem is identical to the saddle value of the other, proving that the two analyses are completely equivalent. In the case of pencil glide, although the constraint functions are nonlinear and neither convex nor concave, the equivalence of the optimal values to the saddle value of the Lagrangian (which is again identical for both problems) is still valid. WHEN a crystal deforms plastically by crystallographic shear, five independent shears are generally required to accommodate five independent strain components specifying the deformation. Assuming slip as the only shear mechanism, Taylor1 in 1938 analyzed the deformation in terms of a minimum work criterion. He hypothesized that of all combinations of five slip systems which are capable of accommodating the deformation, the active combination is that one for which the internal work C is a minimum, where 1 TI is the critical resolved shear stress for slip on the 1-th slip system and is the corresponding simple shear. By further assuming equal 72 for all equivalent slip systems and no Bauschinger effect, Taylor re- duced the minimum work problem to one of minimum and applied the analysis to the case of axisym- Metric flow by {111}(110) slip in fcc crystals. However, he did not consider the question of whether the resolved shear stress has in fact attained the critical value for slip on the newly found active systems without exceeding it on the inactive systems. In 1951 Bishop and ill' put forth the maximum work analysis in which slip is again assumed as the only deformation mechanism. In this analysis, the work o1 done in a given strain ij by a stress ujj not violating the yield condition is maximized. In addition, the analysis takes into account the possibility that the critical resolved shear stress for slip may not be equal among the slip systems and that the slip behavior may exhibit the Bauschinger effect. As with Taylor, a single set of slip systems—{111)(110) — was analyzed numerically. It thus appears that the Bishop and Hill treatment is on a more sound physical basis than the Taylor treatment. However, Bishop and Hill showed that where there is equal hardening among all slip systems and when there is no Bauschinger effect, Eq. [11 ] of Ref. 2, as assumed by Taylor, the results of their maximum work analysis are the same as those of Taylor's minimum work analysis. Hence at least under those conditions there is an implication that the Taylor analysis does lead to a critical resolved shear stress for slip on the predicted active systems without violating the yield condition on the inactive systems. Recently, the Taylor analysis was applied for numerical solutions of the axisymmetric flow problem, for slip on {110}(111), {112}(111). {123)(111) systems as well as a mixture of all three sets of svstems."1 Computational techniques based on the optimization theories of linear and nonlinear programming4 were employed in these solutions. The same techniques were employed in the solutions of an axisymmetric flow problem of deformation by slip on (111) (110) systems and twinning on (111)(112) systems5 which had been considered theoretically from a modified Taylor approach. The utilization of these techniques has led to the realization that the solutions of Taylor's minimum work problem imply the solutions of Bishop and Hill's maximum work problem. The two problems turn out to be dual problems in the well known sense of mathematical programming. It is thus the purpose of this paper to first generalize the minimum and maximum work analyses to include crystallographic slip on several sets of slip systems, twinning mixed with slip, and slip by (non-crystallographic) pencil glide, as well as the possibility of a Bauschinger effect and of unequal hardening

Jan 1, 1970

By F. N. Kissell, W. E. Schroeder, N. I. Jayaraman

This article discusses methods of improving airborne capture by conventional dust suppression systems. The work was performed in a full-scale plywood model of a mine entry with a wooden mockup of a Joy 12 CM continuous miner. Dust was introduced at the front of the entry, and water sprays knocked down part of the dust cloud. Sampling in the return air determined the percentage dust reduction. Based on the results, guidelines are offered for continuous miner spray systems to control dust. Empirical relations have been developed to determine the effect of changes in water spray pressure on air movement, dust knockdown, and pollutant concentrations in the face return. Nozzle selection and arrangement are discussed for high dust knockdown efficiencies.

Jan 1, 1986

By R. L. Moment

The anisotropic elastic behavior of rolled beryllium sheet has been measured, using a pulse echo technique, and compared with X-ray diffraction data. Calculated elastic stiffness constants compared favorably with published values for beryllium single crystals which were attributed to the strong (0002) rolling plane texture. Variations of Young's modulus in the yolling plane could be associated with the velative distribution of (0002) planes out of their ideal position in the rollitzg pkule. WHEN a metal is subjected to cold working such as drawing, forming, or rolling, a crystallographic texture develops which can significantly alter its physical properties. One method for detecting this texture is X-ray diffraction, but Alers and Liu' have recently pointed out how the prediction of anisotropic physical properties from pole figures alone is not always accurate due to differences in interpretation. Variations in Young's modulus with orientation or, more completely, the values of the effective elastic constants of the worked metal, also serve to indicate the presence of a texture. In fact, as Alers and Liu' pointed out, calculated variations in Young's modulus for assumed orientations, when compared with experimental data, can be used to eliminate some of the uncertainty in interpretation of X-ray pole figures. Thus, elasticity measurements can serve not only to clarify any unusual elastic behavior of worked metal, but also to detect and in part determine the nature of its texture. X-ray determination of the texture of rolled beryllium has been reported by Smigelskas and Barrett,2 who found a strong texture of (0002) in the rolling plane with (1070) planes normal to the rolling direction. In the case of metal rolled at room temperature, they reported that [1010] directions also appeared at positions 60 and 120 deg from the rolling direction in the rolling plane, while in more recent work Keeler3 found these directions were also tilted towards the rolling plane. The texture for beryllium rolled at 80O0C, however, only showed (1010) planes normal to the rolling direction and the spread of (0002) planes out of the rolling plane was less. In looking for elastic anisotropy one might consider unidirectional rolling of a metal as introducing an or-thorhombic symmetry through reorientation of the grains, since the three deformations, compression, extension in the rolling direction, and extension in the cross direction, are orthogonal to each other and unequal in magnitude. Thus the rolled sheet could be treated like an orthorhombic single crystal and the nine stiffness constants of the elasticity tensor used to calculate the anisotropy of Young's modulus, the shear modulus and Poisson's ratio. In this case we could write: which is symmetric about its diagonal. Borik and Alers4 have recently used this approach on rolled die steel with very good results. They found, however, that instead of displaying orthorhombic elastic symmetry their specimens could be considered tetragonal in which case Cr1 = c22, c13 = Ca, and c44 =cjj. This conclusion was made solely on the basis of the measured tensor elements, and serves to point out the advantage of this method for studying the anisotropy of rolled metals. Their calculated values for Young's modulus as a function of angle in the rolling plane also checked very well with direct measurements made on different specimens using the resonance technique. In the present study, cross-rolled beryllium was used which had been unidirectionally rolled about 11 pct for the final reduction. This imparted a slight anisotropy in the rolling plane which was detected both by X-ray techniques and elasticity measurements. For purposes of discussion in this paper, the rolling direction is that direction in which the most reduction passes were made and cross direction is the normal to the rolling direction in the rolling plane. It was also decided to consider the rolled sheet as displaying orthorhombic symmetry for the purpose of obtaining elasticity samples with the direction defined as in Table I. Any change in the final symmetry attributed to the sheet would then be made on the basis of the measured elastic stiffnesses. The final data would then be compared with that expected from the X-ray study and that reported for beryllium single crystals. EXPERIMENTAL PROCEDURE Rolling Schedule. The samples used in this study were taken from a large sheet which, because of its size, had to be unidirectionally rolled for the final reduction. The resulting texture was that of cross-rolled metal with a slight unidirectional texture superimposed. A cast beryllium ingot, 9.500 in. sq by 3.325 in. thick, was cross-rolled to 81 pct reduction followed by unidirectional rolling for an additional 11 pct to give a total reduction of 92 pct. The thickness of the final sheet ranged from 0.265 to 0.280 in. Reduction up to 67 pct was done at 980°C and the final 25 pct at 870°C. Analysis for metallic impurities showed aluminum 0.06 pct, iron 0.19 pct, and silicon 0.11 pct, giving a beryllium purity of 99.64 pct.

Jan 1, 1968

By Charles R. Stewart, William W. Owens

Reservoir performance predictions based on laboratory core test data assume that fluid flow is laminar for the laboratory test. A study has been made to determine the validity of this assumption for laboratory tests on various types of porosity found in producing limestone formation. Data are presented which show that turbulence and slippage can occur during laboratory tests on hetero-geneous-type porosity limestones, thus causing serious errors in measured single-phase permeabilities and two-phase relative permeability characteristics. In single-phase flow tests it is possible to eliminate turbulence and correct for slippage or to eliminate both factors by controlling test conditions. It is not always possible to control test conditions and thereby eliminate turbulence and slippage in two-phase .flow tests. A correction method is presented which can be used to calculate the true two-phase laminar flow relative permeability characteristic even though furbulence and slippage exist. .INTRODUCTION It is customary to make use of Darcy's law and modifications of this law, together with laboratory data on formation core samples to predict the performance of producing reservoirs. Such predictions are based on an assumption that fluid flow is in the laminar or streamline region for the laboratory test. It was the purpose of this inves- tigation to determine the extent to which turbulent flow may occur in laboratory fluid flow tests on hetero-geneous porosity limestones. Considering that turbulent flow conditions might exist in some laboratory fluid flow tests, additional emphasis was placed on the development of a method to correct for turbulence when laminar flow conditions could not be attained. FLUID FLOW CONCEPTS FOR POROUS MEDIA The Influence of Pore Geometry on Fluid Flow One of the more important factors influencing fluid flow in porous media is the geometry of the pore space which includes such characteristics of the pores as size, shape, distribution, roughness, uniformity, etc. In general, oil- and gas-producing formations can be divided into two broad types on the basis of pore geometry. One has been called sandstone-type porosity media, which is characterized by a small range in pore size, uniformity in shape of the pores, smooth pore surfaces and a regular and uniform distribution of pores. The other type has been called heterogeneous porosity media and is usually limited to the dolomites and limestones. This type is characterized by a wide variation in the size, shape, and distribution of the pores and rough, irregular pore surfaces. It is therefore apparent that conditions are much more favorable for turbulent flow* in heterogeneous-type porosity media than in sandstone-type porosity media. Interrelationship Between Turbulence and Gas Slippage In studying the problem of turbulent flow in laboratory tests on porous media, it is necessary to be aware of the interrelationship between slippage and turbulence for gas flow. As a result of slippage or the Klinken-berg effecta, apparent perrneabilities to gas are greater than the true value because there is no stationary layer of gas in contact with the walls of the flow channels. Gas slippage decreases as the mean free path of the gas molecules decreases. Since the mean free path of any gas decreases with increasing density, increases in static pressure result in lower apparent gas permeabilities. However, a reduction in gas permeability can also be due to turbulence. Therefore, in studying only turbulent flow in porous media, it is necessary to hold gas density, and slippage, constant or to reduce slippage to a negligible value by operating at high static pressures. Presentation of Laminar and Turbulent Flow Data A graphical relationship between permeability and a pseudo-Reynolds number, N,, will be used to show the two types of fluid flow, i.e., laminar and turbulent. The usual graphical method for such a description has been the use of friction factor-Reynolds number charts4. On such a logarithmic diagram, the laminar region appears as a straight line having a slope of 45 degrees. As the friction factor decreases and the Reynolds number increases, the turbulent region is reached and appears as a deviation from the 45-degree slope line. However, in petroleum engineering literature resistance of por-

By H. J. Ramey

A general solution is presented for the transient temperature distribution caused by radial movement of a cylindrical heat source through a homogeneous medium of infinite extent. This problem represents a highly simplified model of the movement of a combustion front during the thermal recovery of oil. Numerical solutions are presented for a heat source moving at a constant velocity, and a velocity inversely proportional to the radial location of the heat source. Numerical solutions are presented for both finite and infinite vertical thickness of the heat source, i.e., with and without vertical heat losses. The numerical solutions are used to estimate the fuel concentration needed to maintain a combustion front during the thermal recovery of oil. The sensible heat carried to the combustion front by the injected gas stream is discussed in regard to the fuel concentration. Numerical solutions for a heat source (combustion front) of finite vertical thickness indicate the field conditions that may be necessary to sustain the temperature of a combustion front above the ignition temperature of the fuel over considerable distances from an injection well. The results also have implications in regard to the quantity of heat required during ignition. The results of this computation apply to a highly idealized model of a thermal recovery process. But the results may be used as a guide in engineering consideration of the thermal recovery process. INTRODUCTIO N The thermal recovery of crude oil has received considerable attention since the publication by Kuhn and Kochl of laboratory and field tests of this method of oil recovery. Other publications2,3,4,5,6, since that time have dealt with various features of the thermal recovery of oil and have illustrated the extremely complex nature of this process. Vogel and Krueger7 described an electric analog computer designed to solve a moving heat source problem which was analogous to movement of a combustion front during the thermal recovery of oil. The heat source was maintained at constant temperature and was a cylindrical source of infinite height (no vertical heat loss) moving radially at either constant velocity, or at a velocity inversely proportional to the radial location of the heat source. Jenkins and Ramey2 presented analytical solutions to a similar heat conduction problem and pointed out the possible importance of vertical heat losses. The following presents a general solution to the transient heat conduction problem introduced by Jenkins and Ramey4 as well as numerical solutions for a wide range of conditions possibly similar to those that might exist during field operation of the analogous thermal recovery process. DESCRIPTION OF PROBLEM The problem considered in this paper is determination of transient temperature distribution caused by a cylindrical heat source of infinite or finite vertical height, moving radially through an isotropic medium of infinite extent. It is further assumed that the heat flux generated at the surface of the moving heat source may be a function of time such that the fuel concentration required may be constant, or permitted to vary, and that heat is generated at the surface of the heat source only. In regard to the thermal recovery of oil the assumption that heat is generated at the surface of the front is equivalent to assuming that the thickness of the combustion zone is infinitely small, or that the reaction rate between fuel and oxidant gas is infinite. The last assumption does not appear stringent. See laboratory information published by Martin, et al3 Under these assumptions the temperature distribution throughout an infinite medium caused by a heat source moving radially with an arbitrary velocity may be described by The heat generation function or heat flux at the surface of the source, +(t), is Note that the fuel concentration may be defined as a function of time. Eq. 1 is derived in the Appendix.

By E. A. Gulbransen, K. F. Andrew

Weight loss measurements were made using a sensitive microbalance operating in a high vacuum system. The Langmuir equation was used to Calculate the vapor pressures of the several metals and alloys. Aluminum lowered the vapor pressure of i9,on in a 5.4 Al-94.6 Fe alloy, and of iron and chromium in a 4.8 Al-21.5 Cr-73.7 Fe alloy. The influence of oxide films formed on the alloys on the effective vapor pressures of the alloys was also studied. These studies were used to aid in the inte9,pretation of the high temperature oxidation of alloys based on iron, chromium, and aluminum. RECENT studies on the oxidation of chromium1 and the heat resistant alloy 5 Al-22 Cr-73 Fe showed transitions in the rate of oxidation at 900o and 1050oC respectively. The transition for chromium occurred at a temperature where the rate of evaporation of chromium from oxide-free surfaces equalled the rate of chromium atoms reacting with oxygen. This paper presents new vapor pressure studies on iron, chromium, and the specially prepared alloys 5.4 Al-94.6 Fe, 21.9 Cr-78.1 Fe, and 4.8 Al-21.5 Cr-73.7 Fe. The ternary alloy is similar in basic composition to the patented heat-resistant commercial alloy known as Kanthal. The binary alloys were studied to show the individual effects of aluminum on the vapor pressure of iron and of iron on the vapor pressure of chromium. The purpose of the proposed studies was to test the influence of metal volatility on the oxidation behavior of heat-resistant alloys based on iron, chromium, and aluminum. Since oxide films formed on the metal may limit metal transfer to the surface, the influence of oxide films on metal volatility was also studied. LITERATURE A) Vapor Pressure. The vapor pressure of iron has been studied by Jones, Langmuir, and Mackay,' Dornte and Nrton. Edwards. Johnston. and Ditmar. and McCabe, Hudson, and Paton. Stull and SinkeT have calculated a heat of sublimation at 298° K of 99.83 kcalper g atom. The vapor pressure of chromium has been reported by Bauer and Brunner,' Speiser, Johnston, and Blakburn, Gulbransen and Andrew,'' Vintaikin,o McCabe, Hudson, and Paton, and Kubaschewski and Heymer.12 Good agreement was obtained except for the early work of Bauer and Brunner.' Stull and Sinke7 calculate a heat of sublimation at 298°K of 95.0 kcal per g atom. Vapor pressure measurements on aluminum have been made by Brewer and Searcy,13 Bauer and Brunner,aand Farkas.I4 Stull and Sinke,7 giving Brewer and Seary's' data the most weight, derive a heat of sublimation of 77.5 kcal per g atom at 298 oK. B) Method. Two methods15 are used for measuring the vapor pressure of metals: 1) The Langhuir free evaporation method, and 2) the Knudsen effusion method, In the Langmuir method, used in the present work, the vapor pressure of the metal P is given by the equation: Here M is the molecular weight of the vapor species. T is the absolute temperature. (dw/dt) is the rate of sublimation in g per sq cm per sec and a is the condensation coefficient. a is a measure of the efficiency of condensation of molecules striking the metal surface from the vapor. If condensation results from each collision a is unity. This coefficient has been found to be unity,a, le for metals. These results establish the general validity of the Langmuir free evaporation method as applied to metals. EXPERIMENTAL A) Vacuum Microbalance Method. A modification of the Langmuir free evaporation method was used. Strip specimens were suspended from a sensitive microbalance operating inside of a high vacuum system. For alloys, microbalance methods17 must be used on large area samples at relatively low temperatures to minimize composition changes. The

Jan 1, 1962

By C. H. Wilson, J. W. Daly, J. J. Jakosky

Meteor Crater, Arizona, is a natural wonder which for years has been the subject of considerable discussion and study as to its origin and age (Fig. 1). Of the various theories advanced regarding the origin, two have been considered as most probable by a majority of investigators. The first theory is that the crater was formed by a meteorite or a swarm of meteoric material striking the earth at a high velocity and burying itself. The second theory is that the crater was formed by a "steam explosion" from hot solutions or gases coming from underneath the present sedimentary beds which overlie the area. In order to obtain more detailed information regarding the conditions existing at the crater, the present management decided to have a detailed geophysical examination made. The origin of the crater is of considerable scientific importance and interest. More or less active mining operations have been carried out during the past 20 years by various groups who believed that the crater was of meteoric origin, and who hoped that sufficient quantities of such material had been buried to be commercially valuable. Fragments of meteoric material found on the surface of the ground in the debris around the rim gave analysis of approximately 92 per cent iron and 8 per cent nickel, with some platinum and iridium. The material having a gross value of approximately $50 per ton, mining would be commercially profitable if sufficient quantities of the material existed in condition for mining. In the early work numerous drill holes and shallow shafts were put down, mainly in the northern portion of the crater. Pulverized material and debris caused considerable difficulty, and in an effort to overcome this the present management started a 1500-ft., two-compartment shaft on the south rim of the crater and outside of the fill zone. The intention was to sink the shaft in the solid rock outside of the crater, and then to crosscut toward the crater at a depth of 1500 ft. The work was finally stopped when a depth of 650 ft. was reached, because water poured into the shaft through large cracks and fissures in the rock. The enormous force responsible for the crater had also broken

Jan 1, 1932

By F. Bonaccorso, G. Venturello, C. Antonione

A study of the effect of small amounts of interstitial impurities on recovery and re crystallization in high-purity iron (99.995 pct) has been undertaken. This paper gives results on the effect of carbon, introduced in small-dosed amounts @om 0.0005 to 0.0086 wt pct) by heating iron specimens of high purity m a static atmosphere of CO + ,. The materials prepared in this way, cold-worked 80 pct and subjected to a series of isochronal and isothermal annealings, were submitted to examinations by X-rays, micrographs, and hardness tests. It was observed that effect of carbon is remarkable in the sense of blocking recovery of mechanical properties up to the temperature at which re crystallization begins. On the contrary, carbon has a negligible effect on primary re crystallization temperature, when compared with the known effect of substitutional impurities. This is in agreement with the high mobility of interstitials. In effect, only a slight decrease, -20 pct, of the grain-boundary motion rate was noted, due to the interaction between the grain boundary and the carbon atoms. On the other hand, in the samples in which the carbon content is above the solubility limit at the temperature at which re crystallization occurs, a slight increase of nucleation frequency is noted due to the presence of precipitated carbides. ThE effect of purity on recovery and recrystalli-zation phenomena has been known for a long time; however, recently, new attention has been given to the problem since new methods for obtaining metals with extremely low impurity contents have become available. Most of this research work essentially concerns the effects of impurities which cause precipitates or give rise to substitutional solid solutions. The works of Bolling and winegard1 and Aust and utter' on lead, and Vandermeer and orddon' on aluminum, 01sen4 on nickel, and Abrahamson and Blakeney on iron should in particular be referred to. On the effect of this type of impurities a quantitative theory has been formulated by Detert and Lucke6 and has lately been discussed critically by Cahn7 and Gordon and vandermeer.' Interstitial solid solutions in iron have not as yet been studied. As the study of the effect of interstitials is of great interest both from a theoretical and practical standpoint, it was deemed useful to examine the effects of carbon, nitrogen, and boron on recovery and primary recrystallization of iron. There already is some work by Chaudron et al.,1-" on the effect of interstitials on iron; their work, however, mainly concerns secondary recrystallization. The present paper refers in particular to the effect of carbon. EXPERIMENTAL PART Preparation of Materials. The pure iron used for this research was obtained from FeCIS recrystal-lized and purified by extraction with isopropylic ether. From the ferric chloride purified in this manner, hydroxide was precipitated with a solution of very pure ammonia, and, by calcination in pure sintered alumina crucibles, oxide was obtained. Reduction of the oxide to iron sponge was performed in sintered alumina tubes with very pure hydrogen at 650°C; at the end of the operation temperature was increased to 900°C. Specimens for the experiments were obtained by sintering the sponge at 1480°C in pure Hz after a pre sinter ing treatment at 900° C. It is important to note that the treatment at 1480°C in Ha produces a further purification from more volatile elements such as zinc, cadmium, arsenic, lead, and tin. Details on the preparation and characteristics of this type of very pure iron are given in a previous work.'' Only the complete analysis performed by neutron-activation methodsz3 is given here, Table I. Some of the specimens prepared in this way were carburized in the 0 region with very low amounts of carbon by treating them at 700°C in a static atmosphere of Ha containing a definite amount of CO. The set-up used is described in Fig. 1. A gas-tight quartz tube containing the specimen to be carburized and an internal friction control specimen, after being evacuated, was filled with Hz

Jan 1, 1963

By Alan J. Strauss

The Group IV elements germanium, tin, and lead form nine 1:1 compounds with the Group VI elements sulfur, selenium, and tellurium. This paper reviews the properties of the pseudobinary solid solutions formed by these compounds, including the extent of mutual solid solubility, temperature-composition phase diagrams, transport properties, deviations from stoichiometry, optical properties, and energy band structure. Particular emphasis is placed on the Pbl-xSnxTe and Pb1-xSn,Se alloys with rocksalt structure, because of current interest in these malerials for generating and detecting infrared ,radiation. ThE Group IV elements germanium, tin, and lead form nine 1:1 compounds with the Group VI elements sulfur, selenium, and tellurium. Some of the physical and electronic properties of these compounds, including their melting points1-8 and energy gaps,4'9-12 are listed in Table I. Four compounds (SnTe, PbS, PbSe, and PbTe) have the cubic rocksalt (Bl) structure. At room temperature GeTe has a rhombohedra1 structure closely related to the B1 structure, into which it is transformed at about 400°C. Four compounds (GeS, GeSe, SnS, and SnSe) have the orthorhombic B29 structure. In samples which have not been intentionally doped with impurities, the electrical conductivity is due primarily to electrons or holes produced by the ionization of donor or acceptor lattice defects associated with deviations from stoichiometry. Undoped samples of PbS, PbSe, and PbTe may be either n type or p type, depending on whether they contain excess lead or an excess of the Group VI element, respectively, but only p-type samples of the other compounds have been reported. This paper will review the properties of the pseudo-binary solid solutions formed by the nine 1:l compounds. The topics to be considered include the extent of mutual solid solubility, temperature-composition phase diagrams, transport properties, deviations from stoichiometry, optical properties, and band structure. Particular emphasis will be placed on the Pb]-xSnxTe and Pbl-xSnxSe alloys with B1 structure, which are promising materials for generating and detecting infrared radiation in the 8 to 14 µm atmospheric window and beyond. MUTUAL SOLID SOLUBILITY The extent of mutual solid solubility in the pseudo-binary systems has been investigated for fourteen of the eighteen ternary systems (in which the two terminal compounds have a common element) and for nine of the eighteen quaternary systems. In most cases, X-ray diffraction measurements made at room temperature were used to determine the structure and lattice parameter(s) of the phase(s) present in samples of various compositions prepared by freezing from the melt. In many investigations, including the extensive studies of Krebs and co-workers,13'14 the samples were annealed at elevated temperatures before X-ray measurements were made. In some cases, metallographic examination and thermal analysis have also been employed. The results for the ternary and quaternary systems are summarized in Tables II and III, respectively. (The original references should be consulted for the annealing temperatures.) The solubility of one compound in the other is at least 5 mol pct in all cases, and is often much larger. Complete solid solubility has been observed in all systems so far investigated in which the terminal compounds have the same structure, although in the PbS-PbTe system complete solubility is limited to elevated temperatures. Thus complete solid solubility occurs in five systems where both compounds have the B1 structure (SnTe-PbTe, PbS-PbSe, PbS-PbTe at temperatures above 805°C, PbSe-PbTe, and SnTe-PbSe) and in three where both compounds have the B29 structure (GeS-SnS. GeSe-~n~e-, and SnS-SnSe), as well as in two involving GeTe and a compound with B1 structure (GeTe-SnTe and GeTe-PbTe). On the basis of thermal analysis data, complete solid solubility at sufficiently high temperatures has also been reported for the GeSe-GeTe2 and SnS- pbS22 systems. This seems unlikely, however,

Jan 1, 1969

By H. M. Woodrooffe

Three of the world's largest phosphate deposits are located in the USSR. These have an estimated reserve of 2,600 million short tons of elemental phosphorus. The best known lies in the Khibiny Massif, an alkalaic rock formation some 600 square miles in area with mountain peaks rising to an altitude of 3300 ft above the surrounding terrain. The Khibiny Massif is north of the Arctic Circle in the Kola Peninsula at a latitude of 67° 30' north and a longitude of 35° 30' east. It is connected by rail to Murmansk, 100 miles northwest, and Leningrad, 900 rail miles south.

Jan 12, 1972

By R. B. Sudderth

The recovery of copper by various combinations of processes centered around liquid ion exchange has certainly been accepted as a technically superior entity when compared to the widely applied classical method of cementation. While the technical aspect of liquid ion exchange processes has received much attention and acclaim through actual operating data, in particular the ease of handling solutions, percentage recovery, reduction of pollution, etc., the corresponding economic data has either gone unpublished or has been held confidential. Much of this economic data has now been made available for publication; and with this data, we at General Mills Chemicals, Inc., will attempt to present an accurate economic analysis of some liquid ion exchange plants and an economic comparison between these plants and the typical cementation facility.

Jan 1, 1974

By Major N. R. Junner

Although manganese-garnets are rarely seen in hand specimens of the manganese ores and associated rocks, they have been detected by the writer in thin sections of all the ores and wall rocks that have been sliced, with the single exception of the Insuta deposit. Rhodonite, and an almost colorless amphibole, which appears to be manganiferous, are less common. Types of Occurrence (a) Rocks Composed of Manganese-garnet and Quartz The best known examples of this type are from near Bole and Wa in the Northern Territories of the Gold Coast, and from near Chama in the Gold Coast Colony. Normally these are fine and even-grained rocks, in which the garnets cannot be recognized in hand specimens. They are often stained with secondary manganese oxide. Thin sections examined microscopically are seen to be composed of aggregates of minute crystals of garnet in a matrix of mosaic quartz. Some manganese oxide and braumite (?) is often present, but other minerals are exceptional. The relative proportions of garnet and quartz in these rocks vary greatly and there is a complete gradation from quartz-rock to garnet-rock. (6) Rocks Composed Mainly of Manganese-garnet and Manganese Oxides Rocks of this type are known from Kalimbi Hill near Siripe, Pabini Hill near Bangwon, near Odumase and Konongo, and in the Ya Ya creek, about 5 miles east of Dunkwa. These are usually fine-grained banded dark-colored rocks, composed of abundant minute garnets in a matrix of black manganese oxide and braunite (?). A little microcrystalline is often present, and occasionally a little colorless amphibole. quartz (c) Rocks Composed of Rhodonite and Manganese-garnet Rocks of this type are not common. They are known from near Chama, and from the road between Senyong and Mangkuma. Some quartz is usually present. All these rocks resemble closely in texture and composition certain rocks from India, described by Fermorl as "gondite;" and they also resemble certain rocks described by Derby2 from the Queluz district,

Jan 1, 1927

By R. H. Richman, F. Borik

A 0.3 pct C-12 pct Cr-6 pct Ni steel was rolled to 93 pct reduclion in area as austenite at 510°C, and then partially transformed as desired to ~rlartensite by qnenching to - 196°C. Pole figures for the austenitic matrix and for the martensitic product were separately determined by an X-ray transmission method. The deforitration texture of' the warm-worked austenite is characlerized by (110)(225) components, and is thus closely similar to those produced in a brasses. The pole jigure of the martensite in partially transformed material agrees well with that which can be constructed by transfortnation of the {110)(225) orientations according to either the Kuvdjuniov- Sacks or the Nishi-yatuu relatiotship. Howeuer, an important result of this construction is that me-third of the predicted orientations are missing. A graphical analysis can then be used to show that in deformed austenite certain crystallographic variants of martensite (related to the most probable austenite slip systems) are suppressed, resulting in this preferred transformation. The evidence for preferred transformation is corroborated by the measured elastic anisotropy of warm-rolled and fully transformed H-11 steel. EXTENSIVE plastic deformation of a polycrystal-line aggregate in a manner that causes flow predominantly in one direction results in a preferred orientation of the constituent crystallites. The particular orientations that are produced depend upon the crystal structure and composition of the material, as well as upon the temperature, mode, and degree of deformation; in any case, the preferred crystallo-graphic orientations, or textures, are reflected in directionality of mechanical properties. Although such anisotropy may be exploited in certain specialized applications, it is more commonly diminished or eliminated by heat treatment lest it interfere undesirably in subsequent forming operations or in structural design. In the recently developed thermomechanical treatments that significantly enhance the strength of some steels,1,2 considerable deformation of the metastable austenite prior to the martensite transformation is essential to the strengthening process. If the austenite is textured by the deformation, and if the transformation to martensite proceeds according to one of the relationships established for transformation in annealed austenite, then it must be expected that the martensite will also possess a preferred orientation even though the multiplicity of martensite orientations possible in a given austen- ite crystal will tend to restore some degree of randomness. The existence of a residual anisotropy, both mechanical 3-6 and crystallographic,' has been substantiated. In the latter crystallographic investigation, preferred orientations were determined for the martensitic structure of an SAE 4340 steel rolled 72 pct as austenite at 833°C and then quenched. However, the choice of a composition that transformed almost completely to martensite during the quench to room temperature did not permit direct measurement of the prior austenitic texture. In fact, when the "ideal orientations'' associated with well-known fcc rolling textures were converted, alone or in combination, to martensite according to the Kur-djumov-Sachs (K-s)' or Nishiyama8 relations, the agreement obtained with the observed martensite texture was only fair at best. Recently a pertinent aspect of the austenite to martensite transformation was reported by Bokros and parker,10 who found that certain habit-plane variants of martensite were suppressed by tensile deformation of Fe-31.7 Ni single crystals prior to the necessary subzero cooling. It might be anticipated that the consequences of such preferred transformation are sustained during the formation of martensite in warm-worked austenite that has a well-developed deformation texture. The present investigation was undertaken first to establish more firmly the relation between preferred orientations in plastically deformed austenite and in the resulting martensite, and second to examine the textures for evidence of deformation-induced preferred transformation. EXPERIMENTAL PROCEDURES An alloy containing 0.3 pct C, 12 pct Cr, 6 pct Ni, and the balance iron, was selected because the mar-tensite-start temperature (M,) of about -100°C allowed convenient experimental manipulation of either austenite or martensite at room temperature. Furthermore, this composition can be readily deformed as metastable austenite at moderately elevated temperatures without intervention of appreciable isothermal or athermal decomposition products. The alloy was austenitized at 1150°C, aircooled to 510°C, rolled unidirectionally at this temperature to 93 pct reduction of cross-sectional area, and finally oil-quenched to room temperature. Partial transformation to martensite was accomplished by quenching to -196°C as needed. The rolled stock was reduced in thickness from 0.067 to 0.010 in. by etching in a solution of 5 pct HC1, 45 pct HNO3, and 50 pct water, and further thinned by careful mechanical polishing to maintain the two sides of the sheet parallel within 0.0003 in. After mechanical polishing to 0.005 in., electropolishing in 1:9 perchloric-acetic acid solution produced a final thickness of 0.002 in. The preferred orientations were determined from

Jan 1, 1968

By G. Abowitz, R. B. Gordon

The velocity of sound, mass density, and thermal expansion coefficient were measured in dilute lead, tin, zinc, bismuth, cadmium, indium, and potassium amalgams as a function of composition. Measurements in the Hg-TI system extended ovev the wider range of composition, 0 to 40 at. pct TI. The sonic velocity was measured as a function of temperature. The calculuted adiabatic compressibility in all cases was found to decrease much more rapidly than would be expected on the basis of a linear variation with composition. An interpretation of the results is given in terms of changes in the Debye characteris -tic temperature of mercury on alloying and a possible change in the bonding of mercury itself. If a liquid is considered as an isotropic elastic medium, the phase velocity of a longitudinal elastic wave traveling in the liquid is given by where K is the bulk modulus, G the shear modulus, and d the mass density. G can be neglected for liquids of ordinary viscosity, i.e., approximately 1CT2 poise, because the shear stress component is relaxed essentially at once. That is, for liquids of ordinary viscosity and at frequencies of the applied stress up to several hundred megacycles, the relaxation time of the shearing process is much shorter than the period of the applied stress. Eq. [I] then reduces to The adiabatic rather than the isothermal compressibility is measured at the megacycle frequencies used in ordinary ultrasonic velocity measurements. The sonic velocity and compressibility have been determined for many organic liquids and aqueous eletroltes,' and the results have been interpreted in terms of existing theories for these fluids. The velocity of sound and compressibility are important parameters in theories of the liquid state. For example, equations of state have been proposed relating the free volume of liquids to the acoustic velocity.3 Comparatively very little work has been done on the acoustic properties of liquid metals and alloys. leppa,' orddon, and several other investigators

Jan 1, 1963