Note: this leads to a change how to set boundary conditions for the momentum balance. Inflow now becomes Dirichlet and Outflow becomes Neumann. This is consistent with the new notation.