Search Documents

By C. H. Li, R. J. Stokes, T. L. Johnson

The effect of alumina particles on the nucleation and growth of cracks through a silver- chloride matrzx has been investigated. It has been found possible to promote fibrous cracking in dispersion-strengthened silver, chloride under notch-impact conditions at temperatures at which silver chloride alone cleaves brittlely The modification) of fracture beIzavzor is thought to be due to the relaxation of hydrostatrc stress beneath a notch by the nucleation of cavities near alumina particles. In recent years, composite or dispersion-strengthened materials have been studied primarily to understand their high resistance to plastic flow particularly at elevated temperatures. Dislocation models have been developed with which it is possible to deduce with fair success the effects of interparticle distance, particle size, temperature, upon yielding and creep behavior.l-4 Much less attention has been paid to the fracture behavior of these materials (with the notable exception of common structural steels) and little is known experimentally about the manner in which inclusions affect the nucleation and growth of cracks through a matrix. Nevertheless a beginning has been made in connection with fibrous cracking in ductile matrices where inclusions appear to play an essential ro1e.5-7 During the severe localized plastic deformation which accompanies necking in a tensile test, cavities are believed to develop at inclusions; these cavities subsequently grow and coalesce by plastic flow until separation is complete. It is of interest to consider whether inclusions can affect fracture behavior under loading conditions which restrict the plasticity of the matrix itself (for example, cleavage under conditions of a high imposed strain rate at low temperatures). It is particularly interesting to study these effects in a solid which shows a spectrum of behavior ranging from fully ductile to semi-brittle behavior. Such a solid is silver chloride whose mechanical behavior depends sensitively upon temperature and strain rate.'," The present paper is concerned with a study of the influence of inclusions (in the form of alumina particles) on the fracture behavior of silver chloride loaded uniaxially at low strain rates at room temper- ature and also under notch impact conditions over a wide range of temperature. In particular, it will be shown that the alumina particles can exert a startling effect on the ductile-brittle transition temperature of notched silver chloride and that the magnitude and nature of the effect depends upon both the quantity of alumina and the shape of the alumina particles. 1. EXPERIMENTAL PROCEDURE 1.1 Materials Used. Silver chloride powder of analytical reagent (AR) quality having an average particle size 6 was supplied by the Mallinckrodt Chemical Works (st. Louis, MO.). Acid washed 900 mesh alumina powder, designated A38-900, was supplied by the Norton Co. (Cambridge, Mass.). This powder was added to silver chloride in two forms: a) the as-received condition in which the individual particles were of random irregular shape; their statistical average size was 7; b) in a condition in which each particle was spherically shaped by a fusion technique. In this case, the statistical average particle size determined with the optical microscope was approximately the same (about 5p) but electron micrographic evidence indicated that many ultra fine particles were present in the spheroidized powder. 1.2 Preparation of Composite Materials. Silver chloride-alumina composites containing 2.5, 5, and 15 pct by volume of alumina were produced by the extrusion of mechanically mixed powders blended in a ball mill for 24 hr at room temperature. The mixtures were compacted at 50,000 psi at room temperature in the form of billets 3/4 in. in diameter and 1 in. long which were then extruded with a 16:l reduction ratio at 370°C through a radius-type steel die having a 5 deg lead-in angle. An extrusion temperature of 370°C was selected to ensure that all composites had sufficient plasticity to be extruded. Apart from this general requirement, the choice was arbitrary. 1.3 Microstructure of Composites. Attempts were made to check the distribution of alumina metallo-graphically by polishing transverse and longitudinal sections of the extruded rod. Specimens were wet-ground to 600 emery paper and lapped successively with 5 and 0.25-p grades of diamond paste. They were etched for 10 sec in 10 pct sodium thiosulfate solution and lightly polished in concentrated ammonium hydroxide. The most effective way to render the alumina particles and grain boundaries visible was to radiate the surface with intense white light to decorate the grain boundaries and the particle-matrix interfaces photolytically.

Jan 1, 1962

By Arne Bergholm

Australian rutile was chlorinated in the presence of CO or carbon. The chlorination velocity in CO was found to be strongly influenced by temperature and proportional to the CO concentration, but independent of the Cl, concentration. In the presence of carbon, the reaction velocity is much higher. The reactivity of the carbon and the distance between the carbon and the rutile surfaces are important variables. The reaction velocity is approximately proportional to the Cl, concentration and independent of the CO concentration of the surrounding atmosphere. Experiments with fluidized-bed chlorination of carbon-rutile mixtures indicate that the motion of the bed has little influence on the reaction velocity. At low temperatures, the chlorination velocity of dense tablets is much greater than that of TiO, coke mixtures suitable for fluidization. The reaction mechanism is discussed. In the industrial production of TiCl, rutile is chlorinated in the presence of carbon. Disregarding intermediate steps, the reaction may be expressed by the following equations: The ultimate object of this study was to find out whether the reaction velocity was higher in fluidized bed operation compared with chlorination of pelle-tized carbon-rutile mixtures. From a literature survey and preliminary experiments it was learned that some basic knowledge about the reaction mechanism was needed for a good experimental design. Therefore the following sets of experiments were carried out: 1) Studies on the reaction velocity in the chlorination of rutile with CO as the only reducing agent. 2) Chlorination of separate rutile-carbon tablets. 3) Chlorination of rutile-carbon tablets at different temperatures with various kinds of carbon, various grain sizes, and various tablet-making techniques. This series of experiments was carried out not only with pure chlorine but also with mixtures of chlorine with argon, CO and CO,. 4) Chlorination of rutile-carbon tablets made in a strictly standardized manner. 5) Chlorination of static-bed rutile-carbon mixtures. 6) Fluidized-bed chlorination of the same mixture. The experiments 4 to 6 formed the final and direct test of the main question: Static bed vs fluidized bedo. Although there are many patents and papers dealing with the general aspects of chlorination, only few experiments from which detailed information can be obtained have been reported in the literature. Pamfilov and coworkersl3 have studied chlorination of TiO, with CO or carbon as the reducing agent. They found that the weight decrease per hour at 600o to 800°C amounted to 11 to 14 pct with CO and 45 to 51 pct with carbon. They suggest that phosgene might be an intermediate in the chlorination of TiO,. Takimoto and Hattori have chlorinated reduced titanium oxide (TiO) and found a very high rate of TiC1, -production. They proved that the composition of the gas from the chlorination of rutile-carbon mixtures contained mainly CO. They reported for instance 74.7 pct CO, and 5.6 pct CO at 800°C at which temperature the Boudouard equilibrium composition is 12 pct CO and 88 pct CO. Seligman and Segerchano6 have studied the chlorination of TiO. They proved that Ti0 and chlorine react rather rapidly at temperatures as low as 300°C. Above 400°C, the reaction was complete. At 500°C the velocity of this reaction was much higher than was chlorination of TiO, + C. McTaggart,7 Nishimura, et al. and Wilskam have chlorinated mixtures of carbon and various kinds of rutile or beneficated ilmenite. Nishimura, et al. also report the results from reduction of TiO, with H2, CO, or carbon. At 900°C only 1 pct Ti O is formed in 2 hrs. At 1100°C 23.1 pct Ti, 0, was formed if elementary carbon was present. No carbide formation occurred below 1400°C. McIntosh and Cofferll have observed that the CO, content of the exit gases from chlorination of rutile and calcined petroleum coke is appreciably higher than found in the Boudouard equilibrium. At 900°C the ratio (CO, /CO + CO) is about 80 pct, whereas the equilibrium value is about 2 pct. W. E. Dunn12 has studied the chlorination rates of several TiO,-bearing minerals with CO + Cl, or COCl . Chlorinations were carried out either in a fluidized reactor or a fixed-bed reactor, both having a 30-mm diameter. The results obtained in both reactors were comparable. It was proved that benefi-ciated ilmenite (i.e., ilmenite from which the iron oxide had been removed by chlorination) was chlorinated 10 times faster thanrutile. Sore1 slag showed an intermediate rate. When phosgene is used, the

Jan 1, 1962

By N. A. Gokcen

