impedance boundary condition
To prevent reflection of the pressure wave, the rate of momentum change in the control
volume and the incoming momentum flux have to balance to zero:





Vy: this component is computed from the non-equilibrium bounce-back in vertical direction:


