A grid of the channel with 61 × 41 solution points and the corresponding numerical result of the density. (a) The 61 × 41 mesh. (b) Final flow field.