The construction and operation of a simple and accurate vacuum-fusion apparatus are described in detail. Absolute accuracy of the oxygen analysis has been determined by the reduction of oxides weighed to 20.01 mg. The Bureau of Standards steel samples have been analyzed repeatedly, and their inadequacy for the acceptable range of oxygen is presented. The limits of interference of manganese, aluminum, and titanium, and the effects of tin and water-cooling of the furnace tube have been investigated in detail. ANALYSIS by the vacuum-fusion method consists of melting a sample in a degassed graphite crucible under high vacuum and extracting oxygen as carbon monoxide, and nitrogen and hydrogen as gaseous elements. The resulting gas mixture is analyzed for its constituents, from which oxygen, nitrogen, and frequently hydrogen may be determined within 1/2 hr or less. The history and summary of methods1-4 and extensive reviews of literature'.' for the determination of gases in metals have been published in detail. The foundation of the modern vacuum-fusion method of analysis was laid by Jordan and Eckman,11 and Oberhoffer and his associates.1, 10 Diergarten,11, 12 Meyer,13 Thanheiser et al.,14, 15 Ericson and Benedicks,16 Ziegler,17, 18 Sloman,16-23 Thompson et al.,24 and later many others contributed to the improvement, accuracy, and limitations of the method. In some investigations, however, the limitations of the procedure and the errors involved in the analysis require further critical examination and evaluation. In many apparatuses, lunnecessarily complicated and cumbersome features and elaborate precautionary measures do not have conclusive advantages. The purpose of this invest:igation was, therefore, a) to construct a very simple and accurate apparatus, b) to determine the absolute accuracy of results, c) to establish the limits of interference of manganese. aluminum, and titanium, d) to reexamine critically the reduction and recovery of oxygen from oxide powders of various sizes, and e) to evaluate critically the eight steels of the Bureau of Standards."' Apparatus The apparatus used in this investigation is shown in Fig. 1. It has been used for over three years, during which several modifications were made. A brief description of the apparatus is as follows. The transparent silica or Vycor tube, E, is joined to the Pyrex head, B, by mean of a face to face ground joint sealed with mercury. A similarly sealed ball and socket joint would also have been satisfactory. A standard tapered ground joint, successfully used by some investigators, was tried, but in the absence of grease or vacuum cement it was very difficult to disassemble the joint. The graphite crucible, G, 15/16 or l 1/8 in. OD, 2 1/2 or 3 in. high, 1/10 in. wall, shown enlarged in Fig. 3a, is made from the rods containing less than 0.08 pet ash. Two baffles, overlapping sufficiently to prevent spattering, are inserted in two opposite slots cut on a band saw. Two small holes in the baffles permit the temperature measurement of the melt. A second type of crucible, Fig. 2, G, shown enlarged in Fig. 3b, is also satisfactory, though not as convenient as the crucible shown in Fig. 3a. Somewhat similar but more elaborate crucibles with baffles were first tried by Ericson and Benedicks,'" but were abandoned in favor of a crucible with a stopper and a peripheral graphite powder filter. Various versions of their crucible with a stopper have later been used by other analysts, 10-27, 27-30 The crucible is packed directly in E, Figs. 1 and 2, with —35 + 48 graphite powder for shielding and insulation. A minimum layer of 7 mm of powder is necessary for keeping E cool and minimizing the blank. The first use of graphite powder was made by Sloman,10-23 Who found that —200 mesh was the most satisfactory. The author tried various size powders from +20 to —100 mesh and found that a) particles of 28 mesh or larger were not sufficiently insulating and were heated by the induction current, b) the powder finer than 60 mesh packed too much and did not readily permit the escape of gases even when the crucible was heated very slowly; hence, the powder was occasionally blown out of the side in the crucible and in the mercury pump. The optimum size powder is, therefore, 35 to 48 mesh size or 0.4 mm in average diameter. The furnace is designed to eliminate excessive amounts of

Jan 1, 1959

By Thomas L. Garrard, James A. Sprague, Rex B. McLellan, Samuel J. Horowitz

A simple rnodel for interstitial solid solutions has been devised in which each solute atom interacts with the solzlent lattice in such a way as to exclude an integral number of nearest-neighbor sites from being occupied by solute atoms. The vemaining interstitial sites ave not affected. In calculating the entropy of such a solution the increase in the number of interstitial sites available to solute atoms due to the ozlerlap of blocked sites is taken into account. The ther-modynamic functions and the solute activity fov the solution rnodel are calculated. It is shown that considevation of the overlap sues vise to a correction term of the order 8' ( where 0 is the ratio of solute to solvent atoms) in the thevmodynamic functions of the solution and the solute activity. The theo,uetical equation for the solute activity was used to analyze the available data for the activity of carbon in austenite at three different temperatures. The use of a computer enabled many comparisons to be made between the experimental activity data and those predicted by both the complex blocking model and the simple, widely used blocking model in which the ozierlapping of blocked sites is not taken into consideration. Even though at the highest concentvations oj carbon the value of 0 is small it is shown that the best fit to the experimental data is found with the complex blocking model. RECENTLY it has been shown by analysis of gas-solid equilibriums that many interstitial solid solutions can be described by the quasi-regular solution model in which the partial energy of solution E, is independent of composition and, apart from a composition-independent term due to nonconfigurational factors, the partial entropy is that of an ideal solution. Dilute solutions of nitrogen in both bcc and fcc iron,' dilute solutions of oxygen and silver,' carbon in a Fe,' and dilute solutions of carbon in ? Fe2 are all quasi-regular. In addition many hydrogen-metal solutions, both dilute and concentrated (with the notable exceptions of the H-Ni and H-Pd systems), can be described by this model.3 However, it has been shown that, at various temperatures, solutions of carbon in ? Fe begin to depart strongly from quasi-regular behavior at low concentrations. A plot of the activity of carbon in ? Fe relative to pure graphite vs atom fraction of carbon (Cc) calculated from Smith's data4 begins to depart from the (nonlinear) plot for a quasi-regular solution at about Cc = 0.02. Several statistical models have been proposed to account for the thermodynamic behavior of carbon in y Fe. Darken and smith4 set up a model in which a carbon atom in its octahedral site in the fcc lattice has either one or no neighboring carbon atoms in the twelve nearest-neighbor octahedral sites. It is assumed that the number of carbon atoms having more than one nearest neighbor is negligibly small and that a solute atom does not interact with another solute atom located at a greater distance than the first shell of interstitial sites. Thus this model contains two characteristic interaction energies. Darken and Smith showed that this model gives reasonable agreement with the measured variation of carbon activity ac with composition. It is noteworthy that Darken and Smith concluded from their analysis of the activity data using this model that there was a repulsive force between carbon atoms in austenite which reduces the concentration of C-C pairs below that corresponding to random mixing. Recently Aaronson, Domian, and pound5 have considered in detail the statistical model developed by Lachere and by Fowler and Guggenheim7 and have shown that it is compatible with the activity data for carbon in austenite. In Darken's model the energy of solution of carbon in austenite is dependent on temperature and concentration; the departure from Henry's law arises from the concentration dependence of the energy of solution. In Aaronson, Domian, and Pound's discussion of the Lacher, Fowler, Guggenheim model the interaction energy wy between adjacent carbon atoms in austenite is strongly temperature-dependent. Aaronson et al . suggested that this temperature dependence of wy and hence the energy of solution of carbon could be explained either by the formation of clusters of carbon atoms (despite the repulsive nature of wy) or by the dual occupancy of the octahedral and tetrahedral sites by the carbon atoms. However, a temperature-dependent heat of solution

Jan 1, 1968

By S. Nakajima, H. Okazaki

Quantitatiue studies of the deformation texture in drawn tungsten wives were made by the X-vay dif-fractottletetr. Experimental results show that the diffraction Intensities are equal to tilose pvedicted from the (1 10). fiber lexlure but the angxla), spreads of. diffraction peaks in the pole distribution curres are different for different diffraction planes and directions. For this reason a modified (110) fiber lextuve model, in which a kind of anisotropy is assumed, is proposed to explain the results. According to this model the poles lying on a line directing front the (110) to the (110) poles in the (1 10) standard stereograpllic projection should show spreads which are different from those lyitlg on a line directing from the (001) to the (001) poles, which is confirmed by the experiments. The anisolvopy and the spveads of the pole positions are large at the outer part of the wires and decrease gradually lowards the inside of the wire. The possibilily of occurrence of such anisolropy in irrelals with fcc stvuctures is discltssed. THE deformation texture of drawn tungsten wires has been assumed by different investigators to be the simple ( 110) fiber texture.' Recently, however, Leber2,3 has shown that a swaged tungsten rod has a cylindrical texture. It changes gradually to the (110) fiber texture by drawing through dies. However, even after drawing to 0.25 mm in diam, the cylindrical texture can still be found in wires together with the (110) fiber texture. This was deduced from the pole figures obtained from the longitudinal section of these wires. Use was made also of quantitative measurements of the pole distribution curves. Leber stated that the angular spread of the pole distribution curves (henceforward called dispersions) are quite different for (400) 45 deg and (400) 90 deg: the former is always larger than the latter. This inequality is accompanied by deviations of the diffraction intensities from the theoretical values for the ( 110) fiber texture. Bhandary and cullity4 have reported similar results on iron wire and explained them by assuming a cylindrical texture. Both Leber3 and Bhandary4 used only the results of the (400) reflection for the determination of the dispersion. The pole figures found by Leber3 and by Rieck5 are largely different. The model given by Leber to explain the effects is in the authors' opinion in some respects unsatisfactory, especially if one looks at other than the (400) reflections. Intensities and dispersions of diffraction peaks are conclusive factors for the determination of the fine structure in wire textures. For this reason we studied them extensively to come to a model which is more suitable to fit the facts. In the following, after giving the experimental set-up, we report about measurements of X-ray diffraction on drawn tungsten wires. Different models to describe the experimental results will be discussed. EXPERIMENTAL GO-SiO2-A12O3 doped tungsten wires drawn to 0.18 mm in diam were used for the measurements. The wires were chemically etched to various diameters down to 0.03 mm. Measurements were carried out for the different wires in order to determine the dependence of the texture on the radius. The wires were cut to pieces of 10 mm length and fixed with paste closely against each other on a flat, polished glass plate. Parallelism of the wires with the surface of the glass plate should be adequate. For the diffraction studies three different X-ray sources were applied, respectively, giving the CuK,, FeK,, and FeKp emission. The measurements were carried out with a diffrac-tometer with a GM counter. The latter was fixed to a certain diffraction angle 20hkl and the diffraction intensity was recorded as a function of the angle of rotation of the specimen around the axis, lying in the specimen surface and perpendicular to the wire axis, as shown in Fig. 1. Measurements were also done with the detector at angles slightly deviating from the diffraction maxima The measured intensities in this case were taken to be equal to the background level. The deviations were chosen as small as possible but large enough to eliminate the influence of the diffraction maxima. The useful range of the rotation angle x of the specimen is generally limited by the wavelength of the X-rays. We have: where and cp is the angle between the wire axis and the normal of the diffraction plane. Intensity measurements were made to find the necessary corrections for counting loss of the GM counter and for distortion resulting from such effects as absorption of X-rays and from inclination of the reflection plane under study with respect to the surface of the specimen. The counting loss was esti-

