Part XI – November 1968 - Papers - Numerical Solutions to the Finite, Diffusion-Controlled, Two-Phase, Moving-Interface Problem (with Planar, Cylindrical, and Spherical Interfaces)

Treatment of the finite, diffusion-controlled, two-phase, moving-interface homogenization problem has been carried out by numerical methods for a wide range of input parameters. Solutions have been obtained for planar, cylindrical, and spherical diffusion geometries. The effects of parameters such as the interdiffusion coefficients (concentration-independent) of- both phases, the concentrations (time-independent) of each phase at the interlace, mean composition, initial size of the dissolving- phase, and initial composition of- the two phases are considered. The results are displayed to show both the position of the interjace and the degree of homogenization (expressed in terms of concentration at the matrix phase symmetry boundary) as a function of diffusion time. Typical exaw~ples to which the results may be applied are discussed. A limited number of exact, closed-form solutions to the multiphase, diffusion-controlled, moving-interface problem exists in the literature, e.g., Refs. 1 to 7. The exact, closed-form solutions generally consider the outermost phases to be infinite in extent and also assume an interfacial velocity function. It is through use of such boundary conditions that closed-form solutions are readily obtained. An extensive literature also exists on numerical solutions to the mathematically analogous solidification problem, e.g., Refs. 8 to 12. In addition, a numerical solution involving mass transfer with spherical interfaces and infinite boundaries has recently been carried out.13 A practically important class of diffusion problems exists, however, in which the boundaries of each phase away from the interface are finite. For example the finite boundary condition arises when one sets up a model for analyzing the homogenization behavior of two-component powder compacts in two-phase systems or for studying the solution behavior of precipitates. Although exact, closed-form solutions have as yet not been found to this class of problems, solutions are possible using numerical methods. The purpose of this paper is twofold: 1) to describe a numerical method for treating the two-phase, diffusion-controlled, moving-interface problem involving finite boundaries with planar, cylindrical, and spherical interfaces; 2) to present the results of a parametric study which identifies the role played by each variable appearing in the analytical model. R. A. TANZlLLI,Student MemberAIME,onleaveofabsence from the General Electric Co., Re-entry Systems Department and A. W. HECKEL, Member AIME, are Graduate Student in Materials Engineering, and Professor, Department of Metallurgical Engineering, respectively, Drexel Institute of Technology, Philadelphia, Pa. Manuscript submitted May 6. 1968. IMD MATHEMATICAL MODEL The three geometries analyzed are shown in Fig. 1. For each configuration a symmetry element is shown, across whose boundaries the mass transfer is zero. Mean composition, C, and the characteristic solute-rich phase half-thickness (or radius), 1/2, are the independent variables while the characteristic symmetry element half-thickness, L/2 (measured from solute-rich phase center to the matrix phase plane of zero mass transfer), is a dependent variable. The symmetry element for a planar interface accurately represents a real diffusion situation involving the stacking of symmetry elements of constant size. Because cylindrical and spherical symmetry elements of constant size cannot be stacked to fill space, the meaning of the dimension, L/2, for these situations may be considered in terms of the following approximations: and where NA is the number of parallel cylindrical elements per unit area on a surface normal to the cylinder axes and Nv is the number of spherical elements per unit volume. The model consists of a solute-rich phase matrix phase composite of limited mutual solubility. The concentrations at the interface are assumed constant with time; no assumption of local equilibrium at the interface is necessary except for the case where the