Removed flux BC from LinearSolver.

The flux options is removed and the user should apply Neumann_rate
instead. Because, it is not clear for the solver to use which area for
the flux. User should convert flux into the rate and then using
Neumann_rate_single and Neumann_rate_group will be easy.