Jan 1, 1968

By R. C. Ku, A. J. McEvily, T. L. Johnston

The microplastic response of a series ofas-quenched Fe-Ni-C martensites has been measured at 77°K. At strains less than JO'3 the flow stress is governed primarily by the transformation-induced dislocation structure of the martensite. Only at strains in excess of 10-3 is the influence of carbon manifested in the flow stress. At these macroscopic strains, typically 10-2, the solid-solution hardening is proportional to (wt pct C)1/3, and, in an alloy containing 0.39 wt pct C, amounts to 50 pct of the flow stress. THE technological significance of high-strength ferrous martensite has stimulated many investigations of its structure and properties. Although our knowledge of the characteristics of martensite has increased immensely, especially with the advent of high-resolution techniques, an understanding of the basic strengthening mechanism still remains elusive. The purpose of the present paper is to consider certain aspects of micro-plastic behavior of Fe-Ni-C martensite which we feel can help to resolve this important problem. Such alloys are particularly suitable for experimental investigation because their compositions can be adjusted to reduce the M, to a temperature low enough essentially to eliminate the diffusion of carbon in the freshly formed martensite.1 The mechanical properties in this condition are of interest inasmuch as they reflect a state that is free of the important but complicating influence of precipitation processes. In this virgin martensite the carbon is distributed as it was inherited from the parent austenite; i.e., it is present interstitially, and gives rise to tetragonality through strain-induced ordering.' In order to determine the source of strength of such alloys, Winchell and Cohen1 investigated the low-temperature macroscopic stress-strain behavior of a series of virgin martensites of increasing carbon content but of common M, temperature (-35°C). They found that the flow stress increased rapidly with carbon content up to 0.4 wt pct; beyond this point the flow stress increased at a much slower rate. It was concluded that martensite is inherently strong. To account quantitatively for the strength of virgin or as- quenched martensite in terms of the role of carbon, Winchell and cohen3 suggested that the carbon atoms, trapped in their original positions by the diffusionless martensite transformation, interfere with dislocation motion according to a model akin to that of Mott and Nabarro. 4 In this treatment, individual carbon atoms are considered to constitute centers of elastic strain and thereby generate an average stress resisting the motion of dislocations throughout the lattice. The additional stress necessary to move dislocations, over and above that necessary for motion in a carbon-free martensite, is given by where L is an effective length of dislocation capable of motion. L was assumed to be limited to the spacing between the twins that are an essential structural element of Fe-Ni-C martensites. They assumtd the spacing to be invariant and of the order of 100A. However, recent work5 has shown that L is variable and can be in excess of 1000Å, so that the assignment of an appropriate value of L is not straightforward. In contrast to the above conclusion that there is an intrinsically high resistance to plastic flow, it has been suggested by Polakowski6 that freshly quenched martensite is in fact "soft" in the sense that dislocations are initially free to move upon application of stress. The high indentation hardness and macroscopic yield stress of ferrous martensites are then a consequence of rapid strain hardening that depends upon carbon in solution. Consistent with this point of view are the results of Beau lieu and Dubé who measured the rate of recovery of internal friction as a function of aging (tempering) temperature in a freshly quenched steel containing 0.90 wt pct C, 0.37 wt pct Mn, 0.1 wt pct Cr, and 0.07 wt pct Ni. The kinetics were clearly consistent with the idea that many dislocations are unpinned in the as-quenched state and that during aging they become progressively pinned by carbon at a rate controlled by carbon diffusion in the body-centered martensite lattice. In order to provide a basis upon which to distinguish between the "hard" and "soft" interpretations indicated above, we have made studies of the initial stages of plastic deformation in Fe-Ni-C martensites similar to those'used by Winchell and Cohen. It will be shown that the results support the contention that dislocation segments in as-quenched material are indeed

Jan 1, 1967

By A. L. Mincher, W. F. Sheely

Mechanical properties of columbium have been studied over the temperature range of -196 to 1093oC. The decreased strengthening influence of cold-work at temperatures below ambient has been interpreted in terms of the Peierls-Nabarro effect. Maxima in the rate of strain hardening observed during tensile testing in the range 250-600°C. have been correlated with interstitial impurities to indicate the temperature ranges at which carbon, oxygen, and nitrogen, respectively, are responsible for strain aging. THE growing need for structural materials for use above the useful service temperatures of the iron-, nickel-, or cobalt-base alloys has caused the refractory metals to be considered as potential engineering materials. These metals, which include columbium, tantalum, molybdenum, and tungsten, are called refractory because the lowest melting point among them,that of columbium, is about 1000°C higher than the average melting temperatures of conventional high-temperature alloys. They are all body-centered cubic transition metals and, as such, their mechanical properties have basic characteristics which distinguish them from the face-centered cubic metals. For example, all show a much steeper rise in strength with decreasing temperature below room temperature than do the face-centered cubic metals, and their mechanical properties are strongly influenced by interstitially dissolved impurities. In order that these new metals may be used efficiently, it is necessary that their characteristics of behavior be fully known. In this paper, the mechanical properties of columbium will be examined over a wide range of temperatures. In particular, the influences of cold-work and individual species of interstitial impurity atoms on mechanical properties will be described, and basic mechanisms which may control the observed characteristics will be explored. EXPERIMENTAL The material used in this investigation was Union Carbide Metals Co. columbium roundels consolidated to four 4-in. diam ingots, three by consumable-electrode arc melting and one ingot by electron beam melting. Impurity contents of the ingots and methods of ingot conversion and treatment are summarized in Table I. The only metallic impurity occurring in any significant quantity was tantalum at about 0.1 pct. Iron, silicon, titanium, and zirconium were each less than 0.015 pct; boron was 1 ppm or less. This should have no appreciable influence on properties. The electron beam melted material, being the purest, will be used as the basis for comparison in the discussions to follow. Tensile tests were conducted from-196 to 1093oC, on both cold-worked and fully recrystallized arc-melted and electron-beam melted columbium using standard 1/4-in. diam, 1-in. long gage length test specimens. A strain-rate of 0.005 in. per in. per min was employed until the 0.2 pct yield strength was achieved and then the strain-rate was increased to 0.05 in. per in. per min for the balance of the test. Samples were protected in an inert atmosphere at tests above 300°C. The tensile properties obtained on the electron-beam melted columbium, E, in both the cold-swaged and recrystallized conditions are given in Fig. 1. The yield strength data of Dyson, et al.,' obtained on recrystallized electron beam melted columbium and the tensile strength data reported by Tottle2 on powder metallurgy columbium are included in Fig. 1. The material used by Tottle had been purified by vacuum sintering. There is excellent agreement between Dyson's data and those obtained in the present investigation. The tensile strengths obtained by Tottle were slightly greater than those obtained in this investigation on electron-beam melted columbium but varied with temperature in a similar manner. Tottle's data showed a maximum in tensile strength near 500°C, as did our data on electron-beam melted material, and also showed a small maximum at 300°C. The significance of these maxima will become evident later in the discussion. The tensile properties of cold-swaged and recrys-tallized arc melted columbium are plotted in Fig. 2. It was found that the properties of the recrystallized arc-melted columbium from all three heats showed very close agreement except at temperatures between about 500" and 800°C. A reason for this range of disagreement will be suggested in the discussion. The generally good agreement, however, attests to the ability of cold-working and subsequent recrystal-lization to erase the effects of the three different primary breakdown procedures and to produce nearly equivalent structures in the samples derived from the three different heats. wesse13 reported tensile data on columbium having interstitial impurity contents between those of the

Jan 1, 1962

By H. I. Aaronson, C. Wells

