Dense underflow into a stationary reservoir computed. High flow solutions computed and compared with exact solutions. Lower bound on flow rate computed for a number of inflow angles. Non-uniqueness in the solution domain is found near the minimal flow.