Institute of Metals Division - On the Solution of Diffusion Problems Involving Concentration-Dependent Diffusion Coefficients

This paper contains solutions of the differential equation of diffusion in binary alloys if the diffusion coefficient is an exponential function of the concentration of one of the components. THE general differential equation for one-dimensional diffusion in a binary system can be written as: where c is the concentration of one of the two components in any suitable units; x, the distance from a reference plane; t, time; and D, the diffusion coefficient. Hartley and Crank&apos; have pointed out that eq 1 with conventional length units for measuring the distance x holds only if the alloy can be regarded as a rigid framework. In general, this is not the case, especially if the lattice parameter varies appreciably with composition. If a modified coordinate system is used instead of an orthodox system based on the centimeter a:: length unit, eq 52 is obtained, which has the same general form as eq 1 and is generally valid. The necessary transformations are presented in the appendix. If&apos; D is independent of concentration, analytical solutions of eq 1 are known for numerous boundary conditions. If D depends on concentration, solutions based on known analytical functions are not available. Among others, Boltzmann2 and Matano3 have shown that the concentration c depends only on (x/t1/2) in the case of diffusion between two semi- infinite spaces with uniform initial concentrations: c = c, for x < 0, t = 0, i.e., x/tl/2 = — x [21] c = c2 for x > 0, t = 0, i.e., x/tl/2 = + a Moreover, Boltzmann and Matano have shown how D may be calculated as a function of concentration over the whole range of composition from c -c, to c = c2 if c as a function of (x/t1/2) for the initial conditions stated in eq 2 has been determined experimentally. Thus, for instance, in principle, the diffusion coefficient in silver-gold alloys may be calculated as a function of concentration for the whole range of composition from only one diffusion experiment with a couple of pure silver and pure gold welded together. If the diffusion coefficient D as a function of concentration c is known, a solution of eq 1 for boundary conditions as stated in eq 2 with arbitrarily chosen values of c, and C2 may be wanted. For instance, if the diffusion coefficient in an alloy A-B as a function of concentration has been derived from an experiment with a diffusion couple involving 0 and 30 pct B as initial. concentrations, the calculation of the concentration distribution in a diffusion couple involving other initial concentrations, e.g., 5 and 20 pct B, may be desired. This problem may be solved by numerical or graphical integration of eq l. The solution of an ordinary differential equation of second order, which results from using (x/t1/2) as an independent variable, can be computed quite rapidly if two conditions for the same value of the independent variable are given. In eq 2, however, two conditions are given for two