For example, a simple 1D porous flow problem in a reservoir. On the left of the model, there is a production well. So how should I deal with this production well? Adding a Neumann boundary condition at the left? Will it change the physical problem that I am trying to simulate? since if there is a production or injection well, it should be a source term in the mathematical model.