Secondary flow within a river bed and contaminant transport

We have developed a numerical method to simulate the transport of non-sorbing contaminants within the sediment layer of a stream and the leaching of these contaminants in the steam. Typical stream bottom surfaces are uneven with triangularly shaped undulation forms. The flow of the water above such triangular surfaces causes external pressure changes that result in a “pumping effect” and a secondary flow within the sediment. The latter causes a significant contaminant advection within the sediment layer. The flow field in the porous sediment layer is obtained by solving numerically Darcy’s equations. The unsteady mass transfer equation is solved by using a finite-difference method with an up-wind scheme. The effects of parameters, such as channel slope, hydraulic head and dispersion, are studied by quantitatively comparing the numerical results of the total mass flow rate from the contaminant source, the concentration front propagation, and the contaminant mass flow rate into the water column. The “pumping effect,” increases the flow in the vertical direction and, thus, enhances the vertical advective mass transport of the contaminant. This bedform-shape induced flow is largely responsible for the mass transfer of contaminants into the water column. The numerical results also show that the mechanical dispersion inside the sediment bed will significantly increase the contaminant mass flow rate from the source.