Part IV – April 1969 - Papers - A Numerical Method To Describe the Diffusion-Controlled Growth of Particles When the Diffusion Coefficient Is Composition-Dependent

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