All of the simulations have been
carried out using the numerical software Flux-Expert (Simulog,
60 rue Lavoisier, 38 330 Montbonnot, France) based on the finite
element method. A nonlinear steady-state algorithm was used for
hydrodynamic calculations, and a linear steady-state one was used
for the solving of the diffusion-convection equation. Both
algorithms are based on the Gauss inversion method. In this work,
the equations are solved on a two-dimensional Cartesian geometry
(i.e. rectangular channel with infinite depth) and in a dimensional
form on a Silicon Graphics Indigo 2 workstation.