I'm currently modelling a simple unconfined aquifer system with a river flowing through the center of a square domain of 2000x2000m. I'm attempting to model the surface-river interactions by visualising the drawdown cone and flux data from the river boundary conditions (Cauchy) by varing an abstraction bore from the river.

My question is that I seem to be getting a strange result from FeFLow. I've run a series of models and plotted the flux sourced from the river against time. [CAUCHY_BC(FLUX_IN(+)]/[WELL_BC(FLUX_OUT(-))]*100 from the results file. The problem arises when vary the pumping rate at the same distance I get the same normalised distribution regradless of the pumping rate, which I find a bit odd.

For example:
Say the pumping bore is 100m from the river.
I've got 3 different sets of data for different pumping rates (10 l/s, 5 l/s , 2ls/)
On a plot of water sourced from river against time the pumping rates are identical

Could anyone point out any possible errors conceptially that I'm doing?

Big thanks in advance!