Configurations of ferrite crystals have been found in a plain carbon steel which appear to have resulted from the nucleation of new ferrite crystals at the interphase boundaries of previously formed crystals despite the high carbon concentrations which necessarily develop at these boundaries. This phenomenon has been termed sympathetic nucleation. An attempt has been made to reconcile the occurrence of sympathetic nu-cleation with current nucleation theory. THIS investigation is one of a series on the formation of proeutectoid ferrite from austenite. From the viewpoint of chemical composition, this reaction consists of the nucleation and diffusional growth of crystals of carbon-poor ferrite within a matrix of carbon-rich austenite. The austenite adjacent to the austenite-ferrite boundaries will be greatly enriched in carbon, approximately to the value of the y/(a + y) equilibrium curve or its metastable extrapolation at the temperature of transformation. Those areas of austenite appreciably farther removed from the growing ferrite, on the other hand, will be relatively unaltered in composition, especially at the earlier stages of transformation. Since rates of nucleation are considered to decrease exponentially with decreasing supersaturation,' the frequency with which ferrite nuclei appear at austenite-ferrite boundaries should be negligible in relation to that at which they form in other regions of the austenite. During this investigation, however, many groupings of ferrite crystals have been found which appear to have resulted from the nucleation of ferrite at austenite-ferrite boundaries. This phenomenon has been given the name of sympathetic 71.1tcleation. A number of micrographs of morphological configurations caused by sympathetic nucleation will be presented, after which an explanation for this reaction will be proposed in terms of current nucleation theory. Some of the structures to be considered are composed of bainite, an aggregate of ferrite and carbide, rather than of ferrite. Since ferrite and bainite differ only in that bainite forms under conditions which result in the nucleation of carbides behind the advancing austenite-ferrite boundaries,' it will usually be unnecessary, for the purpose of this paper, to distinguish between the two reaction products. All studies were performed on an electric furnace steel (obtained from the Vanadium Alloy Steel Co.) containing 0.29 pct C, 0.76 pct Mn, 0.25 pct Si, 0.005 pct P, and 0.007 pct S. The alloy was cast as a 150 Ib, 7x7 in. cross section ingot and forged into bars 2x2 in. in cross section. These bars were homogenized for 48 hr at 1250°C in an Endo-Gas atmosphere. The depth to which decarburization penetrated during this heat treatment was determined by chemical and microscopic analyses and the affected metal was removed by machining. Specimens for isothermal transformation studies were cut from the remaining material; most of these specimens were 1/2x1/4X1/16 in., though some with a thickness of 1/32 in. were prepared for use at the shorter reaction times and lower reaction temperatures. Specimens were austenitized for 30 min at 1300°C, isothermally reacted for various times at temperatures ranging from 775" to 475 "C, and then quenched in iced water. The austenite grain sizes within individual specimens ranged from ASTM Nos. 1 through —4. A commercial heat-treating salt which was continuously deoxidized by an immersed graphite crucible served to minimize the loss of carbon during austenitizing; thick covers of powdered graphite and immersed graphite rods effectively prevented decarburization in the lead pots employed for the isothermal reaction treatments. The heat-treated specimens were sectioned and mounted in Bakelite. Following the completion of standard grinding and mechanical polishing procedures, the specimens were electrolytically polished with a Buehler-Waisman apparatus and etched in 2 pct nital. Experimental Results Rules of Evidence for Sympathetic Nucleation—On the basis of observations made on a single plane of polish, one precipitate crystal may be considered to have been sympathetically nucleated at the inter-phase boundary of another precipitate crystal when the following conditions are fulfilled: 1) The sympathetically nucleated crystal is not in contact with a grain boundary or a subboundary in the matrix phase. 2) The shape, size, and location of the crystal at whose boundary sympathetic nucleation occurred (hereafter termed the base crystal) and the crystal formed by sympathetic nucleation substantially pre-clude the possibility that the plane of polish em-ployed may have concealed the fact that both crys-tals actually nucleated at a grain boundary or a sub-boundary in the matrix phase.

Jan 1, 1957

By R. A. Wolfe, H. W. Paxton

Self-diffilsion coefficients for cr51 and Fe55 in 12 pct Cr-Fe and 17 pct Cr-Fe for Fe55 in chromium, and for Cr51 in vanadium have been measured. The results are compared with other values for the Fe-Cr system, and with the various theories of diffusion in hcc metals. Some empirical correlations are discussed between Do and Q in hcc systems, or, expressed differently, the constancy of ?G*/T solidus for seveval bcc metals and alloys is noted. It appears very probable that a vacancy mechanism is operative in bcc metals, hut this cannot he stated with certainty. THE great bulk of work on diffusion in metals, both experimental and theoretical, was for many years concentrated on those with close-packed and, in particular, fcc lattices.1,2 There appears to be little doubt that the mechanism of diffusion in these solids is vacancy migration, leading to mass transfer and in substitutional solid solutions to a Kirken-dall effect.3,4 For bcc metals, the picture is much less clear. The Kirkendall effect certainly occurs in several alloys.5-10 However, attempts to understand the factors contributing to the pre-exponential in the usual expression for the diffusion coefficient D =D, exp {-Q/RT) by extension of ideas useful in close-packed lattices have not always been successful. Zener,11 Leclaire,12 and Pound, Paxton, and Bitlerl3 have suggested that various forms of ring diffusion may be important in some bcc metals. For close-packed metals, Do is usually about 1 sq cm per sec and Q - 35Tm kcal per mole (Tm = melting temperature in OK). The theory of Pound et al. suggests for ring diffusion that Do may be about 10-4 and Q, although difficult to calculate with any precision, would be significantly less than 35 T,. The experimental results on self and solute diffusion in ? uranium14,15 and ß zirconium,10 and for solutes in 0 titanium,17 and possibly for self-diffu- sion in chromium below about 0.75 T,," gave some credence to this theory. However, not all bcc materials display low values of DO and Q, and the exceptions were not predicted by any theory. Furthermore, it has recently become apparent that, in bcc materials, log D is not always linear with T-l if a sufficiently wide range of temperature is studied.16,18 This variation may be such that Q may increase18,19 or decrease20 with increasing temperature. The present work was undertaken in an attempt to provide further diffusion data on bcc metals, and to try to understand the factors which contribute to differences in behavior between the various elements. For part of this work, the Fe-Cr system was chosen since it is of considerable technological importance, and data on 12 pct Cr and 17 pct Cr alloys appeared well worthwhile to supplement that existing for the remainder of the stern.18,22 The diffusion of Fe55 in chromium was studied as an example of a more or less "normal" tracer element in a possibly abnormal host lattice. Finally, no data were available for vanadium, the neighbor of chromium in the periodic table, because of lack of a suitable isotope so cr55 was used as a tracer in a few preliminary experiments. For convenience, we shall refer to elements whose Do and Q are low compared to those predicted by Zener's theory as "anomalous". PROCEDURE This investigation determined self-diffusion rates by means of radioactive tracers and the integral-activity method first utilized by Gruzin.23 In this method a thin layer of radioisotope of the diffusing element is plated or coated onto a planar surface of the diffusion sample, which is then given an isothermal-diffusion annealing treatment. The determination of an activity-penetration curve involves measuring the residual activity of the specimen after each successive layer or section has been removed parallel to the original planar surface. The method used here is essentially the same as that used by Gondolf18 and Kunitake.21 Two radioactive tracers, cr51 and Fe55, were used in this investigation. Diffusion coefficients were determined for the diffusion of one or both of these tracers in four different materials, viz., Fe-12 wt pct Cr alloy, Fe-17 wt pct Cr alloy, chromium, and vanadium. The diffusion samples had nominal dimensions of 1.5 cm diameter and 0.5 cm thickness. The grain size was several millimeters for the Fe-Cr alloys and at least 1 mm for the chromium and vanadium samples. Accurately planar surfaces

Jan 1, 1964

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 C. M. Schwartz, A. E. Austin, W. B. Wilson

