IllbO DISSOLUTION OF GRAIN BOUNDARY ALLOTRIOMORPHS IN THE Al-Cu AND Al-Ag SYSTEMS BY ARISTEDES PASPARAKIS B.Sc, University of Pretoria, 1964 M.Sc. (Tech.), University of Sheffield, 1967 A THESIS SUBMITTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY . in the Department of METALLURGY We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA February, 1972 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the Head of my Department or by his representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of Metallurgy The University of British Columbia Vancouver 8, Canada Date April 4, 1972 - ii -ABSTRACT The dissolution behaviour of grain boundary allotriomorphs has been studied in the systems Al-Cu and Al-Ag. Using electron probe microanalysis, isoconcentration contours around dissolving allotriomorphs in an effectively infinite matrix have been determined. The influences of allotriomorph shape, grain boundary misorientation, solute super-saturation in the matrix, volume and grain boundary diffusion, and homologous temperature on the shape of these contours have been examined. It was established that the shape of the isoconcentration contours is dependent only on homologous temperature (TR) for both systems studied. Above T„ = 0.92, no grain boundary diffusion contribution hi to the dissolution process is observed. The grain boundary diffusion contribution was- found to increase with decreasing T and to completely H dominate the dissolution process below T = 0.72. H Using the back scattered electron image on the electron probe microanalyzer, the dissolution rate of indivTatial grain boundary allotriomorphs was determined for various values of T , under conditions in which there was impingement of diffusion fields from adjacent precipitates. An exponential relationship between axial length or width and dissolution time was found to adequately describe the observed dissolution kinetics. A change in axial ratio accompanied the dissolution of the grain boundary allotriomorphs. At high T an increase in axial ratio with dissolution time was observed, whereas at low T„, a decrease in axial ratio (i.e., spheroidization) was observed, rl A model has been proposed to account for this behaviour. - iii -TABLE OF CONTENTS Page 1. INTRODUCTION 1 1.1 Grain Boundary Precipitates 1 1.2 Growth of Precipitates 4 1.2.1 Thickening of Grain Boundary Allotriomorphs. 5 1.2.2 Lengthening of Grain Boundary Allotriomorphs 5 1.2.3 The "Grain Boundary Collector Plate" Model . 6 1.3 Dissolution Studies of Precipitates 7 1.3.1 Observations of the Dissolution of Precipitates by Electron Microscopy 12 1.3.2 Dissolution Studies Using the Electron Probe Mlcroanalyzer 1.4 Dissolution Mechanisms Derived from Precipitate Growth ig 1.5 Objective of Present Investigation 17 2. EXPERIMENTAL PROCEDURE 20 2.1 Reasons for Systems Chosen 22.2 Preparation of Master Alloys 1 2.3 Procedure Used to Grow Large Grain Boundary Allotriomorphs 24 2.4 Dissolution Heat Treatment Procedure 25 2.5 Diffusion Couple Work 28 2.6 Electron Probe Microanalysis 29 2.6.1 Isoconcentration Contours 22.6.1.1 Resolution of Experimental Technique 32 2.6.2 Kinetic Studies 36 - iv -Page 2.6.3 Diffusion Couple Measurements 38 2.7 Validity of Experimental Technique 3^ 2.7.1 Isoconcentration Contour Measurements 38 2.7.2 Kinetic Measurements 41 3. ISOCONCENTRATION CONTOUR RESULTS AND DISCUSSION 44 3.1 Introduction 43.1.1 Diffusion Couple Results 45 3.1.2 The Microstructure of Equilibrated Alloys .. 49 3.1.3 Dissolution Isoconcentration Contours 50 3.2 Rate Controlling Mechanism 54 3.3 The Effect of Precipitate Shape T.S /R 1 56 3.4 Homologous Temperature [T^] 60 3.4.1 Calculations of Dy and Cj 71 3.5 The Effect of Supersaturation (k) and Volume . Diffusion Coefficient (Dy) 77 3.6 The Effect of Grain Boundary Misorientation (R) ... 81 3.6.1 Evaluation of the Grain Boundary Diffusion Coefficient(D , ) 84 g.b. 3.7 Semi-Quantitative Analysis of Isoconcentration Contours 89 3.8 The Al-15.75 wt.% Ag System 94 3.9 Applicability of Contours 9 - v -Page 4. KINETIC RESULTS AND DISCUSSION 101 4.1 Introduction 104.2 Application of the Existing Models for Dissolution Kinetics 102 4.3 An Exponential Model for Planar Dissolution Kinetics 111 4.4 Precipitate Shape Change During Dissolution 122 4.4.1 High Homologous Temperatures 131 4.4.2 Low Homologous Temperatures 137 4.4.3 Intermediate Homologous Temperatures 138 5. CONCLUSIONS: 142 BIBLIOGRAPHY 4 APPENDIX I Correction Procedures for Quantitative Electron Probe Microanalysis 148 1.1 The Al-Cu System1.2 The Al-Ag System 153 APPENDIX II Electron Spot Size Determination 156 APPENDIX III Calculation of the Flux Line Divergence for Different Sections of a Two Dimensional Elliptical geometry 159 APPENDIX IV Dv Calculation from the Diffusion Couple Results 162 APPENDIX V Calculation of Standard Deviation for the Isoconcentration Contours 164 - vi -Page APPENDIX VI Calculation of ij> 166 APPENDIX VII Calculations of the Impingement Resulting from the Interaction of the Confocal Spheroids About an Allotriomorph and a Surrounding Spherical Field 170 - vii -LIST OF FIGURES Figure Page 1 A schematic diagram of the "Grain Boundary Collector Plate" model showing the three steps that take place in the growth of grain boundary allotriomorphs 8 2 (a) Section of an equilibrium phase diagram showing Tg and T^ for an alloy of composition Cn. (b) Dissolution composition profile for an alloy of composition Cn 10 3 A comparison of composition contours for (a) growth, (b) dissolution 18 4 (a) The aluminum-copper equilibrium phase diagram (b) The aluminum rich end of the aluminum-copper equilibrium phase, diagram 22 5 The aluminum-silver equilibrium phase diagram 23 6 Optical micrograph of a typical alloy microstructure for the Al-2.81 wt.% Cu alloy 26 7 Schematic representation of the directions of movement of a specimen relative to the electron beam, in the electron probe microanalyzer 30 8 Optical micrograph of a typical grain boundary allotriomorph, with contamination marks showing the full extent of scans across the precipitate and grain boundary 31 - viii -Figure Page 9 The (a) specimen current, (b) background concentration, and (c) Cu^Kct^ concentration, for the blank profile (1) Through the centre of the precipitate, (2) at the allotriomorph tip and (3) 5 u away from the tip across the grain boundary 34 10 Schematic diagram depicting the effect of spot size on the evaluation of concentration at a point x^ away from the precipitate interface 35 11 A schematic depiction of the dissolution flux lines for; (a) an allotriomorph sectioned above the minor axis b, (b) an allotriomorph sectioned at the minor axis and (c) an allotriomorph sectioned below the minor axis 35 12 Optical micrograph showing (a) electron probe traces across the bonded interface of an Al-CuA^ diffusion couple, (b) microstructure of large CuA^ grain in an a + 6 matrix 39 13 Dissolution isoconcentration contours of a grain boundary allotriomorph on the surface and 3u below the surface (A and B) and a second allotriomorph in the bulk material (C). At T„ = 0.97 for alloy A ... 40 H 14 Population distribution curve for allotriomorphs of varying axial ratio SQ/Ro for alloy A 42 15 Composition contour for Al-CuA^ diffusion couple at 500°C 46 - ix -Figure Page 16 Probability plot of atomic % Cu versus diffusion distance (x) over the solid solubility range of the a phase for the 500°C diffusion couple 47 17 Arrhenius plot of versus 1/T°K ''"for the diffusion couple results 48 18 Size distribution (So) of precipitates as a function of S /R for alloy A 51 o o J 19 Isoconcentration contours at T = 0.91 for alloy B.. 52 rl 20 Concentration contour along the minor axis of a grain boundary allotriomorph at 545°C for alloy A .. 55 21 Isoconcentration contours for three allotriomorphs having axial ratios SQ/Ro of 1) 18, 2) 3.5, 3) 1.7 for alloy A at T • = 0.93 57 22 Isoconcentration contours of two Widmanstatten plates for alloy A at T„ = 0.93 59 rl 23 Isoconcentration contours for alloy A at T = 0.97.. 62 H 24 Isoconcentration contours for alloy B at T = 0.94.. 63 H 25 Isoconcentration contours for alloy Bat T = 0.92.. 64 ri 26 Isoconcentration contours for alloy B at T = 0.905. 65 ti 27 Isoconcentration contours for alloy B at T = 0.86.. 66 H 28 Isoconcentration contours for alloy C at T = 0.80.. 67 ti 29 Isoconcentration contours for alloy C at T = 0.77.. 68 30 Isoconcentration contour for alloy D at T = 0.72... 70 rl 31 Interface concentration C^. versus vTTt" for alloy B at TR = 0.92 73 - X -Figure Page 32 C versus erf (—— ) for alloy B at T„ = 0.92 74 2vfit~ H C(x t) ~ CM 33 Probability plot of [—~2— ] versus diffusion 4 ~ LM distance (x) for alloy B at T„ = 0.92 75 n 34 Isoconcentration contours for alloy A at T = 0.92.. 79 H 35 Isoconcentration contours of two allotriomorph shaped precipitates situated on a low angle grain bounary for the Al-Ag alloy at T„ = 0.84 82 36 Isoconcentration contours for the Al-Ag alloy at TR = 0.84 83 37 Schematic diagram for Shewmon surface diffusion model (after P.G. Shewmon) 85 38 Evaluation of D , for alloy B at T, = 0.86 86 g.b. H m, " ^. flux via grain boundary diffusion ,,s 39 The function — :— TTF^—. W flux via volume diffusion versus temperature T 91 40 A comparison of theoretically derived versus homologous temperature T with the experimental rl , _. _ solute via grain boundary diffusion observation of —;— ; =—: ,. — .. 92 solute via volume diffusion 41 Schematic diagram of method used to calculate 1) solute via grain boundary diffusion, and 2) solute via volume diffusion 94 42 Isoconcentration contours for the Al-Ag alloy at Tu = 0.885 95 rl 43 Isoconcentration contours for the Al-Ag alloy at TR = 0.93 6 - xi -Figure Page 44 Isoconcentration contour for alloy B at T = 0.83... 97 rl 45 Isoconcentration contours for alloy B at T = 0.885. 98 rl 46 Arrhenius plot of Dy versus l/T°K-1for the Al-15.75 wt.% Ag alloy 100 2 47 A plot of R versus t and R versus /t for alloy B at T„ = 0.92 104 ri 48 Log [S/Sq] versus t and log [R/Rq] versus t for alloy A at T = 0.97 106 rl 49 Log [S/Sq] versus t and log [R/Rq] versus t for alloy D at T„ = 0.86 107 rl 50 Log [S/Sq] versus t and log [R/Rq] versus t for alloy B at T„ = 0.86 108 rl 51 Log [S/Sq] versus-t and log [R/Rq] versus t for alloy D at T • = 0.84 109 n 52 Log [S/SO] versus t and log [R/Rq] versus t for alloy D at Tu = 0.76 110 rl •53 (a) Geometrical model of a particle with half-thickness Rq in a matrix with interparticle spacing 111 I. (b) An appropriate phase diagram with growth of precipitates taking place at and dissolution at T^* (c) Dissolution profile prior to impingement. 112 54 (a) The concentration distributions during dissolu tion after impingement takes place, (b) Schematic diagram of model used with a constant concentration C„, at x = 1 113 - xii -Figure Page 55 R versus t for equations 11 and 12 117 56 Log [R/Rq] versus t for dissolution of precipitates in an Al-4 wt.% Cu system at 520°C 119 57 Log [R/RQ] versus t for the dissolution of planar precipitate in a finite system with k = 0.22 120 58 Arrhenius plot of versus 1/T°K "'"for the kinetic results 121 59 Back scattered electron image micrographs,used to observe the dissolution kinetics of an allotriomorph at TR = 0.97 for alloy A 123 60 Back scattered electron image micrographs,used to observe the dissolution kinetics of an allotriomorph at TR = 0.97 for alloy A 124 61 Back scattered electron image micrographs,used to observe the dissolution kinetics of an allotriomorph at T„ = 0.88 for the Al-Ag alloy 125 ri 62 Back scattered electron image micrographs, used to observe the dissolution kinetics of an allotriomorph at TH » 0.88 for alloy B 126 63 Back scattered electron image micrographs,used to observe the dissolution kinetics of an allotriomorph at TD = 0.84 for the Al-Ag alloy 127 n 64 Fraction of precipitates showing an increase in axial ratio greater than 10% (F) versus super-saturation (k) 129 - xiii -Figure Page 65 Fraction of precipitates showing an increase in axial ratio greater than 10% (F) versus T 130 n 66 (a) Symmetrical distribution of precipitates with [S/R ] = 4 131 o o (b) Schematic representation of concentration fields around a symmetrical array of precipitates... 132 (c) Diagram of spherically symmetrical field surrounding a grain boundary allotriomorph 132 67 The resulting impingement between the confocal spheroid.; isoconcentration contours about an allotriomorph and the contours of the surrounding spherical field aftex t = 500 seconds for various [C/C^.] values 134 68 (a) [C/Cj] versus x for the minor axis Ob ^35 (b) [C/C ] versus x for the equatorial plane Oa 136 S ,R 69 Percentage increase in [_, ] versus [S /R ] 139 S /R o o o o 70 (a) Log [R/RQ] versus t (b) [S/R] versus t 140 1-1 Cn versus K_ 152 Cu Cu I- 2 CA versus K. 5 Ag Ag II- l Schematic diagram of the electron spot size 157 III- l Two dimensional sections of an oblate spheroid 159 IV- 1 Composition profile of diffusion couple at 500°C .... 162 VII-1 Schematic diagram comparing actual contour to confocal spheroid 170 -xiv-Figure Page VII-2 Position of error function profile 171 VII-3 Pertaining to Equation 4 172 VII-4 Schematic diagram of the summing up of two in dependent profiles 174 - XV -LIST OF TABLES Table Page I Morphology Classification and Relationship to Boundary Misorientation 3 II (a) Precipitate Growth Equations 13 (b) Precipitate Dissolution Equations 14 III Average Sq for Allotriomorphs with SQ/R0 = 3.5 and Average Width of Precipitate Free Zone for the Different Alloys 50 IV Contour Experiments 61 V Interface Concentration (C^.) and Volume Diffusion Coefficient (Dy) for Contour Experiments 75 VI Comparison of the Effects of k, T and T on the rl Shape of Isoconcentration Contours 80 VII Grain Boundary Diffusion Coefficients (D^ ^ ) Determined from Contours 88 VIII Kinetic Experiments 103 II- l E , A, Z and p Values for Ag Lot, , Al Ka and Cu Ka^ Radiation 156 III- l A9° for Different Values of x/b 160 IV- 1 Dv Values of the Diffusion Couple Results 163 V- l Diffusion Distances, x, a and — 164 2 2 ^ VII-1 Values of (r^ - B ) 2 and n^ of Confocal Spheroids for Different C/Cj. Ratios 172 VII-2 Values of r for Concentric Spheres for Different C/CT Ratios 173 - xvi -ACKNOWLEDGEMENTS The author would like to thank Dr. L.C. Brown for invaluable help and encouragement throughout the duration of this project. Thanks are also due to Dr. D.E. Coates for numerous helpful discussions and guidance during the preparation of this thesis. The criticisms and suggestions of Dr. E.B. Hawbolt, Dr. F. Weinb and R.V. Krishnan of the final draft of the thesis are sincerely appreciated. The assistance of the technical staff and in particular that of Mr. A. Lacis is greatly appreciated. Thanks are extended to Mr. A. Bakas for his help with the figures and drawings for the thesis. Financial assitance provided by the National Research Council under grant number A 2459, and the graduate fellowship awarded by the University of British Columbia are gratefully appreciated. 1. INTRODUCTION The growth of precipitates in solid-state systems has been examined extensively and the results have been analysed theoretically on the basis of diffusion theory. Reasonably good correlation between experiment and theory has been obtained for grain boundary allotrio morph growth, by incorporating volume and grain boundary diffusion in the analysis. Much less experimental effort has been expended on the dissolution of precipitates and, as yet, no conclusive experimental observations have been made on the dissolution of grain boundary allotriomorphs. It has been suggested that the same analysis as used for allotriomorph growth is applicable to allotriomorph dissolution. For a number of reasons, which will be discussed later, this analysis is considered to be at best only a first approximation when applied to dissolution. The theory pertinent to the growth and the dissolution of precipitates will follow a short discussion of the morphology of grain boundary precipitates. 1.1 Grain Boundary Precipitates The morphology of grain boundary precipitates as a function of misorientation across the boundary on which they are situated was - 2 -1 originally studied and classified by Dube et al. This classification 2 3 was amended slightly by Aaronson. In subsequent work Aaronson has shown that grain boundary misorientation can be characterized by the parameter R where: and X, Y, and Z represent the minimum rotations about the three ortho gonal axes required to bring the two adjacent lattices into coincidence. This classification of precipitate morphology as a function of R is given in Table I. The fact that precipitate shape is a function of L 4 only R has been verified by the work of Toney and Aaronson in the Fe-Si system, Clark"* in the Al-18 wt.% Ag system and Hawbolt and Brown6 in the Ag-5.64 wt.% Al system. Grain boundary allotriomorphs, the most common type of grain boundary precipitates, are confined to grain boundary misorientations of R > 15°. These grain boundaries are considered to be high energy disordered boundaries and comprise a large proportion of the boundaries in a random polycrystal.^ The high energy of the boundaries is 'instrumental in: (i) grain boundary allotriomorphs being the first precipitates.to nucleate and hence appear during growth,and (ii) alio-triomorphs forming in preference to Widmanstatten plates at low degrees of supersaturation. In general9allotriomorphs have an oblate spheroidal shape, being circular in the plane of the grain boundary and having a lens shaped R Morphology Classification and Relationship to Boundary Misorientation Grain Boundary Allotriomorphs Primary Sideplate Widmanstatten Secondary Sideplates Primary and Secondary Sawteeth Idiomorphs Morphology ill A Grain Boundary Misorientation Associated with large angle or disordered— type grain boundaries i.e., >15° misorientation » Associated with small angle or dislocation type boundaries and grain interiors <15°. Generally associated with intermediate Disorientations 15"- 25°. Generally associated with intermediate misorientations 15*- 25°. Usually an intragranular morphology although were observed for low angle grain boundaries in the Pe-l.55 Si system. (Clark , Dube and Aaronson ) - 4 -(elliptical) cross section. In the Al-18 wt.% Ag and the Al-4 wt.% g Cu systems, the allotriomorph-matrix interface is generally found to be incoherent and the allotriomorphs develop equally in both 4 8 grains. It has been noted, ' however, that nucleation often takes place in only one grain. This results in a planar, semi-coherent interface with a definite crystallographic relationship between the allotriomorph and the matrix grain. Initial growth of the allotrio morph takes place in the adjacent matrix grain. At a later stage of growth the semi-coherent interface becomes mobile followed by growth into both matrix grains. In an examination of grain boundary ferrite precipitation from 9 austenite in a Co-20% Fe alloy, Ryder and Pitsch and P.L. Ryder et al.^ found that grain boundary precipitates on high angle grain boundaries always have an orientation relationship with at least one and sometimes both of the austenite grains. However, they did not make sufficient measurements to determine whether the austenite-austenite relationships are completely random. Thus some doubt is cast on the observations of an orientation relationship between the precipitate and both grains. 1.2 Growth of Precipitates Theoretical growth rate calculations for the diffusion-controlled growth of planar, cylindrical and spherical precipitates were originally 11 12 made by Zener and Frank. These solutions were extended to cover 13 14 the growth of spheroids and ellipsoids by Ham. ' The rate of growth of the major and minor axes of oblate spheroids, for various • - 5 -axial ratios, was numerically evaluated by Horvay and Cahn.^ These solutions do not include impingement effects and are therefore only valid during the initial period of growth prior to the interaction of adjacent diffusion fields. The experimental observations and analysis of the growth and dissolution of allotriomorphs reported in the literature is summarized below. First, growth is considered in terms of thickening and lengthening kinetics. 1.2.1 Thickening of Grain Boundary Allotriomorphs The first measurements of growth kinetics of grain boundary 16 allotriomorphs were made in commercial steels by Mazanec and Cadek 17 and Hickley and Woodhead. Both reported a parabolic thickening rate for the growth of ferrite precipitates in bulk steel samples, which is consistent with diffusion controlled growth. Hawbolt and Brown, using a statistical approach, found that at high homologous temperatures (T > 0.90) the thickness of g phase rl grain boundary allotriomorphs in the Ag-5.64 wt.% Al system increases 18 parabolically with time. Aaron and Aaronson examined the thickening of eCCuA^) grain boundary allotriomorphs in an Al-4 wt.% Cu alloy for homologous temperatures in the range T = 0.54-0.69. They assumed the true thickness to be that of the thickest allotriomorphs in a series of foils reacted for successively increasing times at each temperature. Their diffusion coefficients, calculated using Zeners^ planar model, were found to be several orders of magnitude higher than the literature volume diffusion coefficients. They conclude - 6 -that at lower homologous temperatures grain boundary diffusion is operative and contributes to the enhanced growth rate. 1.2.2 Lengthening of Grain Boundary Allotriomorphs Less data on the lengthening rate has been reported than on the thickening rate of grain boundary allotriomorphs. Dube et al.''" and 3 Aaronson measured the largest ferrite precipitate resulting from the Fe-C pro-eutectoid reaction and reported a linear lengthening rate up to impingement. Aaronson's results for the Fe-C system showed reasonable agreement with the lengthening rate equation proposed by 19 6 Hillert. Hawbolt and Brown found that the lengthening of grain boundary allotriomorphs increases parabolically with time in the Ag-5.64 wt.% Al system at high homologous temperatures Tu > 0.90. 18 At low homologous temperatures, Aaron and Aaronson found that the lengthening rate of grain boundary allotriomorphs in the Al-4% Cu 19 system, cannot be adequately accounted for by the Hillert analysis This analysis gave values several orders of magnitude higher than the literature value. 1.2.3 The "Grain Boundary Collector Plate" Model The "grain boundary collector plate" model was formulated to account for the greatly enhanced growth kinetics observed by Aaron 18 and Aaronson. This model makes use of the high diffusivity paths provided by the grain boundary and the disordered allotriomorph-matrix interface. The model assumes D , /DTI -»- °° i.e. the volume diffusion g.b.' V contribution is insignificant^and grain boundary diffusion completely dominates the growth process. Growth of grain boundary allotriomorphs - 7 -according to this model occurs by: 1. volume diffusion of solute to the grain boundary; 2. solute transport along the grain boundary collector plate to the allotriomorph (a finite grain boundary area of half length (R^-R) is associated with each allotriomorph); 3. interphase boundary diffusion and deposition of solute over the surface of the growing allotriomorph as shown in Figures 1(a) and 1(b). 20 18 Brailsford and Aaron modified the Aaron and Aaronson Collector Plate Model to accommodate situations in which the volume diffusion contribution is significant. The "Grain Boundary Collector Plate" model predicts a thickening rate proportional to t^*"* and a lengthening rate proportional to t^"2"' for the experimental 18 conditions used in the work of Aaron and Aaronson. These predicted rates are in reasonable agreement with the latter's experimental 0 34 observations of thickening proportional to t " and lengthening 0.27 21 proportional to t . Goldman et al., using the Brailsford and , 20 Aaron analysis for allotriomorph growth rates in the temperature range T = 0.69-0.78, found that volume diffusion directly to the rl allotriomorph contributes significantly to growth only at the highest homologous temperature (TR = 0.78). Below this temperature the major contribution to growth is via grain boundary diffusion. 1.3 Dissolution Studies of Precipitates Dissolution of precipitates occurs when a two-phase material is heated or cooled to a temperature where the precipitate phase becomes unstable. This results in mass transport from the unstable second phase Figure 1. A schematic diagram of the "Grain Boundary Collector Plate" model showing the three steps that take place in the growth of grain boundary allotriomorphs. - 9 -to the stable matrix phase, as shown in Figure 2(a), where alloy Cn may be present as either (a+6) or a, depending on the temperature. The dissolution composition profile for an alloy of composition Cn, after instantaneous heating from T^ to T^, is depicted in Figure 2(b). The comenclature is: Cn = equilibrium alloy concentration; = concentration in the matrix at the precipitate-matrix interface; C^ = matrix concentration; Cp = composition of precipitate; Tg = precipitate growth temperature; Tp = precipitate dissolution temperature. 22 Nolfi et al. developed analytical expressions for the growth and dissolution of spherical precipitates for diffusion controlled, interface reaction controlled,and mixed control mass transport. In subsequent experimental work on the dissolution of Fe^C in ferrite, 23 24 Nolfi et al. ' found that local equilibrium with respect to the iron is;not maintained at the Fe^C-ferrite interface. They conclude the dissolution of Fe^C-ferrite is controlled entirely by an interfacial re action. 25 26 , Using a numerical method, Tanzili and Heckel ' treated: the diffusion controlled dissolution, from an unstable spherical, cylindrical or planar second<phase to a stable finite matrix phase,as a function of , C^, Dg, Cn and R for complete precipitate solution and subsequent overall homogenisation of the matrix. Because the present investigation is primarily concerned with diffusion controlled dissolution, it is appropriate that this topic be discussed in greater detail. For a spherical precipitate the - 10 -Z o I— < \— z uu U z o UJ ? T, D < ? T, c c c CONCENTRATION C, -'O -M •R(t)-R, DISTANCE <r> Figure 2. (a) Section of an equilibrium phase diagram showing T and T for an alloy of composition C^. (b) Dissolution composition profile for an alloy of composition C 0' -11-diffusion field for diffusion controlled dissolution must satisfy the following field equation: DV2C = £ 1 dt where C = C(r,t) is the concentration of solute in the matrix surrounding the precipitate, r (defined as in Figure 2(b))is the distance from the centre of tie precipitate,and it is assumed D f D(C) . The following boundary conditions apply: C(r = R,t) = C. 0 < t < C(r, t=0) = C. M r > R C(r = oo}t) = C. M 0 < t « oo where R is the precipitate radius. C^. is the concentration in the matrix at the precipitate-matrix interface for t > 0 and is the initial matrix composition. A further boundary condition is the interface flux balance: f C - C ) ^ = D lLP V dt D 9r r=R where it is assumed f (r,t). This assumption is justified, as the change in composition of CuAl^ over a range of 548°C is less than 0.5 wt.% Al (see Figure 4a). - 12 -Several approximations have been used to treat the diffusion controlled dissolution of an isolated precipitate in an infinite matrix. These models include the Reversed Growth model, the Laplace 2 dR model (V C= 0), the Stationary Interface model (-j-^- = 0), and the Linear ized Gradients model. Of these,the Stationary Interface approximation appears to be the most accurate. In a theoretical comparison of the different mathematical models for diffusion limited growth and 27 dissolution Aaron et al. showed that the stationary interface approximation is the only one valid for both growth and dissolution of spherical and planar precipitates (see Table II). 1.3.1 Observations of the Dissolution of Precipitates by Electron Microscopy Using electron microscopy at high homologous temperatures, 2 8 Thomas and Whelan (1960) observed precipitate dissolution in an Al-A wt.% Cu alloy. They found that in the final stages of dissolution three out of twelve precipitates satisfy the following equation: = - kD 4 dt They suggest that equation (4) should adequately describe solution kinetics for small k values where: 2<ci-y s is the supersaturation. Equation (5) is derived by an approximation cr°M using the steady state diffusion theory, for small (7;—pr-) ratios. T> I - 13 -TABLE no PRECIPITATE GROWTH SPHERE R = \.(Dt)"s" PLANE S = X.(Dt)1 j Exact Solution Invariant Size Approx. Invariant Field Approx. Linearized Gradients Approx. X, = 2X ^eX2[e~x2-\ erfc \] = - 7 4 lim X1 - Xz k -> 0 lim \1 = x3 k - 0 k < 0: Xx > X2 > X3 lim Xx = XA', k ^ -2j Xx > XA k - 0 x2 = - k lim X2 - X3 k - 0 /4E f(k) = -4k L 1 + 4k Y l_r„X(k). f(k) X4 - /-k/r X1 = 2X X2 k /TT X e erfc X = - -lim Xx = X2 k - 0 k < 0: \x > X2 |k| > 1.98: \ < XA X2 = -k//Tf NOT DEFINED (cannot satisfy the far field condition) 2(1+ § )* lim \A = -k/2 k - 0 Note: k, C , C„, and CT are defined in Section 1.3 p M I The X, j = 1, 2, 3, and 4,correspond to the rate constants for the exact solution, invariant size, invariant field and linearized gradients approximations respectively. - 14 -Exact Solution Invariant Size Approx. Invariant Field Approx. TABLE lib PRECIPITATE DISSOLUTION SPHERE Linearized Gradients Approx. (y + 2p ./> y + t) arctan /T _Y_ + p where' R/R , _ aft "R 2 2 _ kD, p2 - kAn and reduces to y2 = i - T in the limit as p (or k) goes to zero 2 - Tj2 or y2 = 1 - kDt o Sufficiently complex to negate its usefulness. 27 PLANE S = SQ - X.(Dt)2 >.x = 2X X eXZer£c(-\) = k/2 lim X1 = A.2 > X4 k - 0 k > 0: \2 > 0 ^ k ^ 1.5: ^ s ^i > X* X2 = k//ff NOT DEFINED (cp - cMF (Cp 2(1 + | )S lim XA - ~ k - 0 ^ - 15 -The fact that only a small fraction of the precipitates observed obey equation (4) is attributed to surface diffusion enhancement of dissolution. In most alloy systems of interest |k| < 0.3,and values of |k| < 0.1 are quite common. For this small k range, R is a slowly varying function of time,and the solution to equation (1) reduces to that 29 obtained by assuming a stationary interface. Whelan used the stationary interface approximation to obtain the following equation for the dissolution rate: dR - k , fD , dt" = 2 [D/R + 1 6 The D/R term arises from the steady state part of the diffusion field (Equation 4)yand the term arises from the transient part of the field. In the case of planar precipitates (the one dimensional case), there is only a transient term in the expression for dissolution rate^ and one obtains: R = R - k ,F^~ 7 . O y 7T For spherical precipitates (the three dimensional case)9 both terms in equation (6) contribute to the dissolution rate. In the final stages of dissolution of spherical precipitates, the transient term becomes negligible and equation (6) reduces to equation (4). A good physical description of the diffusion fields around dissolving spheres is 30 given by Readey and Cooper. - 16 -1.3.2 Dissolution Studies Using the Electron Probe Microanalyzer. 31 Hall and Hayworth used electron probe microanalysis to determine the concentration profiles around dissolving Widmanstatten plates and needles in an Al-5 wt.% Cu alloy. They found that at high homologous temperatures (T = 0.91-0.95), the plates dissolve via a volume n diffusion mechanism and interface equilibrium is maintained. Eifert 32 et al. used the probe to observe the g (B.C.C.) -* a(F.C.C) transformation in the Cu-12.5 wt.% Al system. They reported a diffusion-controlled lattice transformation in which interface equilibrium is not maintained. 1.4 Dissolution Mechanisms Derived from Precipitate Growth 33 34 Both Tanzili and Heckel and Aaron have proposed an experimental procedure, based on theoretical considerations, whereby dissolution experiments may be used to determine the growth mode of precipitates. That the observed growth mode cannot be used to arrive at dissolution 27 kinetics has been shown by Aaron et al. In their treatment of dissolution at the reverse of growth, they show that to treat dissolu tion as the mathematical inverse of growth is at best a crude approximation, as the boundary conditions are intrinsically different 27 in the two cases. Table II summarizes the relationship between supersaturation (k) and transformation rate constant (^j)> f°r tne different mathematical models used to determine growth and dissolution kinetics. 22 It has further been shown by Nolfi et al. that curvature and interface reaction kinetics play an important part in the early stages of growth and that both tend to decrease the growth velocity. - 17 -Curvature effects on the other hand, are only important in the later stages of dissolution and tend to increase the dissolution velocity. The inference that research on precipitate growth precludes the need for extensive direct investigation of precipitate dissolution, may 35 be shown to be incorrect by the following reasoning. In the case of growth,the formation of a critical nucleus, after some incubation time, is followed by the growth of the precipitate, causing depletion of solute in the matrix immediately ahead of the advancing interface (as shown in Figure 3(a)). At any position r in the matrix phase (r > R(t)), the solute concentration is a monotonically decreasing function of time. During dissolution, on the other hand, there is no incubation or nucleation period. The precipitate has a finite size (RQ) and decreases in size by mass transfer of solute into the matrix behind the receding interface (as shown in Figure 3(b)). The variation of solute concentration with time for dissolution is different and more complex than in growth. At a position well away from the precipitate (» > r >> R(t)), the solute concentration increases with time; close to the precipitate (r > R(t)), the concentration decreases with time; and at intermediate positions, the solute concentration will increase, decrease or remain unchanged. Furthermore, the extent of the regions corresponding to 'near', 'intermediate' and 'far' vary with time. 1.5 Objective of Present Investigation The purpose of the present investigation is to provide experimental data on the dissolution of grain boundary allotriomorphs and thereby - 19 n ellucidate the dissolution process. An attempt to arrive at a satisfactory theoretical analysis to account for the results will be made. The atomic processes involved in dissolution of precipitates have thus far been indirectly deduced from statistical and other re lated kinetic measurements. The objective of this investigation is to determine the mechanisms of the dissolution phenomenon in a more direct and fundamental way by the use of electron probe microanalysis. - 20 -2. EXPERIMENTAL PROCEDURE The electron probe microanalyzer was used as an instrument for both quantitative and qualitative analysis. Quantitative electron probe microanalysis was used to determine the concentration profiles around dissolving grain boundary allotriomorphs. A technique was established whereby the concentration profiles around a single precipitate were measured as a function of dissolution time. Kinetic studies of the dissolution rate of grain boundary allotriomorphs were made using the back-scattered electron image on the electron probe to determine the size of a precipitate as a function of dissolution time. Here once again, the dissolution of individual precipiatates was observed. 2.1 Reasons for Systems Chosen The aluminium-copper system was chosen for the following reasons: 1. the solidus composition changes rapidly with temperature, giving rise to a significant variation in concentration around a 36 37 dissolving precipitate, see the phase diagram (Figure 4); ' 2. the low vapour pressure of the alloy at the experimental temperatures; 3. the existence of reliable phase diagram data; - 21 -4. the well-established relationship between precipitate shape and grain boundary misorientations; and another consideration is 5. the high degree of accuracy that can be achieved in the quantitative analysis of these alloys using the electron probe microanalyzer. This accuracy is due to the small correction required to convert measured Cu X-ray intensity to weight percent of Cu present in an alloy. 36 The aluminium-silver system ' (Figure 5), was chosen for supplementary experiments as it also fulfills all the above criteria. 2.2 Preparation of Master Alloys Four Al-Cu master alloys(A,B,C,and D) of composition Al-4.83wt.% Cu, Al-2.81wt.% Cu,Al-1.57wt.% Cu, and Al-1.22wt.% Cu were prepared by melting 99.99% pure copper and 99.99% pure aluminium in a graphite crucible at 720°C. This was followed by casting in a split graphite mold to produce a slab of dimensions 6"x 2 3/8" x 1/2". An aluminium-15.75 wt.% silver alloy was produced in the same way using silver of 99.95% purity. Following casting, the slabs were homogenized for seven days at 550°C, then cold-rolled in several passes to give a reduction in thickness of 50%. This was followed by a two hour anneal at 550°C and further cold rolling to give a final sheet thickness of 0.05". The sheet was then given a final anneal at 550°C for five days. The resultant sheet had an equiaxed grain structure, the grain size varying from 2.5 to 4 mm. The grain boundaries were straight and extended from the top to the bottom surface of the sheet. The majority of the grain boundaries were found to be stable on subsequent heat treatment. Figure 4. °C 800 (a) The aluminum-copper equilibrium phase diagram °F 36 Atomic Percentage Copper / 2 3 700 600 500 400 300 200 100 Al 'DC 2 B3 4 A5 6 7 8 Weight Percentage Copper i i | L i 1 -660° is i t | K K + L 1 548" i 1 -y. — I i l! 1 1 1 I 1 l | 1 1 --I-1 1 1 1 ' 1 1 1 j I 1 i i -1400 1200 1000 800 600 400 10 Figure 4. (b) The aluminum rich end of the aluminum-copper equilibrium phase diagram. - 23 -Figure 5. The aluminum-silver equilibrium phase diagram. - 24 -The copper content of the aluminium-copper alloys, as well as the silver content of the aluminium-silver alloy, were determined using two different techniques: (the analyses were performed by Warnock-Hersey International Ltd. and Cantest Ltd.). (a) The Atomic Absorption technique, which has an accuracy of + 0.02% up to 2% Cu and + 0.05% for up to 5% Cu. (b) The electroplating of copper on a platinum electrode,which has an accuracy of 0.02% for the range of interest. The analyses were in all cases within + 0.08% of the cast alloy compositions based on the weight of components used. Microscopic sections of specimens were scanned in the electron probe microanalyzer to evaluate the homogeneity of the material. No significant microsegregation was detected in a series of samples taken across the width of the sheet. 2.3 Procedure Used to Grow Large Grain Boundary Allotriomorphs To study the concentration profiles associated with individual precipitates, a heat treating procedure had to be established which would produce: (1) a minimum number of precipitates, (2) precipitates large enough in size for accurate probe analysis (thickness of precipitates > 4y), and (3) large precipitate free zones adjacent to the grain boundaries, as well as large interparticle spacings on the grain boundaries to avoid the impingement of neighbouring diffusion fields at short dissolution times. The desired microstructure was produced in the following way. Ten discs of 1 1/4" diameter were cut from the final sheet of a particular alloy. These discs were tightly stacked together, tied with - 25 -chromel wire>and a hole 1/16" in diameter was made through half of the discs to accommodate a temperature controlling chromel-alumel thermocouple. The assembly was inserted into a 2" diameter vertical tube furnace with a 4" long maximum temperature zone at 550°C. The temperature was controlled to within + 1°C for prolonged periods using a deviation amplifier. After three days the temperature of the furnace was reduced slowly to a temperature approximately 10°C above, the phase diagram solvus temperature. Further cooling at 3°C per day was used to nucleate and grow precipitates. The slow cooling was continued until the precipitates had grown to a suitable size for analysis. The sample was then allowed to equilibrate for seven days to remove any concentration gradients introduced into the matrix by the precipitate growth. The samples were finally quenched in ice water from the equilibration temperatures (T ) (as depicted in Figures 4 and 5). A typical equilibrated microstructure is shown in Figure 6. Small test samples were quenched and optically inspected at various stages of the heat treatment. It was established that nucleation occurred at temperatures within + 3°C of the literature phase diagram solvus temperature. The number of grain boundary allotriomorphs present just below the solvus was comparable to that observed in the equilibrated structure. 2.4 Dissolution Heat Treatment Procedure Specimens approximately 3/8" x 1/4" were cut from the discs (in some cases spark machined), cold mounted and polished to 1 u diamond. Figure 6. Optical micrograph of a typical alloy microstructure for the Al-2.81 wt.% Cu alloy. - 27 -Mechanical polishing,was used in preference to electro-polishing,as this resulted in a surface finish more suitable for microprobe analysis. Polishing produced a slight etching of the surface in the case of the Al-4.83 wt.% Cu, the Al-2.83 wt.% Cu and the Al-Ag alloy, facilitating the location of appropriate grain boundary allotrio morphs. For the other alloys suitable allotriomorphs were selected using the absorbed electron image on the microprobe, thus eliminating the need for etching the specimen to outline the structure. The same precipitates could be relocated after several heat treatments by the probe contamination marks and by micro-hardness indentations made at strategic points on the surface. For dissolution heat treatment,the specimens were immersed in a 60% KNO.j-40% NaNO^ salt bath controlled to + 1°C by a steady state power input unit. Surface contamination by the salt was avoided by wrapping the specimens in 0.0005" thick, high purity tantalum foil. The salt did not penetrate the wrapping foil because of the high viscosity of the salt at the temperatures used. The immersion of the small mass of the sample produced a negligible temperature change in the salt (< 1°C). Temperature stabilization occurred in less than 30 seconds after immersion of the specimen in the salt. It is estimated that the specimen reached bath temperature within two seconds of immersion in the salt. After heat treatment the specimens were quenched in water, the transfer time from the salt to the beaker being less than 1.5 seconds. Using this procedure, as many as twenty successive heat treatments - 28 -were carried out on some specimens without any observable trace of surface contamination by the salt. 2.5 Diffusion Couple Work Diffusion couples were prepared by clamping an aluminium bi-crystal, with a grain boundary angle R = 17°, to a specimen with a structure containing large grains of CuA^ surrounded by a fine eutectic of a + 0 (see Figure 12(b)). The aluminium bi-crystal was used to permit the measurement of both grain boundary and volume diffusion coefficients. The specimen containing the CuA^ grains was produced by melting stoichiometric amounts of copper (53 wt.%) and aluminium in a vacuum sealed quartz tube at 700°C and quenching in brine. Several attempts were made to produce pure CuA^ without success. The bonding surfaces of both the aluminium and the CuA^ were polished to 6 y diamond and ultrasonically cleaned prior to bonding. The couples were diffused at a fixed temperature in a cracked ammonia atmosphere. The temperature-controlling thermocouple was placed just above the couple and the maximum temperature fluctuation registered was < 3°C. Unfortunately, bonding was achieved only in . four couples at temperatures of 480°C, 500°C, 520°C, and 535°C. Repeated attempts at lower temperatures were unsuccessful, probably due to the presence of the oxide film on the aluminium. - 29 -2.6 Electron Probe Microanalysis 2.6.1 Isoconcentration Contours A JEOLCO (JXA-3A) electron probe microanalyzer with a step-scan attachment was used to determine the composition contours around dissolving grain boundary allotriomorphs. Figure 7 gives a schematic representation of the specimen in the microprobe specimen holder and the micrometer screws, as well as the possible directions of movement of the specimen. The step-scan attachment is an electronic device that may be attached to either of the micrometer screws to automatically displace the specimen relative to the electron beam by steps of 1.25 u, 5 u or 20 u. Most of the measurements in this investigation used 1.25 u step scanning. Each displacement is followed by a stationary spot count of 10 seconds and the resultant counts are registered automatically. Grain boundary allotriomorphs situated on straight grain boundaries were selected for investigation and were rotated into a position per-v pendicular to one micrometer. Step scanning was then carried out across the width of the precipitate into the matrix. The scan was continued until the equilibrium matrix concentration was reached. This procedure was repeated along the length of the precipitate and the adjacent grain boundary. The scanning procedure is depicted in Figure 8. The dark contamination streaks mark the trace of each scan. An initial trace was taken on each precipitate before any dissolution heat treatment was carried out to give a blank profile. In the analysis of the Al-Cu alloys, one spectrometer was set to coincide with the CuK^, peak and the second was set at 30' above the 0 1 Specimen 2 Rotating Specimen Holder 3 Electron Beam 4 Spectrometers A&B Micrometer Screws X-Rays Figure 7. Schematic representation of the directions of movement of a specimen relative to the electron beam, in the electron probe microanalyzer. - 31 -Figure 8. Optical micrograph of a typical grain boundary allotriomorph, with contamination marks showing the full extent of scans across the precipitate and grain boundary. CuKa^ peak to measure the background radiation intensity. These settings were chosen in preference to the Al La^ radiation because the maximum variation of the aluminium concentration in these alloys is between 95 wt.% and 98.5 wt.% Al,and therefore cannot be picked up with any degree of accuracy. In the Al-Ag alloy, the Ag Lct-^ peak and a background at 30' above the peak were registered. Appendix I gives the details of the correction procedure used to convert measured Cu Ka^ and Ag Lcc^ intensities to weight percent coppi and weight percent silver, respectively. 2.6.1.1 Resolution of Experimental Technique The accelerating potential used in the evaluation of the concen-38 tration profiles was 25 K.V. This satisfies the criterion which has been established for optimum accuracy in quantitative analysis; K where: V = accelerating potential V„ = excitation potential for any particular radiation. The principal resolution.limitation in quantitative electron probe microanalysis , is the finite spot size (i.e., the volume from which X-rays are generated for any one spot count). The experimental conditions in this work resulted in a theoretically-calculated spot size of about 6 u in diameter in the aluminium-rich matrix,and 3u in diameter in the copper-rich CuA^. In the Al-Ag alloy, the spot size is approximately 5 u in the aluminium-rich matrix and 2.5 u in the silver-rich phase (see Appendix II for calculations). The spot - 33 -size was also determined experimentally from the measured Cu concen tration in the matrix adjacent to a precipitate interface. A copper concentration above the equilibrium matrix concentration was found to extend over a distance of approximately 5 u in the blank profile (Figure 9). This increase at a distance 5 y away from the determined interface was never greater than 1-2% above the matrix composition. Since 1.25 y step scanning was used,the effective spot size will be > 3.75 y in radius,which is in good agreement with the theoretically calculated value. That step-scanning in 1.25 y steps gives meaningful values of concentration, even though the spot size is - 3.75 y in radius, may be appreciated if one considers that an error function profile (erf x) of C versus x is very close to a linear profile over the range x=0tox=0.8. Accordingly for a linear composition profile, adjacent to the allotriomorph interface, the mean composition of the circular area with centre x^ as in Figure 10, is equal to the C^ the composition at the point x^. In the present work the spot size limits the X-ray resolution at the precipitate-matrix interface and the grain boundary. The position of the precipitate-matrix interface was determined from the blank profile by a large decrease in specimen current and a sudden increase in background count, associated with going from the aluminium-rich matrix to the copper-rich precipitate (see Figure 9). The grain boundary position was roughly established from the position of adjacent precipitates,to the precipitate being observed,along the grain boundary (see Figure 6). it was found that the maximum solute concentration on traversing a grain boundary - 34 -E X c 0 k_ u c a) E U 'u cu 10 a OS o "* o u O 00 150 h 100h 501-0^ D u 2000 |r 1000 CM 0 Figure 9 cp I 4000 30001-The (a) specimen current, (b) background concentration, and (c) Cu-Kc^ concentration,for the blank profile (1) through the centre of the precipitate, (2) at the allotriomorph tip and (3) 5u away from the tip across the grain boundary. - 35 -Figure 10. Schematic diagram depicting the effect of spot size on the evaluation of concentration at a point x^ away from the precipitate interface. Figure 11. A schematic depiction of the dissolution flux lines for; (a) an allotriomorph sectioned above the minor axis b, (b) an allo triomorph sectioned at the minor axis,and (c) an allotriomorph sectioned below the minor axis. - 36 -invariably coincided with the position of the grain boundary as established from the back scattered electron image. Concentration gradients on either side of the maximum concentration were also very similar, providing further verification of the grain boundary position. 2.6.2 Kinetic Studies Attempts were made to study precipitate dissolution by heat-treating the specimen, etching and examining the micro-structure using an optical microscope. As it was necessary to polish the specimen each time prior to etching, a different position on the precipitate was observed after each heat treatment. Thus the resulting data is not useful for a kinetic analysis. The back-scattered electron image (b.s.e.i.) of the microprobe . proved much more satisfactory since it did not require etching of the specimen. The experimental procedure developed approximates very well the direct observation of precipitate dissolution. The procedure involves photographing the back-scattered electron image of individual grain boundary allotriomorphs after successive heat treatments in the salt. A polaroid camera,attached to a monitor on the microprobe console, was used with polaroid N52 type film. The largest allotriomorphs found in any one specimen were used for kinetic studies. This was done to ensure that the allotriomorphs studied were sectioned as close to the maximum diameter as possible. As shown in Figure 11, precipitates sectioned above the maximum half length will dissolve more slowly than precipitates sectioned through the - 37 -center, whereas precipitates sectioned below the maximum half length will dissolve faster. For the two-dimensional elliptical geometry depicted in Figure 11, the divergence of the flux lines at x^,(Figure 11(b)),may be related to the change in the angle of the tangent at x^, relative to the tangent at the maximum half thickness b. For an ellipse with an axial ratio of 3.5 to 1,the change in angle A0 is 4° and 10° for x^ equal to 0.25 a and 0.5 a respectively (see Appendix III for calculations). As the precaution was taken to study the dissolution of only the largest grain boundary allotriomorphs in any one sample, tha effect of the flux line divergence, which as calculated is quite small, is considered negligible in this work. Photographic observations were made at magnifications of x500, xlOOO, x2084, and x4168; at higher magnifications the resolution was found to be impaired. f Measurements of the length of precipitates were made directly from the photograph to within 0.001" using vernier callipers. The thickness measurement was made at a position halfway between the tips of the precipitate. The lengthy heat treatments in this part of the work increased the possibility of contamination by the salt, of grain boundary migration and of extensive thermal grooving of the grain boundaries. If there was any detectable evidence of contamination of the specimen surface, the specimen was rejected. - 38 -2.6.3 Diffusion Couple Measurements An etchant of composition 47% HCl, 6% HI and 47% H20, when freshly prepared, delineated the copper composition in diffusion couples. The etchant was used to outline the diffusion field and thus to verify the extent of bonding which had taken place. For each successfully prepared diffusion couple, at least two composition profiles (in some cases three) were taken across the centre of the two largest CuA^ grains at the bonded interface (see Figure 12(a)). The step scan device was used, and in most cases 1.25 u steps were taken. 2.7 Validity of Experimental Technique 2.7.1 Isoconcentration Contour Measurements The validity of the.experimental technique depends on whether the dissolution contours measured on the sample surface are representative of those present about precipitates in the interior of the sample. Experiments were carried out to establish that the surface measurements in this investigation are representative of the bulk material. This was done by measuring dissolution contours as follows: (A) around a given precipitate at the surface; (B) around the same precipitate after a layer of material 2-3 u thick was carefully removed from the specimen surface; and (C) around a second precipitate in the bulk material. The resultant dissolution contours are shown in Figure 13. Within reasonable scatter results(A),(B)and (C) coincide, ruling out any - 39 -Figure 12(b). Microstructure of large CuAl2 grain in an a + 9 matrix. TH=0.97 Alloy-A Figure 13. Dissolution isoconcentration contours of a grain boundary allotriomorph on the surface and 3u below the surface (A and B) and a second allotriomorph in the bulk material (C). At T = 0.97 for H alloy A. - 41 -surface diffusion effect. This also indicates that the effect of flux divergence is small relative to the experimental error, an observa tion that is important for the kinetic studies. A further important question was whether impingement from adjacent diffusion fields just below the surface was affecting the measured results. To detect this possibility: (1) A scan was obtained along the direction of the precipitate minor axis from the bulk matrix on one side to the bulk matrix on the other side. (2) A parallel scan was also obtained just beyond the tips on either side of the precipitate continuing to a distance, on both sides the grain boundary, associated with an equilibrium matrix composition. The total diffusion field around the precipitate could then be examined, and any anomalies in the diffusion field became apparent. 2.7.2 Kinetic Measurements Individual grain boundary precipitates were examined during dissolution. No statistical analysis of dissolution kinetics was carried out in the present work, as is generally done, because there were a relatively small number of precipitates on the grain boundaries, with a large range of precipitate shapes (s /R varying o o from > 10/1 to 2/1) (see Figure 14), With the small number of precipitates in any one specimen, the population distribution curves are inadequate for statistical analysis. 20 18 16 z «/> 14 | 0 12 'a <D 10 Q_ M— 0 8 m O Z 6 --• - >10 i 3.5 4 [So/Ro] 8 10 Figure 14. Population distribution curve for allotriomorphsof varying axial ratio S^/R^ for alloy A. - 43 -To determine if surface diffusion increased the dissolution rate in the kinetic studies, four grain boundary allotriomorphs were examined after a 12 minute dissolution treatment at 527°C. This dissolution heat treatment resulted in a decrease in the initial half-thickness R of about 40%. The precipitate size was observed o using the back-scattered electron image. A very thin surface layer was removed, the specimen was etched with a 5% NaOH solution5and the precipitate measurements were taken again optically. For all four precipitates the optical measurements were in good agreement with measurements obtained using the b.s.e.i., thus eliminating the possibility of surface diffusion contribution to the dissolution. The use of the back-scattered electron image to measure the 'true' precipitate size was tested by measuring the precipitate size optically and comparing the results with those obtained by the b.s.e.i. and the absorbed electron image (a.e.i.). The absorbed electron image was found to be unreliable, as a change in contrast could give a drastically different precipitate size. The b.s.e.i., however, could be picked up at all settings without distorting the image. All the settings gave measurements of precipitate size equal to those obtained optically. - 44 -3. ISOCONCENTRATION CONTOUR RESULTS AND DISCUSSION 3.1 Introduction Experimental determinations of the isoconcentration contours around dissolving grain boundary allotriomorphs, in an effectively infinite matrix, will be presented and analyzed in this chapter. One anticipates that isoconcentration contours associated with dissolving grain boundary precipitates would be influenced by the following parameters: (1) Dissolution rate controlling mechanism, i.e., whether the dissolution process is diffusion controlled or interface reaction controlled. (2) The precipitate shape (S0/RQ). (3) Homologous temperature T . (4) The supersaturation (k) of the alloy. (5) The volume diffusion coefficient (D^) of the alloy. (6) The grain boundary misorientation (R) and the grain boundary diffusion coefficient D , g.b. These parameters have each been considered experimentally. The results are presented, following a discussion of: (a) The diffusion couple results,to establish whether D^ is independent of composition in the a phase; - 45 -(b) The microstructure of the equilibrated alloys,and (c) The nature of isoconcentration contours to establish the reproducibility of the technique used. 3.1.1 Diffusion Couple Results The concentration profiles of dissolving precipitates are strongly dependent on the value of the volume diffusion coefficient of a given material. The analysis of the concentration profiles is greatly simplified if D^ can be considered as independent of alloy composition. To establish to what extent D^ is related to composition, composition measurements were made on a series of Al-CuA^ diffusion couples annealed at 480°C, 500°C, 520°C, and 535°C for 163 hours. The results for the 500°C diffusion couple are shown in Figures 15 and 16. In Figure 15 three runs are shown which are essentially coincident, indicating reproducible measurements. The anomalous rise in concentration near the interface is believed due to a two hour dissolution heat treatment at 535°C. This was done to dissolve CuA^ precipitates formed in the diffusion zone upon furnace cooling of the couple. Representative values of concentration from Figure 15 are plotted on probability paper, versus distance from the couple interface. The result shows a linear dependence, clearly indicating that D^ is independent of composition over the complete range of solid solution in the a phase,( see Figure 16). (The complete diffusion couple data is listed in Appendix IV). The temperature dependence of D^ from the present results is shown in an Arrhenius plot in Figure 17, along 39 with values of D for an Al-0.5 wt.% Cu alloy published by Murphy. 3OO0f-400(-100C[-J 30* COMPOSITION PROFILES - 500 C •O Run-I • • Run-2 Run-3 550 600 tso Pittance, (H) Figure 15. Composition contour for Al-CuAl„ diffusion couple at 500°C. - 47 -Diffusion Couple-500C X x 107 Cm. Figure 16. Probability plot of atomic % Cu versus diffusion distance (x) over the solid solubility range of the a phase for the 500°C diffusion couple. - 48 -VTxlO4,^"1 Figure 17. Arrhenius plot of D versus l/T°K_1for the diffusion couple results. - 49 -The present results are observed to follow an Arrhenius plot and agree closely with the values of Murphy. Since the results reported by Murphy were obtained for a wider temperature range than in the present case, his values of will be used in the analysis of the present results. Grain boundary diffusion coefficients could not be determined, as the pressure applied to achieve bonding caused a certain:.amount of recrystallization in the bi-crystal. 3.1.2 The Microstructure of Equilibrated Alloys The shape of typical grain boundary allotriomorphs in this investigation can be observed in Figure 6 (on page 26 ) which is a representative micrograph of alloy B. The precipitate shape is taken as oblate spheroidal and is defined by the length to width ratio SQ/RQ. The grain boundary allotriomorphs are widely spaced with 3 < SQ/RQ < 6. A number of "snakes" of high axial ratio are also observed in the micrographs, which are believed to result from the impingement of two or more adjacent allotriomorphs. Widmanstatten plates can be observed, however, they are not present in the areas immediately adjacent to the grain boundaries. The presence of a precipitate-free zone adjacent to the grain boundaries is a result of the depletion of vacancies in this area by the grain boundary precitates which form first. The vacancy concentration is 40 41 as a result too low for the nucleation of a second phase. ' - 50 -The size distribution (Sq) of the grain boundary precipitates obtained in Alloy A is plotted in Figure 18 as a function of the axial ratio (S /R ) for the precipitate distribution plotted in o o Figure 14. The average allotriomorph size and the width of the depleted zone is shown in Table III, for the different alloy compositions studied. Table III. Average Sq for Allotriomorphs with S /R =3.5 and Average o o 6 Width of Precipitate Free Zone for the Different Alloys. Alloy Half length of Width of Precipitate allotriomorphs(SQ)y free zone y A 25 90 B 15 70 C 10 60 D 10 30 3.1.3 Dissolution Isoconcentration Contours Five separate grain boundary allotriomorphs were examined by electron probe microanalysis and concentration contours determined for both sides of each precipitate. The results are shown in Figure 19. (One side of one precipitate was not used because of early impingement from an adjacent allotriomorph.) The top half of the figure shows dissolution contours after 90 seconds from the start of dissolution, the bottom after 360 seconds. The five allotriomorphs were in an alloy B specimen (Al-2.81 wt.% Cu) which was homogenized at 400°C (T ) and then upquenched to 520°C (T ). The maximum solubility of Cu 60 -50 ->10 3.5 4 [SQ/RQ] 8 10 Figure 18. Size distribution (S ) of precipitates as a function of S /R for alloy A. ° o o Figure 19. Isoconcentration contours at T = 0.91 for alloy B. - 53 -at 400°C in a is 1.9 wt.%; at 520°C it is 5.2 wt.%. Accordingly the precipitates will dissolve on upquenching until they ultimately disappear. The concentration measurements were confined to the early part of the dissolution process to avoid interaction effects due to solute fields from adjacent allotriomorphs. All composition measure ments were made on scans perpendicular to the allotriomorph equatorial axis and grain boundary. In general, concentration measure ments could not be made at distances less than 6 u from the CuA^ precipitate, since within this distance Cu in the precipitate could contribute to the measured concentration. In Figure 19 two values of concentration are plotted, 2.1 and 2.5 wt.% Cu. Higher values of concentration occurred within 6 u of the precipitate surface and were therefore considered unreliable. Lower values were approaching the matrix composition (1.9 wt.% Cu), thus making composition differences comparable to the experimental uncertainty of the measurements. Examining the values of concentra tion plotted in Figure 19, it is observed that the position of points from sides I and II of a given precipitate are in close proximity indicating a high degree of symmetry of the isoconcentration contours about a given precipitate. However, points of equal concentration for the five precipitates are not in exact coincidence, as might be expected. The points for 2.5 wt.% Cu can essentially be enclosed in an envelope, shown by the two solid lines closest to the precipitate. The 2.1 wt.% Cu points are contained in the envelope furthest from tfie precipitate. To a first approximation, lines drawn through the centre of each zone can be taken as the isoconcentration contours for - 54 -the composition being considered. The standard deviation (a) for these profiles varies from 0.68 u at a mean distance (x) of 6.11 u from the precipitate interface, to a standard deviation (a) of 1.17 u at a mean distance (x) of 17.59 u from the precipitate interface. (See Appendix V for calculations). In the first case a is 11.12% of the mean distance x whereas in the second case a is only 6.6% of the mean distance (x)(calculations in Appendix V). For the evaluation of and C^. values, the longer time dissolution experiments with diffusion distance greater than 20 u from the interface were used, wherever possible, to minimize the effect of scatter just established as inherent in this experimental procedure. 3.2 Rate Controlling Mechanism To establish that the dissolution process is diffusion controlled, values of concentrations at various time intervals were plotted as a function of x//t , where x is the distance along the minor axis of the precipitate starting at the precipitate surface,and t is the time. The average result for two precipitates in the Al-4.83 wt.% Cu alloy at a dissolution temperature of 545°C (T = 0.97) is shown in Figure 20. The points are observed to be effectively coincident, indicating that the composition C(x,t) = f( )• The profile vTTt obtained extrapolates at =0, to an interface concentration CT = /t 1 5.5. This is in good agreement with the phase diagram indicating local equilibrium at the interface and therefore no interface reaction control. A similar analysis of data for the range of homologous temperatures between T = 0.72 and 0.97 gave similar results for diffusion distances 5.5 5-0 D 4.5 V 4.0 3.5 3-0 h • A AO TH= 0.97 Alloy -A. T = 545C Ttme(t) Seconds • 100 • 225 • 400 A 900 • 1600 ± 0.2 0-4 0-6 0.8 1.0 x x 107 Cm. t 1.2 1.4 1.6 Figure 20. Concentration profile along the minor axis of a grain broundary allotriomorph at 545°C fo alloy A. - 56 -along the minor axis of the allotriomorphs. 3.3 The Effect of the Precipitate Shape (S /R ) , . , £ . E 0 0 To determine the effect of precipitate shape on the concentration contours,three precipitates were examined with different axial ratios: (1) a typical allotriomorph with S^/R^ = 3.5, (2) a snake with S^/R^ = 18 and (3) a nearly spherical precipitate having SQ/RQ =1.7 The resultant concentration contours are shown in Figure 21 for concentrations of 3.25 and 3.50 wt.% Cu and dissolution times of 30 and 120 seconds. The dissolution temperature was the same for all three precipitates. The results shown in Figure 21 clearly indicate that the composition contours for the three precipitates are essentially coincident,and therefore that under these experimental conditions the distance from the precipitate surface of the isoconcentration contours is independent of the axial ratio. This is not meant to imply that the overall shape of isoconcentration contours about a sphere and a snake are the same, as may be misconstrued from Figure 21. It should be pointed out that in Figure 21 as in all the figures of isoconcentra tion contours the points at the extreme left hand side denote average diffusion distances from the interface at the precipitate centre. The isoconcentration contours attempt to give an accurate representation of the isoconcentration lines, about the tip of an allotriomorph, and along the adjacent grain boundary, with only an average diffusional distance away from the broad allotriomorph faces. In general, two effects have to be considered for dissolution of precipitates, of different axial ratios: (a) the Gibbs-Thomson effect Figure 21. Isoconcentration contours for three allotriomorphs having axial ratios S /R of 1) 3,5, 2) 18 3) 1.7 for alloy A at T = 0.93. - 58 -and (b) the point effect of diffusion. (a) The Gibbs-Thomson effect is the variation with curvature of the equilibrium solute concentration at the precipitate-matrix 42 interface. Aaron and Kotler have calculated the influence of this effect on the dissolution kinetics of Al-Cu alloys and conclude that it is negligible for the total dissolution time of a precipitate. Therefore the influence of interface curvature on the isoconcentration con tours is expected to be minimal. (b) Because of the point effect of diffusion, the flux from the tip of an allotriomorph is distributed over a greater solid angle than the flux perpendicular to the side of an allotriomorph. Thus the point effect tends to round off the isoconcentration contours at the tip of an allotriomorph. As a further verification that the Gibbs-Thomson effect does not influence the position of the isoconcentration curves around a dissolving precipitate for short dissolution times, concentration measurements were made around two Widmanstatten plates located in the centre of a grain under the same conditions as those for Figure 21. The average results are shown in Figure 22 and show comparable diffusion distances, along the minor axis, to those for the grain boundary allotriomorphs in Figure 21. Hence, the shape of the precipitate does not influence the position of the isoconcentration curve relative to the precipitate surface. In addition, Widmanstatten plates, like grain boundary allotriomorphs, are found to dissolve by a diffusion mechanism in 31 agreement with the results of Hall and Hayworth in an Al-5 wt.% Cu alloy. Figure 22. Isoconcentration contours of two Widmanstatten plates for alloy A at T = 0.93. - 60 -3.4 Homologous Temperature To determine the effect of T„ on the position and shape of the n isoconcentration lines around an allotriomorph, concentration measurements were made on a series of alloys, over a range of dissolution times, for 0.72 < T < 0.97. The specific experiments H are tabulated in Table IV. Some of the corresponding concentration contours are given in Figures 23-29 for the times and concentrations indicated. No scatter bars have been placed on the contour plots; the scatter as established in Section 3.1.3 is expected to be the same for these results. Comparing the isoconcentration contours for T = 0.97 (Figure 23) H with T = 0.77 (Figure 29), it is observed that the isoconcentration rl lines in the former are inclined at approximately 90° to the grain boundary plane, whereas in the latter the corresponding angle is about 20°. The decreasing of the angle of inclination is clear evidence that enhanced grain boundary diffusion has occurred at the lower homologous temperature. Figures 23-29 show a progressive decrease in the angle of inclination with decreasing T and therefore a rl progressively greater contribution from grain boundary diffusion to the dissolution process. The contours for the highest homologous temperature T = 0.97 (Figure 23) exhibit essentially no grain rl boundary contribution to dissolution. The slight rounding-off of the contours at the precipitate tip is thought to be due to the point effect of diffusion. At T = 0.92 there is a very small amount of grain boundary ri diffusion and the diffusion distance along the grain boundary at the tip of the precipitate is slightly greater than the diffusion distance - 61 -Table IV. Contour Experiments T H alloy No. of allotriomorphs Dissolution times (seconds) Figure No. .97 A 2 25,100,225,400,900,1600 23 .94 B 1 400,900,1600 24 .93 A * 5 30,120 21 and 22 .93 Al-Ag 3 40,160 43 .92 A ** 3 100,400,1600 34 .92 B 2 400,900,1600,3200 25 .91 B 9 90,360 19 .905 B 1 400,900,1600 26 .885 B 2 500,1600 45 .885 Al-Ag 2 50,100,200 42 .86 B 1 400,900,1600 27 .84 Al-Ag *** 3 70,150,225,400 35 and 36 .83 B 1 100,400,1600,6400 44 .80 C 1 1600,6400 28 .77 C 2 1600,3600,6400 29 .72 C 2 1600,3600 30 2 Widmanstatten plates, plus 3 different shape grain boundary precipitates. ** 1 precipitate only at 1600 seconds *** 2 precpitates on low angle grain boundary Figure 23. Isoconcentration contours for alloy A at T = 0.97. Figure 25. Isoconcentration contours for alloy B at T = 0.92. Figure 26. Isoconcentration contours for alloy B at T = 0.905. Figure 27. Isoconcentration contours for alloy B at T = 0.86. Figure 28. Isoconcentration contours for alloy C at T = 0.80. TH=Q.77 Alby-C wt./£Cu oUO • UO A 1.60 Figure 29. Isoconcentration contours for alloy C at T = 0.77. - 69 -perpendicular to its flat surface. As the homologous temperature decreases, this effect becomes more pronounced until at the lowest temperature studied T = 0.72 (Figure 30), the diffusion fields along the grain boundaries impinge with those of adjacent allotriomorphs before there is enough volume diffusion from the precipitate to be resolved by the electron probe. The grain boundary diffusion contribution to growth was estimated 21 experimentally using statistical techniques, by Goldman et al. Their results indicate that below T = 0.78 solute transport should be due solely to grain boundary diffusion. No experiments were done above T„ = 0.78. However, between T„ = 0.79 and 0.90, they postulate volume diffusion should become progressively more dominant, and above T = 0.91 they predict transport should be entirely by rl volume diffusion. Hawbolt and Brown,^ in their studies of the growth of grain boundary allotriomorphs in Ag-5.64 wt.% Al, find that at TR = 0.92, growth is entirely volume diffusion controlled. At T = 0.90, however, rl there is some evidence in their results for grain boundary diffusion enhancing the growth rate of precipitates. The present results show quite conclusively that volume diffusion completely dominates the dissolution process above a homologous temperature of between 0.92 and 0.93. It is further shown that at TR = 0.77, the effect of volume diffusion on dissolution is negligible, and that there is essentially total grain boundary diffusion control at T = 0.72. Thus the present allotriomorph dissolution results are in excellent 21 6 agreement with the work of both Goldman et al. and Hawbolt and Brown TH^0.72 Alloy-D p 0.45wt.% Cu • 0-45 wt. L Cu Precipitate 1 Precipitate 2 o o • o • o • o a o 0 O ° O DO o o o • a • o o 10u f o = 3600 Seconds Q O o • o n O OD o o o a • a • • o o o o o p o Figure 30. Isoconcentration contour for alloy D at T = 0.72. - 71 -regarding the transition from complete volume diffusion control to complete grain boundary diffusion control for the growth of allotrio morphs. The relative importance of volume diffusion to grain boundary diffusion,in the transition region between complete volume diffusion control and complete grain boundary diffusion control, is also well illustrated for the first time experimentally in these results. In the following section the volume diffusion coefficient Dy is calculated using the concentration profiles parallel to the minor axis of the grain boundary allotriomorph,where enhancement by grain boundary diffusion can be neglected. Values of C^., the allotriomorph interface composition, are also determined. 3.4.1 Calculations of. and C^ Evaluations of and C^ are made to test the accuracy of the concentration contour experiments. The following assumptions were made in determining and C^.. (1) is independent of concentration. This is based on the diffusion couple results (Section 3.1.1) which indicate that the volume diffusion coefficient is independent of concentration over the complete solid solubility range of the a-phase. (2) The allotriomorph-matrix interface is essentially flat at the minor axis,so that a one-dimensional solution to the diffusion equation can be used with distances measured perpendicular to the centre of the flat surface of the allotriomorph. (Figure 20 indicates an error function profile in this direction). (3) There is very little movement of the interface relative to - 72 -the diffusional distance. The precipitate interface is considered stationary relative to the diffusion distances considered during dissolution. A calculation for the prevailing experimental conditions, 31 by Hall and Hayworth, indicates a diffusion distance which is approximately two orders of magnitude greater than the displacement of the interface. Subject to these assumptions, the solution to Fickls law reduces to: could not be measured experimentally due to the finite size of the electron beam in the electron probe, therefore, both C^ and have unknown values in equation 1. Values of and C^. were obtained from equation 1 using the following graphical analysis. (1) For three values of x (x = 10 y, 15 y, 20 y), the corres ponding experimental values of C, . were obtained from the iso-(x,t) concentration lines in Figure 25 and were used in conjunction with 1/2 a range of (D^t) values to solve equation 1. The graph of C^. vs. 1/2 (D^t) (Figure 31) consists of three curves which cross at a point 1/2 corresponding to approximate values of C^. and (D^t) (2) A plot of C, •> versus erf. [ — T77] was made using 1/2 *KUyt) the value of. (Dyt) obtained from Figure 31. This is a linear plot which is consistent with equation 1. When extrapolated to erf. [ r-TTT-] = 0 a more accurate value of CT is obtained as shown 2(Dyt)i/Z - 73 -- 74 -TH = 0.92 Alloy- B O Side • Side 0.2 0.4 0.6 0.8 Erf (^) Figure 32. C versus erf ( ) for alloy B at TR = 0.92 1.0 2>/Dt - 76 -in Figure 32. C(x t)_CM (3) Accurate values of D were obtained by plotting •' V VSl versus x on probability paper. can be calculated from the i slope of the line in Figure 33, which is equal to 2(Dt)1/2 ' (4) Iteration of (2) and (3) could be used to determine even more accurate values of and C^. This was, however, considered unnecessary considering the good agreement between the and values calculated and the literature values. Table V gives a comparison of the experimentally determined values with the 36 37 equilibrium phase diagram values. ' There is, in most cases, excellent agreement between the experimental and the literature values, which is consistent with the assumption that the dissolution process is diffusion controlled. Table V. Interface Concentration (C^) and Volume Diffusion Coefficients (Dy) for Contour Experiments. Alloy Experimental values Cj wt.% Cu Equilibrium diagram values 5.60 4.76 5.62 4.97 4.57 4.27 3.60 2.85 2.40 2.15 5.6 4.6 5.7 5.1 4.6 4.3 3.4 2.8 2.5 1.9 T H .97 .92 .940 .92 .90 .885 .86 .83 .80 .77 Experimental Dy xlO cm se 12.0 12.6 13.3 11.6 11.6 8.2 5.1 4.8 4.5 2.3 Murphy's values of D^ 13.2 9.5 13.2 10.5 9.5 6.1 3.0 1.9 1.5 0.6 Temperature °C 545 517 545 532 517 500 480 448 442 418 - 77 -The probability plots for the evaluation of were all found to be linear and thus indicated a concentration-independent diffusion coefficient, confirming the diffusion couple results. Table V also gives the calculated values of D^. The volume diffusion coefficients determined from the profile studies are in good agreement with the literature values at the high experimental temperatures. The experimental conditions at the low temperatures, resulted in: (1) smaller concentration gradients in the matrix, (2) shorter diffusion distances,and, (3) a greater possibility of a surface diffusion contribution to dissolution. This effect would enhance the dissolution rate but would leave the shape of the profiles unaffected. All these factors may be instrumental in the less accurate Dy values obtained at low temperatures. The activation energy for volume diffusion has not been evaluated from the results, as the good correlation between individual volume diffusion coefficients and literature values, as well as the excellent agreement between calculated C^. values and equilibrium phase diagram values, is adequate demonstration of the consistency and accuracy of the concentration contour experiments. 3.5 The Effect of Supersaturation (k) and Volume Diffusion Coefficient In any one alloy, changing the temperature changes both the volume diffusion coefficient (D.r) and the supersaturation (k). This makes - 78 -it difficult to estimate the individual contributions of these two variables. They will therefore be considered in one section,and an attempt will be made to show the effect of each of these factors on the shape of the isoconcentration contours. In order to determine the effect of k on the shape of the iso concentration lines, measurements were carried out at T = 0.92 for ri alloy A with k = 0.035. The results are shown in Figure 34,which can be compared to Figure 25,in which the homologous temperature is the same but k is appreciably higher (k = 0.094). The dissolution temperatures were different in the two cases (517°C for Figure 34 and 532°C for Figure 25). Comparing the two figures,it is evident that although the actual concentration values and diffusion distances are quite different the overall shape of the isoconcentration contours is the same. In the above experiment both the temperature and the super-saturation were different. An evaluation of the influence of these two variables will be attempted by making a comparison between pairs of experiments in which either T, k or T are kept constant. Five sets H of such experiments are listed in Table (VI), three sets have T constant, one set has T constant and one has k approximately constant. All three sets of isoconcentration lines having a constant T have ri the same shape. In both the cases of a constant volume diffusion coefficient and constant k, but different T , there is no similarity ri in the shape of the isoconcentration contour pairs. One may safely conclude that the influence of homologous temperature on the shape of the isoconcentration contours is independent of both the volume - 80 -Table VI. Comparison of the Effects of k, T and T on the Shape of n Isoconcentration Contours. Alloy Temp. °K H D xlO10 2 -1 cm sec S/sat. k Shape of Contours Fig. No. A B 790 805 .92 .92 7.3. 10.7 .035 .094 same 34 25 A B 790 790 .92 .905 7.3 7.3 .035 .078 different 34 26 Al-Ag A 823 800 .93 .93 50 9.4 .597 .045 same 43 21 Al-Ag B 786 777 .885 .885 22 5.4 .341 .0653 same 42 45 A B 818 753 .97 .86 7.3 2.8 0.515 0.511 different 23 27 - 81 -diffusion coefficient and the supersaturation. 3.6 The Effect of Grain Boundary Misorientation (R) The Al-15.75 wt.% Ag alloy had a micro-structure with a significant number of low angle grain boundaries. These were identified using the 4 Clark classification of precipitate morphology in this system. On one grain boundary it was found that nearly all the precipitates were in the form of primary side plates and sawteeth (R ^ 9.5°). However, two precipitates on this boundary resembled allotriomorphs., Composition contours were measured around these and were compared with composition contours obtained from an allotriomorph on a high angle grain boundary at the same temperature. In Figure 35 the profiles clearly show that grain boundary diffusion effects are not important for low angle grain boundaries, whereas in Figure 36 there is a significant grain boundary contribution to dissolution for a high angle grain boundary. Investigations of the effect of grain boundary misorientation 43 44 on the grain boundary diffusion coefficient ' have shown that beyond a misorientation of R 'v 20° there is very little significant variation of.D , . Since grain boundary, allotriomorphs generally g.b. form on grain boundaries with a misorientation R > 17°, the grain boundary contribution would not be expected to vary significantly from allotriomorph to allotriomorph. The experimental results appear to confirm this because the diffusion distances along the grain boundary vary by no more than 15% from precipitate to precipitate in any one system, for any one set of experimental conditions. Figure 35. Isoconcentration contours of two allotriomorph shaped precipitates situated on a low angle grain boundary for the Al-Ag alloy at T = 0.84. Figure 36. Isoconcentration contours for the Al-Ag alloy at T = 0.84. - 84 -3.6.1 Evaluation of the Grain Boundary Diffusion Coefficient (D ^ ) 45 A model developed by Shewmon for surface diffusion was used to evaluate D , in the present work. In Shewmon's model, solute is g.b. considered to quickly diffuse out along a surface from a cylindrical source. From the surface the solute then diffuses into the bulk of the material by a volume diffusion process. Shewmon's model is the 46 cylindrical analogue of Fisher's model for one-dimensional surface diffusion (see Figure 37). It suffers from the same approximations as used by Fisher; viSi: that diffusion occurs in the bulk only in the direction perpendicular to the surface and that steady state conditions on the surface are established early in the diffusion process. However Turnbull et al.^7in a critical evaluation of the Fisher Model and 48 the more sophisticated Whipple Model, found that D , values calculated from the different models differed only by a few percent. The Shewmon model is immediately applicable to dissolution from circular precipitates by diffusion along the grain boundaries and has 20 been used in this connection by Brailsford and Aaron. In this case the solution takes the form: C(p,z,x)=C (p,x) erfc(z/2v^~F) g. D . V where Cg.b.(p'T)= SAW [K0(a-P)/K0(a.b)] a2 = 2D../D 6(TID..T)1/2 V g.b. V - 85 -T is the annealing time, and b is the source radius. and the boundary conditions on C , (O.T) are g • • Cg.b.(b'T)= Ci g.b. M -H2b h-High Diffusivity Surface Layer Z •Source Low Diffusivity Material Figure 37. Schematic diagram for Shewmon surface diffusion model (after P.G. Shewmon). The concentration contours at low temperatures in the present work were used to determine D , from equations(2) and (3). In g.b. these contours there was little direct volume diffusion contribution to the gradient along the grain boundary. Figure 38 represents the calculations of D , . The theoretical profiles along the grain 1.0 0.8 3 u % 0,6 u 3 -0-41 0.2 r TH=0.86 Alloy-B W \ \W \ x X X \ \ \. \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ . "3 ^ 2 Dg.b.x lQ^cmPsecT^ 1 = 1 2=5 3 = 7.5 4 = 10 • Exptl. Results(t = 400sec.) 10 20 30 40 50 px lof Cm 60 70 Figure 38. Evaluation of D , for alloy B at T = 0.86. g-b. ^ H - 87 -boundary were calculated for various values of D , and the best fit g.b. with the experimental data of the dissolution of alloy B at a homologous temperature T = 0.86 (see Figure 27) was taken as the experimental value of ^ (see Figure 38). Grain boundary diffusion coefficients were determined for several temperatures, the results being shown in Table VII. An approximate grain boundary diffusion coefficient has also been calculated assuming: ^g.b. 2^volume ng.b. r«VOl. r> on 2. D° = Dn = 0.29 cm /sec The activation energy for volume diffusion is taken as Q , = 31.13 OJ ^volume 39 kcal/g.mole. The result is included in Table VII,for T = 480°C. The experimental value of ^ is in reasonable agreement with the value of D , calculated at 480°c. However, the inconsistency g.b. of the values of D , calculated at the other temperatures is not g.b. surprising considering:(1) the approximations made in arriving at values of Q , and D^*^'(2) that the experimental observations show g.b. 0 a diffusional distance along the grain boundary which obeys an almost parabolic relationship with time, whereas the sbewmon model 0 25 predicts a t " relationship, and (3) the fact that the experimental technique used is not one from which accurate evaluations of ^ are to be expected. These D , calculations are intended as peripheral g.b. information based on the experimental results. - 88 -Table VII. Grain Boundary Diffusion Coefficients (D b ) Determined from Contours 8* " T°C D , x 10 cm sec g.b. * 520 7.5 ** 504 9.0 480 0.75 448 7.5 420 5.0 480 0.83 (calculated) ' The grain boundary diffusion coefficient has been determined at these temperatures,even though the relative amount of volume to grain boundary diffusion is very large resulting in a high degree of error in the calculation of D - 89 -3.7 Semi-Quantitative Analysis of Isoconcentration Contours It has generally been found in diffusion studies that no effect due to grain boundary diffusion can be seen at temperatures above 0.6 T^. It is rather surprising therefore that in the present work a grain boundary diffusion contribution persists to a homologous temperature greater than 0.92. One of the essential differences between the two phenomena is statistical in nature. For grain boundary diffusion in a single phase polycrystalline material, one must consider the .mean time that an atom spends in the grain boundary as well as the probability of its getting onto the boundary. In dealing with the dissolution of grain boundary allotriomorphs, this statistical aspect is eliminated since the precipitate Is, in fact, situated on the grain boundary. A theoretical analysis of the relative contributions of grain boundary and volume diffusion to allotriomorph dissolution has been developed by D.E. Coates (see Appendix VI for details). The main features of this analysis are: 2 (1) The use of Laplace's equation (V C=0) to calculate the volume diffusion field about the allotriomorph. This approximation is good for low supersaturations and for longer diffusion times. (2) The use of two extreme conditions: (i) where volume diffusion is the sole contributor to dissolution, as in Figure 22, and (ii) where grain boundary diffusion is the overwhelming contributor to dissolution as in Figure 30. (3) A function ty is evaluated}where ty is the relative contribu tion to the total flux out of an allotriomorph due to grain boundary - 90 -diffusion and volume diffusion, viz: (f J . dA) ^ = —£L_J (Jl \ • dl)vol. If ^ >> 1 grain boundary diffusion dominates. If ty « 1 volume diffusion dominates. As shown in Appendix VI, ty is given by Kl[a V -Q/2RT 2e i 2(TTErc)1/2 f (a/b) (TTDT)1/2 where and Kq are modified Bessel functions of first and zero orde respectively and it is assumed -Q/2RT D /D = e g.b. V where Q is the activation energy for volume diffusion. 1/2 For a given precipitate size a and diffusion distance (TTDT) equation (2) is essentially an equation of the form ty = IJJ(T). ty is plotted in Figure 39 for the experimental conditions corresponding t three of the alloys used. The progressive transformation from grain boundary diffusion control to volume diffusion control with increasing temperature can b - 91 -- 92 -Figure 40. A comparison of theoretically derived i|; versus homologous tempera ture TR with the experimental observation of solute via grain boundary diffusion solute via volume diffusion - 93 -clearly seen. The three plots are drawn as a function of homologous temperature in Figure 40, and are found to coincide giving one curve. This is a fortuitous effect, as the increase in the precipitate size (precipitate half length a) with increasing alloy content compensates 1/2 for the longer experimental dissolution times (TTDT) (see equation 2). It is thus clear why for the present experimental conditions iji is only a function of T , rather than a function of the actual temperature. The dissolution of grain boundary allotriomorphs solely by grain boundary diffusion occurs below T = 0.72, which is in excellent ri agreement with the theoretically predicted value at which grain boundary diffusion dominates. The theoretical calculations confirm the significance of the contribution of grain boundary diffusion to the total flux even at high homologous temperatures. From the experimental data, it is possible to make an approximate calculation of the total flux due to volume diffusion and that due to grain boundary diffusion. The composition contour corresponding to volume diffusion is essentially an oblate spheroid. The extra contribution is assumed to be due to grain boundary diffusion as in Figure 41. By measuring the two areas and integrating over the solid angle corresponding to an oblate spheroid, the relative grain boundary and volume contributions to dissolution can be calculated The experimental points in Figure 40 represent the ratio: area due to g.b. diffusion .. , ^ , , . ^ •, . ... -,-, 5 ° .. ,. rr : , calculated as depicted schematically in area due to volume diffusion Figure 41 for all the profiles measured. The experimental results are in good agreement with the theoretical values. The theoretical analysis (equation 2) predicts no effect of supersaturation and little - 94 -1 2 3 4 G.B. Allotriomorph Solute via Volume Diffusion Solute via G.B. Diffusion Grain Boundary (G.B.) Figure 41. Schematic diagram of method used to calculate 1) solute via grain boundary diffusion, and 2) solute via volume diffusion. effect of axial ratio, both of which are in agreement with the experimental observations (see Appendix VI). 3.8 The Al-15.75 wt.% Ag System Figures 36, 42, and 43 show the concentration contours at three different homologous temperatures in the Al-Ag system. As can be seen by comparing these results with those in the Al-Cu system (viz: Figures 36 and 44, Figure 42 and Figure 23, Figure 43 and 45) there is excellent agreement in the shape of these contours at equivalent homologous temperatures. In both the Al-Cu and the Al-Ag systems, allotriomorphs having a relatively complex structure are dissolving into an f.c.c. matrix. The diffusion coefficients of these alloys differ by a factor of 39 49 five ' at equivalent homologous temperatures. The supersaturation Figure 42. Isoconcentration contours for the Al-Ag alloy at T = 0.885. Figure 43. Isoconcentration contours for the Al-Ag alloy at T = 0.93. TH=i0.83 Alloy-B Figure 44. Isoconcentration contour for alloy B at T = 0.83. - 99 -however may vary by more than an order of magnitude (see Table V). The observation that the isoconcentration contour shapes are comparable for comparable homologous temperatures confirms the singular importance of the homologous temperature; changing the supersaturation only alters the values for the concentration contours, not their shape. Calculations of C in this system are in good agreement with phase diagram data(as in Figure 5). is found to be independent of concentration over the total solid solubility range, as was found by 49 Heuman and Bohmer. Figure 46 compares values of calculated from the contours with values of Dy found in the literature. The agreement is good for both individual values of as well as for the calculation of the activation energy. 3.9 Applicability of Contours It has been found that the dissolution contours in both the Al-Cu and Al-Ag systems are very similar at comparable homologous tempera tures. It thus seems probable that the present results are applicable to any dissolving grain boundary allotriomorph. The diffusion coefficients in Al-Cu and Al-Ag are similar at similar homologous temperatures and so the transformation from grain boundary diffusion to volume diffusion control occurs at the same T in both cases. The H relationship of volume diffusion coefficient to homologous temperature in other systems, will determine the T at which the transformation between grain boundary and volume diffusion control occurs. The basic shapes of the contours will however, be the same and can simply be translated to different T1T values. - 101 -4. KINETIC RESULTS AND DISCUSSION 4.1 Introduction In this chapter, observed dissolution kinetics of grain boundary allotriomorphs are analysed. The results are discussed in terms of the most appropriate existing analytical models for allotriomorph dissolution. The latter include Whelan's spherical and planar models 28 29 50 for precipitate dissolution ' and Aaron's planar model. It is important to note that under the present experimental conditions, impingement of the diffusion fields from adjacent precipitates takes place at a very early stage in the dissolution process (see Table III). On the other hand, the three models just cited are based on the assumption that the diffusion fields of adjacent precipitates do not impinge at all (i.e., these models relate to the dissolution of a single precipitate in an infinite matrix). Accordingly it is anticipated that the correlation between these models and the present experimental results will be, at best, semi-quantitative. Thus an attempt is made to provide a firmer basis for interpretation of the results by formulating a dissolution model which accounts for impingement of the diffusion fields. - 102 -4.2 Analysis of Dissolution Kinetics Using Existing Planar and Spherical Models Table VIII gives a resume of all the kinetic experiments and includes the number of precipitates studied, the dissolution times and temperatures, and the supersaturation values (k). It should be noted that the time for complete dissolution of precipitates varied from 5 minutes to 116 hours, depending on the experimental conditions. 28 29 Whelan's spherical model ' was formulated for application to the later stages of the dissolution of spherical precipitates. Diffusion field impingement is completely ignored in the derivation. For the later stages one can approximate that only the steady state diffusion field contributes to the dissolution rate (see Section 1.3.1, equation (4)). 2 To apply the Whelan spherical model a plot of R versus t was made for each precipitate studied. If the model is valid, such plots should be linear. In fact, these plots were not linear, indicating the model is not consistent with the present data (see Figure 47). 2 It should be noted however that whereas the plots of R versus t are not linear for the full extent of the dissolution heat treatment. There is a tendency at the later stages of the dissolution process for 2 the relationship R versus t to become linear,indicating at least a 29 qualitative agreement with Whelan's experimental observations. 29 50 Whelan's planar and Aaron's planar models are formulated for planar precipitates and for the early stages of dissolution before impingement takes place,and are thus not applicable to the present set _ 103 -Table VIII. Kinetic Experiments. Alloy Number of Precipitates T H Super-saturation (k) T°C Time D 2 0.76 0.014 424 0,10,20,37,65,125,180,125, 307,442,560,686,816,1003 min Al-Ag 6 0.84 0.123 460 116,220,320,457,678,1123 min C 6 0.79 0.036 450 0,55,110,220,440,700 min B 4 0.83 0.032 460 0,1/3,1,3,4,6 1/2,9,12 1/2, 16,20,25,34,44,51,62,89, 116 hrs. D 6 0.84 0.039 489 0,5,8,11,15,24,36,50,69, 82 min D 5 0.86 0.092 500 0,1/2,1 1/2,3,6,12 min B 8 0.88 0.065 500 0,1 1/2,6,24,48,72,96,144 mir B 6 0.86 0.059 489 15,30,45,60,80,100,156,212, 297 min D 5 0.89 0.116 534 0,1/2,1,2,3,5 min Al-Ag 4 0.88 0.328 500 5,17,38,60,90 min B 5 0.90 0.078 513 1 1/2,6,12,18,24,48 min Al-Ag 6 0.92 0.597 534 1/2,1 1/2,3,5,7 1/2 min A 5 0.95 0.048 520 1,4,16,66 min A 9 0.97 0.052 534 5,10,20,30,50 min B 5 0.92 0.09 534 3,6,10,15,20 min TH = 0.92 Alloy-B - 105 -of experimental data. 29 Whelan's planar model (see Section 1.3.1) differs from Aaron's planar model"^ by a factor of /n"/2 in the constant k of equation (7) in Section 1.3.1. A plot of R versus /t for the experimental results should give a linear relationship if this model is valid. The experimental results for some precipitates give a linear relationship for the total amount of dissolution observed^with a significant amount of scatter. A linear fit is observed in the case of a few precipitates for only the earlier dissolution times (see Figure 47). However, no conclusion can be drawn as to the applicability of this model to the early stages (that is before impingement),as the majority of the observations were made after a certain amount of impingement had occurred. In an attempt to find the best correlation between R and t, the data was plotted In the form log R/RQ versus t and, to a good approximation, a linear relationship was obtained. Figures 48 to 52 are typical examples of log R/RQ versus t and log S/Sq versus t plots for the homologous temperature range from T„ = 0.97 to T„ = 0.76. The observed linear relationship indicates H H that a relationship of the form -k t R = R e 1 o is applicable to the change in half thickness and -k2t S •• » S e o 2 - 106 -0 10 20 30 40 50 Time(min) Figure 48. Log [S/S ] versus t and log [R/R ] versus t for alloy A at T = 0.97. u O H - 107 -- 108 -- 109 -- 110 -0 200 400 600 800 1000 Time(min) Figure 52. Log [S/S ] versus t and log [R/R ] versus t for alloy D at 1 = 0.76. - Ill -is valid for the change in half length. In the following section, the observed exponential nature of the dissolution kinetics is rationalized on the basis of an analytical model which accounts for diffusion field impingement. The model assumes planar symmetry. 4.3 An Exponential Model for Planar Dissolution Kinetics It was shown in section 3.4.1 that for the present experimental conditions, impingement of the diffusion fields occurs early in the dissolution process, before there is any significant movement of the interphase boundary. In view of this fact, it was decided to use an approximate solution to the diffusion equation for a finite system, assuming that the movement of the phase boundary can be neglected (i.e., the stationary interface approach is adopted). The initial situation for this geometrical model is shown in Figure 53(a) i Figure 53(a) Geometrical model of a particle with half-thickness R in a matrix with interparticle spacing £. where R^ is the initial particle half thickness and H is the interparticle distance. It is assumed throughout this derivation that RQ << i/2. The approrpriate phase diagram is shown in Figure 53(b). - 112 -Cq CQ Ci Figure 53(b) An appropriate phase diagram with growth of precipitates taking place at.l^ and dissolution at 1^. is the annealing temperature for growth and is tne dissolution temperature. From the lever rule it follows that R 1/2 C -C, o 3 C -C P o since the sample is assumed to be completely equilibrated at temperature (i.e., prior to upquenching to the dissolution temperature T^). Assume that the dissolution process has commenced at 1^. Prior to impingement^ the concentration distributions are as shown in Figure 53(c). Figure 53(c) Dissolution profile prior to impingement. - 113 -During the impingement period, the distributions are as shown in Figure 54(a). Ultimately the alloy will have a uniform composition C . Figure 54(a) The concentration distributions during dissolutions after impingement takes place. Hence the concentration at the centre of the impingement region (&/2) will gradually approach C^, but will attain it only after the second phase has completely disappeared and the concentration gradients in the matrix have been eliminated. There exists no analytical solution to the diffusion equation for the impingement situation. Therefore an approximate solution is used in conjunction with the assumption that there is a constant concentration Figure 54(b) Schematic diagram of model used with a constant concentration C-*3 £L fc X —~ & • - 114 -at x = H, which over a reasonably long time period yields concentra tion contours near the interface which approximate those due to impingement as shown in Figure 54(b). A reasonable value for this concentration is C0 = 2C - C 2 o I as this yields a concentration of Cq (the bulk composition) at i/2 for infinite time. The present model is deficient in that it predicts a finite gradient at ft/2 for all times, whereas there should be no flux across this point. However, it is emphasized that we are concerned with the gradients at x = 0 (the interface) and, therefore, the approximation should involve little error. The solution to the diffusion equation for a plate of thickness £, faces of constant concentration C_ and C_ = 2C -C_, and an initial concentration C„ is 1 2 o I 3 given in Crank"^ (page 47 ). This solution involves a series expansion. For the present problem it is assumed that the leading term in this series dominates. This assumption requires that 2 2 _ n DTT t 2 e % « 1 (n = 2,3 ) 2 Under these conditions the solution is „ 2 _ DTT t C = CT + 2(C -CT)x/£, - - (C -Cjsin ^ e ft2 3 I O I TT O 3 I On differentiating and setting x = 0, one obtains - 115 -9£ 9x = 2 <Co-V (C -C„) -/ ° 3 4 e DTT2t x=0 This gradient is substituted into the interface mass balance (C -C ) ^ % V dt 3x x=0 to give dR dt pr2 -<w * (cQ-c3) £ L (C -CT) + (c-cT) J e £ p "I' p -I' The ratio (C -C )/(C -C_) is f ,, the supersaturation for dissolution I o p I d' (Figure 53(b)). In the second term the ratio C -C- C -C_ 0 3^ o 3 C -CT P I C -C p o since >> C^, Cq (i.e., for low supersaturation) In view of Equation (1), it follows that C -CL R o 3 ^ o Vci 1,2 Using Equation (8) and the definition of f^, Equation 6 becomes - 116 -This is integrated from t = 0 (R = R ) to t(R), to give 2Df,t 8R RQ-R = -f ~f (e ll - 1) 0 TT 2 Assuming 8/IT =1, Equation 10 becomes R = R e ft - 2 —— 11 o ft If one substitutes values of the parameters D, t, ft and f^ which are typical of the present experimental conditions, it is found that the 2Dfdt linear term —-— only becomes important in the later stages of dissolution (see Figure 55). Thus for a considerable (intermediate) portion of the dissolution period the kinetics should follow the equation R = R e ft2 o 2 i.e. ftn R/R = - ]L~ 12 ft2 which is in agreement with the experimental results (Figures 48-52). An attempt was made to generalize the above treatment to spherical symmetry. Unfortunately, in order to obtain an equation of the same form as Equation (12),so many assumptions were required that the derivation was lacking in analytical rigor and was therefore deemed unacceptable. It was felt that it would be worthwhile to check the validity of the exponential relationship with dissolution kinetics obtained by r Time(sec) Figure 55. R versus t for equations 11 and 12. - 118 -52 Baty et al. They used a statistical method to follow the dissolution of C11AI2 precipitates in an Al-4wt.% Cu system at 520°C. Using their data, a plot of log R/RQ versus t is shown in Figure 56 and a reasonably linear fit is obtained. Further support for the exponential model is shown in Figure 57 where the numerically calculated solution 25 by Tanzilli and Heckel for dissolution of planar precipitates in a finite system is plotted in logarithmic form. The supersaturation (k = 0.22) is much greater in this instance than for most of the present experimental work. Equation 12 was derived using a stationary interface assumption, the reliability of which increases with decreasing supersaturation. That a reasonably linear fit was obtained from the numerically derived solution which involves a relatively high supersaturation is very encouraging. The experimental results obey the exponential relationship as established in Equation 12, for the whole range of homologous temperatures observed. The greatest deviation from linearity in the plots of £n R/RQ versus t is in most cases found during the early stages of dissolution, before impingement. This is to be expected as the model does not account for dissolution kinetics prior to impingement. The values of calculated using equation 12 are consistently higher than the literature values as may be seen in Figure 58. All other things being equal, the dissolution rate of a curved surface is faster than that of a planar surface. The surfaces of the allotriomorphs are in fact curved. Thus the application of a planar dissolution model to the dissolution of a curved surface will of necessity yield values which are higher than the actual values. Furthermore the ommision of the Time(sec) Figure 56. Log [R/RQ] versus t for dissolution of precipitates in an Al-4 wt.% Cu system at 520°C. 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Const. Xt Figure 57. Log [R/R0] versus t for the dissolution of planar precipitate in a finite system with k = 0.22. - 121 -- 122 -Dfdt linear term (2 ——) in the final form of the exponential equation I (equation 12) will also result in higher calculated Dv values. Finally,any surface diffusion effect present will increase the observed Dy value. It should be noted, however, that a surface diffusion effect if present would become more important at lower temperatures tending , . , r_, calculated lT. literature, ^. ^ •> , to give higher [D^ /D^ J ratios at low homologous temperatures. If anything the results show the opposite trend. 4.4 Precipitate Shape Change During Dissolution In general it was observed that a constant axial ratio S/R was not maintained during dissolution. Indeed, in certain cases the final axial ratio, S /R , was significantly different from the initial axial r r ratio, S /R . The present section deals with this observed change in o o r & axial ratio. Figures 59-63 are examples of situations in which VRo - VV VRo > VRF AND VRo < VR! It was shown in Section 3.4 that the shape of the concentration profiles during dissolution and hence the mode of dissolution, for a certain homologous temperature is independent of supersaturation. Figure 64 is a plot of supersaturation versus the fraction (F) of precipitates in a given experiment showing an increase in axial ratio of greater than 10%, i.e., the fraction of precipitates for which (VV > 1.1 o o This procedure is crude for several reasons. - 123 -- 124 -Figure 60. Back scattered electron image micrographs,used to observe the dissolution kinetics of an allotriomorph at = 0.97 for alloy A. - 125 -Figure 61. Back scattered electron image microgranhs,used to observe the dissolution kinetics of an allotriomorph at T = 0.88 for ri the Al-Ag alloy. - 126 -Figure 62. Back scattered electron image micrographs, used to o the dissolution kinetics of an allotriomorph at T}{ = for alloy B. - 127 -Figure 63. Back scattered electron image micrographs,used to observe the dissolution kinetics of an allotriomorph at = 0.84 for the Al-Ag alloy. - 128 -1. The use of 10% as the fiducial axial ratio increase is admittedly arbitrary. 2. The comparison is made on the basis of the last axial ratio measurement taken. However, there may have been further changes in axial ratio before the precipitate completely dissolves. This however could not be observed due to the resolution of the technique used. 3. The number of precipitates studied in any experiment was typi cally 3-10, which is insufficient for unquestionably meaningful statistical studies. In spite of these deficiences Figure 64 indicates that at a given T , the change in axial ratio is independent n of supersaturation which is in agreement with the concentration contour studies. In Figure 65 the same function (F), as in Figure 64, is plotted versus the homologous temperature. This plot shows a definite increase in F with increasing homologous temperature. At high homologous temperatures, TR > 0.95, all the precipitates studied show an increase in their axial ratio during dissolution (see Figures 59 and 60). About half of the precipitates examined at 0.95 > T > 0.80 show an ri increase in axial ratio and at low homologous temperatures T < 0.80 ri the precipitates either show no change or else a slight decrease in axial ratio (see Figure 63). These observations can be explained in the following manner. 1.00 0.75 aso a25 Alloy TH Alloy TH Alloy TH Alloy TH • A 0.9/ • B 0.90 • D 0.86' 0 C 0-79 O A 0.95 A • D 0.89 B V B A Al-Ag 0.86 0 D 0-76. 0 B 0.92 A • B 0-88 0.84 C © Al-Ag 0.92. E Al-Ag 0.88. * A B 0.84 A D 0.84. 0 • • Q L 0.02 0.2 D a 0.3 0.4 0.5 0.03 0.04 0.05 0.1 Supersaturation(k) Figure 64. Fraction of precipitates showing an increase in axial ratio greater than 10% (F) versus supersaturation (k). * Dissolution experiment below the solvus, resulted in spheriodi-sation (see Figure 70). A Al-Cu Alloys O Al-Ag Alloys - 131 -A.4.1 High Homologous Temperatures At high homologous temperatures (T > 0.92), the contribution of grain boundary diffusion to the dissolution process may be neglected and the precipitate may be considered to be an oblate spheroid dissolving by volume diffusion into the matrix. Figure 66(a) is a schematic diagram of the distribution of precipitates with an axial ratio of approximately 4 to 1 and an interparticle spacing of L = 10 x R . These values are consistent with most of the experimental work. Figure 66(b) gives a schematic representation of the concentration fields around these precipitates during the initial stages of dissolution. Because only bulk diffusion contributes, these contours are confocal spheroids. In view of the fact that a particular precipitate particle has a symmetrical distribution of other particles about it, it is reasonable to assume that the diffusion field from surrounding precipitates has approximately spherical symmetry as shown in Figure 66(c). o 28 p 0 0 0 0 0 Figure 66(a) Symmetrical distribution of precipitates with [S /R ] = 4. - 132 -I <§ c° (in /®\ nt Figure 66(b) Schematic representation of concentration fields around a symmetrical array, of precipitates. Figure 66(c) Diagram of spherically symmetrical field surrounding a grain boundary allotriomorph. However, even if one assumes a random distribution of precipitates about any one precipitate, the diffusion field from these latter precipitates is on the average radially symmetrical with respect to the precipitate in question. Accordingly, as a crude approximation the concentration contours arising from neighbouring precipitates can be "smeared" into confocal spherical contours with respect to the precipitate in question. Figure 66(c) shows clearly that impingement first - 133 -takes place at the tip rather than the faces of the allotriomorph. The tips are therefore expected to dissolve more slowly resulting in an increase in S/R. Because the precipitates in question are grain boundary allotrio morphs, the assumption of a random distribution of precipitates about a particular precipitate is open to question. However, any asymmetry would be due to the presence of a number of precipitates on the same grain boundary being in closer proximity to the precipitate under observation than the Widmansta'tten plates in the adjacent matrix. Such a distribution will further enhance the likelihood of impingement first taking place in the equatorial plane (i.e., the plane of the grain boundary). The impingement profile is assumed to be a superposition of the confocal spheroids of the central precipitate and the surrounding spherical field (see Appendix VII for details). Figure 67 gives the calculated profiles after 500 seconds. Significant interaction of the diffusion fields in the equatorial plane is evident. In Figure 68(a) it can be seen that on section Ob (from Figure 67), which is along the minor axis, the diffusion field adjacent to the precipitate is hardly affected by diffusion inward from the spherical field (i.e., from the adjacent precipitates). In Figure 68(b), however, the magnitude of the gradients in section Oa, which is the equatorial plane, are reduced significantly by diffusion inward from the spherical surface. The concentration contours will not actually be confocal spheroids but will have somewhat higher axial ratio. This will result in impingement along the equator at earlier times and will tend to accentuate the above effect. - 134 -t= 500 Seconds D=SxlO"10cm2sec1 ®-CA Figure 67. The resulting impingement between the confocal spheroid isoconcentra tion contours about an allotriomorph and the contours of the surrounding spherical field after t = 500 seconds for various [C/Cj]. - 137 -Thus if during the early stages of impingement, the low curvature surfaces are dissolving at a rate such as to maintain the axial ratio, then because of impingement,the high curvature equatorial region is dissolving at a rate slower than that required to maintain the axial ratio. The net result is therefore an increase in axial ratio. This effect will continue even after significant impingement. 4.4.2 Low Homologous Temperatures At low homologous temperatures,the grain boundary diffusion coefficient is so much higher than the volume diffusion coefficient, that the grain boundary is effectively supersaturated with solute before any significant dissolution takes place, (see Figure 30, for concentration profiles at low T ). This results in a solute concentration along the grain boundary which is equal to the interface solute concentration C^.. Thus the tip of the allotriomorph tends to dissolve preferentially, with the corresponding solute moving along the grain boundary and dispersing from here into the matrix. This model for dissolution suggests a spheroidizing process with the tip dissolving much faster than the low curvature faces of the precipitate. However this 53 tendency is probably reduced by interfacial diffusion along the precipitate-matrix interface which would tend to feed solute to the grain boundary from other areas over the precipitate surface. A slight spheroidizing tendency is observed at low homologous in the present work. - 138 -4.4.3 Intermediate Homologous Temperatures At intermediate homologous temperatures the tendencies applicable to both high and low homologous temperatures are probably operative. The relative importance of these tendencies is determined, among other things, by: (1) The interparticle spacing on the grain boundary: For instance, it was qualitatively observed that in cases where a grain boundary precipitate was situated between two other grain boundary precipitates, the rate of change in half length S was reduced significantly. (2) The initial axial ratio (SQ/R ): This ratio is important in determining the degree to which earlier impingement of diffusion fields takes place at the tip as opposed to the sides of the allotriomorph. The model for high homologous temperatures indicates that the higher the initial axial ratio,(S /R ), is,the earlier impingement of ' o o diffusion fields at the tip takes place, and therefore the greater is the expected increase in axial ratio. is plotted against (SO/RQ). There is a definite trend towards a larger change in (S/R) with higher initial axial ratios, a result which is consistent with the preceding discussion. experiment at a temperature below the solvus line for an intermediate homologous temperature (T = 0.83 for alloy B)) in Figure 70(a), shows a good fit for the early stages of dissolution. A plot of [S/R] versus t in Figure 70(b) shows a slight overall decrease in [S/R] for the In Figure 69 the percentage change in axial ratio ( F F/S /R ) o o A plot of log R/R versus t, (for a dissolution o corresponding early stages of dissolution. All four precipitates 7o Increase in SEZEEI So/R< KJ O o - 6CT " - 140 -9 • 7 5 • • # • 5° 3 -2 25 50 t.hrs.75 ,0° 125 3.0 1 • • • •• • 2.5 -LOO£ • • 2.0 • • • 1.6 • Figure 70. (a) Log [R/R ] versus t. (b) [S/R] versus t. - 141 -in this experiment show a small change in [S/R] for the initial stages of dissolution. This is in good agreement with the predicted change in [S/R] with homologous temperature shown in Figure 65. It should be noted that the point of departure of the results from the exponential relationship in Figure 70(a) coincides reasonably well with the point in Figure 70(b) where rapid spheroidization begins to take place. - 142 -5. CONCLUSIONS The results of the present study of the dissolution behaviour of grain boundary allotriomorphs in the Al-Cu and Al-Ag systems yield the following conclusions. 1. Isoconcentration contours around dissolving grain boundary allotriomorphs can be accurately determined by the use of electron probe microanalysis. 2. For allotriomorphs in an effectively infinite matrix, the shape of the isoconcentration contours is dependent only on the homologous temperature, T . H 3. At a particular T , the relative contribution of grain boundary diffusion to volume diffusion can be obtained directly from the isoconcentration contours. 4. The dissolution of grain boundary allotriomorphs is dominated by volume diffusion above Tu = 0.92 and by grain boundary diffusion rl below T = 0.72. For the intermediate range Tu = 0.92-0.72, the grain n H boundary diffusion contribution increases continuously as T decreases.. H 5. Enhanced dissolution due to grain boundary diffusion is limited to precipitates on high angle grain boundaries. Grain boundary allotric-morphs are therefore the only precipitate morphology that will exhibit this phenomenon. 6. The dissolution kinetics of an individual grain boundary allotrio morph can be established from its back scattered electron image on the electron probe microanalyzer. 7. The dissolution kinetics of grain boundary allotriomorphs, under conditions in which impingement of diffusion fields from adjacent - 143 -precipitates takes place early in the dissolution process, are found to satisfy the following equations: _ TT2Dt R - R e £2 o for the change in half thickness and _ v Dt S = S e I1 o for the change in half length. 8. At high homologous temperatures, where volume diffusion dominates, an increase in the axial ratio of the allotriomorphs is observed as dissolution proceeds. On the other hand, at low homologous temperatures where grain boundary diffusion dominates, a slight decrease in axial ratio (i.e., slight spheroidization) is observed. A continuous change in behaviour from one extreme to the other is observed over the range of intermediate homologous temperatures. - 144 -BIBLIOGRAPHY 1. Dube, C.A., Aaronson, H.I. and Mehl, R.F., Rev. Met., 55_, 201 (1958). 2. Aaronson, H.I., The Mechanism of Phase Transformations in Metals, Inst, of Metals Monograph 18, 47 (1956). 3. Aaronson, H.I., Decomposition of Austenite by Diffusional Processes, Interscience Pub. N.Y. (1962), p. 387. 4. Clark, J.B., in High-Temperature, High-Resolution Metallography, Gordon and Beach Pub., N.Y., (1967). 5. Toney, S., and Aaronson, H.I., Trans. A.I.M.E., 221, 909 (1961). 6. Hawbolt, E.B., and Brown, L.C., Trans. A.I.M.E., 239, 1916 (1967). . 7. Read, W.T., Jr., and Shockley, W. , Phys. Rev., 78, 275 (1950). 8. Vanghan, D., Acta. Met., 16, 563 (1968). 9. Ryder, P.L. and Pitsch, W., Acta. Met., 14, 1437 (1966). 10.. Ryder, P.L., Pitsch, W., and Mehl, R.F., Acta. Met., 15, 1431 (1967). 11. Zener, C., J. Appl. Phys., 20, 950 (1949). 12. Frank, F.C., Proc. Roy. Soc, A201, 586 (1950). 13. Ham, F.S., J. Phys. Chem. 335 (1959a). 14. Ham, F.S., J. Appl. Phys. 30, 915 (1959c). 15. Horvay, G., and Cahn, J.W., Acta. Met., 9_, 695 (1961). 16. Mazanec, K., and Cadek, J., Rev. Met., 212, 501 (1958). 17. Hickley, CM., and Woodhead, J.H., J. Iron and Steel. Inst. (London), 176, 129 (1954). 18. Aaron, H.B., and Aaronson, H.I., Acta. Met., 16, 789 (1968). 19. Hillert, M., Jernkontovets Ann., 141, 757 (1957). - 145 -20. Brailsford, A.D., and Aaron, H.B., J. Appl. Phys. 40, 1702 (1969). 21. Goldman, J., Aaronson, H.I., and Aaron, H.B., Met. Trans., JL, 1805 (1970). 22. Nolfi, Jr., F.V., Shewmon, P.G., and Foster, J.S., Trans. A.I.M.E., 245, 1427 (1969). 23. Nolfi, Jr., F.V., Shewmon, P.G., and Foster, J.S., Met. Trans. 1, 789 (1970). 24. Nolfi, Jr., F.V. , Shewmon, P.G., and Foster, J.S., Met. Trans. 1_, 2291 (1970). 25. Tanzilli, R.A., and Heckel, R.W., Trans. A.I.M.E., 242, 2313 (1968). 26. Tanzilli, R.A., and Heckel, R.W., Trans. A.I.M.E., 245, 1363 (1969). 27. Aaron, H.B., Fainstein, D., and Kotler, G.R., Publication Preprint Ford Motor Co. Scientific Laboratory, April, 1970. 28. Thomas, G., and Whelan, M.J., Phil. Mag., i5, 1103 (1961). 29.. whelan, M.J., Metal. Sci. Int. 3, 95 (1969). 30. Readey, D.W., and Cooper, Jr., A.R., Chem. Eng. Science,21, 917 (1966). 31. Hall, M.G., and Haworth, CW. , Acta. Met., 18, 331 (1970). 32. Eifert, J.R., Chatfield, D.A., Powell, G.W., and Spretnack, J.W., Trans. A.I.M.E., 242, 66 (1968). 33. Tanzilli, R.A. and Heckel, R.W. Publication Preprint, Drexel Institute of Technology, Philadelphia, Pa. U.S.A. 34. Aaron, H.B., Acta. Met., JL7, 407 (1969). 35. Aaron, H;.B., and Kotler, G.R., Met. Trans. _2, 393 (1971). 36. Hansen, M., Constitution of Binary Alloys, 2nd ed., McGraw-Hill Book Co., N.Y. (1958), p. 2. 37. Metals Handbook, A.S.M., Cleveland, Ohio (1948), p. 1160. - 146 -38. Brown, L.C., Department of Met., U.B.C., Private Communication. 39. Murphy, J.B., Acta. Met., j), 563 (1961). 40. Embury, J.D., and Nicholson, R.B., Acta. Met., JL3, 403 (1965). 41. Pashley, D.W., Rhodes, J.W., and Sendorek, A., J. Inst, of Metals, 94, 41 (1966). 42. Aaron, H.B., and Kotler, G.R., Publication Preprint, Ford Motor Co., Scientific Laboratory (1970). 43. Achter, M.R., and Smoluchowski, R., Journal of Applied Physics, 22, 1260 (1951). 44. Smoluchowski, R., Physical Review, 8J3, 482 (1952). 45. Shewmon, P.G., J. Appl. Phys., 34, 755 (1963). 46. Fisher, J.C., J. Appl. Phys. 22, 74 (1951). 47. Turnbull, D., and Hoffman, R.E., Acta. Met., 2, 419 (1954). 48. Whipple, R.T.P., Phil. Mag., 45, 1225 (1954). 49. Heumann, Th., and Bohmer, H., J. Phys. Chem. Solids, 29, 237 (1968). 50. Aaron, H.B., Metal. Sci. 2, 192 (1968). 51. Crank, J., The Mathematics of Diffusion, The University Press, Oxford, p. 47 (1957). 52. Baty, D.L., Tanzilli, R.A., and Heckel, R.W., Met. Trans., 1, 1651 (1970). 53. Weinberg, F., Dept. of Met., U.B.C., Private Communication. 54. Reed, S.J.B., and Long, R.V.P., Proc. 3rd. Int. Symp. on X-Ray Optics and Micro-Analysis, Stanford Academic Press, N.Y. (1962). - 147 -55. Heinrich, K.F.J., Proc. of the Symp. on the Electron Microprobe, Washington, John Wiley, New York and London, 296 (1962). 56. Brown, L.C and Thresh, H., Tools and Techniques in Physical Metallurgy, M.Dekker, N.Y. (1970). 57. Duncumb, P., and Reed, S.J.B., Technical Report No. 221, Tube Investments Research Laboratories, Hinxton Hall, Cambridge (May 1967). 58. Bethe, H.A. , Ann. Phys. Lpz. , 5_, 325 (1930). 59. Belk, J.A., 4th Int. Congress on X-Ray Optics and Micro-Analysis Hermann Paris, Sept. 1965. 60. Poole, D.M., Proceedings of Seminar on Quantitative Electron Probe Microanalysis, Nat, Bur. of Stand. Special Pub. 298, Washington, p. 93 (1968). 61. Philibert, J., Proc. 3rd. Int. Symp. on X-Ray Optics and Micro-Analysis, Academic Press, N.Y., Stanford (1962). 62. Duncumb, P., and Shields, P.K., The Electron Microprobe, Wiley, N.Y. (1966). 63. Castaing, R., Advan. Electron. Electron Phys. , 13, 317 (1960). 64. Murdoch, D.C, Analytic Geometry with an Introduction -to Vectors and Matrices, John Wiley and Sons, Inc., N.Y. (1966) 65. Coates, D.E., Dept. of Met., U.B.C, Private Communication. 66. Hawbolt, E.B., Ph.D. Thesis, Dept. of Met., U.B.C. (1967). 67. Brown, L.C, Dept. of Met., U.B.C, unpublished work (1964). 68. Crank, J., Mathematics of Diffusion, The University Press, Oxford, p. 86 (1957). 69. Ham, F.S., Quart. Appl. Math. 17, 137 (1959). - 148 -APPENDIX I. Correction Procedures for Quantitative Electron Probe Microanalysis The measured X-ray intensity of element A in alloy AB may be converted to the true mass concentration (C^) of element A in alloy AB, by applying the following corrections in the order indicated: 1. dead time; 2. background; 3. secondary fluorescent enhancement by characteristic radiation and white radiation; 4. X-ray mass absorption in the sample; 5. the atomic number effect, which is made up of two separate components: (a) the electron stopping power (S); and (b) the electron backscattering effect (R). 1.1 The Al-Cu System Pure Cu was used as a reference standard for CuKa^ radiation. The CuKa^ peak was found at a Bragg angle of 26°34* + 3' using a quartz analyzing crystal to reflect the X-rays to the gas flow proportional counter. The background count was picked up on the second spectrometer at an angle of 27°4'. To convert the measured intensity of CuKa^ radiation to the concentration of Cu in the Al-Cu alloy only the deadtime, background, and atomic number corrections need be applied to the experimentally determined ratio of intensities (K ) where: - 149 -Kc^ Cu Ka1 I n = the Ka., radiation intensity of Cu in an Al-Cu alloy, Ka^ and I = the Ka^ radiation intensity of Cu in pure Cu. A negligible secondary fluorescent enhancement is expected in this 54 system and this was verified by applying the Reed and Long correction for fluorescent excitation. The absorption for this system is also negligible as"'"' Cu Cu yAl = ^Cu where u^j* = 49.6 = (mass absorption coefficient of CuKa^ radiation by Al) and = 53.7 = (mass absorption coefficient of CuKa^ radiation by Cu). (1) The dead time,d, may be calculated from the following 56 relationship ^(bbs.) '(true) 1-C, , xxd (obs.) C. =true counts after the dead time, correction is taken into (true) consideration, C, , . = observed counts, d = dead time = .4 ysec. (obs.) ' for the electron probe used. This simple correction need only be - 150 -applied for high count rates (eg.>2,000 c.p.s.). (2) The background count per ten seconds is registered for each ten second CuKa^ count. The background count is subtracted directly from the CuKa^ count, thereby eliminating the effect of background radiation. (3) The Duncumb and Read"'7 atomic number correction was used to convert the intensity ratio of CuKa^ X-rays produced (KQ ) TO the true mass concentration of Cu (C„ ). This correction has the form: Cu K = C * Ra1~Cu Scu Al-Cu where R = backscatter coefficient; S = electron stopping power; Al-Cu refers to the alloy specimen and Cu refers to the standard. 58 The stopping power is evaluated using the following equation: _ " _ . Z 1 „. (1.66 • E) S = constant * — • — • m — A E where Z = atomic number; A = atomic weight; J = mean ionization potential (tabulated values in reference 57); v E = electron energy ; and the backscatter coefficient, R, is given by: 1 n(W).Q/S.dW FW K 1 - R = E K - 151 -where: Eq = incident electron energy; J E = critical K excitation potential; Q = ionization cross-section; • W = E/E ; S as in equation 2. n (W) .=. backscattered electron energy distribution expressed in integral form. Values of R are tabulated by Duncumb and Reid"*7 as a function of 1/U (U = V/Vv) and Z where: V = accelerating potential; V^ = excitation potential of Ka, characteristic K i radiation. Figure 1-1 is a plot of versus . for the Al-Cu system, after making the corrections described above. 59 The Belk atomic number correction was applied to the Al-Cu system as a check on the validity of the Duncumb and Reed correction used and was found to be in good agreement with it. This correction procedure is very simple to apply and has the form: K - "A A Z -Z 1+ <mr> „ -7 ZACA+ZBCB where Z = CA+CB and Z = atomic number. fin The above has been found by D. M. Poole to be as accurate as any of the more sophisticated correction procedures. - 152 -- 153 -1.2 The Al-Ag System For this system the background, the dead time, the absorption and the atomic number correction, S, were used to convert measured La^ X-ray intensity to the concentration of Ag in the alloy (C. ). The Ag dead time, background and atomic number correction were applied in the same manner as for the Al-Cu system. 61 The following Philibert absorption correction was used F(X) " f(x)- ?(0) = (1 +'x/a) tl + h(l + x-/a>] where f (x) = the fractional transmission after absorption and represents the absorption effect only. F(0) = the atomic number effect only, h = 1.2 A/Z2; A =.atomic weight; X = p/p (sec e); u/p = mass absorption coefficient; 9 = take of angle on the probe = 20°; 0 = modified Leonard's coefficient as calculated by Duncumb 62 and Shields and is given by 2.39 x 105 E1-5 - EJ'5 o K Ag 55 vll = 779.2 pfS = 521.9 Ag - 154 -Figure 1-2 gives a calibration curve of C. versus K. for the Al-Ag Ag Ag 6 system, after applying the corrections described above. - 155 -- 156 -APPENDIX II. Electron Spot Size Determination. The electron spot size for the absorbed, backscattered and topo graphic electron image on the JEOLCO JXA 300 is approximately 0.7 u. This corresponds to the resolution for these images. The X-ray spot 6 3 size however varies according to the relationship: spot size S = 0.033 (E1,7 - E^,7)A/PZ o K K where Eq = operating voltage = 25 K.V. E = excitation potential for the radiation considered; and P = specimen density. The relevent data is given in Table II-l. Table II-l. E , A, Z and p Values for Ag La.. , Al Ka, and Cu Ka.. Radiation. E__ (K.V.) A Z p g cm Ag La^ 4.0 107.868 47 10.49 Al Ka-^ 1.6 26.98 13 2.699 Cu Ka-L 9.0 63.546 29 8.96 A schematic diagram of the electron distribution in the sample is shown in Figure II-l. The total spot size is the sum of the electron beam diameter and the calculated spot size (S). On the basis of the above, the total spot size for the aluminium-rich matrix is - 5.5 u, - 157 -Figure II-l. Schematic diagram of the electron spot size. - 158 -taking into account the fact that the density of the ct-solid solution is slightly larger than that for pure aluminum. The spot size for C11AI2 is approximately 3 u , the spot size in the Ag rich matrix of the Ag-Al system is approximately 5 y, and the spot size in Ag^Al -2.25 u. - 159 -APPENDIX III. Calculation of the Flux Line Divergence for Different Sections of a Two Dimensional Elliptical Geometry. The general equation for an ellipsoid as in Figure Ill-l(a) is given by: 2 2 2 (a) (b) (c) Figure III-l. Two dimensional sections of an oblate spheroid. If a = d > b one obtains an oblate spheroid, whereas is a = b < d one obtains a prolate spheroid. Assume a grain boundary allotriomorph can be approximated by an oblate spheroid and assume the grain boundary lies on the y = 0 plane (i.e. x-z plane),(Figure Ill-l(b)). A cross section perpendicular to the grain boundary, the z=0 plane (i.e. x-y plane),is shown in Figure III-l(c) which is an ellipse of equation 2 2 1 If c is the focal length of the ellipse, the following relationships , 64 exist - 160 -e = c/a = eccentricity, and If c -+ a, e = 1, and b -* 0, therefore the spheroid reduces to a disc. If c + 0, e = 0, and a = b, therefore the spheroid reduces to a sphere. In the Al-Cu system the most common allotriomorphs present have an axial ratio of 3.5. With an axial ratio, a/b = 3.5, Eq (1) is used to determine the angle between the tangent to the ellipse at x = 0, y = b, and at arbitrary values of x and y. Results are given in Table III-l. Table III-l A0° for Different Values of x/b. x/b y/b A6° 0 1 0 0.875 0.97 4 1.75 0.867 10 2.625 0.658 20 3.20 0.407 39 - 161 -APPENDIX IV. Dv Calculation from the Diffusion Couple Results As shown in Figure 16, a probability plot of C^ ^ versus x within the supersaturated a phase gives a straight line for the diffusion couple annealed at 500°C. Consequently a volume diffusion coefficient Dy, which is independent of concentration, is indicated. All the diffusion couple results when plotted on probability paper gave diffusion coefficients independent of concentration. In Figure IV-1, the concentration profile at 500°C is used to evaluate D^ from the equation cv - IF where Z is the point at which the tangent to the concentration gradient at the interface cuts the concentration line C=0. This equation gives accurate values of D^ when D^ £ f (concentration), and the interface is essentially stationary. From Figure IV-1 _ _ (2.80 x 10"2)2 -10 2, D - -r-r = 4.26 x 10 cm /sec V TT x 163 x 3,600 Dy values were also determined in an alternate manner using the probability plots for each of the four diffusion couples. Using a standard diffusion coefficient, D = 5 x 10 ^" cm2sec \ and the appropriate value of time, a standard line was calculated and drawn on each probability plot. Then D^ could be calculated from the relationship: Distance in (fJ) Figure IV-1. Composition profile of diffusion couple at 500°C. (?'2 xv where Xg and x^. are the diffusion distances in the aluminum corresponding to the same copper concentration for the standard and experimental lines respectively. The value of D^. calculated in this way at 500°C differs only in the second decimal place from the D^ value calculated using equation 1. Table IV-1. Dv Values of the Diffusion Couple Results. x_(50-l%) = 560 y x 2 = 31.36 x 104 y2 - 163 -D.. T°C T°K 1/T°K-Ixl04 vy (50-1%) 2 in"4 2 xv xlO y Dvxl0 480 753 13.28 368 13.54 2.16 500 773 12.94 514 26.40 4.21 520 793 12.60 700 49.00 7.80 535 808 12.38 745 55.50 8.84 - 164 -APPENDIX V. Calculation of Standard Deviation for the Isoconcentration Contours. The mean (X) of a set of data X^ is given by: Mean = X n = E X. i=l 1 n The standard deviation (a) is given by: Standard deviation = a = - 2 E (X.-X)Z 1=1 1 n-1 Table V-l gives the mean (X) and standard deviation (a) for the dissolution results of 5 precipitates in alloy B at T = 0.91 shown H in Figure 19. Table V-l. Diffusion Distance, X, a and "=f Precipitate number At centre 90 sec At centre 360 sec 2.1% 2.5% 2.1% 2.5% 1 10.62 6.87 20 13.2 2 9.37 6.87 17.5 15.0 1 10.00 5.63 17.5 10.0 2 8.75 6.25 13.5 1 10.00 6.87 17.5 11.25 2 8.75 5.0 16.25 11.25 1 8.75 5.63 17.5 11.25 1 10.00 6.25 11.25 8.75 5.63 16.87 11.25 Table continued Table V-l Continued - 165 -2.1% 2.5% 2.1% _2-5% X 9.44 6.11 17.59 11.88 o 0.73 0.73 1.17 1.47 a/X 0.08 0.11 0.066 0.11 - 166 -APPENDIX VI. Calculation of \p Coates^^ has attempted to calculate, as a function of temperature, the relative contributions of volume and grain boundary diffusion to dissolution of grain boundary allotriomorphs. The following is an outline of his calcu-tion. Consider the grain boundary to lie on the x-y plane. The allotriomorph is approximated by an oblate spheroid whose major axis is in the x-y plane (i.e., an ellipse in the x-z plane is rotated about its minor axis, the z axis). The foci are at x = + c and the semi-major and semi-minor axes are a and b respectively (i.e., the spheroid intersects the x, y and z axes at + a, + a and + b respectively). The first step in the calculation is to determine the distribution of solute (the diffusion field) about the precipitate assuming only volume diffusion. Provided the supersaturation is low, the diffusion 2 field can be approximated by a solution of Laplace's equation V C=0. Physically this is equivalent to assuming a steady state distribution of solute. Using this approximation the diffusion field about the allotriomorph is determined. Then the interface gradient of this field V^C is computed to give in turn the interface flux J^. = -D^V^C. Finally the total interface outflow of material from the allotriomorph, J« dA; is calculated to give 1 where CT and C., are the interface and bulk matrix concentrations I M respectively,and a/b is the axial ratio of the spheroid. - 167 -The next principal step is the calculation of the diffusion field about the allotriomorph assuming the material leaves via the grain boundary only. Thus one is effectively considering diffusion from a cylindrical source (the precipitate) of radius a and height 6 46 47 (the grain boundary thickness). Shewmon has modified the Fisher analysis of grain boundary diffusion so that it is applicable to the same cylindrical symmetry as is involved here. Using Shewmon's solution for the distribution of solute in the grain boundary, Coates calculates the interface flux J* = -DCVTC. The total interface outflow of material from the allotriomorph is simply Jj -oA = 27ra6J , i.e., cf 3yaA)_ u = 2Tra6D0a(CT-Cj K^Caa) g.b. -M-V^g^j. R (33) where a = , .. 6(TTDVT) I DS/DV x is time and KQ and are modified Bessel Functions of zeroth and first order respectively. Now a function \\i is defined such that (f 3V2A)O K Jj_ I g.b. (f J-3A) . I vol. which is to be used as an index of the relative contributions of volume and grain boundary diffusion. - 168 -If ij; >> 1,grain boundary diffusion dominates. If 1/ip » 1,volume diffusion dominates. After substituting equations (1), (2) and (3) into (4), Coates obtains, after considerable rearrangement: DQ/D„6 K1(aa) * " / ' TT^a) f (a/b) 9, _ .1/2 Vaa; 2(TTDVT) where f (a/b) - a^2 ~ ' ^[{U/b)2-!}^2] /(a/b)2-] Notice that the supersaturation (C^-C^) has cancelled out in Equation (6). For a sphere (a/b = 1) and a flat disk (a/b = °°)} the corresponding values of f (a/b) are 1 and TT/2 respectively. Clearly therefore, axial ratio has little influence since f(a/b) is the only te in equation (5) in which axial ratio appears. The factor /TTDT is approximately the characteristic diffusion distance for the matrix. Notice that D appears only in the ratio D /D (cf. equations o o V —0 /RT (3) and (5)). Assume that = D^e , where Q is the activation energy for volume diffusion. If the activation energy for grain boundary diffusion is % Q/2 and the pre-exponential remains unchanged, then VDV > eQ/2RT which can be substituted into equations (3) and (5) to give Coates' final result: - 169 -^2e-Q/2RT 6 0TDVT)1/2 eQ/2RT6 * _____ f(a/b) 2(,D x)1/2 v . /2e-Q/2RT . V K [a / JJ] For a given precipitate size a and diffusion distance (irD T)^2 equation (8) is essentially an equation of the form ty =. KT) i.e. a function describing the relative contributions of volume and grain boundary diffusion as a function of temperature. - 170 -APPENDIX VII. Calculations of the Impingement Resulting from the Interaction of the Confocal Spheroids about an Allo triomorph and a Surrounding Concentric Spherical Field. For the field around the ellipsoid A in Figure 66(c) assume that the diffusion distances are shorter than the sphere radius so that we can consider the precipitate is dissolving in an infinite medium. Consider that the concentration contours are a series of confocal spheroids; in the case of a dissolving precipitate the profiles will be less spherical than this at short times as may be seen in Figure VII-1. As discussed in the text this will make the axial ratio increase faster on dissolution. The equation for confocal oblate Actual Point effect negligable Confocal Spheroid Point effect important Figure VII-1. Schematic diagram comparing actual contour to confocal spheroid. spheroids is given by 69 2 2 + = 1 In two dimensions the contour shapes are defined by: - 171 -The value of g (the focal length) may be obtained directly for the 2 2 case of growth from a plot of the square of the minor axis^.^ -0 )versus the nl axial ratio ^—2 1/2— ^or var-"-ous values of the supersaturation (nl "r?7) (f) by L.C. Brown.' From analogy with growth at f = 0.1 and an axial ratio of 4, g = 1.095 and Is found to vary from 1.13 to 2.82. The actual concentrations are calculated assuming that the profile at the centre of the flat face of the oblate spheroid is an error function as in Figure VTI-2. c=erfc m Figure VII-2. Position of error function profile. Distances are —— = [ (n., 2-g2)1/2 - 0.283] where: B = 1.095, n, = 1.13; . = 4; and (n, '-0 ) = 0.283. The (ni 6 } -10 -1 confocal spheroids are calculated for D - 5 x 10 cm sec and t = 500 seconds. 2/67 = 2 /25 x 10 = 10 y 2 2 1/2 The length of the semi-minor axis will therefore be 10 x (n^ -g ) and the semi-major axis will be 10n^. These values are listed in Table VII-1. - 172 -Table VII-1. Values of (n1 - B )2 and n1 of Confocal Spheroids for Different C/C Ratios. [C/CT] - x 10(n.2-B2)1/2 10n. 1 2/57 11 -1 0 0 2.83 11.3 0.75 0.23 2.3 5.1 12.1 0.5 0.48 4.8 7.6 13.5 0.25 0.81 • 8.1 10.9 15.5 0.1 1.16 11.6 14.4 18.1 0.05 1.39 13.9 16.7 20.0 For diffusion inward from the sphere, again consider that diffusion distances are shorter than the sphere radius so that the presence of the oblate spheroid can be neglected. The profile is of the form oo C/Cj. = a/r I [erfc C2n+1)a~r - erfc (2n+l)a+r ] 4 n=0 2/nt 2/Dt Figure VII-3. Pertaining to Equation 4 - 173 -Calculations for the concentric spheres are made using a plot of C 2 68 versus r/a for different values of Dt/a in Crank (page 86, Figure 6.1). In this' case jDtj = 5 x 10-10 x 500 m 0Q319 a (28 x 10 ) Table VII-2 gives the locations of the concentric spheres for various [C/CT] ratios. Table VII-2. Values of r for Concentric Spheres for Different C/CT Ratios. [C/Cj] r/a 1 1 28.3 0.75 0.89 24.9 0.50 0.79 22.2 0.25 0.66 18.5 0.10 0.51 14.3 0.05 0.42 11.8 The actual profile can be found by summing these two independent profiles as in Figure VII-4. The final profiles are given in Figures 67 and 68 on pages 134-136. - 174 -Figure VII-4. Schematic diagram of the summing up of two independent profiles.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Dissolution of grain boundary allotriomorphs in the...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Dissolution of grain boundary allotriomorphs in the Al-Cu and Al-Ag systems Pasparakis, Aristedes 1972-03-23
pdf
Page Metadata
Item Metadata
Title | Dissolution of grain boundary allotriomorphs in the Al-Cu and Al-Ag systems |
Creator |
Pasparakis, Aristedes |
Publisher | University of British Columbia |
Date Issued | 1972 |
Description | The dissolution behaviour of grain boundary allotriomorphs has been studied in the systems Al-Cu and Al-Ag. Using electron probe microanalysis, isoconcentration contours around dissolving allotriomorphs in an effectively infinite matrix have been determined. The influences of allotriomorph shape, grain boundary misorientation, solute super-saturation in the matrix, volume and grain boundary diffusion, and homologous temperature on the shape of these contours have been examined. It was established that the shape of the isoconcentration contours is dependent only on homologous temperature (T[sub H]) for both systems studied. Above T[sub H] = 0.92, no grain boundary diffusion contribution to the dissolution process is observed. The grain boundary diffusion contribution was found to increase with decreasing T[sub H] and to completely dominate the dissolution process below T[sub H] = 0.72. Using the back scattered electron image on the electron probe microanalyzer, the dissolution rate of individual grain boundary allotriomorphs was determined for various values of T[sub H], under conditions in which there was impingement of diffusion fields from adjacent precipitates. An exponential relationship between axial length or width and dissolution time was found to adequately describe the observed dissolution kinetics. A change in axial ratio accompanied the dissolution of the grain boundary allotriomorphs. At high T[sub H] an increase in axial ratio with dissolution time was observed, whereas at low T[sub H], a decrease in axial ratio (i.e., spheroidization) was observed. A model has been proposed to account for this behaviour. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2011-03-23 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0079022 |
URI | http://hdl.handle.net/2429/32830 |
Degree |
Doctor of Philosophy - PhD |
Program |
Metals and Materials Engineering |
Affiliation |
Applied Science, Faculty of Materials Engineering, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-UBC_1972_A1 P38.pdf [ 12.53MB ]
- Metadata
- JSON: 831-1.0079022.json
- JSON-LD: 831-1.0079022-ld.json
- RDF/XML (Pretty): 831-1.0079022-rdf.xml
- RDF/JSON: 831-1.0079022-rdf.json
- Turtle: 831-1.0079022-turtle.txt
- N-Triples: 831-1.0079022-rdf-ntriples.txt
- Original Record: 831-1.0079022-source.json
- Full Text
- 831-1.0079022-fulltext.txt
- Citation
- 831-1.0079022.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0079022/manifest