<p>In this study, based on an existing analytical solution for the two-dimensional transport of contaminants in a saturated soil layer, for the pulse source, the R program code was developed. The simulation is used to obtain the profiles of contaminant concentration, for a steady groundwater velocity, as a function of distance from the source and time. The problem is analytically solved by leveraging the similarity between the Gaussian (normal) distribution and contaminant concentration distribution, the development of the analytical model (i.e., solution of partial differential equation) by using this similarity is explained step by step since it may be of interest to researchers. Contaminant propagation is modeled using R software, which helps to understand how contaminants migrate through a saturated soil layer. This approach aids in comprehending the mechanisms and spatial dynamics of contaminant dispersion, facilitating the prediction and management of soil and groundwater contamination. The provided R code can be altered for different parameters and time intervals.</p>
Read full abstract