High-temperature X-ray diffraction studies were conducted with Th-U alloys with up to 10 wt pet U. The solid solubility of uranium in thorium as a function of temperature was determined by the method of lattice parameters. Thorium will dissolve up to 2.5 wt pet U at 950°, 4.5 wt pet U at 1150°, and 7.5 wt pet U at 1250°C. Determinations were made of the temperature of the transition of thorium and of the effect of uranium upon the transition. The a to ß transition for thorium was observed to occur at 1330' ±20°C. Mean coefficients of expansion were calculated for thorium and two alloys, and for ThO2 in contact with the thorium, using X-ray lattice-parameter data. Values obtained at 950° C for thorium and Tho2 were 12.1 and 9.40 X 10-6 per OC, respectively. Impurities obtained during the X-ray exposure were identified by diffraction and were essentially Tho2 and ThC, with two additional unknown phases being detected. The effect of the impurities upon the results is discussed. DIRECT investigation (i.e., high-temperature X-ray diffraction studies) of the phase diagram of thorium-rich uranium alloys has been shown to be necessary since recent work1 - has disclosed the presence of an allotropic transformation near 1400°C in pure thorium with the room temperature face-centered-cubic phase transforming to a body-cen-tered-cubic structure at the elevated temperature. The effect upon the transition of the addition of uranium to thorium and of the solubility of uranium in thorium at high temperatures remained unknown, yet was of interest in understanding fabrication procedures and elevated temperature use. The present work was undertaken to provide information in this area by determining the transition temperature, the effect of uranium on the transition temperature, and the solubility of uranium in thorium as a function of temperature. Experimental Work The high-temperature diffraction data were obtained using a camera especially designed for the purpose3 and capable of reaching temperatures in excess of 2000°C at pressures as low as 1 x 10." mm Hg. Temperature regulation was provided by regulating the power input to ±0.1 pet variation, and by regulating the water flow through the camera jacket to provide a constant thermal load. The X-ray sample was a rod nominally 80 mil diam, which was further turned down to 20 mil and then etched to 18 mil over 1/2 in. of one end. This was placed in the sample holder and mounted on the camera so that the smaller part was surrounded by a cylindrical tantalum-sheet radiation-type heating element. Diffraction from the sample was recorded on film after passing through a slot in the heating element and radiation-baffle shield and through beryllium vacuum windows. The X-ray film mounting was of the Straumanis type4 with a camera diameter of 114.59 mm. Since previous work of Chiotti' indicated that impurities considerably alter the transition temperature, chemical analysis of the arc melted iodide crystal-bar thorium samples was obtained prior to testing. The analysis disclosed the material to contain as low as 0.001 ±0,0002 wt pet H and 0.007 ±0.001 wt pet 0. Carbon was 0.003 wt pet and nitrogen less than 0.002 wt pet. This material was sealed in mild steel in an inert atmosphere and subsequently hot rolled to 3/8-in. diam rods at a temperature of 732°C. Following removal of the jacket, the material was pickled and cold swaged to 1/8-in. rods, from which the diffraction samples were prepared. The alloys were similarly prepared, with the uranium being added during arc melting. The uranium analyses of the alloys prepared appear in Table I. The experimental procedure for diffraction examination of the three samples of high-purity thorium differed from those of the Th-U alloys. The original practice, later modified, consisted of pumping down the camera with the diffusion pump on and then admitting liquid nitrogen to the cold trap of the system. This was modified for the Th-U alloys by maintaining liquid nitrogen in the cold trap at all times before and while the diffusion pump was heated. This minor change produced a reduction in the amount of carbon pickup by the sample during exposure to the diffusion-pump vapors. The sample was brought to the desired test temperature and exposed for 21/2 hr at pressures which were usually 2 x 10-6 mm Hg, or lower. Exposures were made at

Jan 1, 1959

By E. Buehler, H. J. Levinstein

The temperature ranges in which the three inter-metallic phases in the Nb-Sn system form have been determined and the composition and structure of two of the three phases has been established. The kinetics of the formation of Nb3Sn in cored wire samples has been studied in the temperature range of 800° to 1050°C. From 800°to 950°C the rate of formation increases by four orders of magnitude. The rate-controlling step for the formation process in this temperature range appears to be the diffilsion of tin through NbSn. At higher temperatu~es a change occurs in the mechanism of the formation process such that up to a temperature of 1050°C the rate of formation of Nb3Sn does not increase above the rate observed at 950°C. For temperatures helow 950°C the current-carrying capacity of the wire increases with increased percent reaction reaching a maximum value when the formation process is 90 to 95 pct complete. The maximum current-carrying capacity obtainable in this temperature range is independent of the temperature. Above 950°C tlze current-carrying capacity obtainable in the wire decreases with increasing temperature of formation. A model is proposed which accounts for the ohserved behavior. RECENTLY, Buehler et a1.l reported the results of an investigation of the process variables which influence the superconducting properties of Nb3Sn-cored wire. These results indicated that at least four variables affect the properties of the manufactured wire. These include composition, particle size of the starting powder mix, temperature of heat treatment, and time of heat treatment. In order to understand completely the role of these variables, it is necessary to have an accurate knowledge of the phase equilibria in the Nb-Sn system. At the present time, phase-equilibrium diagrams for the Nb-Sn system have been published by a number of investigators.2-5 The diagrams differ as to the number of phases present, the composition of the phases, and the temperature range of stability of the phases. The present investigation was undertaken in order to resolve these differences. Since the investigation of Buehler et al. demon- strated that the length of time at the temperature of heat treatment affected the superconducting properties of Nb3Sn, it is apparent that it is necessary to understand the kinetics of the formation process as well as the equilibrium conditions before a complete understanding of the system is possible. As a result, the kinetics of formation of the various phases in the system were also studied in this investigation. EXPEFUMENTAL PROCEDURE Diffusion couples and sintered powdered compacts were employed in the phase-diagram investigation. The diffusion couples were made by filling 1/8-in.-ID monel-sheathed niobium tubes with tin. The monel sheath was employed to facilitate drawing.' The tubes were then drawn to a tin-core diameter of 32 mils. Samples approximately 3 in. long were then cut from the drawn composite. The tin was drilled out of the ends to a depth of 1/4 in. and niobium-wire plugs were inserted into the ends and peened over. The monel was removed by etching in concentrated nitric acid, after which the samples were sealed in evacuated quartz bulbs and heat-treated in a resistance-wound tube furnace. The samples were quenched into ice water upon removal from the furnace. The diffusion couple samples were examined metallographically employing a chemical etching solution consisting of 10 ml of saturated chromic acid per g of NaF. In addition, two anodizing solutions were used for phase-identification purposes. The first was the picklesimer7 solution; the second consisted of equal parts by volume of 30 pct H2O2 and concentrated NH4OH to which 1 g of NaF was added per 25 ml of solution. The anodizing conditions for the second solution were 2 v and 100 ma with a tin cathode. The powdered compacts were made by pressing previously mixed powders of 99.9 pct pure Sn and 99.6 pct pure Nb supplied by the United Mineral Co. into cylinders 3/8 in. in diameter by 1/2 in. long. The cylinders were then sealed in quartz tubes and heat-treated in the same manner as the diffusion couples. The samples were examined metallographically and by X-ray diffraction techniques. Since it was desirable to be able to correlate the kinetic data with current-carrying capacity, the type of specimen chosen for this part of the investigation had to be a compromise between the optimum system for studying kinetics and one which was suitable for making current-carrying capacity

Jan 1, 1964

By Y. A. Rocher, J. E. Dorn, L. A. Shepard

