ABSTRACT
Many systems are often better described by finite layer geometry. The ability to easily specify arbitrary initial concentration profiles is another important advantage of using the finite layer system models over semi-infinite domain models. The fate and transport mechanisms of diffusion, adsorption, and first-order reactive decay are considered with a variety of surface boundary conditions. A higher number of terms is required for predicting results at short time intervals. Truncation of the series without a sufficient number of terms often gives rise to oscillations in the predicted solution. A simple bisection algorithm is used to find the roots to the transcendental function. Although a number of more efficient approaches can be used to find the roots, the bisection can be reliably applied to all of the root-finding needs in this volume without fears of missing roots.
