Institute of Metals Division - Analysis of Interstitial Diffusion Using Activity Methods

Thermodynamic activity rather than chemical composition is basic to the analysis of diffusion. This is the essential conclusion reached by Darken1-3 and by Birchenall and Mehl.4 If so, it is reasonable to expect that diffusion calculations should be carried out using activities. This paper will illustrate a method of employing activities in analyzing interstitial diffusion. For interstitial diffusion (of carbon in iron and probably generally) the diffusion of the solvent atoms, volume changes, and other such complexities not explicitly provided for in the usual diffusion equations cause a negligible error of the order of 1 pct. Therefore it is relatively easy to test the validity of an analysis based on activities in the instance of interstitial diffusion. The mathematical treatment of the more complex substitutional diffusion could be simplified once the adequacy of the activity method had been established for substitutional diffusion. Thus the present work has the purpose not only of providing a superior analysis of interstitial diffusion but also of indicating a possible approach to the more important question of an adequate treatment of substitutional diffusion. Development of the Diffusion Equation Although activities are convenient quantities to determine experimentally, they are derived from the more fundamental fugacities.9 Possible confusion in using activities can be avoided by basing equations involving their use on the relationship, dP = -D1d/ax dt [I] where, dP = the number of grams of solute crossing one cm2 in the time dt sec. D, = the diffusion constant for use with fugacities; it is assumed to be constant at a given temperature, and has units of secs. f = the fugacity of the solute in the solid solution; the units are those of pressure, dynes per cm2. x = the distance in cm. Eq 1 can be considered to be the basic form of the first Fick Law for one-dimensional diffusion. In order to convert Eq 1 into an equation involving activities, a choice of a standard or reference fugacity, f°, must be made: then. a = f [2] where a is the activity corresponding to the fugacity f. However, the value of a is also determined by the relation, a = a\c [3] where a is the activity coefficient and c is the concentration. For a given value of a (such as a = 1) it is evident that the units used for expressing concentration will affect the numerical value of a corresponding to a given amount of carbon dissolved in y-iron. Since weight per cent concentration is so widely used, this unit will be adopted as the standard in this paper. When the value off given by Eq 2 is substituted in Eq 1, the first Fick Law becomes, dP= -D,f°dl [4] Since D,f° is constant (at a given temperature), the second form of Fick's Law can be obtained in the usual manner, -dt D'f° dx [5] where c' is the concentration in g per cm3. In terms of weight per cent concentration, c, Eq 5 becomes, dc _ 100 d2a Tt = T /f w [6] where p is the density of the solid solution and will be assumed to be constant. In order to simplify Eq6, it is necessary to consider the standard state, f°. Standard states are chosen for convenience. There appears to be one especially convenient choice for the analysis of interstitial diffusion. If the activity coefficient, a, is set equal to unity at infinite dilution of the solute, 100 then — 100D1fl° is approximately equal P to the usual diffusion constant in dilute solutions. This is true since the activity, a, in Eq 6 can be replaced by the concentration, c, with little error in an infinitely dilute solution.* Here f1° is the standard fugacity necessary to achieve this standard state. It is proposed that a diffusion constant Da' be defined, Da1=100/pD1f1° [7] where A is the metal whose diffusion is * The second derivative of a is equal to the second derivative of c in dilute solution only if -da/dr- is. zero. Although theory does not require this condition, in the systems C in Fe, and Zn in Cu, dc is in fact found to be essentially zero.