Creep activation energies for single aluminum crystals favorably oriented for shear by (010) [101] glide were detemined over the temperature range from 78" to 900°K. Observations of slip bands on the specimen surface were made in conjunction with the investigation. From 78" to 780°K, the activation energies obtained in this imestigation agreed closely with those previously found for creep by (111) [101] slip. Between 78" and 140°K, the activation energy was identified with the Peierls process, while between 260°and 780°K the activation energy was close to that for cross-slip. The coarse wavy slip bands nominally parallel to the (010) plane observed above 260°K were attributed to fine cross-slip. From 800" to 900°K, unusually high apparent activation energies ranging from 28,000 to 54,000 cal per mole were obtained. These apparent activation energies were attributed to re crystallization. AS illustrated in Fig. 1, a recent investigation1 has shown that creep of aluminum single crystals by the (111) [i01] mechanism is controlled by three unique processes, each of which is characterized by a single activation energy which is independent of the applied stress and the creep strain. A comparison of the observed activation energies with theoretically calculated values permits a fairly clear identification of the three operative creep processes. Below 450°K, where the activation energy for creep is 3,400 cal per mole, the deformation is controlled by the Peierls process, the activation energy for creep agreeing well with that calculated by seeger2 for the energy required to nucleate the motion of a dislocation loop against the atomic forces of the lattice. Between 590° and 750°K, the observed activation energy for creep of about 28,000 cal per mole agrees well with the energy necessary to induce cross-slip. Seeger and schoeck3 estimate that the activation energy is about 24,000 cal per mole whereas Friedel4 recently calculated this activation energy to be 28,000 cal per mole. Above 800°K the activation energy of 35,500 cal per mole that was observed for creep agrees well with that estimated for self-diffusion in aluminum.= In this range the operative rate-controlling slip process has been clearly identified as that arising from the climb of edge dislocations. The objective of this investigation is to ascertain whether a single crystal of aluminum favorably oriented for simple shear in the [loll direction on the (010) plane might exhibit uniquely different activation energies for creep from those obtained previously for (111) [101] slip. Whereas the exis- tence of such unique activation energies would constitute incontrover table evidence for new mechanisms of slip, the absence of any new activation energies might suggest that slip of aluminum is confined to the (111) [loll mechanism. Several factors prompted the selection of the (010) [101] orientation for study. First, there are more reported observations of (010) [loll slip than of any other nonoctahe-dral mechanism.8-10Secondly, Chalmers and Martius1l have concluded from considerations of the energies of dislocations that (010) slip is the second most favored mechanism in face-centered-cubic metals. Finally, favorable orientations for simple shear by the (010) [loll mechanism provide the least favored orientation for slip by the (111) [101] mechanism. EXPE-RIMENTAL PRO-CEDURE The high-purity aluminum stock, specimen preparation, shear fixture, extensometry, and experimental technique used in this investigation were the same as those previously reported.' Single-crystal spheres grown from the melt of 99.995 pct pure Al* were _ *The high-purity aluminum used in this investigation was graciously given by the Aluminum Company of America. oriented, carefully machined into dumbbell-shaped shear specimens, annealed, and chemically polished. The finished specimen had a central reduced section 0.190 in. wide and 0.590 in. in diam and 1/4-in. grip sections at both sides, 0.690 in. in diameter. The specimen was oriented in the stainless steel grips of the shear fixture with the (010) plane perpendicular to the dumbbell axis and the [loll direction parallel to the stress axis within 2 deg. Creep activation energies were calculated in the previously described manner1 from determinations of the instantaneous change in shear strain rate produced by an abrupt 15 to 20 deg increase or decrease in test temperature. If is the instantaneous strain rate at strain y and temperature T1, and ?2 the instantaneous rate at y and T2,

Jan 1, 1960

By E. Miller, J. M. Eldridge, K. L. Komarek

The liquidus curve of the Mg-Pb system was accurately redetermined. The compound Mg2Pb decomposes peritectically at 538.2° ± 0.3°C to liquid and to a compound p' which melts congruently at 35.0 at. pct Pb and 549.0° ± 0.3°C. The solidus curve of ß' was determined. X-ray diffraction studies indicate that 4' has an orthorhombic structure. Activity values of magnesium calculated from the phase diagram agree with those published in the literature. EXPERIMENTAL thermodynamic properties of binary metallic systems have to be consistent with values calculated from the phase diagram. In systems forming intermetallic compounds the shape of the liquidus curve near a compound is determined by the thermodynamic properties of the coexisting solid and liquid phases. Hauffe and Wagner' neglected the temperature dependence of the chemical potentials and obtained the potential differences of the components of the liquid alloys, relative to stoichiometric liquid. Their calculations were based on the liquidus curve and on the heat of fusion of the compound, and were only valid near the congruent melting point. Steiner, Miller, and Komarek2 developed equations which account for the temperature dependence and obtained the chemical potentials of liquid Mg-Sn alloys over the entire phase diagram from the liquidus and solidus curves and from enthalpy values with the pure components as the standard states. The Mg-Pb phase diagram has been studied by several investigators whose results have been compiled and critically evaluated by Hansen.3 Although the liquidus curve was poorly defined, the general features of the diagram, i.e., one congruent melting compound, Mg2Pb, of essentially stoichiometric composition, two eutectics, and limited terminal solid solubilities, seemed to be suitable for a similar thermodynamic analysis. A careful redeter-mination of the liquidus by thermal analysis revealed, however, the existence of another compound. The liquidus curve between the two eutectics was precisely delineated and the structure and solidus curve of the new compound were investigated. The revised phase diagram was thermodynamic ally analyzed to evaluate the activity of magnesium in the liquid alloys. EXPERIMENTAL PROCEDURE The magnesium metal (Dominion Magnesium Ltd., Toronto, Canada) had a purity of 99.99+ pct; lead (American Smelting and Refining Co.) contained 99.999 pct Pb. Most experiments were carried out in graphite crucibles. Several experiments were made in high-purity alumina (Triangle R.R., Mor-ganite, Inc.) and in Armco iron crucibles to test the inertness of the graphite crucibles. Chemical analysis of magnesium and detailed description of the procedure for thermal analysis have been given previously. For the determination of the solidus curve of the compounds, specimens of initial composition Mg2Pb were equilibrated in a closed isothermal system with magnesium vapor. The source of the magnesium vapor was an alloy which had a gross composition lying in the 0' + L field at the temperature of equilibration. As equilibrium was approached, the specimens lost magnesium to the two-phase reservoir thereby lowering the activity of magnesium in the specimens until activity and composition equaled that of the ß'/ß' + L boundary. Crucibles (1.9 cm ID by 2.2 cm OD by 4.1 cm high) and tightly fitting lids were machined from a molybdenum rod; small, shallow trays were fashioned from thin (0.005 in.) molybdenum sheet, and all the molybdenum components were degreased in hot carbon tetrachloride and then dried. The pieces were then degassed in vacuum at 950°C for about 6 hr. The two-phase alloy was placed at the bottom of the crucible and small specimens of the Mg2Pb compound, weighed on an analytical balance, were placed in two molybdenum trays above the two-phase alloy. The crucible was closed by forcing its lid on and then inserted in a titanium crucible. This crucible was evacuated, flushed twice with argon, and welded under argon. The specimens were equilibrated for about 1 week in a resistance furnace regulated by a Celectray controller, and the runs were terminated by water quenching. The specimens were again weighed and the equilibrium compositions were calculated on the basis that the weight losses were solely due to a loss of magnesium to the two-phase alloy. The structure of the B' phase was investigated by the Debye-Scherrer X-ray diffraction technique. Selected ingots from thermal-analysis experiments containing about 35 at. pct Pb were re-melted, slowly cooled, and crushed in an argon-filled glovebox until the entire ingot passed through a 50-mesh sieve. The powder was thoroughly

Jan 1, 1965

By R. L. Ash, R. R. Rollins, C. J. Konya

An investigation was performed to determine conditions for optimizing the spacing of simultaneously initiated multiple explosive columns. This was done by using models of mortar, dolomite, and Plexiglas with 10-grain mild detonating fuse as the explosive charge. It was desired to simulate blastholes with multiple primers initiated by detonating fuse or when high-velocity explosives are used in low-velocity materials. It was found that optimum spacing between multiple charges was strongly influenced by charge length. At less than optimum charge length, the spacing at which complete shearing was possible between adjacent charges decreased exponentially with a subsequent loss of broken material volume. For charges fired simultaneously, larger burdens and spacings were possible as compared to those necessary for single-crater charges. For each material studied, there was a characteristic optimum charge length and a maximum attainable spacing at any given burden. Proper selection of the spacing distance between charges is fundamental to successful blasting. Its value directly affects the cost of drilling and explosives used per unit of broken material. In addition, the choice of a spacing that is Compatible with a given set of blasting conditions aids in the control of fragmentation sizing, ground vibrations, overbreak, and throw which in turn, influence other production costs. For example, normally loaded blastholes that are spaced too closely invariably promote overbreak and usually give coarse fragmentation. Unless care is taken, airblast and violent flyrock will occur and under certain conditions cutoffs and misfires may result. Too large a spacing, on the other hand, frequently leads to conditions that form bootlegs or toes. The choice of a particular spacicg to use, however, is largely a matter of individual experience and judgment, usually based on trial and error. Very little is known or can be found in the literature with regard to how the spacing between charges is related to field conditions and charge geometry. As a general rule, the firing time sequence of adjacent charges and properties of a material are thought to have the most significant influence on the spacing distance best suited for any given field condition. For example, delayed initiation of adjacent charges usually always requires a closer spacing than when charges are fired at the same time. This should be expected if one considers that the energy normally dissipated and lost in the surrounding ground from charges fired independently would be captured and utilized for breaking material between charges when they are initiated together. Spacing can be extended also when charges are aligned with structural planes of a material, such as jointing, along which shearing is relatively easy. It is customary to relate the spacing (S) between charges to their common burden (B) in the form of a spacing ratio, or SIB. The burden normally is considered as the optimum depth or distance from any single charge perpendicular to the nearest free or open face at which the desired fragmentation and maximum crater yield are obtained. For production blasting, value of the ratio is generally considered to vary from 1 to 2, depending on conditions.1-6 When adjacent charges are fired independent of one another, the value varies from 1 to about 1.4, the closer amount being employed to square corners or produce craters having the ideal 90" apex angle. The larger ratio is the geometric balance value for craters having an apex angle of 135". The basic ideal crater forms in the plane of the charge diameter for charges fired independently are shown in Fig. 1. In the event charges are fired simultaneously, geometric balance in the plane of their charge diameters suggests that a spacing ratio near 2 would be appropriate, as illustrated by Fig. 2. In practice, however, some compromise ratio value must be selected to conform with the specific ground conditions. An example would be where the jointing planes tend to produce 60° or 120° crater angles, the appropriate geometrically balanced charge arrangement being given by Fig. 3. In this condition, the spacing ratio is 1.15, not 1 or 1.4 as suggested for the 90° cratering of independently fired adjacent charges. In view of the foregoing, it would seem logical to assume that whenever charges all having the same burden are fired at the same time, spacing distances always can be greater than those permitted by charges fired independently. In practice this is not the case, however.

Jan 1, 1970

By Charles H. Grant

Some of the limiting factors relative to explosive crushing of rock and ways to overcome a few of these problems are presented. Relationships between borehole diameters, bench heights, and spacings, along with a review of the influence geometry has on energy as these are changed, are discussed. Efficiency in use of explosives and the decay of energy as it moves through rock and is absorbed and dissipated, is described, along with fragmentation as a function of spacings and energy zoning, etc. Communications are one of the major problems encountered. In an effort to provide a better understanding of the use of explosives, it is necessary to take a little different view of what explosives are, how to look at them as tools to fragment rock, and some of the problems encountered in doing so. First, take the explosive: although there are many factors involved, consider these as being reduced to only two — shock-strain imparted to the rock by the high early development of energy, and the gas effect which is a combination of heat, moles of gas formed, rate of formation of these gases which develop pressures, etc. First, consider shock energy by itself and assume there is no gas effect in the reaction. Fig. 1 illustrates a block or cube of rock, in the center of which is detonated an explosive charge which is 100% shock energy. Tensile slabbing would be seen on the surface and probably the cube of rock would generally hang together even though microcracks were formed. If the situation is reversed and an explosive whch has no shock energy and only gas effect (Fig. 2) is considered, the cube of rock would act as a pressure vessel and contain the pressure from the gas effect until it exceeded the rock-vessel strength; then the rock would break in a few large pieces. If these two kinds of energy are put together and the area of shock-strain around the explosive (Fig. 3) is considered, the two energies will be seen working together to furnish broken rock. The gas effect applies pressure to the microcracks formed from the shock energy to weaken the rock-pressure vessel and propagate these cracks to break the rock apart. It not only will be broken more finely, but will break apart at a lower pressure than the gaseffect case, since the shock energy has first weakened the rock vessel. Although tensile spalling from the shock-strain imparts momentum to the rock, the main source of displacement comes from the gas effect. The term "rock" is being used to mean any material to be blasted. These energies are absorbed by the rock in different ways. First, classify rock into two main categories: "elastic" and "plastic-acting." Elastic rock should be thought of as rock which can transmit a shock wave and is high in compressive strength, such as granite or quartzite. Since this elastic rock transmits a shock wave well, it makes good use of the shock energy from the explosive-forming cracks, etc., for the gas effect to work on. Plastic-acting rocks are rock masses which are relatively low in compressive strength and absorb shock energy at a much faster rate, thereby making poor use of the shock energy by not developing as extensive a cracked zone for the gas effect to work on. Rocks of this type are generally softer materials such as some limestones, sandstones, and porphyries. For the most part, the shockenergy part of the explosive reaction is wasted in plastic-acting rock, leaving most of the work to the gas effect. Since the ratio of gas effect to shock energy is different in different explosives, it is easy to understand why some explosives perform well in elastic rock and poorly in plastic-acting rock, and vice versa. Some of the most difficult blasting situations arise when mixtures of plastic-acting and elastic rock are encountered (Fig. 4). Fig. 4 shows an example of granite boulders cemented together with something like a decomposed quartz monzonite which is plastic-acting. The elastic granite boulders will transmit the shock-strain within itself, but when this shock tries to move through the monzonite to the next boulder, its intensity is absorbed by the monzonite and little shock-strain is placed on the adjoining boulder. In addition to this loss by absorbtion, shock reflection at the surface of the boulder will effect tensile spalling. The net effect is poor breakage of the boulders which do not have drillholes in them as they simply will be popped out with the muck. The same is true (Fig. 5) when layers and joints make

Jan 1, 1970

By S. R. Seagle, R. Schuhmann, N. A. Parlee

Rates of CO evolution and CO absorption were measured for liquid-iron alloys containing from 0.15 to 4.4 pet C, using a modified Sieverts apparatus. The alloys were held in alumina crucibles, so that both crucible-metal and gas-metal reactions occurred simultaneously. The data are interpreted on the hypothesis that the C-O reaction rate is controlled by oxygen diffusion across the boundary layers at the gas-metal and crucible-metal surfaces. CARBON-monoxide evolution from liquid Fe-C-0 alloys is a key reaction in steelmaking processes, both in the steelmaking furnace and in the ingot mold. Also, this reaction probably is an essential step in the desulfurization of iron under blast-furnace conditions, for which carbon is the principal reducing agent or deoxidizer present in the iron. The thermodynamic properties of the solutions of C and 0 in Fe and the equilibria of these solutions with gaseous CO and CO, have been investigated in some detail and are reasonably well understood. However, while these equilibrium data have accumulated in the laboratory, other data have accumulated to show that equilibrium conditions are often not achieved under plant operating conditions. Thus, an understanding of the rate and mechanism of the carbon-oxygen reaction not only has theoretical interest, but ultimately may assume considerable practical importance. Reaction Mechanisms and Rate Equations The over-all reaction under consideration is C (in liquid Fe) + 0 (in liquid Fe) + CO (gas) [I] Although this appears to be a straightforward heterogeneous reaction, involving just two phases (liquid metal and gas), a wide variety of reaction mechanisms has been proposed. Different authors have expressed divergent views as to the nature of the rate-controlling steps. Accordingly, a brief discussion of the possible reaction steps will be given. 1) Homogeneous Reaction Within Liquid Iron— Feild' and Jette2 have applied chemical reaction-rate theory as if the C-0 reaction were a homogeneous, second-order reaction occurring within the liquid-metal phase. This mechanism necessarily produces CO "molecules" dissolved in liquid iron. In order to yield CO gas, the dissolved CO must nucleate CO gas bubbles, and the dissolved CO also must be transported to the bubble surface where it enters the gas phase. However, this mechanism is inconsistent with present concepts of the structure of liquid-metal solutions. 2) Reaction at Gas-Metal Interface—At the surface of contact between gas and metal, a second-order reaction between dissolved C and 0 may occur to produce a CO molecule which enters the gas phase. This reaction presumably involves an "activated complex" structure in the interface. Some such surface mechanism appears essential because the carbon and oxygen do not occur in the same species in the two phases and therefore must react in some way at the surface. At liquid-iron temperatures, calculations based on reaction-rate theory indicate that the surface reaction must proceed extremely rapidly and thus cannot be a rate-determining step in the C-O reaction under ordinary conditions.3 ccordingly when the over-all reaction is measurably slow, the surface reaction may be considered at equilibrium; that is Cc* Co* = m'Pco [2] in which Co* and Co* are the concentrations of carbon and oxygen (weight per unit volume of metal), respectively, in the liquid-metal phase at the interface, Poo* is partial pressure of CO in the gas phase at the surface (atmospheres), and m' is the mass-action constant for Reaction [I]. Since activities of C and 0 are not strictly proportional to concentrations, m' is not a true thermodynamic equilibrium constant but varies with carbon content.' 3) Mass Transport—If the reaction is to proceed at the gas-metal surface, as described above, dissolved carbon and oxygen in the liquid iron must come to the surface and gaseous carbon monoxide must move away from the surface. Details of the mass-transport mechanisms depend on the kind of system under consideration. For a stirred metal bath in contact with a CO gas phase, the transport of C and 0 to the surface can be described by 1 Do Mols CO evolved per sec = — 1/16— Do/do Agm (CO-Co*) = 1/12 Do/do Agm *(Co-Co*) [3] In this equation, Do and D, are the diffusion constants for oxygen and carbon in liquid iron, 6, and 6,. are the boundary-layer thicknesses, A,,,, is the

Jan 1, 1959

By K. Tanaka, J. D. Mote, J. E. Dorn

It) an attempt to ulLcoce.lP the operative strain-rate-contl-olliy: dislocation nieclzanistns, specially oviented sizgle clystals of the intel-nzediate 1zexagonal phase containing Ag plus 33 at. pct A1 were tested in tension over a wide range of temperatures. Slip was observed to take place by the {0001} <1120> {l100} mechani fracture took place across the(i100) plane and winning occurred by the (i01Z) ?lechanisn. Basal slip exhibited a strong yield point over the -alzge from 77 to 450°K, the upper ,esolved shear st]-ess having the exceptionally high value of 10,500 psi over this entire ?-a?zge of tenzpei,atuves. The critical 9-esolved shear stress for prismatic slip decreased f7-om 48,000 psi at 4.3"K to 23,000 psi at 170°K (Region 1) follozcirg zt:lzich it decl>eased sloz&ly to 21,500 psi at 475°K (Res&apos;on II); from 475" to 575°K (Regioz III), the c7-itical esolced shear stress dec&apos;-eased precipitously to 2000 psi; and from 575" to 750°K (Region IV) it decreased less afi&apos;dly to a low value of about 500 psi. Pvistintic slip in Region I was pobably controlled by the tliel-nally activated riecharzisui of nucleation and g,-ozcth of kinks in dislocations lying in Peierls potential troughs. In Region II for prismatic slip the critical 1-esolved shear stress was slzocn to be deteemined by sh0l.t-range 01-dering, Overall the forgiorz fo basal slip, 7.c.lre1-e a Strong yield-point phenorlienu ia7as observed, the critical vesolved slzea?-stress was shoztn to be determined by n conibirzation 0-f Szizuki locking and short-range-order Izavderzizg, The precipitous decrease in the critical resolved shear stress with increase in ter,/pe7-atrir-e over Region HI was tentatively ascribed to a decrease in the degree of slort-)ange 07-del;iqq (0)- clusteing) and also the effect of fluctuations the degree of o?der, It is at pgreser2t zrtzce)taitz as to 1t1hethe1- these or other possi1)le effects are also ,esponsible. fo- the data obsel-ved 172 Region IV. 1NTEREST in inter metallic compounds stems not only from their role in dispersion hardening of polyphase alloy ystems but equally from their potentialities for high strength, hardness, and stability not only at atmospheric temperatures but especially at elevated temperatures. As summarized in a re- cent symposium of the Electrochemical Society on "Mechanical Properties of Inter metallic Compound", most of the experimental evidence regarding the mechanical behavior of intermetallic compounds centers about the effect of temperature on the hardness and ductility of polycrystalline specimens. The available data reveal that the plastic behavior of intermetallic compounds might be rationalized in terms of the usual dislocation mechanisms appropriate to a solid solutions providing the additional complexities arising from crystal structure, long-range ordering, short-range ordering, and defect lattices are taken into consideration. It is apparent, however, in terms of the history on a solid solutions, that a complete detailed mechanistic rationalization of dislocation processes may not be possible until the deformation processes are studied in single crystals of intermetallic compounds. The present paper contains a preliminary report on the plastic behavior of single crystals of the hexagonal Ag-A1 intermetallic phase over a wide range of temperatures. The results confirm the thesis that single crystal data provide a most effective method of identifying operative dislocation mechanisms in intermetallic compounds. EXPERIMENTAL TECHNIQUES Several factors prompted the selection of the hexagonal Ag-A1 intermetallic phase for this preliminar investigation on the plastic properties of single crystals of intermetallic compounds: 1) This phase has a wide solubility range5 which would permit future investigations on the effect of composition and axial ratios on slip mechanisms. 2) Although it undoubtedly exhibits short-range ordering (or clustering) this intermetallic phase is free from complexities arising from long-range ordering.6 3) Since the atomic radii of aluminum and silver are practically identical, the possible complications due to Cottrell locking are minimized. 4)Whereas the dislocations on the basal planes are expected to dissociate into Shockley partials and are thus susceptible to Suzuki locking, those on the prismatic planes probably remain complete. 5) The axial ratio, being 1.61, is almost ideal, suggesting that short-range ordering may be almost spherically symmetrical. The present investigation was conducted exclusively with the hexagonal Ag-A1 alloy containing 33 at. pct Al. Preliminary investigations revealed that this alloy undergoes basal slip by the (0001)

Jan 1, 1962

By E. Miller, J. M. Eldridge, K. L. Komarek

The thermodynamic properties of liquid Mg-Ge alloys have been determined between 1000°and 1500°K by an isopiestic method. Germanium specimens, heated in a temperature gradient and contained in covered graphite crucibles of special geometry, were equilibrated with magrtesium vapor in closed titanium tubes. The crucible design allowed free access of magnesium vapor to the samples during the equilibration to form alloys of magnesium and germanium, but prevented magnesium losses from the crucibles on quenching the titaniuin tubes to terminate the experimental runs, thus preserving the equilibrium alloy compositions. The activities and partial molar enthalpies of magnesium and the integral thermodynamic properties of the system were calculated from the experimental data. THE Mg-Ge phase diagram&apos; shows one congruent melting compound, Mg2Ge, of essentially stoichio-metric composition, two eutectics, and very limited terminal solid solubilities. Very little information is available on the thermodynamic properties of the Mg-Ge system. The free energy of formation of Mg,Ge was recently deter-mined2 by a Knudsen cell technique in the temperature range 610° to 760°C. The standard enthalpy of formation of Mg,Ge was measured calorimetrically by Bever and coworkers.3 The present study was undertaken as part of a general investigation of the thermodynamic properties of the homologous series of Mg-Group IVB systems, i.e., Mg-Pb,4 Mg-Sn,5 Mg-Ge, and Mg-Si. An isopiestic technique was used which was developed by the authors5 for investigating the thermodynamic properties of liquid Mg-Sn alloys. Specimens of the nonvolatile component, contained in covered graphite crucibles, are heated in a temperature gradient in an evacuated and sealed titanium reaction tube, and equilibrated with magnesium vapor of known pressure. The method employs crucibles of special geometry which preserve the high-temperature equilibrium composition of liquid alloys having a highly volatile component such as magnesium on termination of the experimental runs by quenching the crucibles to room temperature. EXPERIMENTAL PROCEDURE First reduction germanium of 99.999+ pct purity (Eagle-Pitcher Co., Cincinnati, Ohio) and 99.99+ pct magnesium metal (Dominion Magnesium Ltd., Toronto, Canada) were used. The graphite crucibles were machined from high-density (1.92 g per cu cm) graphite rods (Basic Carbon Corp., Sanborn, N.Y.) which had a maximum ash content of less than 0.04 pct. The non-reactivity of graphite with germanium at the temperatures used in this study had been previously established by Scace and Sleck.6 The experimental procedure has been previously described in detail.5 The selection of a particular crucible geometry for a run was determined by a combination of imposed experimental conditions, the principle being that more tightly covered crucibles were required to preserve alloy compositions during quenching when higher magnesium pressures and higher specimen temperatures were used. Depending upon the composition range of the equilibrated alloys the source of the magnesium vapor was either pure magnesium or a two-phase mixture of Mg2Ge + Ge-rich liquid of known magnesium pressure. The experimental runs can be divided into the following three groups on the basis of crucible geometry and magnesium source material. Crucibles with Small Holes and Pure Magnesium Reservoirs. The crucible dimensions were identical to those of the Mg-Sn investigation5 except that the hole diameters were reduced to 0.010 in. because of the higher temperatures and higher magnesium pressures involved in the Mg-Ge system. During an equilibration run, magnesium vapor diffused from the reservoir to each specimen through the small holes, one drilled through the crucible lid and two others drilled through graphite baffles positioned vertically inside the crucible between the lid hole and the specimen. Since the magnesium pressure was high, i.e., in the range 117 to 277 Torr, during the equilibration time of approximately 24 hr, equilibration was not impeded by these holes. A specimen composition at equilibrium was fixed by the relative temperatures of the specimen and the reservoir, and by the thermodynamic properties of the system. Upon brine quenching the titanium reaction tube to end a run the vapor pressure of magnesium above the liquid alloys decreased exponentially with decreasing temperature, and the small cross-sectional areas of the holes (4.9 x 10"* sq cm) drastically reduced magnesium losses from the crucibles. Because of its low vapor pressure, germanium losses from crucibles during a run were at most 0.2 mg for pure germanium and correspondingly less for the alloys. This crucible geometry satisfactorily retained the equilibrium alloy compositions on quenching for magnesium-rich (from 3 to 33 at. pct Ge) alloys provided their temperatures were below the melting

Jan 1, 1967

By C. Atkinson

A method is described for the numerical solution of the diffusion equation with a composition-dependent diffusion coefficient and applied to the radial growth of a cylinder; the radial growth of a sphere, and the symmetric growth of an ellipsoid. Sample applications of the method are made to the growth of particles of proeutectoid ferrite into austenite. RECENTLY&apos; we described a method for numerical solution of the diffusion equation with a composition-dependent diffusion coefficient for the case of the growth of a planar interface. In this paper we extend this method to describe the radial growth of a cylinder, the radial growth of a sphere, and the symmetric growth of an ellipsoid. In the latter case, limiting values of the axial ratios of the ellipsoid reduces the problem to one of a cylinder, a sphere, or a plane depending on the axial ratio. A check on these limiting values is made in the results section. In all of these cases we consider growth from zero size. A natural consequence of this assumption as applied to the sphere, for example, is that the radius of the sphere is proportional to the square root of the time. This is consistent with the condition that the radius is zero initially, i.e., grows from zero size. It may be argued that it is more realistic to consider particles which grow from a nucleus of finite initial size; even in this case the analysis of this paper is likely to be applicable. This can be seen if a comparison is made of the work of Cable and Evans,2 who consider a sphere of initially finite size growing by diffusion in a matrix with a constant diffusion coefficient, with the results of Scriven3 for growth from zero size. This comparison shows that the rates of growth in each case differ trivially by the time the particle has grown to about five times its initial size." This investigation is a generalization of those of Zener,4 Ham,5 and Horvay and cahn6 to the situation often encountered experimentally, in which the diffusion coefficient varies with concentration. First let us consider each of the cases separately. I) GROWTH OF SPHERICAL PARTICLES FROM ZERO SIZE In this case the differential equation in the matrix depends only on R, the radius in spherical coordinates, and can be written: ? 1 <^\ ^13D . , dt U\dRz + R 3Rj + dR dR [ J where C is the composition, t is the time, and D is the diffusion coefficient which depends on c. The boundary conditions will be: c = c, at the moving interface in the matrix, c = c, at infinity in the matrix (and at t = 0, everywhere in the matrix), c = X, is the composition in the spherical particle. Each of the above compositions is assumed constant. In addition there is the flu condition at the moving interface which can be written: , dR0 ~/3c dt \dR/H =Ra where R,, which is a function of t, is the position of the moving interface. We make the substitution q = RI~ in [I] reducing this equation to: & - m - *ws) »i where we have written D = D,F(c) or simply D,F, and Do = D(c,). Thus F[c(q0)] = 1 where q, = ~,/a is the value of the dimensionless parameter q evaluated at the interface. Multiplying Eq. [2] by dq/dc and integrating, we find: where the lower limit of the integral has been chosen so that dc/dq — 0 as c — c,, thereby satisfying the boundary condition at infinity. We require, then, to solve Eq. [3] subject to the condition c = c, when q = q, (this follows from putting R = R, at the interface) together with the flux condition which can be rewritten in terms of q as: Eqs. [3] and [4] together with the condition c = c, at q = q0 enable us to find 77, and the concentration profile c = c(q). Numerical Method. We treat Eq. [3] in the same way as we did the corresponding equation for the planar interface problem&apos; i.e., by dividing the interval c, to c, into n equal steps so that: cr = ca -rbc [5] where r takes the values 0, 1, ... n and we call no,, q1, ... nn the values of n corresponding to the compositions c,, c,, ... c,.

Jan 1, 1970