Search Documents

By G. W. Thomas

Jan 1, 1967

By R. W. Vieth, S. A. Bradford

THE orientations of tetragonal tin crystals are commonly examined by either the transmission Laue method or the back-reflection Laue method. The predominant planes and zone axes of the pattern are plotted on a Wulff net and indexed by measuring the angles between planes and directions. The problem of identifying the Laue spots has been approached in several different ways. Nicholas 1 tabulated interplanar angles for ten planes but did not relate the angles to permutation of indices or sign. A table of angles between adjacent planes and of angles between adjacent zones was published by Chandrasekhar and veal.' These tables are very helpful, provided the film shows all the low-order zones that occur within about 20 deg of each other. Frounfelker and Hirthe's compilation3 of angles between the (001) plane and twenty-nine other planes is useful for the indexing, only if the (001) spot can be found on the film. For all but the simplest orientations, however, more information is needed to index the diffraction spots. Standard stereographic projections, such as those published by Bicelli, 4 provide part of the information required. By their nature, these projections show the positions of only a few zone axes. In the tetragonal system many low-order zone axes are not normal to the planes having the same indices. Consequently, planes and directions can have entirely different positions on a stereographic projection. For example, the [011] zone axis is 33 deg from the normal to the (011) plane. The crystallographic orientation usually has to be deduced from a combination of the predominant zone axes and a few low-order planes that appear on the film. In addition to a standard stereographic projection of the planes, a stereographic projection of the directions is needed. Such a projection can be constructed from the angles given in Table I. To provide more information about angles between planes and directions, we have calculated interplanar angles up to (332) and angles between planes and directions up to [332] on an IBM 1620 computer. The angles are differentiated according to permutations of sign and indices. Those doing research on tin single crystals are welcome to request a copy of these tables from the authors.

Jan 1, 1967

By J. T. McCartney

OPTICAL PROPERTIES of U.S. coals of differing rank have been investigated by the Bureau of Mines, for about 20 years. These investigations have been remade to further the development of an objective basis for petrographic analysis of coal and to add to knowledge of the structure and composition of coal and coal components by determining their optical constants in various regions of the electromagnetic spectrum. Various experimental techniques have been used in these studies, including determination of reflectance from bulk specimens and transmittance of thin and ultrathin sections, interferometry, electron microscopy, and election diffraction. Similar studies leave also been made of graphite, the end product of coal metamorphism. The bureau of Mines conducted the first U.S. investigation of reflectance properties of coal components, laying the foundation for petrographic analysis based on arbitrary reflectance classes and for automated reflectance analysis of coals. Optical constants of coals of differing rank and of graphite have been determined in the ultraviolet-visible-near infrared spectrum. Maxima in the optical constants of graphite were observed in the ultraviolet region and were interpreted to be related to transitions of electrons between the valence and conduction bands. Similar maxima for coals, especially high-rank coals, appeared to be characteristic of the presence of graphitelike layers. A special application of optical theory to absorbing, anisotropic materials like graphite and high-rank coals has been developed. Electron microscopic and electron diffraction studies of coals have revealed various ultrafine structures in coal components and have shown the evolution of an oriented, crystalline graphite structure in high-rank coals.

Jan 1, 1967

By J. Grosso, D. J. Nagel

In many metallurgical studies, it is necessary to reveal the composition and structural properties of materials. Usually this is done by etching, but some substances do not etch well or easily. Anodizing is known to provide an alternate treatment in many cases. In general, this technique is simple, often more sensitive to composition differences than etching, has excellent resolution, and provides variable contrast in those elements and alloys which anodize. Here we report the use of anodizing in diffusion studies in the TiCb system. Fig. 1 shows a cross section of the diffusion zone in a TiCb couple. The anodizing has clearly revealed both composition and crystallographic features. The maximum bulk penetration of the titanium into the columbium is marked by a sharp line because of the steep concentration gradient. The titanium grain size and orientation differences are also evident. Variation in the grain size reflects the gradually decreasing columbium concentration.

Jan 1, 1967

By Maurice Deul

