An analytical solution for one-dimensional contaminant diffusion through multi-layered system and its applications

An analytical solution for one-dimensional contaminant diffusion through multi-layered media is derived regarding the change of the concentration of contaminants at the top boundary with time. The model accounts for the arbitrary initial conditions and the conditions of zero concentration and zero mass flux on the bottom boundary. The average degree of diffusion of the layered system is introduced on the basis of the solution. The results obtained by the presented analytical solutions agree well with those obtained by the numerical methods presented in the literature papers. The application of the analytical solution to the problem of landfill liner design is illustrated by considering a composite liner consisting of geomembrane and compacted clay liner. The results show that the 100-year mass flux of benzene at the bottom of the composite liner is 45 times higher than that of acetone for the same composite liner. The half-life of the contaminant has a great influence on the solute flux of benzene diffused into the underlying aquifer. Results also indicates that an additional 2.9-5.0 m of the conventional (untreated) compacted clay liner under the geomembrane is required to achieve the same level of protection as provided by 0.60 m of the Hexadecyltrimethylammonium (HDTMA)-treated compacted clay liners in conjunction with the geomembrane. Applications of the solution are also presented in the context of a contaminated two-layered media to demonstrate that different boundary and initial conditions can greatly affect the decontamination rate of the problem. The method is relatively simple to apply and can be used for performing equivalency analysis of landfill liners, preliminary design of groundwater remediation system, evaluating experimental results, and verifying more complex numerical models.