On the Transient Atomic/Heat Diffusion in Cylinders and Spheres with Phase Change: A Method to Derive Closed-Form Solutions

Analytical solutions for the transient single-phase and two-phase inward solid-state diffusion and solidification applied to the radial cylindrical and spherical geometries are proposed. Both solutions are developed from the differential equation that treats these phenomena in Cartesian coordinates, which are modified by geometric correlations and suitable changes of variables to achieve closed-form solutions. The modified differential equations are solved by using two well-known closed-form solutions based on the error function, and then equations are obtained to analyze the diffusion interface position as a function of time and position in cylinders and spheres. These exact correlations are inserted into the closed-form solutions for the slab and are used to update the roots for each radial position of the moving boundary interface. The predictions provided by the proposed analytical solutions are validated against the results of a numerical approach. Finally, a comparative study of diffusion in slabs, cylinders, and spheres is also presented for single-phase and two-phase solid-state diffusion and solidification, which shows the importance of the effects imposed by the radial cylindrical and spherical curvatures with respect to the Cartesian coordinate system in the process kinetics. The analytical model is proved to be general, as far as, a semi-infinite solution for diffusion problems with phase change exists in the form of the error function, which enables exact closed-form analytical solutions to be derived, by updating the root at each radial position the moving boundary interface.