A special short-term research project was initiated to develop plant design parameters for treating mine waters to yield an effluent acceptable for discharge into streams and a sludge amenable to effective and economic disposal. Neutralization was investigated to provide simple and direct treatment of ferruginous acid waters using lime; coarse limestone, and limestone followed by lime. Limestone treatment was accomplished using a small cement mixer as a simple reactor to provide abrasive agitation of limestone and mine water. Results from tests with nine mine waters encompassing a wide range of iron and acid concentrations showed that coarse limestone, one of the cheapest neutralizing agents known, is potentially useful for treating mine water discharges. However, process variables must be defined more completely before the practical applications and limitations of the process can be fully established. In all the tests a pH of >7.2 was obtained ultimately with limestone. Reaction time to reduce the iron content below 7 ppm in the water was dependent on the original ferrous iron concentration. For waters with high ferrous iron content, treatment with limestone followed by lime rapidly produced an acceptable treated water, Sludges with good settling and handling characteristics resulted when limestone was used for neutralization either alone or with subsequent treatment. Reaction rates are expressed graphically as pH change and iron depletion curves; characteristic solids settling curves are shown. The temperature dependency of ferrous iron oxidation was verified.

Jan 1, 1967

By H. I. Aaronson, H. A. Domian, G. M. Pound

