**5.2. Numerical model**

The "pressure-based" solver that is used has a SIMPLE-C algorithm implemented (see [33]), together with a multi grid technique for increasing the rate of convergence of stationary flow problems (Ansys Fluent). The spatial discretization is of second order accuracy, with a under-relaxation coefficient of 0.5 for both the pressure and the momentum. The solution has converged when the global <sup>2</sup>*L* - norm of the pressure and the velocity residuals is lower than 10-8.

The Reynolds number sets were selected such that the flow in channel is laminar and stationary.
