Computational methods for multiphase flow and reactive transport problems arising in subsurface contaminant remediation

A mathematical formulation and some numerical approximation techniques are described for a system of coupled partial differential and algebraic equations describing multiphase flow, transport and interactions of chemical species in the subsurface. A parallel simulator PARSIM has been developed based on these approximation techniques and is being used to study contaminant remediation strategies. Numerical results for a highly complex geochemistry problem involving strontium disposal in a pit at Oak Ridge National Laboratory are presented.