Zener's two-parameter theory of the y a reaction in Fe-X alloys is extended to encornpass austenite-stabilizing as well as fewite-stabilizing elements, and is then cottzbitzed with statistical thermodynamic theories of Fe-C alloys to permit development of a ther-modynamic description of the proeutectoid fewite reaction in Fe-C-X alloys. Equilibrium tielines are calculated for tile a + y rep'on from experimental data and frort extrapolations of the y/y +a equilibrium surface. Sectzovs of the y/y + 0 surface are calculaled for the tzetastable sstuation LYL u-kich alloyrtg eleitzents do not partition betzseen azrstenite andfer-rite. The finding that the metastable y/a + y curves Lie close to their equilibriuwz counterparts when X = Si, .Vo, Co, Al, and Cu but well below them when X = A/ln or Ari pr-elides a thermodunamic explanation for the expcri1?7e?ztally obserted absence of partition during the proeutec toid ferrite reartlon in the fortner alloys and for the occurrence of partition at hzgher ternperatures in the latter alloys. Cotarison ofno- patitzoa free-energy changes and y/y i a curr,es in Fe-C-X alloys with the equilibrium 11al1es of these quantities in Fe-C alloys furnishes additional qualitative insight into the irqluence of allojling elements upon the kinetics of the proeutectoicl ferrite reaction. THE introduction of a substitutional alloying element, X, appreciably complicates calculation of the thermodynamics of the formation of proeutectoid ferrite from austenite. As noted in the preceding paper,' the positional entropy of the interstitial species is the principal component of the free energy of an interstitial solid solution with which theory has so far been able to deal. One would expect, however, that other components of the free energy of this type of solid solution may be significantly altered by the addition of a substitutional alloying element. Even if the assumptions are made, by analogy to the case of Fe-C alloys,' that changes in the positional entropy of carbon represent a significant part of the thermodynamic effects of an alloying element, and that the remaining effects can be taken into account by fitting the equations developed on this basis to experimental data on Fe-C-X phase diagrams or on the activities of carbon in alloyed austenite and ferrite, the experimental information available on either of these quantities is not yet sufficiently accurate or extensive, respectively, to make such an approach useful. An attempt made by Zener on the former basis to explain the effects of alloying elements on the thermodynamics of Fe-C alloys in terms of a temperature-independent "free-energy change" (actually enthalpy hane') required to transfer 1 mole of an alloying element from austenite to ferrite, in which the "free-energy change" was determined by fitting the theoretical relationships developed to Fe-C-X phase diagrams, thus proved inadequate in part because of deficiencies in the available ternary-phase-diagram data.= Other difficulties of a more fundamental nature, however, also indicated the desirability of a different approach to the problem.3y4y6 zener6 subsequently proposed that the free-energy change associated with y — a transformation in pure iron can be decomposed into magnetic and nonmagnetic components. Alloying elements were assumed to affect those components separately. One parameter was used to describe the quantitative effect exerted on each component. These parameters were evaluated from Ms (martensite-start) temperature data and from other experimental information usually either readily available or measurable with acceptable accuracy. Zener applied this treatment only to the calculation of Fe-X phase diagrams in which X is a ferrite-stabilizing element. In the present study, this treatment is extended to include austenite-stabilizing elements, and then combined with treatments previously considered for Fe-C alloys1 to permit calculation of the thermodynamic quantities of interest in the austenite - proeutectoid ferrite transformation in Fe-C-X alloys, where a number of representative, and commonly used alloying elements are chosen for X, including Si, Mn, Co, Mo, Al, Cr, and Cu. The results are used to explain important features of the partition of alloying elements between austenite and proeutectoid ferrite, as reported in a companion paper7 on the basis of electron-probe analysis, and to provide some additional understanding of the influence of alloying elements upon the kinetics of nucleation and growth of proeutectoid ferrite. Portions of this treatment6 have been clarified in the course of a review by Kaufman and Cohen.3 The consolidated summary of these developments with which this section is begun provides a basis for further clarification of the treatment, from which extension to include austenite- as well as ferrite-stabilizing elements is a natural consequence. The division of the free-energy change associated with the 7 - a transformation in pure iron, into two independent components is formally stated as:

Jan 1, 1967

By Manuel Gomez

Published experimental coke data from a pilot plant slot-type oven were examined by statistical regression analysis. Quantitative relationships were found to exist between the coking rate, bulk density, and fine grinding of the coal, sole-oven expansion, volatile matter, moisture, ash, and carbon con-tent of the coal and coke parameters indicative of size and strength. Regression equations developed from coke data were shown to yield computed values that were in good agreement with experimental results for 11 coke properties. Prediction equations developed were evaluated with an independent set of data. Present results indicate that the variables cited undergo much interaction and the sum of these interactions determines coke size and strength.

Jan 1, 1967

By E. A. Sondreal

The Bureau of Mines investigated the use of lignite char for producing carbon disulfide; CS2 production rates were as high as when wood charcoal was used as raw material. The maximum space velocity for CS2 produced in the bench-scale study was 9,280 hour-1, reacting sulfur and lignite char at 792° C. The reaction temperature and the temperature at which the char was carbonized were the two most important factors in determining CS2 production rates. Increasing the carbonization temperature in nitrogen from 400° to 900° C advanced the temperature at which carbon disulfide was first produced from 550° to 850° C. The average activation energy for reacting chars with sulfur was 44,000 cal/g mole. The reaction was first order with respect to the partial pressure of sulfur. Increasing the rate of sulfur admission caused a decrease in rate of production at otherwise similar conditions. Rates were not affected by either char size or sodium carbonate catalysis. High space velocities at moderate temperatures using inexpensive lignite chars could offer an economic advantage in the manufacture of carbon disulfide provided that yields were attractive.

Jan 1, 1967

Jan 1, 1967

The Ringelmann Smoke Chart fulfills an important need in smoke abatement work and in certain problems in the combustion of fuels. A knowledge of its history and method of preparation is, therefore, of interest to many. Since instructions on its use are not shown on the recent edition of the chart, those included in this revision of the previous Bureau of Mines publication now are a necessary complement to the chart, More detail regarding the use of the chart is included than was given in the earlier version.

Jan 1, 1967

By Richard L. Ash

Usually stresses generated by explosive charges during blasting are considered to propagate through materials in a spherical manner. Although conical stress-waves had been suggested by a few investigators they have never been clearly defined or critically observed. An experimental technique using models and lengths of Mild Detonating Fuse was developed to produce conical stress-waves and study their effects. Several materials with different properties were tested to compare the influence of those properties on wave shape and crater formation. The investigations showed that whenever column charges were end-initiated, conical stress-forms were produced if the explosive's velocity exceeded the energy propagation velocity of the surrounding material. The stress-cone angle significantly affected fragment size and shape as well as the geometry of craters produced.

Jan 1, 1967

By E. F. Foerster

The Bureau of Mines used a diphenyl ether calorimeter to determine the heat contents of natural and synthetic blast-furnace slags. For solid slags, enthalpy (HT-H298) at 1,448° to 1,680° K ranged from 282.1 to 348.5 cal/g (±3 percent); liquid slag values were 400.1 to 514 cal/g (±3 percent) in the temperature range between 1,626° and 1,890° K. Experimental results for liquid slags were 14.3 to 16.2 percent higher than the enthalpies calculated from thermodynamic tables; values for solid slags ranged from 5.3 percent lower to 3.9 percent higher than calculated values. Experimental results for liquid synthetic slags ranged from 14.9 to 24.3 percent higher than the calculated enthalpies; values for solid synthetic slags were from 1.0 to 3.4 percent higher.

Jan 1, 1967

By E. A. Breitenbach

A computer program has been developed to afford rapid and complete quantitative log analysis for exploration and production decisions. The computation consists of automatic selection of tops and bottoms of porous intervals from the digitized data, and then point-by-point calculations within each selected interval. Nearly all log types can be analyzed. This paper presents the calculation techniques found to be appticable to machine evaluation and gives examples of their use. INTRODUCTION As the quantitative interpretation of well logs entails repetitive use of charts and equations, it is natural that digital computer programs would be written. A Universal Log Interpretation Computer Program (ULICP) has been developed to afiord rapid and complete analysis for exploration and production decisions. Application of this program enables the log analyst to: (I) apply rapid quantitative log analysis for exploration and production decisions on either a field, a well or on a single zone basis; (2) apply concepts of interpretation requiring detailed numerical analysis; (3) analyze all porous intervals on each log, rather than a few selected zones, and complete the analysis in less time than previously required. Computations are reported for every digitized point in each zone; (4) use empirical techniques applicable to a given area as an integral part of the computation; (5) experiment with the empirical coefficients and exponents in the interpretation equations to find the best possible solution; and (6) automatically plot both the original and computed data to scales analogous to the field prints. The primary utility of a log interpretation program stems from its ability to do an overwhelming amount of work with very little man-power. If a calculation procedure or thought process can be formalized to the extent that step-by-step logic can be written, a computer program can be developed that follows this logic. The major question is one of economics. A feasibility study of the costs for such a program indicates that digital processing is economical primarily as a means of increasing the productivity of the log analyst. In effect, the probability of missing productive intervals in any well, because of a lack of time to do detailed calculations, is greatly reduced. The nucleus of ULICP is programmed to compute large sections of a suite of logs by selecting zones automatically and then performing all pertinent computations on a data point by data point basis within each zone.'.' An entire suite of logs can be processed in this manner with very little manual intervention. Sufficient programming logic is available so that each log analyst can request computations pertinent to his area. These requests are made by simple additions or deletions of information on the input header cards. Hence, the log analyst is always in complete control of the computation process. The evaluation of a suite of logs requires pre-editing, digitization, computation, presentation of results and interpretation. Work by the log analyst is reqilired only in pre-editing and interpretation. Thus, he is allowed more time for comprehensive interpretation, rather than calculation. For continuity, the discussion of ULICP is organized sequentially: pre-editing, digitization, computation and presentation of results. DISCUSSION PRE-EDITING The extent of pre-editing prior to computation is dependent on the format of the original data. For analog prints, it requires inspection, correlation and editing of the logs, plus the entering of required data an special forms. For digitized data such as magnetic tape field recordings, only the special forms are necessary. The process for analog prints will be given here. The inspection and correlation process involves the selection of sections for digitization and the correlation of the traces to a common depth. To decrease the cost of digitization, traces that exhibit large variations in shale formations can be redrawn to a non-zero baseline. The next step is to enter all pertinent data on the input header cards. The header cards presently used for ULICP are presented as Figs. 1 through 6. Cards 1, 2, 3 and 4 (Fig. 1) are defined as the Main Header Cards. They describe the particular well for output identification and give basic information. Cards 5 through 15 (Figs. 2 through 6) are defined as the Block Header Cards. As such, they define the log types and the interpretation parameters for the block of data immediately following. Cards 1 through 10 are required for every computer run. Cards 11 through 14 pertain only to nuclear logs. Card 11 is required if any nuclear log is supplied, and cards 12 through 14 are required only when gamma-ray, density or neutron log data, respectively, are supplied. Card 15 (Fig. 6) must be supplied when cross-plot calculations are required. Either two- or three-component

Jan 1, 1967

By M. A. Jones

Performance of a pilot flood in eastern Kansas indicates improved recovery and accelerated production resulting from mobility ratio control obtained by adding a high moleculur weight polymer to injected water. Ultimate pilot flood recovery of the 75-cp reservoir oil is anticipated to he 350 bbllacre-ft during the six-year flood life. Cumulative recovery to date totals 242 bbl/acre-ft after 30 months of injection. Reservoir parameters and flood performance are disc~rssed in detail to provide a case history relating to the feasibility of polymer flooding. INTRODUCTION The Vernon Polymer flood of the Brazos Oil & Gas Co., located in Sec. 2, 24 S., 16 E., Woodson County, Kans., was begun in late Oct., 1963. Injection of a 500 ppm solution of Pusher* polymer was initiated in a 15-acre pilot to determine feasibility of improving oil recovery by mobility ratio control. Performance to date and the anticipated recovery of 115,000 bbl from the 329 acre-ft pilot indicate efficient displacement of the viscous reservoir oil. Dependence of flood recovery on the displacing fluid-displaced oil mobility ratio has long been recognized. Displacement efficiency of the reservoir oil decreases with increasing mobility ratio. Therefore, economic reduction of this parameter is an attractive approach to improving oil recovery. particularly from viscous oil reservoirs where mobility ratios are often quite high. Waterflood mobility ratios can be controlled by dilute solutions of certain polymers'; which provide a displacing fluid with considerably reduced mobility. Reservoir parameters and flood performance are emphasized in this paper to provide data relating to polymer flood applicability. THEORY The polymer used in the Vernon pilot exhibits the property of unusually high resistance to flow through porous media.'.' A small concentration in injected water effects a mobility reduction considerably beyond that predicted from increases in bulk viscosity as measured in a capillary or rotating cup viscosimeter. Water mobility in the subject flood is reduced by a factor of 7.55 with the addition of 500 ppm Pusher chemical to injected water while the bulk viscosity is increased from 1.03 to only 1.40 cp at 75F. The exact flow mechanism causing resistance to flow has not been established; however, it is observed only in tortuous flow systems and appears to be influenced by brine properties, rock properties and oil saturation. Consequently, mobility reduction is estimated from flow studies for a given sample of reservoir rock and is expressed in terms of a "resistance factor".' This parameter is defined as the ratio of brine mobility to polymer solution mobility measured at residual oil saturation: RF= ?m/?p= (km/µm)/(kp/µp)M.......(1) where RF = resistance factor (dimensionless) k*. and A, = mobility of water and polymer solution, respectively k,,, and k,, = water and polymer permeability, respectively (darcies) µie and µp = viscosity of water and polymer, respectively (poise). Eq. 1 is similar to the expression commonly used to relate the relative mobility of the water bank at residual oil saturation to oil-bank mobility at irreducible water saturation: where M =mobility ratio of water to oil (dimension-less) A, = mobility of oil k, = permeability to oil (darcies) )i* = oil viscosity (poise). From Eqs. 1 and 2 it is apparent that the polymer-solution-to-oil-mobility ratio can be approximated by dividing the brine-to-oil ratio by the resistance factor RF: Mpn=M10/RF...........(3) where M,-, = mobility ratio of polymer solution to oil (dimensionless). This simplified presentation assumes no dilution of the polymer solution by connate water. Also, the presence of a connate water bank separating the polymer solution from the displaced oil is not considered. Displacement of oil by waterflooding is increasingly dominated by viscous forces as the viscosity of the oil displaced increases. When the oil mobility is less than that of the injected water, the oil-water interface is not a piston-like front.3,' Instead, instabilities occur at the flood front in the form of viscous fingers which tend to penetrate the less mobile oil bank. The magnitude of these protrusions, or the degree of channeling, increases

