Numerical Solution of Fick's Equation with Concentration-Dependent Diffusion Coefficients

Abstract

A method is developed for obtaining numerical solutions to Fick's diffusion equation when the diffusion coefficient is of the form D=D₀Cⁿ, n being an arbitrary positive parameter. The method is applicable to two different types of boundary condition corresponding to frequently encountered physical situations. Some numerical results are presented.