This thesis investigated the dissolution of dense non aqueous phase liquids (DNAPL) source zones in the subsurface and the effectiveness of enhanced bioremediation for the treatment of chlorinated ethene DNAPLs, using numerical modeling. For this purpose, an existing multiphase numerical model was extended to include comprehensive models for the processes of dissolution and reaction.
The first part of the thesis examined DNAPL dissolution. First, a thermodynamic-based dissolution model was validated using experimental data from two complex heterogeneous DNAPL releases. Model predictions for DNAPL spatial distribution and effluent concentrations agreed well with experimental measurements, without requiring calibration. This is the first successful application of a predictive dissolution model in the literature. Model results showed the important effects of relative permeability and interfacial areas on dissolution rates. Then, the thermodynamic dissolution model was compared to simpler models typically used in the literature. Five Sherwood-Gilland (SG) empirical correlations were evaluated and their limitations were illustrated. A new dissolution model was proposed that combined the predictive ability of the thermodynamic model and the simplicity of SG models, and is applicable for complex source zones. Lastly, the relationship between the DNAPL source architecture and downstream concentrations was investigated, focusing on multistage concentration profiles. A new upscaled model was proposed that is able to capture such complex behavior.
In the second part of this thesis the thermodynamic dissolution model was combined with a model for reductive dechlorination of chlorinated ethenes to simulate DNAPL bioremediation. Simulations were conducted for simple DNAPL source zones to investigate the impact of dissolution-related processes on bioremediation effectiveness. Dissolution kinetics and back-partitioning of daughter products in the DNAPL were shown to affect dechlorination. Then, the investigation was extended to DNAPL source zones of complex architectures in heterogeneous domains, illustrating the importance of the source zone architecture for the effectiveness of DNAPL bioremediation.
Overall, this thesis presents a comprehensive numerical model that will be an important research tool for evaluating the effectiveness of in-situ bioremediation for DNAPL source zones, and will provide the means for a better understanding and control of the critical factors affecting this technology in the field.Ph