Jan 1, 1967

By H. O. McLeod, J. M. Campbell

Data analysis of pentane and hexane adsorption from natural gas in a fixed bed of silica gel shows that constant length mass transfer zones form, the curvature of the adsorption isotherm controls the growth of the mass transfer zone and surface diffusion of molecules inside the silica gel particle controls the mass transfer rate. Curvature of the hexane isothem is more than the curvature of the pentane isotherm. Because of this curvature the hexane adsorption zones reached a constant length. In contrast, the pentane adsorption zones were always increasing in length during each run. A procedure was developed to obtain correct mass transfer coefficients using effluent curve slopes. These transfer roefficients increase with the amount of hydrocarbon adsorbed on the silica gel particle. The characteristic shape of the hexane effluent curves also show that molecular diffusion inside the silica gel particle controls the adsorption rate of pentane and hexane. INTRODUCTION The purpose of this study was to determine the mechanisms that control the dynamic adsorption of hydrocarbons from a natural gas onto silica gel. Before one can deal effectively with multicomponent adsorption, the transfer mechanisms by which a single hydrocarbon component is adsorbed from the gas stream must be defined. Two principal investigations of this system have been published ' and indicate that diffusion through the gas around the particle controls the adsorption rate. Some of the experimental observations in each study either do not support this transfer mechanism or are inconsistent with the mathematical model used in analysis. In this study surface diffusion of molecules inside the particle controls the mass transfer rate of pentane and hexane. This mechanism is indicated by the effluent curve shape for a constant length transfer zone and by the variation of the mass transfer coefficient with concentration of the adsorbed hydrocarbon. THEORY AND DEFINITIONS — MATHEMATICAL MODELS Mathematical solutions for the isothermal adsorption of a trace component from a carrier gas are derived from three relationships: the mass balance or continuity equation, the equilibrium relationship between the gas and solid phases, and a mass transfer rate equation. The transfer rate is proportional to the adsorbate concentration gradient within either the gas or solid phase. Mathematical solutions of these equations usually give the adsorbate concentration as a function of time and distance from the bed inlet. That part of the bed in which the adsorbate concentration changes from a maximum to a minimum value is called the transfer zone. This transfer zone is directly related to a plot of the effluent concentration vs time which has a characteristic S-shape. This general shape is determined by the continuity equation and occurs in many processes of diffusional transfer.

