Solute dispersion from continuous point sources in ice-covered turbulent flows with bed absorption

A numerical investigation is conducted to study the steady-state concentration field
when a solute is released from multiple continuous line sources in an ice-covered
channel with an absorbing bed under turbulent flow conditions. The governing
equations are solved using the Crank-Nicolson scheme by adopting a two-power law
velocity and a quartic eddy diffusivity profile, which is influenced by the roughness of
the bed layer and the ice cover. Validation against earlier numerical results for a
specific case reveals strong consistency in the concentration profiles. The findings
highlight how the roughness of the boundaries affects the solute concentration. It
further demonstrates the effect of the bed absorption parameter in the early mixing
stages when solute is released near the bed. For zero bed absorption, the solute
concentration asymptotically attains a uniform far-field value of unity, while any non-
zero bed absorption leads to complete depletion of solute downstream.