Jan 1, 1967

By K. P. V. Lei

The Bureau of Mines investigated molten salt electrorefining 'n helium atmosphere cells as a means for the production of ductile vanadium. 'Electrolytes composed of NaCl-VCl2, LiCl-NaCl-VCl2, BaCl2-KCl-NaCl-VCl2 KCl-LiCl-VCl2, and CaCl2-NaCl-VCl2 were tested. Two electrolytes, KCl-LiCl-VCl2 and CaCl2-NaCl-VC12, were selected for laboratory-scale production tests. A molten salt electrorefining process was developed for the production of ductile 99.6 percent vanadium from a commercial 90-percent grade. Using the CaCl2-NaCl-VC12 electrolyte, an 87-percent recovery of ductile vanadium was made; an 85-percent recovery was made with the KCl-LiCl-VCl2 electrolyte.

Jan 1, 1967

The 63rd Regular Meeting of the Rocky Mountain Coal Mining Institute convened at the Stanley Hotel, Estes Park, Colorado, and was in session at 10:00 a.m. The first session was called to order by Past President George A. Nugent, Vice President and General Manager, Big Horn Coal Company, Sheridan, Wyoming. He thanked and offered his congratulations to the Program Committee for a very wonderful job, namely, John Reeves, Chairman, E. G. Fox, Lee Jensen, O. F. Lewis, J. Peperakis, R. M. von Storch, J. E. Ward and R. R. Williams, Jr. He thanked the suppliers for a fine cocktail party the evening before. Mr. Nugent's Opening Address and introduction of President James Westfield, Assistant Director, U.S. Bureau of Mines, Washington, D.C., follows: It is a real privilege to call this 63rd Regular Meeting of The Rocky Mountain Coal Mining Institute to order and to say welcome.

Jan 1, 1967

By J. D. Marchant

The Bureau of Mines studied solid-state electromigration to determine the value of the technique as a method of purifying cerium metal. Cerium metal in bars 6 inches long and 1/2 inch in diameter was refined at 600° C by solid-state electrolysis to remove iron, copper, and manganese. Iron and manganese were electrotransported to the anode and copper was electrotransported to the cathode. The concentration of iron at the anode of a cerium bar was 800 times greater than the concentration at the cathode. The concentration of manganese at the anode was twice the value at the cathode. The concentration of copper was 170 times greater at the cathode than at the anode. A radiotracer technique was used to follow the iron migration in the cerium metal. Molybdenum, silicon, and aluminum impurities were not electrotransported.

Jan 1, 1967

By Frank E. Williams

This Bureau of Mines report deals with the urbanization of Tucson, Ariz., a southwestern city with an arid environment. Zoning regulations concerning known aggregate sources are described, and information relating to projected needs for the urban area is furnished. History of growth is outlined not only for the local aggregate industry, but also for the entire urban area. Past conditions are pointed out and background material for future conditions is provided.

Jan 1, 1967

By S. G. Gibbs, A. B. Neely

Jan 1, 1967